Expokit: A software package for computing matrix exponentials

Expokit provides a set of routines aimed at computing matrix exponentials. More precisely, it computes either a small matrix exponential in full, the action of a large sparse matrix exponential on an operand vector, or the solution of a system of linear ODEs with constant inhomogeneity. The backbone of the sparse routines consists of matrix-free Krylov subspace projection methods (Arnoldi and Lanczos processes), and that is why the toolkit is capable of coping with sparse matrices of large dimension. The software handles real and complex matrices and provides specific routines for symmetric and Hermitian matrices. The computation of matrix exponentials is a numerical issue of critical importance in the area of Markov chains and furthermore, the computed solution is subject to probabilistic constraints. In addition to addressing general matrix exponentials, a distinct attention is assigned to the computation of transient states of Markov chains.

Erosion/productivity modelling of maize farming in the Philippine uplands part II: Simulation of alternative farming methods

A version of the Agricultural Production Systems Simulator (APSIM) capable of simulating the key agronomic aspects of intercropping maize between legume shrub hedgerows was described and parameterised in the first paper of this series (Nelson et al., this issue). In this paper, APSIM is used to simulate maize yields and soil erosion from traditional open-field farming and hedgerow intercropping in the Philippine uplands. Two variants of open-field farming were simulated using APSIM, continuous and fallow, for comparison with intercropping maize between leguminous shrub hedgerows. Continuous open-field maize farming was predicted to be unsustainable in the long term, while fallow open-field farming was predicted to slow productivity decline by spreading the effect of erosion over a larger cropping area. Hedgerow intercropping was predicted to reduce erosion by maintaining soil surface cover during periods of intense rainfall, contributing to sustainable production of maize in the long term. In the third paper in this series, Nelson et al. (this issue) use cost-benefit analysis to compare the economic viability of hedgerow intercropping relative to traditional open-field farming of maize in relatively inaccessible upland areas. (C) 1998 Elsevier Science Ltd. All rights reserved.

Erosion/productivity modelling of maize farming in the Philippine uplands part III: Economic analysis of alternative farming methods

Two previous papers in this series (Nelson et al., this issue) described the use of the Agricultural Production Systems Simulator (APSIM) to simulate the effect of erosion on maize yields from open-field farming and hedgerow intercropping in the Philippine uplands. In this paper, maize yields simulated with APSIM are used to compare the economic viability of intercropping maize between leguminous shrub hedgerows with that of continuous and fallow open-field farming of maize. The analysis focuses on the economic incentives of upland farmers to adopt hedgerow intercropping, discussing farmers' planning horizons, access to credit and security of land tenure, as well as maize pricing in the Philippines. Insecure land tenure has limited the planning horizons of upland farmers, and high establishment costs reduce the economic viability of hedgerow intercropping relative to continuous and fallow open-field farming in the short term, In the long term, high discount rates and share-tenancy arrangements in which landlords do not contribute to establishment costs reduce the economic viability of hedgerow intercropping relative to fallow open-field farming, (C) 1998 Elsevier Science Ltd. All rights reserved.

The Heisenberg antiferromagnet on an anisotropic triangular lattice: linear spin-wave theory

We consider the effect of quantum spin fluctuations on the ground-state properties of the Heisenberg antiferromagnet on an anisotropic triangular lattice using linear spin-wave (LSW) theory. This model should describe the magnetic properties of the insulating phase of the kappa-(BEDT-TTF)(2)X family of superconducting molecular crystals. The ground-state energy, the staggered magnetization, magnon excitation spectra, and spin-wave velocities are computed as functions of the ratio of the antiferromagnetic exchange between the second and first neighbours, J(2)/J(1). We find that near J(2)/J(1) = 0.5, i.e., in the region where the classical spin configuration changes from a Neel-ordered phase to a spiral phase, the staggered magnetization vanishes, suggesting the possibility of a quantum disordered state. in this region, the quantum correction to the magnetization is large but finite. This is in contrast to the case for the frustrated Heisenberg model on a square lattice, for which the quantum correction diverges logarithmically at the transition from the Neel to the collinear phase. For large J(2)/J(1), the model becomes a set of chains with frustrated interchain coupling. For J(2) > 4J(1), the quantum correction to the magnetization, within LSW theory, becomes comparable to the classical magnetization, suggesting the possibility of a quantum disordered state. We show that, in this regime, the quantum fluctuations are much larger than for a set of weakly coupled chains with non-frustrated interchain coupling.

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].

Social identity and self-categorization processes in organizational contexts

Although aspects of social identity theory are familiar to organizational psychologists, its elaboration, through self-categorization theory, of how social categorization and prototype-based depersonalization actually produce social identity effects is less well known. We describe these processes, relate self-categorization theory to social identity theory, describe new theoretical developments in detail, and show how these developments can address a: range of organizational phenomena. We discuss cohesion and deviance, leadership, subgroup and sociodemographic structure, and mergers and acquisitions.

Effect of density on spatial distribution, pit formation and pit diameter of Myrmeleon acer Walker, (Neuroptera : Myrmeleontidae): patterns and processes

The effect of increasing population density on the formation of pits, their size and spatial distribution, and on levels of mortality was examined in the antlion Myrmeleon acer Walker. Antlions were kept at densities ranging from 0.4 to 12.8 individuals per 100 cm(2). The distribution of pits was regular or uniform across all densities, but antlions constructed proportionally fewer and smaller pits as density increased. Mortality through cannibalism was very low and only occurred at densities greater than five individuals per 100 cm(2). Antlions in artificially crowded situations frequently relocated their pits and when more space became available, individuals became more dispersed with time. Redistribution of this species results from active avoidance of other antlions and sand throwing associated with pit construction and maintenance, rather than any attempt to optimise prey capture per se.

Decoherence and coherent population transfer between two coupled systems

We show that an arbitrary system described by two dipole moments exhibits coherent superpositions of internal states that can be completely decoupled fi om the dissipative interactions (responsible for decoherence) and an external driving laser field. These superpositions, known as dark or trapping states, can he completely stable or can coherently interact with the remaining states. We examine the master equation describing the dissipative evolution of the system and identify conditions for population trapping and also classify processes that can transfer the population to these undriven and nondecaying states. It is shown that coherent transfers are possible only if the two systems are nonidentical, that is the transitions have different frequencies and/or decay rates. in particular, we find that the trapping conditions can involve both coherent and dissipative interactions, and depending on the energy level structure of the system, the population can be trapped in a linear superposition of two or more bare states, a dressed state corresponding to an eigenstate of the system plus external fields or, in some cases. in one of the excited states of the system. A comprehensive analysis is presented of the different processes that are responsible for population trapping, and we illustrate these ideas with three examples of two coupled systems: single V- and Lambda-type three-level atoms and two nonidentical tao-level atoms, which are known to exhibit dark states. We show that the effect of population trapping does not necessarily require decoupling of the antisymmetric superposition from the dissipative interactions. We also find that the vacuum-induced coherent coupling between the systems could be easily observed in Lambda-type atoms. Our analysis of the population trapping in two nonidentical atoms shows that the atoms can be driven into a maximally entangled state which is completely decoupled from the dissipative interaction.

Kinematic and non-linear analysis of foldable barrel vaults

Kinematic analysis is conducted to derive the geometric constraints for the geometric design of foldable barrel vaults (FBV) composed of polar or angulated scissor units. Non-linear structural analysis is followed to determine the structural response of FBVs in the fully deployed configuration under static loading. Two load cases are considered: cross wind and longitudinal wind. The effect of varying member sizes, depth-to-span ratio and geometric imperfections is examined. (C) 2000 Elsevier Science Ltd. All rights reserved.

Quantum measurement and stochastic processes in mesoscopic conductors

In quantum measurement theory it is necessary to show how a, quantum source conditions a classical stochastic record of measured results. We discuss mesoscopic conductance using quantum stochastic calculus to elucidate the quantum nature of the measurement taking place in these systems. To illustrate the method we derive the current fluctuations in a two terminal mesoscopic circuit with two tunnel barriers containing a single quasi bound state on the well. The method enables us to focus on either the incoming/ outgoing Fermi fields in the leads, or on the irreversible dynamics of the well state itself. We show an equivalence between the approach of Buttiker and the Fermi quantum stochastic calculus for mesoscopic systems.

Dual processes in recognition: Does a focus on measurement operations provide a sufficient foundation?

Current theoretical thinking about dual processes in recognition relies heavily on the measurement operations embodied within the process dissociation procedure. We critically evaluate the ability of this procedure to support this theoretical enterprise. We show that there are alternative processes that would produce a rough invariance in familiarity (a key prediction of the dual-processing approach) and that the process dissociation procedure does not have the power to differentiate between these alternative possibilities. We also show that attempts to relate parameters estimated by the process dissociation procedure to subjective reports (remember-know judgments) cannot differentiate between alternative dual-processing models and that there are problems with some of the historical evidence and with obtaining converging evidence. Our conclusion is that more specific theories incorporating ideas about representation and process are required.

The relationship between radiation-induced chemical processes and transverse relaxation times in polymer gel dosimeters

The effects of ionizing radiation in different compositions of polymer gel dosimeters are investigated using FT-Raman spectroscopy and NMR T-2 relaxation times. The dosimeters are manufactured from different concentrations of comonomers (acrylamide and N,N'-methylene-bis-acrylamide) dispersed in different concentrations of an aqueous gelatin matrix. Results are analysed using a model of fast exchange of magnetization between three proton pools. The fraction of protons in each pool is determined using the known chemical composition of the dosimeter and FT-Raman spectroscopy. Based on these results, the physical and chemical processes in interplay in the dosimeters are examined in view of their effect on the changes in T-2 The precipitation of growing macroradicals and the scavenging of free radicals by gelatin are used to explain the rate of polymerization. The model describes the changes in T-2 as a function of the absorbed dose up to 50 Gy for the different compositions. This is expected to aid the theoretical design of new, more efficient dosimeters, since it was demonstrated that the optimum dosimeter (i.e, with the lowest dose resolution) must have a range of relaxation times which match the range of T-2 values which can be determined with the lowest uncertainty using an MRI scanner.