Evolutionary dynamics on rugged fitness landscapes: Exact dynamics and information theoretical aspects

The parallel mutation-selection evolutionary dynamics, in which mutation and replication are independent events, is solved exactly in the case that the Malthusian fitnesses associated to the genomes are described by the random energy model (REM) and by a ferromagnetic version of the REM. The solution method uses the mapping of the evolutionary dynamics into a quantum Ising chain in a transverse field and the Suzuki-Trotter formalism to calculate the transition probabilities between configurations at different times. We find that in the case of the REM landscape the dynamics can exhibit three distinct regimes: pure diffusion or stasis for short times, depending on the fitness of the initial configuration, and a spin-glass regime for large times. The dynamic transition between these dynamical regimes is marked by discontinuities in the mean-fitness as well as in the overlap with the initial reference sequence. The relaxation to equilibrium is described by an inverse time decay. In the ferromagnetic REM, we find in addition to these three regimes, a ferromagnetic regime where the overlap and the mean-fitness are frozen. In this case, the system relaxes to equilibrium in a finite time. The relevance of our results to information processing aspects of evolution is discussed.

Lower Bounds on the Exchange-Correlation Energy in Reduced Dimensions

Bounds on the exchange-correlation energy of many-electron systems are derived and tested. By using universal scaling properties of the electron-electron interaction, we obtain the exponent of the bounds in three, two, one, and quasione dimensions. From the properties of the electron gas in the dilute regime, the tightest estimate to date is given for the numerical prefactor of the bound, which is crucial in practical applications. Numerical tests on various low-dimensional systems are in line with the bounds obtained and give evidence of an interesting dimensional crossover between two and one dimensions.

Dissipation as a mechanism of energy gain

Some properties of the annular billiard under the presence of weak dissipation are studied. We show, in a dissipative system, that the average energy of a particle acquires higher values than its average energy of the conservative case. The creation of attractors, associated with a chaotic dynamics in the conservative regime, both in appropriated regions of the phase space, constitute a generic mechanism to increase the average energy of dynamical systems.

Gene Expression Complex Networks: Synthesis, Identification, and Analysis

Thanks to recent advances in molecular biology, allied to an ever increasing amount of experimental data, the functional state of thousands of genes can now be extracted simultaneously by using methods such as cDNA microarrays and RNA-Seq. Particularly important related investigations are the modeling and identification of gene regulatory networks from expression data sets. Such a knowledge is fundamental for many applications, such as disease treatment, therapeutic intervention strategies and drugs design, as well as for planning high-throughput new experiments. Methods have been developed for gene networks modeling and identification from expression profiles. However, an important open problem regards how to validate such approaches and its results. This work presents an objective approach for validation of gene network modeling and identification which comprises the following three main aspects: (1) Artificial Gene Networks (AGNs) model generation through theoretical models of complex networks, which is used to simulate temporal expression data; (2) a computational method for gene network identification from the simulated data, which is founded on a feature selection approach where a target gene is fixed and the expression profile is observed for all other genes in order to identify a relevant subset of predictors; and (3) validation of the identified AGN-based network through comparison with the original network. The proposed framework allows several types of AGNs to be generated and used in order to simulate temporal expression data. The results of the network identification method can then be compared to the original network in order to estimate its properties and accuracy. Some of the most important theoretical models of complex networks have been assessed: the uniformly-random Erdos-Renyi (ER), the small-world Watts-Strogatz (WS), the scale-free Barabasi-Albert (BA), and geographical networks (GG). The experimental results indicate that the inference method was sensitive to average degree k variation, decreasing its network recovery rate with the increase of k. The signal size was important for the inference method to get better accuracy in the network identification rate, presenting very good results with small expression profiles. However, the adopted inference method was not sensible to recognize distinct structures of interaction among genes, presenting a similar behavior when applied to different network topologies. In summary, the proposed framework, though simple, was adequate for the validation of the inferred networks by identifying some properties of the evaluated method, which can be extended to other inference methods.

Molecular dynamics simulation of liquid trimethylphosphine

Structural and dynamical properties of liquid trimethylphosphine (TMP), (CH(3))(3)P, as a function of temperature is investigated by molecular dynamics (MD) simulations. The force field used in the MD simulations, which has been proposed from molecular mechanics and quantum chemistry calculations, is able to reproduce the experimental density of liquid TMP at room temperature. Equilibrium structure is investigated by the usual radial distribution function, g(r), and also in the reciprocal space by the static structure factor, S(k). On the basis of center of mass distances, liquid TMP behaves like a simple liquid of almost spherical particles, but orientational correlation due to dipole-dipole interactions is revealed at short-range distances. Single particle and collective dynamics are investigated by several time correlation functions. At high temperatures, diffusion and reorientation occur at the same time range as relaxation of the liquid structure. Decoupling of these dynamic properties starts below ca. 220 K, when rattling dynamics of a given TMP molecules due to the cage effect of neighbouring molecules becomes important. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3624408]

Full-dimensional analytical ab initio potential energy surface of the ground state of HOI

Extensive ab initio calculations using a complete active space second-order perturbation theory wavefunction, including scalar and spin-orbit relativistic effects with a quadruple-zeta quality basis set were used to construct an analytical potential energy surface (PES) of the ground state of the [H, O, I] system. A total of 5344 points were fit to a three-dimensional function of the internuclear distances, with a global root-mean-square error of 1.26 kcal mol(-1). The resulting PES describes accurately the main features of this system: the HOI and HIO isomers, the transition state between them, and all dissociation asymptotes. After a small adjustment, using a scaling factor on the internal coordinates of HOI, the frequencies calculated in this work agree with the experimental data available within 10 cm(-1). (C) 2011 American Institute of Physics. [doi: 10.1063/1.3615545]

Land-Water interactions in the amazon

Biogeochemistry is hosting this special thematic issue devoted to studies of land-water interactions, as part of the Large-scale Biosphere-Atmosphere Experiment in Amaznia (LBA). This compilation of papers covers a broad range of topics with a common theme of coupling land and water processes, across pristine and impacted systems. Findings highlighted that hydrologic flowpaths are clearly important across basin size and structure in determining how water and solutes reach streams. Land-use changes have pronounced impacts on flowpaths, and subsequently, on stream chemistry, from small streams to large rivers. Carbon is produced and transformed across a broad array of fluvial environments and wetlands. Surface waters are not only driven by, but provide feedback to, the atmosphere.

Small rivers in the southwestern Amazon and their role in CO(2) outgassing

A recent estimate of CO(2) outgassing from Amazonian wetlands suggests that an order of magnitude more CO(2) leaves rivers through gas exchange with the atmosphere than is exported to the ocean as organic plus inorganic carbon. However, the contribution of smaller rivers is still poorly understood, mainly because of limitations in mapping their spatial extent. Considering that the largest extension of the Amazon River network is composed of small rivers, the authors` objective was to elucidate their role in air-water CO(2) exchange by developing a geographic information system ( GIS)- based model to calculate the surface area covered by rivers with channels less than 100 m wide, combined with estimated CO(2) outgassing rates at the Ji-Parana River basin, in the western Amazon. Estimated CO(2) outgassing was the main carbon export pathway for this river basin, totaling 289 Gg C yr(-1), about 2.4 times the amount of carbon exported as dissolved inorganic carbon ( 121 Gg C yr(-1)) and 1.6 times the dissolved organic carbon export ( 185 Gg C yr(-1)). The relationships established here between drainage area and channel width provide a new model for determining small river surface area, allowing regional extrapolations of air - water gas exchange. Applying this model to the entire Amazon River network of channels less than 100 m wide ( third to fifth order), the authors calculate that the surface area of small rivers is 0.3 +/- 0.05 million km(2), and it is potentially evading to the atmosphere 170 +/- 42 Tg C yr(-1) as CO(2). Therefore, these ecosystems play an important role in the regional carbon balance.

Cadmium telluride nanocrystals as luminescent sensitizers in flow analysis

A fully automated multipumping flow system (MPFS) using water-soluble CdTe quantum dots (QD) as sensitizers is proposed for the chemiluminometric determination of the anti-diabetic drugs gliclazide and glipizide in pharmaceutical formulations. The nanocrystals acted as enhancers of the weak CL emission produced upon oxidation of sulphite by Ce(IV) in acidic medium, thus improving sensitivity and expanding the dynamical analytical concentration range. By interacting with the QD, the two analytes prevented their sensitizing effect yielding a chemiluminescence quenching of the Ce(IV)-SO(3)(2-)CdTe QD system. The pulsed flow inherent to MPFS assured a fast and efficient mixing of all solutions inside the flow cell, circumventing the need for a reaction coil and facilitating the monitoring of the short-lived generated chemiluminescent species. QD crystal size, concentration and spectral region for measurement were investigated. (C) 2011 Elsevier B.V. All rights reserved.

Low dimensional behavior in three-dimensional coupled map lattices

The analysis of one-, two-, and three-dimensional coupled map lattices is here developed under a statistical and dynamical perspective. We show that the three-dimensional CML exhibits low dimensional behavior with long range correlation and the power spectrum follows 1/f noise. This approach leads to an integrated understanding of the most important properties of these universal models of spatiotemporal chaos. We perform a complete time series analysis of the model and investigate the dependence of the signal properties by change of dimension. (c) 2008 Elsevier Ltd. All rights reserved.

Multifractality in the random parameter model for multivariate time series

The Random Parameter model was proposed to explain the structure of the covariance matrix in problems where most, but not all, of the eigenvalues of the covariance matrix can be explained by Random Matrix Theory. In this article, we explore the scaling properties of the model, as observed in the multifractal structure of the simulated time series. We use the Wavelet Transform Modulus Maxima technique to obtain the multifractal spectrum dependence with the parameters of the model. The model shows a scaling structure compatible with the stylized facts for a reasonable choice of the parameter values. (C) 2009 Elsevier B.V. All rights reserved.

Exploring the emission intensities of ICP OES aided by chemometrics in the geographical discrimination of mineral waters

A rapid method for classification of mineral waters is proposed. The discrimination power was evaluated by a novel combination of chemometric data analysis and qualitative multi-elemental fingerprints of mineral water samples acquired from different regions of the Brazilian territory. The classification of mineral waters was assessed using only the wavelength emission intensities obtained by inductively coupled plasma optical emission spectrometry (ICP OES), monitoring different lines of Al, B, Ba, Ca, Cl, Cu, Co, Cr, Fe, K, Mg, Mn, Na, Ni, P, Pb, S, Sb, Si, Sr, Ti, V, and Zn, and Be, Dy, Gd, In, La, Sc and Y as internal standards. Data acquisition was done under robust (RC) and non-robust (NRC) conditions. Also, the combination of signal intensities of two or more emission lines for each element were evaluated instead of the individual lines. The performance of two classification-k-nearest neighbor (kNN) and soft independent modeling of class analogy (SIMCA)-and preprocessing algorithms, autoscaling and Pareto scaling, were evaluated for the ability to differentiate between the various samples in each approach tested (combination of robust or non-robust conditions with use of individual lines or sum of the intensities of emission lines). It was shown that qualitative ICP OES fingerprinting in combination with multivariate analysis is a promising analytical tool that has potential to become a recognized procedure for rapid authenticity and adulteration testing of mineral water samples or other material whose physicochemical properties (or origin) are directly related to mineral content.

Influence of Granulometry and Organic Treatment of a Brazilian Montmorillonite on the Properties of Poly(styrene-co-n-butyl acrylate)/Layered Silicate Nanocomposites Prepared by Miniemulsion Polymerization

The influence of granulometry and organic treatment of a Brazilian montmorillonite (MMT) clay on the synthesis and properties of poly(styrene-co-n-butyl acrylate)/layered silicate nanocomposites was studied. Hybrid latexes of poly(styrene-co-butyl acrylate)/MMT were synthesized via miniemulsion polymerization using either sodium or organically modified MMT. Five clay granulometries ranging from clay particles smaller than 75 mu m to colloidal size were selected. The size of the clay particles was evaluated by Specific surface area measurements (BET). Cetyl trimethyl ammonium chloride was used as an organic modifier to enhance the clay compatibility with the monomer phase before polymerization and to improve the clav distribution and dispersion within the polymeric matrix after polymerization. The sodium and organically modified natural clays as well as the composites were characterized by X-ray diffraction analysis. The latexes were characterized by dynamic light scattering. The mechanical, thermal, and rheological properties of the composites obtained were characterized by dynamical-mechanical analysis, thermogravimetry, and small amplitude oscillatory, shear tests, respectively. The results showed that smaller the size of the organically modified MMT, the higher the degree of exfoliation of nanoplatelets. Hybrid latexes in presence of Na-MMT resulted in materials with intercalated structures. (C) 2009 Wiley, Periodicals, Inc. J Appl Polym Sci 112: 1949-1958, 2009

SYNCHRONIZATION OF A CLASS OF SECOND-ORDER NONLINEAR SYSTEMS

The asymptotic behavior of a class of coupled second-order nonlinear dynamical systems is studied in this paper. Using very mild assumptions on the vector-field, conditions on the coupling parameters that guarantee synchronization are provided. The proposed result does not require solutions to be ultimately bounded in order to prove synchronization, therefore it can be used to study coupled systems that do not globally synchronize, including synchronization of unbounded solutions. In this case, estimates of the synchronization region are obtained. Synchronization of two-coupled nonlinear pendulums and two-coupled Duffing systems are studied to illustrate the application of the proposed theory.

Recursive linear estimation for general discrete-time descriptor systems

This paper considers the optimal linear estimates recursion problem for discrete-time linear systems in its more general formulation. The system is allowed to be in descriptor form, rectangular, time-variant, and with the dynamical and measurement noises correlated. We propose a new expression for the filter recursive equations which presents an interesting simple and symmetric structure. Convergence of the associated Riccati recursion and stability properties of the steady-state filter are provided. (C) 2010 Elsevier Ltd. All rights reserved.