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.

Off-diagonal long-range order in a supersymmetric integrable model of correlated electrons

A new model for correlated electrons is presented which is integrable in one-dimension. The symmetry algebra of the model is the Lie superalgebra gl(2\1) which depends on a continuous free parameter. This symmetry algebra contains the eta pairing algebra as a subalgebra which is used to show that the model exhibits Off-Diagonal Long-Range Order in any number of dimensions.

Integrability of a t - J model with impurities

A t - J model for correlated electrons with impurities is proposed. The impurities are introduced in such a way that integrability of the model in one dimension is not violated. The algebraic Bethe ansatz solution of the model is also given and it is shown that the Bethe states are highest weight states with respect to the supersymmetry algebra gl(2/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].

Extended integrability regime for the supersymmetric U model

An extension of the supersymmetric U model for correlated electrons is given and integrability is established by demonstrating that the model can he constructed through the quantum inverse scattering method using an R-matrix without the difference property. Some general symmetry properties of the model are discussed and from the Bethe ansatz solution an expression for the energies is presented.

Synthesis of silica films at the air/water interface: Effect of template chain length and ionic strength

We have grown surfactant-templated silicate films at the air-water interface using n-alkyltrimethylammonium bromide and chloride in an acid synthesis with tetraethyl orthosilicate as the silicate source. The films have been grown with and without added salt (sodium chloride, sodium bromide) and with n-alkyl chain lengths from 12 to 18, the growth process being monitored by X-ray reflectometry. Glassy, hexagonal, and lamellar structures have been produced in ways that are predictable from the pure surfactant-water phase diagrams. The synthesis appears to proceed initially through an induction period characterized by the accumulation of silica-coated spherical micelles near the surface. All syntheses, except those involving C(12)TACl, show a sudden transformation of the spherical micellar phase to a hexagonal phase. This occurs when the gradually increasing ionic strength and/or changing ethanol concentration is sufficient to change the position of boundaries within the phase diagram. A possible mechanism for this to occur may be to induce a sphere to rod transition in the micellar structure. This transformation, as predicted from the surfactant-water phase diagram, can be induced by addition of salts and is slower for chloride than bromide counteranions. The hexagonal materials change in cell dimension as the chain length is changed in a way consistent with theoretical model predictions. All the materials have sufficiently flexible silica frameworks that phase interconversion is observed both from glassy to hexagonal and from hexagonal, to lamellar and vice versa in those surfactant systems where multiple phases are found to exist.

Transfer matrix eigenvalues of the anisotropic multiparametric U model

A multiparametric extension of the anisotropic U model is discussed which maintains integrability. The R-matrix solving the Yang-Baxter equation is obtained through a twisting construction applied to the underlying U-q(sl (2/1)) superalgebraic structure which introduces the additional free parameters that arise in the model. Three forms of Bethe ansatz solution for the transfer matrix eigenvalues are given which we show to be equivalent.

A director theory for visco-elastic folding instabilities in multilayered rock

A model for finely layered visco-elastic rock proposed by us in previous papers is revisited and generalized to include couple stresses. We begin with an outline of the governing equations for the standard continuum case and apply a computational simulation scheme suitable for problems involving very large deformations. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered beam under compression. We analyse folding up to 40% shortening. The standard continuum solution becomes unstable for extreme values of the shear/normal viscosity ratio. The instability is a consequence of the neglect of the bending stiffness/viscosity in the standard continuum model. We suggest considering these effects within the framework of a couple stress theory. Couple stress theories involve second order spatial derivatives of the velocities/displacements in the virtual work principle. To avoid C-1 continuity in the finite element formulation we introduce the spin of the cross sections of the individual layers as an independent variable and enforce equality to the spin of the unit normal vector to the layers (-the director of the layer system-) by means of a penalty method. We illustrate the convergence of the penalty method by means of numerical solutions of simple shears of an infinite layer for increasing values of the penalty parameter. For the shear problem we present solutions assuming that the internal layering is oriented orthogonal to the surfaces of the shear layer initially. For high values of the ratio of the normal-to the shear viscosity the deformation concentrates in thin bands around to the layer surfaces. The effect of couple stresses on the evolution of folds in layered structures is also investigated. (C) 2002 Elsevier Science Ltd. All rights reserved.

Large amplitude folding in finely layered viscoelastic rock structures

We analyze folding phenomena in finely layered viscoelastic rock. Fine is meant in the sense that the thickness of each layer is considerably smaller than characteristic structural dimensions. For this purpose we derive constitutive relations and apply a computational simulation scheme (a finite-element based particle advection scheme; see MORESI et al., 2001) suitable for problems involving very large deformations of layered viscous and viscoelastic rocks. An algorithm for the time integration of the governing equations as well as details of the finite-element implementation is also given. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered plate. The domain is compressed along the layer axis by prescribing velocities along the sides. First, for the viscous limit we consider the response to a series of harmonic perturbations of the director orientation. The Fourier spectra of the initial folding velocity are compared for different viscosity ratios. Turning to the nonlinear regime we analyze viscoelastic folding histories up to 40% shortening. The effect of layering manifests itself in that appreciable buckling instabilities are obtained at much lower viscosity ratios (1:10) as is required for the buckling of isotropic plates (1:500). The wavelength induced by the initial harmonic perturbation of the director orientation seems to be persistent. In the section of the parameter space considered here elasticity seems to delay or inhibit the occurrence of a second, larger wavelength. Finally, in a linear instability analysis we undertake a brief excursion into the potential role of couple stresses on the folding process. The linear instability analysis also provides insight into the expected modes of deformation at the onset of instability, and the different regimes of behavior one might expect to observe.

The game to play: expanding the co-opetition proposal through the strategic games matrix

Purpose - Using Brandenburger and Nalebuff`s 1995 co-opetition model as a reference, the purpose of this paper is to seek to develop a tool that, based on the tenets of classical game theory, would enable scholars and managers to identify which games may be played in response to the different conflict of interest situations faced by companies in their business environments. Design/methodology/approach - The literature on game theory and business strategy are reviewed and a conceptual model, the strategic games matrix (SGM), is developed. Two novel games are described and modeled. Findings - The co-opetition model is not sufficient to realistically represent most of the conflict of interest situations faced by companies. It seeks to address this problem through development of the SGM, which expands upon Brandenburger and Nalebuff`s model by providing a broader perspective, through incorporation of an additional dimension (power ratio between players) and three novel, respectively, (rival, individualistic, and associative). Practical implications - This proposed model, based on the concepts of game theory, should be used to train decision- and policy-makers to better understand, interpret and formulate conflict management strategies. Originality/value - A practical and original tool to use game models in conflict of interest situations is generated. Basic classical games, such as Nash, Stackelberg, Pareto, and Minimax, are mapped on the SGM to suggest in which situations they Could be useful. Two innovative games are described to fit four different types of conflict situations that so far have no corresponding game in the literature. A test application of the SGM to a classic Intel Corporation strategic management case, in the complex personal computer industry, shows that the proposed method is able to describe, to interpret, to analyze, and to prescribe optimal competitive and/or cooperative strategies for each conflict of interest situation.

Family Business: how family and ownership shapes business professionalization

This article`s main purpose consists in showing how family and ownership cultures may influence the process of making a ""well-performing"" organization, based on an empirical study in family business in Brazil. The study aimed to find critical moments of company`s history and the focus was to compare critical moments with the three-dimension model of family business development proposed by Davis et al. (1996). Through facts sequence, research was organized so as to find how the process influenced company`s professionalization. The article concludes that family and its value and culture may impact on the evolution, and the first step to organize a company is to organize the family that leads the company.

In the eye of the storm: exploring the introduction of environmental issues in the production function in Brazilian companies

This article aims to analyse the introduction of environmental issues in the context of the production function, which has been referred to as the organisational area to lead corporate environmental management. With that purpose, the theoretical references for corporate environmental management and the necessary alterations in production function have been organised to include environmental aspects, especially in terms of product and process development, quality management, and logistics. Considering that this research field still lacks empirical evidence for Brazilian companies, four case studies were conducted using companies located in the country. The environmental management maturity level of those companies tends to follow the rate with which the environmental issue is introduced in production sub-areas, especially in the product development process. However, in most cases we found that the companies had difficulties in structuring the insertion of the environmental dimension in logistics. The final notes point out the distance observed between what is recommended by international literature and the reality of Brazilian companies in the challenge of making the production function environmentally friendly.

Shell utilization by the land hermit crab Coenobita scaevola (Anomura, Coenobitidae) from Wadi El-Gemal, Red Sea

We conducted a study in order to determine the shell utilization pattern of the land hermit crab Coenobita scaevola (Forskal, 1775), the only species representing the family Coenobitidae in the Red Sea. Hermit crabs were collected during July 2003 and January 2004 along the sandy shores of protected area of Wadi-Elgemal, south Red Sea. Animals were fixed in 10% formalin and transported to the laboratory where they were weighed and measured for cephalothoracic shield length (CSL) and width, left propodus length and height. Gastropod shells species were identified, weighed and measured for shell aperture width and length and shell internal volume. A total of 391 individuals were collected (219 females, 172 males) and were found occupying ten shell species, with clear significant occupation of Nerita undata. A positive relationship was obtained between the size of the shells occupied and the hermit crabs. Analysis of shell internal volume and crab dimensions demonstrated that this shell dimension constitutes mainly the determinant for C. scaevola shell utilization. With respect to the size of the animals and the occupied shell type, Nerita undata was occupied by a wide range of CSL (2.5-8.5mm). Small sized crabs (2.5-3.5mm CSL) occupied Planaxis sulcatus and Nassarius arcularius plicatus while larger specimens (8.5-9.5mm CSL) occupied Turbo radiatits, Polinices milanostomus and Monodonta canilifera. Variations in the shell occupation were also recognized among male and females. Comparisons among populational and shell use features led us to suggest the use of this land hermit crab as key-species in the preserving program of shores and protected areas, since this species is the first organism to disappear from any shore when a new tourist establishment is implemented.

Reward drive and rash impulsiveness as dimensions of impulsivity: Implications for substance misuse

One of the primary personality dimensions or traits that has consistently been linked to substance abuse is impulsivity. However, impulsivity is not a homogenous construct and although many of the measures of impulsivity are correlated, the most recent review of published factor analytic studies has proposed two independent dimensions of impulsivity: reward sensitivity, reflecting one of the primary dimension of J. A. Gray's personality theory, and rash impulsiveness. These two facets of impulsivity derived from the field of personality research parallel recent developments in the neurosciences where changes in the incentive value of rewarding substances has been linked to alterations in neural substrates involved in reward seeking and with a diminished capacity to inhibit behavior due to chronic drug exposure. In this paper, we propose a model that integrates the findings from research into individual differences with recent models of neural substrates implicated in the development of substance misuse. (C) 2004 Elsevier Ltd. All rights reserved.

Dendrimer intercalation in layered double hydroxides

New organic/inorganic (O/I) hybrid assemblies based on Layered Double Hydroxide (LDH) with polyamide amine dendrimer (PAMAM, generation -0.5 and generation +0.5) were prepared by two different routes using either the direct coprecipitation at constant pH or the anion exchange procedure in double surfactant S(+)S(-) phases. The obtained materials were characterized by means of powder X-ray diffraction, thermal gravimetric analysis associated with mass spectrometry, and Fourier-transform infrared spectroscopy. X-ray powder diffraction pattern of the O/I LDH assembly exhibit characteristic profiles of LDH-based materials with basal spacing depending on the nature of the dendrimer. Indeed, for both synthetic procedures, interleaved PAMAM -0.5 gives rise to an interlayer space in agreement with a perpendicular molecular arrangement against the layer of the host structure. For PAMAM+0.5, considering its spherical dimension, a much smaller basal spacing was observed. This observation was interpreted as shrinkage of the molecule to accommodate the interlayer LDH gap, which was rendered possible by the bond angle twisting within PAMAM-0.5. FTIR spectra confirm the presence of both moieties inside both Zn(2)Al/PAMAM G-0.5 and Zn(2)Al/PAMAM G+0.5 assemblies. Finally, thermal analysis associated with mass spectrometry confirm this composition, and in situ temperature XRD data reveal that the highly constrained arrangement for the generation +0.5 is not accompanied by a gain in thermal structural stability; in fact, the assembly prepared from PAMAM -0.5 is more stable. Both O/I PAMAM LDH assemblies constitute well-defined materials which are candidate for catalytic applications.