Metal-insulator transition and charge ordering in the extended Hubbard model at one-quarter filling

We study, with exact diagonalization, the zero temperature properties of the quarter-filled extended Hubbard model on a square lattice. We find that increasing the ratio of the intersite Coulomb repulsion, V, to the bandwidth drives the system from a metal to a charge ordered insulator. The evolution of the optical conductivity spectrum with increasing V is in agreement with the observed optical conductivity of several layered molecular crystals with the theta and beta crystal structures.

Predicting plant leaf area production: shoot assimilate accumulation and partitioning, and leaf area ratio, are stable for a wide range of sorghum population densities

Predicting plant leaf area production is required for modelling carbon balance and tiller dynamics in plant canopies. Plant leaf area production can be studied using a framework based on radiation intercepted, radiation use efficiency (RUE) and leaf area ratio (LAR) (ratio of leaf area to net above-ground biomass). The objective of this study was to test this framework for predicting leaf area production of sorghum during vegetative development by examining the stability of the contributing components over a large range of plant density. Four densities, varying from 2 to 16 plants m(-2), were implemented in a field experiment. Plants were either allowed to tiller or were maintained as uniculm by systematic tiller removal. In all cases, intercepted radiation was recorded daily and leaf area and shoot dry matter partitioning were quantified weekly at individual culm level. Up to anthesis, a unique relationship applied between fraction of intercepted radiation and leaf area index, and between shoot dry weight accumulation and amount of intercepted radiation, regardless of plant density. Partitioning of shoot assimilate between leaf, stem and head was also common across treatments up to anthesis, at both plant and culm levels. The relationship with thermal time (TT) from emergence of specific leaf area (SLA) and LAR of tillering plants did not change with plant density. In contrast, SLA of uniculm plants was appreciably lower under low-density conditions at any given TT from emergence. This was interpreted as a consequence of assimilate surplus arising from the inability of the plant to compensate by increasing the leaf area a culm could produce. It is argued that the stability of the extinction coefficient, RUE and plant LAR of tillering plants observed in these conditions provides a reliable way to predict leaf area production regardless of plant density. Crown Copyright (C) 2002 Published by Elsevier Science B.V. All rights reserved.

Microhabitat use by the rainbowfish Melanotaenia duboulayi in a subtropical Australian stream

Microhabitat use and feeding behavior of the rainbowfish Melanotaenia duboulayi (Castelnau) were investigated in a slow-flowing stream adjacent to riparian forest in south-eastern Queensland, Australia. Fish were more abundant in vegetated areas, but did not enter dense Vallisneria beds, where predators were observed. In sunny conditions shoals of juveniles occurred near the water surface feeding floating material on the surface, but larger fish tended to occur at the bottom near submerged vegetation, often utilizing the overhanging aquatic plant community as a refuge and food source. In the middle of the day, juveniles and small fish seemed to show behavioral thermoregulation at the surface in the warmest site. Under cloudy conditions, however, fish of all sizes preferred deeper water. The present study suggests that in still and sunny pools thermal change caused by sunlight influences the microhabitat choice of small fish. A field experiment using a kingfisher model implies that fish swimming at the surface could escape from aerial predators in sunlit conditions by responding to moving shadows, but could not do so under cloudy conditions.

Habitat-predator association and avoidance in rainbowfish (Melanotaenia spp.)

The ability to recall the location of a predator and later avoid it was tested in nine populations of rainbowfish (Melanotaenia spp.), representing three species from a variety of environments. Following the introduction of a model predator into a particular microhabitat, the model was removed, the arena rotated and the distribution of the fish recorded again. In this manner it could be determined what cues the fish relied on in order to recall the previous location of the predator model. Fish from all populations but one (Dirran Creek) were capable of avoiding the predator by remembering either the location and/or the microhabitat in which the predator was recently observed. Reliance on different types of visual cues appears to vary between populations but the reason for this variation remains elusive. Of the ecological variables tested (flow variability, predator density and habitat complexity), only the level of predation appeared to be correlated with the orientation technique employed by each population. There was no effect of species identity, which suggests that the habitat that each population occupies plays a strong role in the development of both predator avoidance responses and the cues used to track predators in the wild.

A unified model of flames as gasdynamic discontinuities

Viewed on a hydrodynamic scale, flames in experiments are often thin so that they may be described as gasdynamic discontinuities separating the dense cold fresh mixture from the light hot burned products. The original model of a flame as a gasdynamic discontinuity was due to Darrieus and to Landau. In addition to the fluid dynamical equations, the model consists of a flame speed relation describing the evolution of the discontinuity surface, and jump conditions across the surface which relate the fluid variables on the two sides of the surface. The Darrieus-Landau model predicts, in contrast to observations, that a uniformly propagating planar flame is absolutely unstable and that the strength of the instability grows with increasing perturbation wavenumber so that there is no high-wavenumber cutoff of the instability. The model was modified by Markstein to exhibit a high-wavenumber cutoff if a phenomenological constant in the model has an appropriate sign. Both models are postulated, rather than derived from first principles, and both ignore the flame structure, which depends on chemical kinetics and transport processes within the flame. At present, there are two models which have been derived, rather than postulated, and which are valid in two non-overlapping regions of parameter space. Sivashinsky derived a generalization of the Darrieus-Landau model which is valid for Lewis numbers (ratio of thermal diffusivity to mass diffusivity of the deficient reaction component) bounded away from unity. Matalon & Matkowsky derived a model valid for Lewis numbers close to unity. Each model has its own advantages and disadvantages. Under appropriate conditions the Matalon-Matkowsky model exhibits a high-wavenumber cutoff of the Darrieus-Landau instability. However, since the Lewis numbers considered lie too close to unity, the Matalon-Matkowsky model does not capture the pulsating instability. The Sivashinsky model does capture the pulsating instability, but does not exhibit its high-wavenumber cutoff. In this paper, we derive a model consisting of a new flame speed relation and new jump conditions, which is valid for arbitrary Lewis numbers. It captures the pulsating instability and exhibits the high-wavenumber cutoff of all instabilities. The flame speed relation includes the effect of short wavelengths, not previously considered, which leads to stabilizing transverse surface diffusion terms.

Dynamical properties of a strongly correlated model for quarter-filled layered organic molecular crystals

The dynamical properties of an extended Hubbard model, which is relevant to quarter-filled layered organic molecular crystals, are analyzed. We have computed the dynamical charge correlation function, spectral density, and optical conductivity using Lanczos diagonalization and large-N techniques. As the ratio of the nearest-neighbor Coulomb repulsion, V, to the hopping integral, t, increases there is a transition from a metallic phase to a charge-ordered phase. Dynamical properties close to the ordering transition are found to differ from the ones expected in a conventional metal. Large-N calculations display an enhancement of spectral weight at low frequencies as the system is driven closer to the charge-ordering transition in agreement with Lanczos calculations. As V is increased the charge correlation function displays a collective mode which, for wave vectors close to (pi,pi), increases in amplitude and softens as the charge-ordering transition is approached. We propose that inelastic x-ray scattering be used to detect this mode. Large-N calculations predict superconductivity with d(xy) symmetry close to the ordering transition. We find that this is consistent with Lanczos diagonalization calculations, on lattices of 20 sites, which find that the binding energy of two holes becomes negative close to the charge-ordering transition.

Optimal debt, endogenous fertility, and human capital externalities in a model with altruistic bequests

and human capital externalities. Because of such externalities, education investment is too low and fertility is too high. While education subsidies are the conventional means to deal with these problems, we show that the optimal policy also comprises debt even when distortionary taxes are used. The reason is that debt tips the usual trade-off between children's quantity and quality in favor of the latter by increasing the bequest cost of children. The optimal debt-output ratio exceeds 10% for plausible parameterization. (C) 2002 Elsevier B.V. All rights reserved.

On a resampling approach for tests on the number of clusters with mixture model-based clustering of tissue samples

We consider the problem of assessing the number of clusters in a limited number of tissue samples containing gene expressions for possibly several thousands of genes. It is proposed to use a normal mixture model-based approach to the clustering of the tissue samples. One advantage of this approach is that the question on the number of clusters in the data can be formulated in terms of a test on the smallest number of components in the mixture model compatible with the data. This test can be carried out on the basis of the likelihood ratio test statistic, using resampling to assess its null distribution. The effectiveness of this approach is demonstrated on simulated data and on some microarray datasets, as considered previously in the bioinformatics literature. (C) 2004 Elsevier Inc. All rights reserved.

Numerical simulation of abrupt contraction flows using the Double Convected Pom-Pom model

The Double Convected Pom-Pom model was recently introduced to circumvent some numerical and theological defects found in other formulations of the Pom-Pom concept. It is used here for the simulation of a benchmark problem: the flow in an abrupt planar contraction. The predictions are compared with birefringence measurements and show reasonable quantitative agreement with experimental data. A parametric study is also carried out with the aim of analysing the effect of the branching parameter on vortex dynamics and extrudate swell. The results show that the Double Convected Pom-Pom model (DCPP) model is able to discriminate between branched and linear macromolecular structures in accordance with experimental observations. In that respect, the role of the extensional properties in determining complex flow behaviour is stressed. Also, the ratio of the first normal stress difference to the shear stress appears to play a major role in die swell observation. For the time being, the role of the second normal stress difference appears to be less obvious to evaluate in this complex flow. (C) 2004 Elsevier B.V. All rights reserved.

Spatio-temporal scanning and statistical test of the accelerating moment release (AMR) model using Australian earthquake data

The Accelerating Moment Release (AMR) preceding earthquakes with magnitude above 5 in Australia that occurred during the last 20 years was analyzed to test the Critical Point Hypothesis. Twelve earthquakes in the catalog were chosen based on a criterion for the number of nearby events. Results show that seven sequences with numerous events recorded leading up to the main earthquake exhibited accelerating moment release. Two occurred near in time and space to other earthquakes preceded by AM R. The remaining three sequences had very few events in the catalog so the lack of AMR detected in the analysis may be related to catalog incompleteness. Spatio-temporal scanning of AMR parameters shows that 80% of the areas in which AMR occurred experienced large events. In areas of similar background seismicity with no large events, 10 out of 12 cases exhibit no AMR, and two others are false alarms where AMR was observed but no large event followed. The relationship between AMR and Load-Unload Response Ratio (LURR) was studied. Both methods predict similar critical region sizes, however, the critical point time using AMR is slightly earlier than the time of the critical point LURR anomaly.

A simulation model of cereal-legume intercropping systems for semi-arid regions I. Model development

Cereal-legume intercropping plays an important role in subsistence food production in developing countries, especially in situations of limited water resources. Crop simulation can be used to assess risk for intercrop productivity over time and space. In this study, a simple model for intercropping was developed for cereal and legume growth and yield, under semi-arid conditions. The model is based on radiation interception and use, and incorporates a water stress factor. Total dry matter and yield are functions of photosynthetically active radiation (PAR), the fraction of radiation intercepted and radiation use efficiency (RUE). One of two PAR sub-models was used to estimate PAR from solar radiation; either PAR is 50% of solar radiation or the ratio of PAR to solar radiation (PAR/SR) is a function of the clearness index (K-T). The fraction of radiation intercepted was calculated either based on Beer's Law with crop extinction coefficients (K) from field experiments or from previous reports. RUE was calculated as a function of available soil water to a depth of 900 mm (ASW). Either the soil water balance method or the decay curve approach was used to determine ASW. Thus, two alternatives for each of three factors, i.e., PAR/SR, K and ASW, were considered, giving eight possible models (2 methods x 3 factors). The model calibration and validation were carried out with maize-bean intercropping systems using data collected in a semi-arid region (Bloemfontein, Free State, South Africa) during seven growing seasons (1996/1997-2002/2003). The combination of PAR estimated from the clearness index, a crop extinction coefficient from the field experiment and the decay curve model gave the most reasonable and acceptable result. The intercrop model developed in this study is simple, so this modelling approach can be employed to develop other cereal-legume intercrop models for semi-arid regions. (c) 2004 Elsevier B.V. All rights reserved.

A simulation model of cereal-legume intercropping systems for semi-arid regions II. Model application

Smallholder farmers in Africa practice traditional cropping techniques such as intercropping. Intercropping is thought to offer higher productivity and resource milisation than sole cropping. In this study, risk associated with maize-bean intercropping was evaluated by quantifying long-term yield in both intercropping and sole cropping in a semi-arid region of South Africa (Bloemfontein, Free State) with reference to rainfall variability. The crop simulation model was run with different cultural practices (planting date and plant density) for 52 summer crop growing seasons (1950/1951-2001/2002). Eighty-one scenarios, consisted of three levels of initial soil water, planting date, maize population, and bean population, were simulated. From the simulation outputs, the total land equivalent ratio (LER) was greater than one. The intercrop (equivalent to sole maize) had greater energy value (EV) than sole beans, and the intercrop (equivalent to sole beans) had greater monetary value (MV) than sole maize. From these results, it can be concluded that maize-bean intercropping is advantageous for this semi-arid region. Soil water at planting was the most important factor of all scenario factors, followed by planting date. Irrigation application at planting, November/December planting and high plant density of maize for EV and beans for MV can be one of the most effective cultural practices in the study region. With regard to rainfall variability, seasonal (October-April) rainfall positively affected EV and MV, but not LER. There was more intercrop production in La Nina years than in El Nino years. Thus, better cultural practices may be selected to maximize maize-bean intercrop yields for specific seasons in the semi-arid region based on the global seasonal outlook. (c) 2004 Elsevier B.V. All rights reserved.

Predator-specific changes in the morphology and swimming performance of larval Rana lessonae

1. We investigated the morphological responses of larval Rana lessonae to the presence of two predators with substantially different prey-detection and capture techniques; larval dragonflies (Aeshna cyanea) and the Pumpkinseed Sunfish (Lepomis gibossus). 2. We also examined the functional implications of any predator-induced morphological variation on their swimming ability by assessing performance during the initial stages of a startle response. 3. We found the morphological responses of larval R. lessonae were dependent on the specific predator present. Tadpoles raised in the presence of dragonfly larvae preying upon conspecific tadpoles developed total tail heights 5.4% deeper and tail muscles 4.7% shallower than tadpoles raised in a non-predator environment, while tadpoles raised with sunfish possessed tails 2% shallower and tail muscles 2.5% higher than non-predator-exposed tadpoles. 4. Predator-induced morphological variation also significantly influenced swimming performance. Tadpoles raised with sunfish possessed swimming speeds 9.5 and 14.6% higher than non- and dragonfly predator groups, respectively. 5. Thus, the expression of these alternative predator-morphs leads to a functional trade-off in performance between the different environments.

On the optimal ratio of heavy to light chain genes for efficient recombinant antibody production by CHO cells

Monoclonal antibodies (Mab) are heterotetramers consisting of an equimolar ratio of heavy chain (HC) and light chain (LC) polypeptides. Accordingly, most recombinant Mab expression systems utilize an equimolar ratio of heavy chain (he) to light chain (lc) genes encoded on either one or two plasmids. However, there is no evidence to suggest that this gene ratio is optimal for stable or transient production of recombinant Mab. In this study we have determined the optimal ratio of hc:lc genes for production of a recombinant IgG(4) Mab, cB72.3, by Chinese hamster ovary (CHO) cells using both empirical and mathematical modeling approaches. Polyethyleneimine-mediated transient expression of cB72.3 at varying ratios of hc:lc genes encoded on separate plasmids yielded an optimal Mab titer at a hc:lc gene ratio of 3:2; a conclusion confirmed by separate mathematical modeling of the Mab folding and assembly process using transient expression data. On the basis of this information, we hypothesized that utilization of he genes at low hc:lc gene ratios is more efficient. To confirm this, cB72.3 Mab was transiently produced by CHO cells at constant he and varying lc gene dose. Under these conditions, Mab yield was increased with a concomitant increase in lc gene dose. To determine if the above findings also apply to stably transfected CHO cells producing recombinant Mab, we compared the intra- and extracellular ratios of HC and LC polypeptides for three GS-CHO cells lines transfected with a 1:1 ratio of hc:lc genes and selected for stable expression of the same recombinant Mab, cB72.3. Intra- and extracellular HC:LC polypeptide ratios ranged from 1:2 to 1:5, less than that observed on transient expression of the same Mab in parental CHO cells using the same vector. In conclusion, our data suggest that the optimal ratio of hc:lc genes used for transient and stable expression of Mab differ. In the case of the latter, we infer that optimal Mab production by stably transfected cells represents a compromise between HC abundance limiting productivity and the requirement for excess LC to render Mab folding and assembly more efficient.

Unidirectional inversion of ibuprofen in Caco-2 cells: Developing a suitable model for presystemic chiral inversion study

Intestinal chiral inversion of ibuprofen is still lacking direct evidence. In a preliminary experiment, ibuprofen was found to undergo inversion in Caco-2 cells. This investigation was thus conducted to determine the characteristics and influence of some biochemical factors on the chiral inversion of ibuprofen in Caco-2 cells. The effects of substrate concentration (2.5-40 mu g/ml), cell density (0.5-2 x 10(6) cells/ well), content of serum (0-20%), coexistence of S ibuprofen (corresponding doses), sodium azide (10mm), exogenous Coenzyme A (CoA) (0.1 - 0.4 mm),. and palmitic acid (5-25 mu m) on inversion were examined. A stereoselective HPLC method based on the Chromasil-CHI-TBB column was developed for quantitative analysis of the drug in cell culture medium. The inversion ratio (F-i) and elimination rate constant were calculated as the indexes of inversion extent. Inversion of ibuprofen in Caeo-2 cells was found to be both dose and cell density dependent, indicating saturable characteristics. Addition of serum significantly inhibited the inversion, to an extent of 2.7 fold decrease at 20% content. Preexistence of S enantiomer exerted a significant inhibitory effect (p < 0.01 for all tests). Sodium azide decreased the inversion ratio from 0.43 to 0.32 (p < 0.01). Exogenous CoA and palmitic acid significantly promoted the inversion at all tested doses (p < 0.01 for all tests). This research provided strong evidence to the capacity and capability of intestinal chiral inversion. Although long incubation times up to 120 h were required, Caco-2 cells should be a suitable model for chiral inversion research of 2-APAs considering the human-resourced and well-defined characteristics from the present study.