Redundant coordinates in harmonic force-field calculations

Redundancy relations between vibrational coordinates may be linear (as for rectilinear coordinates used in deriving a G matrix), or non-linear (as for curvilinear coordinates used in formulating model force fields). It is shown that geometrically defined internal coordinates are necessarily curvilinear. Hence it is shown that linear force constants can occur in model force field calculations involving redundant coordinates, in disagreement with the recent proposal of Gussoni and Zerbi.

The Klein-Gordon equation on the half line: a Riemann-Hilbert approach

We solve an initial-boundary problem for the Klein-Gordon equation on the half line using the Riemann-Hilbert approach to solving linear boundary value problems advocated by Fokas. The approach we present can be also used to solve more complicated boundary value problems for this equation, such as problems posed on time-dependent domains. Furthermore, it can be extended to treat integrable nonlinearisations of the Klein-Gordon equation. In this respect, we briefly discuss how our results could motivate a novel treatment of the sine-Gordon equation.

Statistical downscaling of river flows

An extensive statistical ‘downscaling’ study is done to relate large-scale climate information from a general circulation model (GCM) to local-scale river flows in SW France for 51 gauging stations ranging from nival (snow-dominated) to pluvial (rainfall-dominated) river-systems. This study helps to select the appropriate statistical method at a given spatial and temporal scale to downscale hydrology for future climate change impact assessment of hydrological resources. The four proposed statistical downscaling models use large-scale predictors (derived from climate model outputs or reanalysis data) that characterize precipitation and evaporation processes in the hydrological cycle to estimate summary flow statistics. The four statistical models used are generalized linear (GLM) and additive (GAM) models, aggregated boosted trees (ABT) and multi-layer perceptron neural networks (ANN). These four models were each applied at two different spatial scales, namely at that of a single flow-gauging station (local downscaling) and that of a group of flow-gauging stations having the same hydrological behaviour (regional downscaling). For each statistical model and each spatial resolution, three temporal resolutions were considered, namely the daily mean flows, the summary statistics of fortnightly flows and a daily ‘integrated approach’. The results show that flow sensitivity to atmospheric factors is significantly different between nival and pluvial hydrological systems which are mainly influenced, respectively, by shortwave solar radiations and atmospheric temperature. The non-linear models (i.e. GAM, ABT and ANN) performed better than the linear GLM when simulating fortnightly flow percentiles. The aggregated boosted trees method showed higher and less variable R2 values to downscale the hydrological variability in both nival and pluvial regimes. Based on GCM cnrm-cm3 and scenarios A2 and A1B, future relative changes of fortnightly median flows were projected based on the regional downscaling approach. The results suggest a global decrease of flow in both pluvial and nival regimes, especially in spring, summer and autumn, whatever the considered scenario. The discussion considers the performance of each statistical method for downscaling flow at different spatial and temporal scales as well as the relationship between atmospheric processes and flow variability.

Estimation of the stoichiometry of volatile fatty acid production in the rumen of lactating cows

The purpose of this study was to improve the prediction of the quantity and type of Volatile Fatty Acids (VFA) produced from fermented substrate in the rumen of lactating cows. A model was formulated that describes the conversion of substrate (soluble carbohydrates, starch, hemi-cellulose, cellulose, and protein) into VFA (acetate, propionate, butyrate, and other VFA). Inputs to the model were observed rates of true rumen digestion of substrates, whereas outputs were observed molar proportions of VFA in rumen fluid. A literature survey generated data of 182 diets (96 roughage and 86 concentrate diets). Coefficient values that define the conversion of a specific substrate into VFA were estimated meta-analytically by regression of the model against observed VFA molar proportions using non-linear regression techniques. Coefficient estimates significantly differed for acetate and propionate production in particular, between different types of substrate and between roughage and concentrate diets. Deviations of fitted from observed VFA molar proportions could be attributed to random error for 100%. In addition to regression against observed data, simulation studies were performed to investigate the potential of the estimation method. Fitted coefficient estimates from simulated data sets appeared accurate, as well as fitted rates of VFA production, although the model accounted for only a small fraction (maximally 45%) of the variation in VFA molar proportions. The simulation results showed that the latter result was merely a consequence of the statistical analysis chosen and should not be interpreted as an indication of inaccuracy of coefficient estimates. Deviations between fitted and observed values corresponded to those obtained in simulations. (c) 2005 Elsevier Ltd. All rights reserved.

Modelling effects of prolonged conditioning on dormancy and germination of Striga hermonthica

The impact of environment on the germination biology of the parasite was studied in the laboratory with seeds conditioned at various water potentials, urea concentrations and at 17.5 to 37.5°C for up to 133 days. Maximum germination was observed at 20 to 25°C. Water stress and urea suppressed maximum germination. The final percentage germination response to period of conditioning showed a non-linear relationship and suggests the release of seeds from dormancy during the initial period and later on dormancy induction. Germination percentage increased with increase in conditioning period to a threshold and remained stable for variable periods followed by a decline with further extension of conditioning time. The decline in germination finally terminated in zero germination in most treatments before the end of experimentation. The investigated factors of temperature, water potential and urea showed clear effects on the expression of dormancy pattern of the parasite. The effects of water potential and urea were viewed as modifying a primary response of seeds to temperature during conditioning. The changes in germinability potential during conditioning were consistent with the hypothesis that dormancy periods are normally distributed within seed populations and that loss of primary dormancy precedes induction of secondary dormancy. Hence an additive mathematical model of loss of primary dormancy and induction of secondary as affected by environment was developed as: G = {[Φ-1 (Kp+ (po+pnN+pwW) (T-Tb) t)]-[Φ-1 (Ks+ ((swW+sa)+sorT)t)]}[Φ-1(aT2+bT+c+cwW)].

Describing long-term trends in precipitation using generalized additive models

With the current concern over climate change, descriptions of how rainfall patterns are changing over time can be useful. Observations of daily rainfall data over the last few decades provide information on these trends. Generalized linear models are typically used to model patterns in the occurrence and intensity of rainfall. These models describe rainfall patterns for an average year but are more limited when describing long-term trends, particularly when these are potentially non-linear. Generalized additive models (GAMS) provide a framework for modelling non-linear relationships by fitting smooth functions to the data. This paper describes how GAMS can extend the flexibility of models to describe seasonal patterns and long-term trends in the occurrence and intensity of daily rainfall using data from Mauritius from 1962 to 2001. Smoothed estimates from the models provide useful graphical descriptions of changing rainfall patterns over the last 40 years at this location. GAMS are particularly helpful when exploring non-linear relationships in the data. Care is needed to ensure the choice of smooth functions is appropriate for the data and modelling objectives. (c) 2008 Elsevier B.V. All rights reserved.

When are ant-attractant devices a worthwhile investment? Vicia faba extrafloral nectaries and Lasius niger ants

Most studies aiming to determine the beneficial effect of ants on plants simply consider the effects of the presence or exclusion of ants on plant yield. This approach is often inadequate, however, as ants interact with both non-tended herbivores and tended Homoptera. Moreover, the interaction with these groups of organisms is dependent on ant density, and these functional relationships are likely to be non-linear. A model is presented here that segregates plant herbivores into two categories depending on the sign of their numerical response to ants (myrmecophiles increase with ants, non-tended herbivores decline). The changes in these two components of herbivores with increasing ant density and the resulting implications for ant-plant mutualisms are considered. It emerges that a wide range of ant densities needs to be considered as the interaction sign (mutualism or parasitism) and strength is likely to change with ant density. The model is used to interpret the results of an experimental study that varied levels of Aphis fabae infestation and Lasius niger ant attendance on Vicia faba bean plants. Increasing ant density consistently reduced plant fitness and thus, in this location, the interaction between the ants and the plant can be considered parasitic. In the Vicia faba system, these costs of ants are unlikely to be offset by other beneficial agents (e.g., parasitoids), which also visit extrafloral nectaries.

GeneMCL in microarray analysis

Accurately and reliably identifying the actual number of clusters present with a dataset of gene expression profiles, when no additional information on cluster structure is available, is a problem addressed by few algorithms. GeneMCL transforms microarray analysis data into a graph consisting of nodes connected by edges, where the nodes represent genes, and the edges represent the similarity in expression of those genes, as given by a proximity measurement. This measurement is taken to be the Pearson correlation coefficient combined with a local non-linear rescaling step. The resulting graph is input to the Markov Cluster (MCL) algorithm, which is an elegant, deterministic, non-specific and scalable method, which models stochastic flow through the graph. The algorithm is inherently affected by any cluster structure present, and rapidly decomposes a graph into cohesive clusters. The potential of the GeneMCL algorithm is demonstrated with a 5730 gene subset (IGS) of the Van't Veer breast cancer database, for which the clusterings are shown to reflect underlying biological mechanisms. (c) 2005 Elsevier Ltd. All rights reserved.

Novel charge transfer supramolecular assemblies with Keggin anions and 2-amino-5-nitropyridine

Several novel compounds with the non-linear optical chromophore 2-amino-5-nitropyridine (2A5NP) and Keggin polyoxoanions (alpha-isomers), having the general formula (2A5NP)(m)H-n[XM12O40]center dot xH(2)O, M = Mo, W, were synthesised. Compounds were obtained with X = P, n = 3, m = 3 and 4 and X = Si, n = m = 4 ( x = 2 - 6). Thus, for each of the anions [PMo12O40](3-) and [PW12O40](3-) two different compounds were obtained, with the same anion and organic counterpart but with a different stoichiometric ratio. These presented different charge transfer properties and thermal stability. All compounds were characterised by spectroscopic and analytical techniques. The single crystal X-ray diffraction structure of (2A5NP)(4)H-3[PMo12O40]center dot 2.5H(2)O center dot 0.5C(2)H(5)OH showed that the water solvent molecules and the organic chromophores are assembled via infinite one-dimensional chains of hydrogen bonds with formation of open channels, which accommodate [ PMo12O40] 3- and ethanol solvent molecules.

The kinetics of the gas-phase decomposition of the 2-methyl-2-butoxyl and 2-methyl-2-pentoxyl radicals

The kinetics of the title reactions have been studied by relative-rate methods as a function of temperature. Relative-rate coefficients for the two decomposition channels of 2-methyl-2-butoxyl have been measured at five different temperatures between 283 and 345 K and the observed temperature dependence is consistent with the results of some previous experimental studies. The kinetics of the two decomposition channels of 2-methyl-2-pentoxyl have also been investigated, as a function of temperature, relative to the estimated rate of isomerisation of this radical. Room-temperature rate coefficient data for the two decomposition channels of both 2-methyl-2-pentoxyl and 2-methyl-2-butxoyl (after combining the relative rate coefficient for this latter with a value for the rate coefficient of the major channel, extrapolated from the data presented by Batt et al., Int. J. Chem. Kinet., 1978, 10, 931) are shown to be consistent with a non-linear kinetic correlation, for alkoxyl radical decomposition rate data, previously presented by this laboratory (Johnson et al., Atmos. Environ., 2004, 38, 1755-1765).

Torque in mooring chain. Part 1: Background and theory

Chain in both its forms - common (or stud-less) and stud-link - has many engineering applications. It is widely used as a component in the moorings of offshore floating systems, where its ruggedness and corrosion resistance make it an attractive choice. Chain exhibits some interesting behaviour in that when straight and subject to an axial load it does not twist or generate any torque, but if twisted or loaded when in a twisted condition it behaves in a highly non-linear manner, with the torque dependent upon the level of twist and axial load. Clearly an understanding of the way in which chains may behave and interact with other mooring components (such as wire rope, which also exhibits coupling between axial load and generated torque) when they are in service is essential. However, the sizes of chain that are in use in offshore moorings (typical bar diameters are 75 mm and greater) are too large to allow easy testing. This paper, which is in two parts, aims to address the issues and considerations relevant to torque in mooring chain. The first part introduces a frictionless theory that predicts the resultant torques and 'lift' in the links as non-dimensionalized functions of the angle of twist. Fortran code is presented in an Appendix, which allows the reader to make use of the analysis. The second part of the paper presents results from experimental work on both stud-less (41 mm) and stud-link (20.5 and 56 mm) chains. Torsional data are presented in both 'constant twist' and 'constant load' forms, as well as considering the lift between the links.

Torque in mooring chain. Part 2: Experimental investigation

Chain is a commonly used component in offshore moorings where its ruggedness and corrosion resistance make it an attractive choice. Another attractive property is that a straight chain is inherently torque balanced. Having said this, if a chain is loaded in a twisted condition, or twisted when under load, it exhibits highly non-linear torsional behaviour. The consequences of this behaviour can cause handling difficulties or may compromise the integrity of the mooring system, and care must be taken to avoid problems for both the chain and any components to which it is connected. Even with knowledge of the potential problems, there will always be occasions where, despite the utmost care, twist is unavoidable. Thus it is important for the engineer to be able to determine the effects. A frictionless theory has been developed in Part 1 of the paper that may be used to predict the resultant torques and movement or 'lift' in the links as non-dimensional functions of the angle of twist. The present part of the paper describes a series of experiments undertaken on both studless and stud-link chain to allow comparison of this theoretical model with experimental data. Results are presented for the torsional response and link lift for 'constant twist' and 'constant load' type tests on chains of three different link sizes.

Extraction and determination of ellagic acid content in chestnut bark and fruit

Chestnuts are an important economic resource in the chestnut growing regions, not only for the fruit, but also for the wood. The content of ellagic acid (EA), a naturally occurring inhibitor of carcinogenesis, was determined in chestnut fruits and bark. EA was extracted with methanol and free ellagic acid was determined by HPLC with UV detection, both in the crude extract and after hydrolysis. The concentration of EA was generally increased after hydrolysis due to the presence of ellagitannins in the crude extract. The concentration varied between 0.71 and 21.6 ing g(-1) (d.w.) in un-hydrolyzed samples, and between 2.83 and 18.4 mg g(-1) (d.w.) ill hydrolyzed samples. In chestnut fruits, traces of EA were present in the seed, with higher concentrations in the pellicle and pericarp. However, all fruit tissues had lower concentrations of EA than had the bark. The concentration of EA in the hydrolyzed samples showed a non-linear correlation with the concentration in the unhydrolyzed extracts. (C) 2008 Elsevier Ltd. All rights reserved.

Effect of minerals on casein micelle stability of cows' milk

The effects of minerals on casein micelle stability of individual cows' milk, throughout a complete lactation, were investigated. Calcium and calcium ions, magnesium, phosphorus, sodium, potassium and citrate contents were analysed, together with the following physical properties of milk; pH, ethanol stability, rennet clotting time and coagulum firmness. There was an inverse non-linear relationship between free calcium ion concentration and ethanol stability (ES; r=0.84). Rennet coagulation time showed a weaker relationship with free calcium ion concentration (r=0.44) but a stronger relationship with pH (r=0.66). In addition, samples containing higher amounts of free calcium ions produced a firmer gel. Citrate in natural samples acts as a stabilizing factor, as it slightly improves milk stability. Potassium, on the other hand, exhibited a negative correlation, but only with rennet clotting time (r=-0.52). Throughout lactation the average values were; free Ca 21 concentration 1.88 mm,, pH 6.63, ES 83.2% and clotting time 13.6 min. The equilibrium relationship between pH and free Ca2+ concentration was investigated by adjusting milk pH from 5.9 to 7.1, using acid and alkali. There was a good inverse linear relationship between pH and log (free Ca 21) for individual milk samples, with a gradient of -0.62 and a standard deviation of 0.042.

State-space model identification of a wireless THz network

We model the large scale fading of wireless THz communications links deployed in a metropolitan area taking into account reception through direct line of sight, ground or wall reflection and diffraction. The movement of the receiver in the three dimensions is modelled by an autonomous dynamic linear system in state-space whereas the geometric relations involved in the attenuation and multi-path propagation of the electric field are described by a static non-linear mapping. A subspace algorithm in conjunction with polynomial regression is used to identify a Wiener model from time-domain measurements of the field intensity.