Experimental identification of chaotic fibers

In a 2D parameter space, by using nine experimental time series of a Clitia`s circuit, we characterized three codimension-1 chaotic fibers parallel to a period-3 window. To show the local preservation of the properties of the chaotic attractors in each fiber, we applied the closed return technique and two distinct topological methods. With the first topological method we calculated the linking, numbers in the sets of unstable periodic orbits, and with the second one we obtained the symbolic planes and the topological entropies by applying symbolic dynamic analysis. (C) 2007 Elsevier Ltd. All rights reserved.

Classical simulation of deposition of thiophene oligomers on TiO(2)-anatase: Relevance of long-range electrostatic interactions

We performed classical molecular dynamics simulations of the vapor-deposition of alpha-T4 oligomers on the TiO(2)-anatase (101) surface, comparing different sets of charges associated with the atoms of the model. The potential energy surfaces for alpha-T4 and TiO(2) were described by re-parametrizations of the Universal force field with charges given by the charge equilibration (QEq) scheme, or with fixed charges obtained by an ab initio method using the Hirshfeld partition. The two sets of charges lead to completely different results for the interface formation, and for the characteristics of the organic film, with a clearly defined alpha-T4 contact layer in the QEq case, and a more homogeneous molecular distribution when using Hirshfeld charges. The main reason for the discrepancy was found to be the incorrect charge assignment given by QEq to the sulfur and alpha-carbon atoms in thiophenes, and highlight the relevance of long-range interactions in the organization of molecular films. (C) 2009 Elsevier B.V. All rights reserved.

Collapse of rho (xx) Ringlike Structures in 2DEGs Under Tilted Magnetic Fields

In the quantum Hall regime, the longitudinal resistivity rho (xx) plotted as a density-magnetic-field (n (2D) -B) diagram displays ringlike structures due to the crossings of two sets of spin split Landau levels from different subbands [see, e.g., Zhang et al., in Phys. Rev. Lett. 95:216801, 2005. For tilted magnetic fields, some of these ringlike structures ""shrink"" as the tilt angle is increased and fully collapse at theta (c) a parts per thousand 6A degrees. Here we theoretically investigate the topology of these structures via a non-interacting model for the 2DEG. We account for the inter Landau-level coupling induced by the tilted magnetic field via perturbation theory. This coupling results in anticrossings of Landau levels with parallel spins. With the new energy spectrum, we calculate the corresponding n (2D) -B diagram of the density of states (DOS) near the Fermi level. We argue that the DOS displays the same topology as rho (xx) in the n (2D) -B diagram. For the ring with filling factor nu=4, we find that the anticrossings make it shrink for increasing tilt angles and collapse at a large enough angle. Using effective parameters to fit the theta=0A degrees data, we find a collapsing angle theta (c) a parts per thousand 3.6A degrees. Despite this factor-of-two discrepancy with the experimental data, our model captures the essential mechanism underlying the ring collapse.

Cross-situational learning of object-word mapping using Neural Modeling Fields

The issue of how children learn the meaning of words is fundamental to developmental psychology. The recent attempts to develop or evolve efficient communication protocols among interacting robots or Virtual agents have brought that issue to a central place in more applied research fields, such as computational linguistics and neural networks, as well. An attractive approach to learning an object-word mapping is the so-called cross-situational learning. This learning scenario is based on the intuitive notion that a learner can determine the meaning of a word by finding something in common across all observed uses of that word. Here we show how the deterministic Neural Modeling Fields (NMF) categorization mechanism can be used by the learner as an efficient algorithm to infer the correct object-word mapping. To achieve that we first reduce the original on-line learning problem to a batch learning problem where the inputs to the NMF mechanism are all possible object-word associations that Could be inferred from the cross-situational learning scenario. Since many of those associations are incorrect, they are considered as clutter or noise and discarded automatically by a clutter detector model included in our NMF implementation. With these two key ingredients - batch learning and clutter detection - the NMF mechanism was capable to infer perfectly the correct object-word mapping. (C) 2009 Elsevier Ltd. All rights reserved.

Spectroscopic investigations of OH- influence on near-infrared fluorescence quenching of Yb3+/Tm3+ co-doped sodium-metaphosphate glasses

Energy transfer processes were studied in two sets of Yb3+ and Tm3+ co-doped sodium-metaphosphate glasses, prepared in air and nitrogen atmospheres. Using Forster, Dexter, and Miyakawa theoretical models, the energy transfer parameters were calculated. The main ion-ion energy transfer processes analyzed were energy migration among Yb3+ ions, cross-relaxations between Yb3+ and Tm3+ ions, and interactions with OH- radicals. The results indicated that Yb -> Tm energy transfer favors 1.8 mu m emissions, and there is no evidence of concentration quenching up to 2% Tm2O3 doping. As expected, samples prepared in nitrogen atmosphere present higher fluorescence quantum efficiency than those prepared in air, and this feature is specially noted in the near-infrared region, where the interaction with the OH- radicals is more pronounced. (c) 2007 Published by Elsevier B.V.

Structural matching of 2D electrophoresis gels using deformed graphs

2D electrophoresis is a well-known method for protein separation which is extremely useful in the field of proteomics. Each spot in the image represents a protein accumulation and the goal is to perform a differential analysis between pairs of images to study changes in protein content. It is thus necessary to register two images by finding spot correspondences. Although it may seem a simple task, generally, the manual processing of this kind of images is very cumbersome, especially when strong variations between corresponding sets of spots are expected (e.g. strong non-linear deformations and outliers). In order to solve this problem, this paper proposes a new quadratic assignment formulation together with a correspondence estimation algorithm based on graph matching which takes into account the structural information between the detected spots. Each image is represented by a graph and the task is to find a maximum common subgraph. Successful experimental results using real data are presented, including an extensive comparative performance evaluation with ground-truth data. (C) 2010 Elsevier B.V. All rights reserved.

A fast and robust statistical test based on likelihood ratio with Bartlett correction to identify Granger causality between gene sets

We propose a likelihood ratio test ( LRT) with Bartlett correction in order to identify Granger causality between sets of time series gene expression data. The performance of the proposed test is compared to a previously published bootstrapbased approach. LRT is shown to be significantly faster and statistically powerful even within non- Normal distributions. An R package named gGranger containing an implementation for both Granger causality identification tests is also provided.

Network genealogy of 195-bp satellite DNA supports the superimposed hybridization hypothesis of Trypanosoma cruzi evolutionary pattern

Trypanosoma cruzi is highly diverse genetically and has been partitioned into six discrete typing units (DTUs), recently re-named T. cruzi I-VI. Although T. cruzi reproduces predominantly by binary division, accumulating evidence indicates that particular DTUs are the result of hybridization events. Two major scenarios for the origin of the hybrid lineages have been proposed. It is accepted widely that the most heterozygous TcV and TcVI DTUs are the result of genetic exchange between TcII and TcIII strains. On the other hand, the participation of a TcI parental in the current genome structure of these hybrid strains is a matter of debate. Here, sequences of the T. cruzi-specific 195-bp satellite DNA of TcI, TcII, Tat, TcV, and TcVI strains have been used for inferring network genealogies. The resulting genealogy showed a high degree of reticulation, which is consistent with more than one event of hybridization between the Tc DTUs. The data also strongly suggest that Tat is a hybrid with two distinct sets of satellite sequences, and that genetic exchange between TcI and TcII parentals occurred within the pedigree of the TcV and TcVI DTUs. Although satellite DNAs belong to the fast-evolving portion of eukaryotic genomes, in >100 satellite units of nine T. cruzi strains we found regions that display 100% identity. No DTU-specific consensus motifs were identified, inferring species-wide conservation. (C) 2010 Elsevier B.V. All rights reserved.

New Relativistic Atomic Natural Orbital Basis Sets for Lanthanide Atoms with Applications to the Ce Diatom and LuF3

New basis sets of the atomic natural orbital (ANO) type have been developed for the lanthanide atoms La-Lu. The ANOs have been obtained from the average density matrix of the ground and lowest excited states of the atom, the positive ions, and the atom in an electric field. Scalar relativistic effects are included through the use of a Douglas-Kroll-Hess Hamiltonian. Multiconfigurational wave functions have been used with dynamic correlation included using second-order perturbation theory (CASSCF/CASPT2). The basis sets are applied in calculations of ionization energies and some excitation energies. Computed ionization energies have an accuracy better than 0.1 eV in most cases. Two molecular applications are inluded as illustration: the cerium diatom and the LuF3 molecule. In both cases it is shown that 4f orbitals are not involved in the chemical bond in contrast to an earlier claim for the latter molecule.

Solvation in Pure Liquids: What Can Be Learned from the Use of Pairs of Indicators?

The solvation of six solvatochromic probes in a large number of solvents (33-68) was examined at 25 degrees C. The probes employed were the following: 2,6-diphenyl-4-(2,4,6-triphenylpyridinium-1-yl) phenolate (RB); 4-[(E)2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePM; 1-methylquinolinium-8-olate, QB; 2-bromo-4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePMBr, 2,6-dichloro-4-(2,4,6-triphenyl pyridinium-1-yl) phenolate (WB); and 2,6-dibromo-4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePMBr(2), respectively. Of these, MePMBr is a novel compound. They can be grouped in three pairs, each with similar pK(a) in water but with different molecular properties, for example, lipophilicity and dipole moment. These pairs are formed by RB and MePM; QB and MePMBr; WB and MePMBr(2), respectively. Theoretical calculations were carried out in order to calculate their physicochemical properties including bond lengths, dihedral angles, dipole moments, and wavelength of absorption of the intramolecular charge-transfer band in four solvents, water, methanol, acetone, and DMSO, respectively. The data calculated were in excellent agreement with available experimental data, for example, bond length and dihedral angles. This gives credence to the use of the calculated properties in explaining the solvatochromic behaviors observed. The dependence of an empirical solvent polarity scale E(T)(probe) in kcal/mol on the physicochemical properties of the solvent (acidity, basicity, and dipolarity/polarizability) and those of the probes (pK(a), and dipole moment) was analyzed by using known multiparameter solvation equations. For each pair of probes, values of E(T)(probe) (for example, E(T)(MePM) versus E(T)(RB)) were found to be linearly correlated with correlation coefficients, r, between 0.9548 and 0.9860. For the mercyanine series, the values of E(T)(probe) also correlated linearly, with (r) of 0.9772 (MePMBr versus MePM) and 0.9919 (MePMBr(2) versus MePM). The response of each pair of probes (of similar pK(a)) to solvent acidity is the same, provided that solute-solvent hydrogen-bonding is not seriously affected by steric crowding (as in case of RB). We show, for the first time, that the response to solvent dipolarity/polarizability is linearly correlated to the dipole moment of the probes. The successive introduction of bromine atoms in MePM (to give MePMBr, then MePMBr(2)) leads to the following linear decrease: pK(a) in water, length of the phenolate oxygen-carbon bond, length of the central ethylenic bond, susceptibility to solvent acidity, and susceptibility to solvent dipolarity/polarizability. Thus studying the solvation of probes whose molecular structures are varied systematically produces a wealth of information on the effect of solute structure on its solvation. The results of solvation of the present probes were employed in order to test the goodness of fit of two independent sets of solvent solvatochromic parameters.

Volsurf Descriptors to Analyse Anti-HCV and Cytotoxic Activities of Sesquiterpene Lactones from Asteraceae Family

Various significant anti-HCV and cytotoxic sesquiterpene lactones (SLs) have been characterized. In this work, the chemometric tool Principal Component Analysis (PCA) was applied to two sets of SLs and the variance of the biological activity was explored. The first principal component accounts for as much of the variability in the data as possible, and each succeeding component accounts for as much of the remaining variability as possible. The calculations were performed using VolSurf program. For anti-HCV activity, PC1 (First Principal Component) explained 30.3% and PC2 (Second Principal Component) explained 26.5% of matrix total variance, while for cytotoxic activity, PC1 explained 30.9% and PC2 explained 15.6% of the total variance. The formalism employed generated good exploratory and predictive results and we identified some structural features, for both sets, important to the suitable biological activity and pharmacokinetic profile.

Time Evolution of the Activation Energy in a Batch Chemical Oscillator

The batch-operated bromate/phosphate/acetone/dual catalyst system was studied at four temperatures between 5 and 35 degrees C. The dynamics was simultaneously followed by potential measurements with platinum and bromide selective electrodes, and spectroscopically at two different wavelengths. By simultaneously recording these four time series it was possible to characterize the dynamics of the sequential oscillations that evolve in time. The existence of three sequential oscillatory patterns at each temperature allowed estimating the activation energies in each case. Along with the activation energy of the induction period, it was possible to trace the time evolution of the overall activation energy at four different stages as the reaction proceeds. The study was carried out for two different sets of initial concentrations and it was observed that the overall activation energy increases as reactants turn into products. This finding was propounded as a result of the decrease in the driving force, or the system`s affinity, of the catalytic oxidative bromination of acetone with acidic bromate, as the closed system evolves toward the thermodynamic equilibrium.

Genetic heterogeneity of residual variance : estimation of variance components using double hierarchical generalized linear models

Background: The sensitivity to microenvironmental changes varies among animals and may be under genetic control. It is essential to take this element into account when aiming at breeding robust farm animals. Here, linear mixed models with genetic effects in the residual variance part of the model can be used. Such models have previously been fitted using EM and MCMC algorithms. Results: We propose the use of double hierarchical generalized linear models (DHGLM), where the squared residuals are assumed to be gamma distributed and the residual variance is fitted using a generalized linear model. The algorithm iterates between two sets of mixed model equations, one on the level of observations and one on the level of variances. The method was validated using simulations and also by re-analyzing a data set on pig litter size that was previously analyzed using a Bayesian approach. The pig litter size data contained 10,060 records from 4,149 sows. The DHGLM was implemented using the ASReml software and the algorithm converged within three minutes on a Linux server. The estimates were similar to those previously obtained using Bayesian methodology, especially the variance components in the residual variance part of the model. Conclusions: We have shown that variance components in the residual variance part of a linear mixed model can be estimated using a DHGLM approach. The method enables analyses of animal models with large numbers of observations. An important future development of the DHGLM methodology is to include the genetic correlation between the random effects in the mean and residual variance parts of the model as a parameter of the DHGLM.

Functional timing analysis of VLSI circuits containing complex gates

The recent advances in CMOS technology have allowed for the fabrication of transistors with submicronic dimensions, making possible the integration of tens of millions devices in a single chip that can be used to build very complex electronic systems. Such increase in complexity of designs has originated a need for more efficient verification tools that could incorporate more appropriate physical and computational models. Timing verification targets at determining whether the timing constraints imposed to the design may be satisfied or not. It can be performed by using circuit simulation or by timing analysis. Although simulation tends to furnish the most accurate estimates, it presents the drawback of being stimuli dependent. Hence, in order to ensure that the critical situation is taken into account, one must exercise all possible input patterns. Obviously, this is not possible to accomplish due to the high complexity of current designs. To circumvent this problem, designers must rely on timing analysis. Timing analysis is an input-independent verification approach that models each combinational block of a circuit as a direct acyclic graph, which is used to estimate the critical delay. First timing analysis tools used only the circuit topology information to estimate circuit delay, thus being referred to as topological timing analyzers. However, such method may result in too pessimistic delay estimates, since the longest paths in the graph may not be able to propagate a transition, that is, may be false. Functional timing analysis, in turn, considers not only circuit topology, but also the temporal and functional relations between circuit elements. Functional timing analysis tools may differ by three aspects: the set of sensitization conditions necessary to declare a path as sensitizable (i.e., the so-called path sensitization criterion), the number of paths simultaneously handled and the method used to determine whether sensitization conditions are satisfiable or not. Currently, the two most efficient approaches test the sensitizability of entire sets of paths at a time: one is based on automatic test pattern generation (ATPG) techniques and the other translates the timing analysis problem into a satisfiability (SAT) problem. Although timing analysis has been exhaustively studied in the last fifteen years, some specific topics have not received the required attention yet. One such topic is the applicability of functional timing analysis to circuits containing complex gates. This is the basic concern of this thesis. In addition, and as a necessary step to settle the scenario, a detailed and systematic study on functional timing analysis is also presented.

The consistency of welfare judgements with a representative consumer

This paper is about economies with a representative consumer. In general a representative consumer need not exist, although there are several well known sets of sufficient conditions under which Qne will. It is common practice, however, to use the representative consumer hypothesis without specifically assuming any of these. We show, firstly, that it is possible for the utility of the representative consumer to increase when every actual consumer is made worse off. This shows a serious shortcoming of welfare judgements based on the representatíve consumer. Secondly, in economies where this does not occur, there exists a social welfare function, which we construct, which is consistent with welfare judgements based on the utility of the representative consumer. Finally we provide a converse to Samuelson' s 1956 representative consumer result, which relates it to Scitovsky's community indifference curves.