A mathematical model for estimating the extent of solute- and water-flux heterogeneity in multiple sample percolation experiments

Multiple sampling is widely used in vadose zone percolation experiments to investigate the extent in which soil structure heterogeneities influence the spatial and temporal distributions of water and solutes. In this note, a simple, robust, mathematical model, based on the beta-statistical distribution, is proposed as a method of quantifying the magnitude of heterogeneity in such experiments. The model relies on fitting two parameters, alpha and zeta to the cumulative elution curves generated in multiple-sample percolation experiments. The model does not require knowledge of the soil structure. A homogeneous or uniform distribution of a solute and/or soil-water is indicated by alpha = zeta = 1, Using these parameters, a heterogeneity index (HI) is defined as root 3 times the ratio of the standard deviation and mean. Uniform or homogeneous flow of water or solutes is indicated by HI = 1 and heterogeneity is indicated by HI > 1. A large value for this index may indicate preferential flow. The heterogeneity index relies only on knowledge of the elution curves generated from multiple sample percolation experiments and is, therefore, easily calculated. The index may also be used to describe and compare the differences in solute and soil-water percolation from different experiments. The use of this index is discussed for several different leaching experiments. (C) 1999 Elsevier Science B.V. All rights reserved.

Tidal effects on sea water intrusion in unconfined aquifers

A variable-density groundwater model is used to analyse the effects of tidal fluctuations on sea-water intrusion in an unconfined aquifer. It is shown that the tidal activity forces the sea-water to intrude further inland and it also creates a thicker interface than would occur without tidal effects. Moreover, the configuration of the interface is radically changed when the tidal fluctuations are included. This is because of the dramatic changes in the flow pattern and velocity of the groundwater near the shoreline. For aquifer depths much larger than tidal amplitudes, the tidal fluctuation does not have much effect on how far the sea-water intrudes into the aquifer; nevertheless, a significant change in the configuration of concentration contours because of the effect of tidal fluctuations is observed. This change is more noticeable at the top of the aquifer, near the water table, than at the bottom of the aquifer, and is caused by the infiltration of salt water into the top of the aquifer at higher tidal levels. A flatter beach slope, therefore, intensifies this phenomenon. The interface configurations do not change noticeably over the course of a tidal cycle. Neglecting tidal fluctuation effects results in an inaccurate evaluation of the water table elevation at the land end of the aquifer, although no distinguishable difference is seen between the water tables near the shoreline. Where the landward boundary condition is a constant head, the effects of tidal fluctuations on sea-water intrusion are more pronounced than for cases where the landward boundary condition is a specified flux. Also it is shown that the effects of tidal fluctuations are more significant for a sloping beach than for a vertical shoreline and the salt water intrudes further inland for the sloping case. (C) 1999 Elsevier Science B.V. All rights reserved.

Integrable Kondo impurities in the one-dimensional supersymmetric U model of strongly correlated electrons

Integrable Kondo impurities in the one-dimensional supersymmetric U model of strongly correlated electrons are studied by means of the boundary graded quantum inverse scattering method. The boundary K-matrices depending on the local magnetic moments of the impurities are presented as non-trivial realizations of the reflection equation algebras in an impurity Hilbert space. Furthermore, the model Hamiltonian is diagonalized and the Bethe ansatz equations are derived. It is interesting to note that our model exhibits a free parameter in the bulk Hamiltonian but no free parameter exists on the boundaries. This is in sharp contrast to the impurity models arising from the supersymmetric t-J and extended Hubbard models where there is no free parameter in the bulk but there is a free parameter on each boundary.

Population balance model of in vivo neutrophil formation following bone marrow rescue therapy

In this paper, we develop a simple four parameter population balance model of in vivo neutrophil formation following bone marrow rescue therapy. The model is used to predict the number and type of neutrophil progenitors required to abrogate the period of severe neutropenia that normally follows a bone marrow transplant. The estimated total number of 5 billion neutrophil progenitors is consistent with the value extrapolated from a human trial. The model provides a basis for designing ex vivo expansion protocols.

V4 region of small subunit rDNA indicates polyphyly of the Fellodistomidae (Digenea) which is supported by morphology and life-cycle data

There is no morphological synapomorphy for the disparate digeneans, the Fellodistomidae Nicoll, 1909. Although all known life-cycles of the group include bivalves as first intermediate hosts, there is no convincing morphological synapomorphy that can be used to unite the group. Sequences from the V4 region of small subunit (18S) rRNA genes were used to infer phylogenetic relationships among 13 species of Fellodistomidae from four subfamilies and eight species from seven other digenean families: Bivesiculidae; Brachylaimidae; Bucephalidae; Gorgoderidae; Gymnophallidae; Opecoelidae; and Zoogonidae. Outgroup comparison was made initially with an aspidogastrean. Various species from the other digenean families were used as outgroups in subsequent analyses. Three methods of analysis indicated polyphyly of the Fellodistomidae and at least two independent radiations of the subfamilies, such that they were more closely associated with other digeneans than to each other. The Tandanicolinae was monophyletic (100% bootstrap support) and was weakly associated with the Gymnophallidae (< 50-55% bootstrap support). Monophyly of the Baccigerinae was supported with 78-87% bootstrap support, and monophyly of the Zoogonidae + Baccigerinae received 77-86% support. The remaining fellodistomid species, Fellodistomum fellis, F. agnotum and Coomera brayi (Fellodistominae) plus Proctoeces maculatus and Complexobursa sp. (Proctoecinae), formed a separate clade with 74-92% bootstrap support. On the basis of molecular, morphological and life-cycle evidence, the subfamilies Baccigerinae and Tandanicolinae are removed from the Fellodistomidae and promoted to familial status. The Baccigerinae is promoted under the senior synonym Faustulidae Poche, 1926, and the Echinobrevicecinae Dronen, Blend & McEachran, 1994 is synonymised with the Faustulidae. Consequently, species that were formerly in the Fellodistomidae are now distributed in three families: Fellodistomidae; Faustulidae (syn. Baccigerinae Yamaguti, 1954); and Tandanicolidae Johnston, 1927. We infer that the use of bivalves as intermediate hosts by this broad range of families indicates multiple host-switching events within the radiation of the Digenea.

Phase diagram of a Heisenberg spin-Peierls model with quantum phonons

Using a new version of the density-matrix renormalization group we determine the phase diagram of a model of an antiferromagnetic Heisenberg spin chain where the spins interact with quantum phonons. A quantum phase transition from a gapless spin-fluid state to a gapped dimerized phase occurs at a nonzero value of the spin-phonon coupling. The transition is in the same universality class as that of a frustrated spin chain, to which the model maps in the diabatic limit. We argue that realistic modeling of known spin-Peierls materials should include the effects of quantum phonons.

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.

Integrable Kondo impurities in the one-dimensional supersymmetric extended Hubbard model

An integrable Kondo problem in the one-dimensional supersymmetric extended Hubbard model is studied by means of the boundary graded quantum inverse scattering method. The boundary K-matrices depending on the local moments of the impurities are presented as a non-trivial realization of the graded reflection equation algebras in a two-dimensional impurity Hilbert space. Further, the model is solved by using the algebraic Bethe ansatz method and the Bethe ansatz equations are obtained.

On modifications to the long-term survival mixture model in the presence of competing risks

A mixture model for long-term survivors has been adopted in various fields such as biostatistics and criminology where some individuals may never experience the type of failure under study. It is directly applicable in situations where the only information available from follow-up on individuals who will never experience this type of failure is in the form of censored observations. In this paper, we consider a modification to the model so that it still applies in the case where during the follow-up period it becomes known that an individual will never experience failure from the cause of interest. Unless a model allows for this additional information, a consistent survival analysis will not be obtained. A partial maximum likelihood (ML) approach is proposed that preserves the simplicity of the long-term survival mixture model and provides consistent estimators of the quantities of interest. Some simulation experiments are performed to assess the efficiency of the partial ML approach relative to the full ML approach for survival in the presence of competing risks.

Using Lindenmayer systems to model morphogenesis in a tropical pasture legume Stylosanthes scabra

This paper summarizes the processes involved in designing a mathematical model of a growing pasture plant, Stylosanthes scabra Vog. cv. Fitzroy. The model is based on the mathematical formalism of Lindenmayer systems and yields realistic computer-generated images of progressive plant geometry through time. The processes involved in attaining growth data, retrieving useful growth rules, and constructing a virtual plant model are outlined. Progressive output morphological data proved useful for predicting total leaf area and allowed for easier quantification of plant canopy size in terms of biomass and total leaf area.

New integrable model of correlated electrons with off-diagonal long-range order from so(5) symmetry

We present a new integrable model for correlated electrons which is based on so(5) symmetry. By using an eta-pairing realization we construct eigenstates of the Hamiltonian with off-diagonal long-range order. It is also shown that these states lie in the ground state sector. We exactly solve the model on a one-dimensional lattice by the Bethe ansatz.

Ameloblast apoptosis and IGF-1 receptor expression in the continuously erupting rat incisor model

Enamel-producing cells (ameloblasts) pass through several phenotypic and functional stages during enamel formation. In the transition between secretory and maturation stages, about one quarter of the ameloblasts suddenly undergo apoptosis. We have studied this phenomenon using the continuously erupting rat incisor model. A special feature of this model is that all stages of ameloblast differentiation are presented within a single longitudinal section of the developing tooth. This permits investigation of the temporal sequence of gene and growth factor receptor expression during ameloblast differentiation and apoptosis. We describe the light and electron microscopic morphology of ameloblast apoptosis and the pattern of insulin-like growth factor-1 receptor expression by ameloblasts in the continuously erupting rat incisor model. In the developing rat incisor, ameloblast apoptosis is associated with downregulated expression of the insulin-like growth factor-1 receptor. These data are consistent with the hypothesis that ameloblasts are hard wired for apoptosis and that insulin-like growth factor-1 receptor expression is required to block the default apoptotic pathway. Possible mechanisms of insulin-like growth factor-1 inhibition of ameloblast apoptosis are presented. The rat incisor model may be useful in studies of physiological apoptosis as it presents apoptosis in a predictable pattern in adult tissues.

Life-cycle preferences over consumption and health: when is cost-effectiveness analysis equivalent to cost-benefit analysis?

This paper studies life-cycle preferences over consumption and health status. We show that cost-effectiveness analysis is consistent with cost-benefit analysis if the Lifetime utility function is additive over time, multiplicative in the utility of consumption and the utility of health status, and if the utility of consumption is constant over time. We derive the conditions under which the lifetime utility function takes this form, both under expected utility theory and under rank-dependent utility theory, which is currently the most important nonexpected utility theory. If cost-effectiveness analysis is consistent with cost-benefit analysis, it is possible to derive tractable expressions for the willingness to pay for quality-adjusted life-years (QALYs). The willingness to pay for QALYs depends on wealth, remaining life expectancy, health status, and the possibilities for intertemporal substitution of consumption. (C) 1999 Elsevier Science B.V. All rights reserved.

RWMODEL II: Computer simulation of the Rescorla-Wagner model of Pavlovian conditioning

RWMODEL II simulates the Rescorla-Wagner model of Pavlovian conditioning. It is written in Delphi and runs under Windows 3.1 and Windows 95. The program was designed for novice and expert users and can be employed in teaching, as well as in research. It is user friendly and requires a minimal level of computer literacy but is sufficiently flexible to permit a wide range of simulations. It allows the display of empirical data, against which predictions from the model can be validated.