Examining random number generators used in stochastic iteration algorithms

This paper investigates random number generators in stochastic iteration algorithms that require infinite uniform sequences. We take a simple model of the general transport equation and solve it with the application of a linear congruential generator, the Mersenne twister, the mother-of-all generators, and a true random number generator based on quantum effects. With this simple model we show that for reasonably contractive operators the theoretically not infinite-uniform sequences perform also well. Finally, we demonstrate the power of stochastic iteration for the solution of the light transport problem.

A matrix iteration for dynamic network summaries

We propose a new algorithm for summarizing properties of large-scale time-evolving networks. This type of data, recording connections that come and go over time, is being generated in many modern applications, including telecommunications and on-line human social behavior. The algorithm computes a dynamic measure of how well pairs of nodes can communicate by taking account of routes through the network that respect the arrow of time. We take the conventional approach of downweighting for length (messages become corrupted as they are passed along) and add the novel feature of downweighting for age (messages go out of date). This allows us to generalize widely used Katz-style centrality measures that have proved popular in network science to the case of dynamic networks sampled at non-uniform points in time. We illustrate the new approach on synthetic and real data.

Preferred orientation in an angled intercalation site of a chloro-substituted lambda- [Ru(TAP)2(dppz)]2+ complex bound to d(TCGGCGCCGA)2

The crystal structure of the ruthenium DNA ‘light-switch’ complex -[Ru(TAP)2(11-Cl-dppz)]2+ (TAP = tetraazaphenanthrene, dppz = dipyrido[3,2-a':2',3'-c]phenazine)) bound to the oligonucleotide duplex d(TCGGCGCCGA)2 is reported. The synthesis of the racemic ruthenium complex is described for the first time, and the racemate was used in this study. The crystal structure, at atomic resolution (1.0 Å), shows one ligand as a wedge in the minor groove, resulting in the 51 kinking of the double helix, as with the parent lambda-[Ru(TAP)2(dppz)]2+. Each complex binds to one duplex by intercalation of the dppz ligand and also by semi-intercalation of one of the orthogonal TAP ligands into a second symmetrically equivalent duplex. The 11-Cl substituent binds with the major component (66%) oriented with the 11-chloro substituent on the purine side of the terminal step of the duplex.

Dynamic load balancing in parallel KD-tree k-means

One among the most influential and popular data mining methods is the k-Means algorithm for cluster analysis. Techniques for improving the efficiency of k-Means have been largely explored in two main directions. The amount of computation can be significantly reduced by adopting geometrical constraints and an efficient data structure, notably a multidimensional binary search tree (KD-Tree). These techniques allow to reduce the number of distance computations the algorithm performs at each iteration. A second direction is parallel processing, where data and computation loads are distributed over many processing nodes. However, little work has been done to provide a parallel formulation of the efficient sequential techniques based on KD-Trees. Such approaches are expected to have an irregular distribution of computation load and can suffer from load imbalance. This issue has so far limited the adoption of these efficient k-Means variants in parallel computing environments. In this work, we provide a parallel formulation of the KD-Tree based k-Means algorithm for distributed memory systems and address its load balancing issue. Three solutions have been developed and tested. Two approaches are based on a static partitioning of the data set and a third solution incorporates a dynamic load balancing policy.

Quantification of physical and biological uncertainty in the simulation of the yield of a tropical crop using present-day and doubled CO2 climates

The impacts of climate change on crop productivity are often assessed using simulations from a numerical climate model as an input to a crop simulation model. The precision of these predictions reflects the uncertainty in both models. We examined how uncertainty in a climate (HadAM3) and crop General Large-Area Model (GLAM) for annual crops model affects the mean and standard deviation of crop yield simulations in present and doubled carbon dioxide (CO2) climates by perturbation of parameters in each model. The climate sensitivity parameter (lambda, the equilibrium response of global mean surface temperature to doubled CO2) was used to define the control climate. Observed 1966-1989 mean yields of groundnut (Arachis hypogaea L.) in India were simulated well by the crop model using the control climate and climates with values of lambda near the control value. The simulations were used to measure the contribution to uncertainty of key crop and climate model parameters. The standard deviation of yield was more affected by perturbation of climate parameters than crop model parameters in both the present-day and doubled CO2 climates. Climate uncertainty was higher in the doubled CO2 climate than in the present-day climate. Crop transpiration efficiency was key to crop model uncertainty in both present-day and doubled CO2 climates. The response of crop development to mean temperature contributed little uncertainty in the present-day simulations but was among the largest contributors under doubled CO2. The ensemble methods used here to quantify physical and biological uncertainty offer a method to improve model estimates of the impacts of climate change.

Percentile growth charts for biomedical studies using a porcine model

Increasing rates of obesity and heart disease are compromising quality of life for a growing number of people. There is much research linking adult disease with the growth and development both in utero and during the first year of life. The pig is an ideal model for studying the origins of developmental programming. The objective of this paper was to construct percentile growth curves for the pig for use in biomedical studies. The body weight (BIN) of pigs was recorded from birth to 150 days of age and their crown-to-rump length was measured over the neonatal period to enable the ponderal index (Pl; kg/m(3)) to be calculated. Data were normalised and percentile curves were constructed using Cole's lambda-mu-sigma (LMS) method for BW and PI. The construction of these percentile charts for use in biomedical research will allow a more detailed and precise tracking of growth and development of individual pigs under experimental conditions.

Nutrition of chicks and layers

When formulating least-cost poultry diets, ME concentration should be optimised by an iterative procedure, not entered as a fixed value. This iteration must calculate profit margins by taking into account the way in which feed intake and saleable outputs vary with ME concentration. In the case of broilers, adjustment of critical amino acid contents in direct proportion to ME concentration does not result in birds of equal fatness. To avoid an increase in fat deposition at higher energy levels, it is proposed that amino acid specifications should be adjusted in proportion to changes in the net energy supplied by the feed. A model is available which will both interpret responses to amino acids in laying trials and give economically optimal estimates of amino acid inputs for practical feed formulation. Flocks coming into lay and flocks nearing the end of the pullet year have bimodal distributions of rates of lay, with the result that calculations of requirement based on mean output will underestimate the optimal amino acid input for the flock. Chick diets containing surplus protein can lead to impaired utilisation of the first-limiting amino acid. This difficulty can be avoided by stating amino acid requirements as a proportion of the protein.

Clonal variation in acetylcholinesterase biomarkers and life history traits following OP exposure in Daphnia magna

Two clones of Daphnia magna (Standard and Ruth) were exposed for 7 days to sub-lethal concentrations of acephate (5.0 and 10.0 mg/L). Survivorship, individual growth, reproduction and the population growth rate (lambda) were evaluated over three weeks. Acetylcholinesterase (AChE) activity was measured on days 2, 7 and 21. Acephate exposure inhibited AChE activity but had no direct effect on life history (LH) traits. There was also no effect of clone on AChE activity, LH and lambda. However, a significant interaction between clone and acephate concentration was found on both fecundity and AChE inhibition at 48 h was associated with a decrease in lambda the Standard clone and an increase in lambda in clone Ruth. Therefore, our findings show that genotypic variation will influence the link between AChE activity and toxic effects at higher levels of biological organisation in D. magna. (c) 2007 Elsevier Inc. All rights reserved.

Synthesis and characterization of fluorescent poly(aromatic amide) dendrimers

The synthesis of a series of poly(aromatic amide) dendrimers up to the second generation is described herein. The AB, building block used throughout the synthesis of the dendrimers was the allyl ester of 3,5-diaminocinnamic acid, which has been synthesized from 3,5-dinitrobenzoic acid in good yield with use of a four-step procedure. Dendron synthesis was achieved via a convergent approach with use of a sequence of deprotection/coupling steps. Two commercially available alcohols, L-menthol and citronellol, were coupled to the AB(2) monomer by using an alkyl diacid spacer and two core units; 1,7-diaminoheptane and tris(2-aminoethyl)amine have been used to produce the final dendrimers. Characterization was carried out by NMR and IR spectroscopies, MALDI-TOF mass spectrometry, GPC, and DSC. The novel monomer and dendritic derivatives exhibited a strong fluorescence emission in the visible region (lambda approximate to 500 nm) of the spectrum and a weak emission in the near-infrared (lambda approximate to 850 nm) upon excitation in the near-UV region. The fluorescence emission characteristics were found to be solvent and dendrimer generation dependent.

Direct detection of dimethylstannylene and tetramethyldistannene in solution and the gas phase by laser flash photolysis of 1,1-dimethylstannacyclopent-3-enes

The photochemistry of 1,1-dimethyl- and 1,1,3,4-tetramethylstannacyclopent-3-ene (4a and 4b,respectively) has been studied in the gas phase and in hexane solution by steady-state and 193-nm laser flash photolysis methods. Photolysis of the two compounds results in the formation of 1,3-butadiene (from 4a) and 2,3-dimethyl-1,3-butadiene (from 4b) as the major products, suggesting that cycloreversion to yield dimethylstannylene (SnMe2) is the main photodecomposition pathway of these molecules. Indeed, the stannylene has been trapped as the Sn-H insertion product upon photolysis of 4a in hexane containing trimethylstannane. Flash photolysis of 4a in the gas phase affords a transient absorbing in the 450-520nm range that is assigned to SnMe2 by comparison of its spectrum and reactivity to those previously reported from other precursors. Flash photolysis of 4b in hexane solution affords results consistent with the initial formation of SnMe2 (lambda(max) approximate to 500 nm), which decays over similar to 10 mu s to form tetramethyldistannene (5b; lambda(max) approximate to 470 nm). The distannene decays over the next ca. 50 mu s to form at least two other longer-lived species, which are assigned to higher SnMe2 oligomers. Time-dependent DFT calculations support the spectral assignments for SnMe2 and Sn2Me4, and calculations examining the variation in bond dissociation energy with substituent (H, Me, and Ph) in disilenes, digermenes, and distannenes rule out the possibility that dimerization of SnMe2 proceeds reversibly. Addition of methanol leads to reversible reaction with SnMe2 to form a transient absorbing at lambda(max) approximate to 360 nm, which is assigned to the Lewis acid-base complex between SnMe2 and the alcohol.

A variational method for the calculation of spin-rovibronic energy levels of any triatomic molecule in an electronic triplet state

The triatomic spin-rovibronic variational code RVIB3 has been extended to include the effect of two uncoupled electrons, for both (3)Sigma(-) and (3)Pi (Renner-Teller) electronic states. The spin-orbital-rotational kinetic energy is included in the usual way, via terms (J+L+S). The phenomenological terms AL.S and lambda 2/3(3S(z)(2)) are introduced to reproduce the 3 spin-orbit and spin-spin splittings, respectively. Calculations are performed to evaluate the spin-rovibronic energy levels of CCO (X) over tilde (3) Sigma(-) and CCO (A) over tilde (3) Pi for which the Born-Oppenheimer potentials are derived from high-accuracy ab initio calculations.

Synthesis, structure, solution chemistry and the electronic effect of para substituents on the vanadium center in a family of mixed-ligand [(VO)-O-V(ONO)(ON)] compexes

Four tridentate dibasic ONO donor hydrazone ligands derived from the condensation of benzoylhydrazine with either 2-hydroxyacetophenone or its para substituted derivatives (H2L1-4, general abbreviation H2L) have been used as primary ligands and 8-hydroxyquinoline (Hhq, a bidentate monobasic ON donor species) has been used as auxiliary ligand. The reaction of [(VO)-O-IV(acac)21 with H2L in methanol followed by the addition of Hhq in equimolar ratio under aerobic condition afforded the mixed-ligand oxovanadium(V) complexes of the type [(VO)-O-V(L)(hq)] (1-4) in excellent yield. The X-ray structure of the compound [(VO)-O-V(L-4)(hq)] (4) indicates that the H2L4 ligand is bonded with vanadium meridionally in a tridentate dinegative fashion through its deprotonated phenolic-O, deprotonated enolic-O and imine-N atoms. The V-O bond length order is: oxo < phenolato < enolato. H-1 NMR spectra of 4 in CDCl3 solution indicates that it's solid-state structure is retained in solution. Complexes are diamagnetic and exhibit only ligand to metal charge transfer (LMCT) transition band near 530 nm in CH2Cl2 solution in addition to intra-ligand pi-pi* transition band near 335 rim and they display quasi-reversible one electron reduction peak near -0.10 V versus SCE in CH2Cl2 solution. lambda(max) (for LMCT transition) and the reduction peak potential (E-p(c)) values of the complexes are found to be linearly related with the Hammett (sigma) constants of the substituents in the aryloxy ring of the hydrazone ligands. lambda(max) and E-p(c) values show large dependence d lambda(max)/d sigma = 32.54 nm and dE(p)(c)/d sigma = 0.19 V, respectively, on the Hammett constant. (c) 2006 Elsevier B.V. All rights reserved.

Higher oxides of chlorine: absorption cross-sections of Cl2O6 and Cl2O4, the decomposition of Cl2O6, and the reactions of OClO with O and O-3

The absorption cross-sections of Cl2O6 and Cl2O4 have been obtained using a fast flow reactor with a diode array spectrometer (DAS) detection system. The absorption cross-sections at the wavelengths of maximum absorption (lambda(max)) determined in this study are those of Cl2O6: (1.47 +/- 0.15) x 10(-17) cm(2) molecule(-1), at lambda(max) = 276 nm and T = 298 K; and Cl2O4: (9.0 +/- 2.0) x 10(-19) cm(2) molecule(-1), at lambda(max) = 234 nm and T = 298 K. Errors quoted are two standard deviations together with estimates of the systematic error. The shapes of the absorption spectra were obtained over the wavelength range 200-450 nm for Cl2O6 and 200-350 nm for Cl2O4, and were normalized to the absolute cross-sections obtained at lambda(max) for each oxide, and are presented at 1 nm intervals. These data are discussed in relation to previous measurements. The reaction of O with OCIO has been investigated with the objective of observing transient spectroscopic absorptions. A transient absorption was seen, and the possibility is explored of identifying the species with the elusive sym-ClO3 or ClO4, both of which have been characterized in matrices, but not in the gas-phase. The photolysis of OCIO was also re-examined, with emphasis being placed on the products of reaction. UV absorptions attributable to one of the isomers of the ClO dimer, chloryl chloride (ClClO2) were observed; some Cl2O4 was also found at long photolysis times, when much of the ClClO2 had itself been photolysed. We suggest that reports of Cl2O6 formation in previous studies could be a consequence of a mistaken identification. At low temperatures, the photolysis of OCIO leads to the formation of Cl2O3 as a result of the addition of the ClO primary product to OCIO. ClClO2 also appears to be one product of the reaction between O-3 and OCIO, especially when the reaction occurs under explosive conditions. We studied the kinetics of the non-explosive process using a stopped-flow technique, and suggest a value for the room-temperature rate coefficient of (4.6 +/- 0.9) x 10(-19) cm(3) molecule(-1) s(-1) (limit quoted is 2sigma random errors). The photochemical and thermal decomposition of Cl2O6 is described in this paper. For photolysis at k = 254 nm, the removal of Cl2O6 is not accompanied by the build up of any other strong absorber. The implications of the results are either that the photolysis of Cl2O6 produces Cl-2 directly, or that the initial photofragments are converted rapidly to Cl-2. In the thermal decomposition of Cl2O6, Cl2O4 was shown to be a product of reaction, although not necessarily the major one. The kinetics of decomposition were investigated using the stopped-flow technique. At relatively high [OCIO] present in the system, the decay kinetics obeyed a first-order law, with a limiting first-order rate coefficient of 0.002 s(-1). (C) 2004 Elsevier B.V. All rights reserved.