Quantum nonlinear dynamics of continuously measured systems

Classical dynamics is formulated as a Hamiltonian flow in phase space, while quantum mechanics is formulated as unitary dynamics in Hilbert space. These different formulations have made it difficult to directly compare quantum and classical nonlinear dynamics. Previous solutions have focused on computing quantities associated with a statistical ensemble such as variance or entropy. However a more diner comparison would compare classical predictions to the quantum predictions for continuous simultaneous measurement of position and momentum of a single system, in this paper we give a theory of such measurement and show that chaotic behavior in classical systems fan be reproduced by continuously measured quantum systems.

Integrable anisotropic spin-ladder model

We present an integrable spin-ladder model, which possesses a free parameter besides the rung coupling J. Wang's system based on the SU(4) symmetry can be obtained as a special case. The model is exactly solvable by means of the Bethe ansatz method. We determine the dependence on the anisotropy parameter of the phase transition between gapped and gapless spin excitations and present the phase diagram. Finally, we show that the model is a special case of a more general Hamiltonian with three free parameters.

Integrable open boundary conditions for the Bariev model of three coupled X Y spin chains

The integrable open-boundary conditions for the Bariev model of three coupled one-dimensional XY spin chains are studied in the framework of the boundary quantum inverse scattering method. Three kinds of diagonal boundary K-matrices leading to nine classes of possible choices of boundary fields are found and the corresponding integrable boundary terms are presented explicitly. The boundary Hamiltonian is solved by using the coordinate Bethe ansatz technique and the Bethe ansatz equations are derived. (C) 2001 Elsevier Science B.V. All rights reserved.

Can increasing the viscosity of formulations be used to reduce the human skin penetration of the sunscreen oxybenzone?

The effect of adding thickening agents on the penetration of a sunscreen benzophenone-3 through epidermal and a high-density polyethylene membrane was studied using both very thick (infinite dose) and thin tin use) applications. Contradictory results were obtained. Thickening agents retard skin penetration, in a manner consistent with a diffusional resistance in the formulation, when applied as an infinite dose. In contrast, when applied as in thin (in use) doses, thickening agents promote penetration, most likely through greater stratum corneum diffusivity arising from an enhanced hydration by the thicker formulations. The two key implications from this work are (i) a recognition of the danger in the potential extrapolation of infinite dosing to in use situations, and (ii) to recognize that thicker formulations may sometimes enhance the penetration of other topical agents when applied in use.

Ligand-field analysis of an Er(III) complex with a heptadentate tripodal N4O3 ligand

Polarized absorption and emission spectra of trigonal single crystals of an Er(III) complex coordinated to a heptadentate tripodal ligand are reported at temperatures between 8 and 298 K. The assigned energy levels below the onset of ligand absorption (< 25 000 cm(-1)) are fitted to a parametrized electronic Hamiltonian. The C-3 site symmetry of the Er(HI) ion requires eight parameters for a full description of the ligand field within a one-electron operator description. This compound shows unusually large splittings of the multiplets, and the fitted parameters imply that this heptadentate ligand imparts the largest ligand field reported for an Er(III) complex. The ligand field was also interpreted within the angular overlap model (AOM). We derive the AOM matrix to include both sigma and anisotropic pi bonding and show that a useful description of the C-3 ligand field can be made using only five parameters. The success of the AOM description is encouraging for applications on isomorphous complexes within the lanthanide series and in describing the ligand field of low-symmetry complexes with less parameters than in the usual spherical harmonic expansion.

Derivation of the probability distribution function for the local density of states of a disordered quantum wire via the replica trick and supersymmetry

We consider the statistical properties of the local density of states of a one-dimensional Dirac equation in the presence of various types of disorder with Gaussian white-noise distribution. It is shown how either the replica trick or supersymmetry can be used to calculate exactly all the moments of the local density of states.' Careful attention is paid to how the results change if the local density of states is averaged over atomic length scales. For both the replica trick and supersymmetry the problem is reduced to finding the ground state of a zero-dimensional Hamiltonian which is written solely in terms of a pair of coupled spins which are elements of u(1, 1). This ground state is explicitly found for the particular case of the Dirac equation corresponding to an infinite metallic quantum wire with a single conduction channel. The calculated moments of the local density of states agree with those found previously by Al'tshuler and Prigodin [Sov. Phys. JETP 68 (1989) 198] using a technique based on recursion relations for Feynman diagrams. (C) 2001 Elsevier Science B.V. All rights reserved.

Lanczos subspace filter diagonalization: Homogeneous recursive filtering and a low-storage method for the calculation of matrix elements

We develop a new iterative filter diagonalization (FD) scheme based on Lanczos subspaces and demonstrate its application to the calculation of bound-state and resonance eigenvalues. The new scheme combines the Lanczos three-term vector recursion for the generation of a tridiagonal representation of the Hamiltonian with a three-term scalar recursion to generate filtered states within the Lanczos representation. Eigenstates in the energy windows of interest can then be obtained by solving a small generalized eigenvalue problem in the subspace spanned by the filtered states. The scalar filtering recursion is based on the homogeneous eigenvalue equation of the tridiagonal representation of the Hamiltonian, and is simpler and more efficient than our previous quasi-minimum-residual filter diagonalization (QMRFD) scheme (H. G. Yu and S. C. Smith, Chem. Phys. Lett., 1998, 283, 69), which was based on solving for the action of the Green operator via an inhomogeneous equation. A low-storage method for the construction of Hamiltonian and overlap matrix elements in the filtered-basis representation is devised, in which contributions to the matrix elements are computed simultaneously as the recursion proceeds, allowing coefficients of the filtered states to be discarded once their contribution has been evaluated. Application to the HO2 system shows that the new scheme is highly efficient and can generate eigenvalues with the same numerical accuracy as the basic Lanczos algorithm.

Subspace wavepacket evolution with Newton polynomials

Time-dependent wavepacket evolution techniques demand the action of the propagator, exp(-iHt/(h)over-bar), on a suitable initial wavepacket. When a complex absorbing potential is added to the Hamiltonian for combating unwanted reflection effects, polynomial expansions of the propagator are selected on their ability to cope with non-Hermiticity. An efficient subspace implementation of the Newton polynomial expansion scheme that requires fewer dense matrix-vector multiplications than its grid-based counterpart has been devised. Performance improvements are illustrated with some benchmark one and two-dimensional examples. (C) 2001 Elsevier Science B.V. All rights reserved.

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.

Solution techniques for transport problems involving steep concentration gradients: application to noncatalytic fluid solid reactions

Some efficient solution techniques for solving models of noncatalytic gas-solid and fluid-solid reactions are presented. These models include those with non-constant diffusivities for which the formulation reduces to that of a convection-diffusion problem. A singular perturbation problem results for such models in the presence of a large Thiele modulus, for which the classical numerical methods can present difficulties. For the convection-diffusion like case, the time-dependent partial differential equations are transformed by a semi-discrete Petrov-Galerkin finite element method into a system of ordinary differential equations of the initial-value type that can be readily solved. In the presence of a constant diffusivity, in slab geometry the convection-like terms are absent, and the combination of a fitted mesh finite difference method with a predictor-corrector method is used to solve the problem. Both the methods are found to converge, and general reaction rate forms can be treated. These methods are simple and highly efficient for arbitrary particle geometry and parameters, including a large Thiele modulus. (C) 2001 Elsevier Science Ltd. All rights reserved.

Dynamical tunnelling of ultracold atoms

The divergence of quantum and classical descriptions of particle motion is clearly apparent in quantum tunnelling(1,2) between two regions of classically stable motion. An archetype of such nonclassical motion is tunnelling through an energy barrier. In the 1980s, a new process, 'dynamical' tunnelling(1-3), was predicted, involving no potential energy barrier; however, a constant of the motion (other than energy) still forbids classically the quantum-allowed motion. This process should occur, for example, in periodically driven, nonlinear hamiltonian systems with one degree of freedom(4-6). Such systems may be chaotic, consisting of regions in phase space of stable, regular motion embedded in a sea of chaos. Previous studies predicted(4) dynamical tunnelling between these stable regions. Here we observe dynamical tunnelling of ultracold atoms from a Bose-Einstein condensate in an amplitude-modulated optical standing wave. Atoms coherently tunnel back and forth between their initial state of oscillatory motion (corresponding to an island of regular motion) and the state oscillating 180 degrees out of phase with the initial state.

A new diffusion and flow theory for activated carbon from low pressure to capillary condensation range

In this paper, we develop a theory for diffusion and flow of pure sub-critical adsorbates in microporous activated carbon over a wide range of pressure, ranging from very low to high pressure, where capillary condensation is occurring. This theory does not require any fitting parameter. The only information needed for the prediction is the complete pore size distribution of activated carbon. The various interesting behaviors of permeability versus loading are observed such as the maximum permeability at high loading (occurred at about 0.8-0.9 relative pressure). The theory is tested with diffusion and flow of benzene through a commercial activated carbon, and the agreement is found to be very good in the light that there is no fitting parameter in the model. (C) 2001 Elsevier Science B.V. All rights reserved.

Nucleation, growth and breakage phenomena in agitated wet granulation processes: a review

Wet agglomeration processes have traditionally been considered an empirical art, with great difficulties in predicting and explaining observed behaviour. Industry has faced a range of problems including large recycle ratios, poor product quality control, surging and even the total failure of scale up from laboratory to full scale production. However, in recent years there has been a rapid advancement in our understanding of the fundamental processes that control granulation behaviour and product properties. This review critically evaluates the current understanding of the three key areas of wet granulation processes: wetting and nucleation, consolidation and growth, and breakage and attrition. Particular emphasis is placed on the fact that there now exist theoretical models which predict or explain the majority of experimentally observed behaviour. Provided that the correct material properties and operating parameters are known, it is now possible to make useful predictions about how a material will granulate. The challenge that now faces us is to transfer these theoretical developments into industrial practice. Standard, reliable methods need to be developed to measure the formulation properties that control granulation behaviour, such as contact angle and dynamic yield strength. There also needs to be a better understanding of the flow patterns, mixing behaviour and impact velocities in different types of granulation equipment. (C) 2001 Elsevier Science B.V. All rights reserved.

A semi-analytical solution for vibration of rectangular plates with abrupt thickness variation

A semi-analytical analysis of free vibration of plates with cross-sectional discontinuities due to abrupt changes in thickness is presented. A basic square element divided into suitable subdomains dependent upon the positions of these abrupt changes is used as the basic building element, Admissible functions that satisfy the essential or geometric boundary conditions are used to define the transverse deflection of each subdomain. Continuities in the displacement, slope, moment and higher derivatives between adjacent subdomains are enforced at the interconnecting edges. The resulting global energy functional from the proper assembly of the coupled strain and kinetic energy contributions of each subdomain is then minimized via the Ritz procedure to extract the frequencies and mode shapes. Contour plots of a range of new mode shapes are presented for the enhancement of understanding the dynamic behavior of this class of plates, (C) 2001 Elsevier Science Ltd, All rights reserved.

Knowledge and markets

An economy is a coordinated system of distributed knowledge. Economic evolution occurs as knowledge grows and the structure of the system changes. This paper is about the role of markets in this process. Traditionally, the theory of markets has not been a central feature of evolutionary economics. This seems to be due to the orthodox view of markets as information-processing mechanisms for finding equilibria. But in economic evolution markets are actually knowledge-structuring mechanisms. What then is the relation between knowledge, information, markets and mechanisms? I argue that an evolutionary theory of markets, in the manner of Loasby (1999), requires a clear formulation of these relations. I suggest that a conception of knowledge and markets in terms of a graphical theory of complex systems furnishes precisely this.