FLOWS WITHOUT WANDERING POINTS ON COMPACT CONNECTED SURFACES

Given a compact 2 dimensional manifold M we classify all continuous flows phi without wandering points on M. This classification is performed by finding finitely many pairwise disjoint open phi-invariant subsets {U(1), U(2), ..., U(n)} of M such that U(i=1)(n) (U(i)) over bar = M and each U(i) is either a suspension of an interval exchange transformation, or a maximal open cylinder made up of closed trajectories of phi.

Multistability and Self-Similarity in the Parameter-Space of a Vibro-Impact System

The dynamics of a dissipative vibro-impact system called impact-pair is investigated. This system is similar to Fermi-Ulam accelerator model and consists of an oscillating one-dimensional box containing a point mass moving freely between successive inelastic collisions with the rigid walls of the box. In our numerical simulations, we observed multistable regimes, for which the corresponding basins of attraction present a quite complicated structure with smooth boundary. In addition, we characterize the system in a two-dimensional parameter space by using the largest Lyapunov exponents, identifying self-similar periodic sets. Copyright (C) 2009 Silvio L.T. de Souza et al.

New constraints on muon-neutrino to electron-neutrino transitions in MINOS

This paper reports results from a search for nu(mu) -> nu(e) transitions by the MINOS experiment based on a 7 x 10(20) protons-on-target exposure. Our observation of 54 candidate nu(e) events in the far detector with a background of 49.1 +/- 7.0(stat) +/- 2.7(syst) events predicted by the measurements in the near detector requires 2sin(2)(2 theta(13))sin(2)theta(23) < 0.12(0.20) at the 90% C.L. for the normal (inverted) mass hierarchy at delta(CP) = 0. The experiment sets the tightest limits to date on the value of theta(13) for nearly all values of delta(CP) for the normal neutrino mass hierarchy and maximal sin(2)(2 theta(23)).

Dynamic polarizability, Cauchy moments, and the optical absorption spectrum of liquid water: A sequential molecular dynamics/quantum mechanical approach

The dynamic polarizability and optical absorption spectrum of liquid water in the 6-15 eV energy range are investigated by a sequential molecular dynamics (MD)/quantum mechanical approach. The MD simulations are based on a polarizable model for liquid water. Calculation of electronic properties relies on time-dependent density functional and equation-of-motion coupled-cluster theories. Results for the dynamic polarizability, Cauchy moments, S(-2), S(-4), S(-6), and dielectric properties of liquid water are reported. The theoretical predictions for the optical absorption spectrum of liquid water are in good agreement with experimental information.

Electronic properties of liquid ammonia: A sequential molecular dynamics/quantum mechanics approach

The electronic properties of liquid ammonia are investigated by a sequential molecular dynamics/quantum mechanics approach. Quantum mechanics calculations for the liquid phase are based on a reparametrized hybrid exchange-correlation functional that reproduces the electronic properties of ammonia clusters [(NH(3))(n); n=1-5]. For these small clusters, electron binding energies based on Green's function or electron propagator theory, coupled cluster with single, double, and perturbative triple excitations, and density functional theory (DFT) are compared. Reparametrized DFT results for the dipole moment, electron binding energies, and electronic density of states of liquid ammonia are reported. The calculated average dipole moment of liquid ammonia (2.05 +/- 0.09 D) corresponds to an increase of 27% compared to the gas phase value and it is 0.23 D above a prediction based on a polarizable model of liquid ammonia [Deng , J. Chem. Phys. 100, 7590 (1994)]. Our estimate for the ionization potential of liquid ammonia is 9.74 +/- 0.73 eV, which is approximately 1.0 eV below the gas phase value for the isolated molecule. The theoretical vertical electron affinity of liquid ammonia is predicted as 0.16 +/- 0.22 eV, in good agreement with the experimental result for the location of the bottom of the conduction band (-V(0)=0.2 eV). Vertical ionization potentials and electron affinities correlate with the total dipole moment of ammonia aggregates. (c) 2008 American Institute of Physics.

Systematic investigation of a family of gradient-dependent functionals for solids

Eleven density functionals are compared with regard to their performance for the lattice constants of solids. We consider standard functionals, such as the local-density approximation and the Perdew-Burke-Ernzerhof (PBE) generalized-gradient approximation (GGA), as well as variations of PBE GGA, such as PBEsol and similar functionals, PBE-type functionals employing a tighter Lieb-Oxford bound, and combinations thereof. On a test set of 60 solids, we perform a system-by-system analysis for selected functionals and a full statistical analysis for all of them. The impact of restoring the gradient expansion and of tightening the Lieb-Oxford bound is discussed, and confronted with previous results obtained from other codes, functionals or test sets. No functional is uniformly good for all investigated systems, but surprisingly, and pleasingly, the simplest possible modifications to PBE turn out to have the most beneficial effect on its performance. The atomization energy of molecules was also considered and on a testing set of six molecules, we found that the PBE functional is clearly the best, the others leading to strong overbinding.

Observation of charge-dependent azimuthal correlations and possible local strong parity violation in heavy-ion collisions

Parity (P)-odd domains, corresponding to nontrivial topological solutions of the QCD vacuum, might be created during relativistic heavy-ion collisions. These domains are predicted to lead to charge separation of quarks along the orbital momentum of the system created in noncentral collisions. To study this effect, we investigate a three-particle mixed-harmonics azimuthal correlator which is a P-even observable, but directly sensitive to the charge-separation effect. We report measurements of this observable using the STAR detector in Au + Au and Cu + Cu collisions at root s(NN) = 200 and 62 GeV. The results are presented as a function of collision centrality, particle separation in rapidity, and particle transverse momentum. A signal consistent with several of the theoretical expectations is detected in all four data sets. We compare our results to the predictions of existing event generators and discuss in detail possible contributions from other effects that are not related to P violation.

The MX2 macromolecular crystallography beamline: a wiggler X-ray source at the LNLS

The Brazilian Synchrotron Light Laboratory [Laboratorio Nacional de Luz Sincrotron (LNLS), Campinas, SP, Brazil] is the first commissioned synchrotron light source in the southern hemisphere. The first wiggler macromolecular crystallography beamline (MX2) at the LNLS has been recently constructed and brought into operation. Here the technical design, experimental set-up, parameters of the beamline and the first experimental results obtained at MX2 are described. The beamline operates on a 2.0 T hybrid 30-pole wiggler, and its optical layout includes collimating mirror, Si( 111) double-crystal monochromator and toroidal bendable mirror. The measured flux density at the sample position at 8.7 eV reaches 4.8 x 10(11) photons s(-1) mm(-2) (100 mA)(-1). The beamline is equipped with a MarResearch Desktop Beamline Goniostat (MarDTB) and 3 x 3 MarMosaic225 CCD detector, and is controlled by a customized version of the Blu-Ice software. A description of the first X-ray diffraction data sets collected at the MX2 LNLS beamline and used for macromolecular crystal structure solution is also provided.

Crystallization and preliminary X-ray diffraction analysis of recombinant chlorocatechol 1,2-dioxygenase from Pseudomonas putida

Chlorocatechol 1,2-dioxygenase from the Gram-negative bacterium Pseudomonas putida (Pp 1,2-CCD) is considered to be an important biotechnological tool owing to its ability to process a broad spectrum of organic pollutants. In the current work, the crystallization, crystallographic characterization and phasing of the recombinant Pp 1,2-CCD enzyme are described. Reddish-brown crystals were obtained in the presence of polyethylene glycol and magnesium acetate by utilizing the vapour-diffusion technique in sitting drops. Crystal dehydration was the key step in obtaining data sets, which were collected on the D03B-MX2 beamline at the CNPEM/MCT - LNLS using a MAR CCD detector. Pp 1,2-CCD crystals belonged to space group P6(1)22 and the crystallographic structure of Pp 1,2-CCD has been solved by the MR-SAD technique using Fe atoms as scattering centres and the coordinates of 3-chlorocatechol 1,2-dioxygenase from Rhodococcus opacus (PDB entry 2boy) as the search model. The initial model, which contains three molecules in the asymmetric unit, has been refined to 3.4 A resolution.

Low-energy electron scattering from methanol and ethanol

Measured and calculated differential cross sections for elastic (rotationally unresolved) electron scattering from two primary alcohols, methanol (CH(3)OH) and ethanol (C(2)H(5)OH), are reported. The measurements are obtained using the relative flow method with helium as the standard gas and a thin aperture as the collimating target gas source. The relative flow method is applied without the restriction imposed by the relative flow pressure conditions on helium and the unknown gas. The experimental data were taken at incident electron energies of 1, 2, 5, 10, 15, 20, 30, 50, and 100 eV and for scattering angles of 5 degrees-130 degrees. There are no previous reports of experimental electron scattering differential cross sections for CH(3)OH and C(2)H(5)OH in the literature. The calculated differential cross sections are obtained using two different implementations of the Schwinger multichannel method, one that takes all electrons into account and is adapted for parallel computers, and another that uses pseudopotentials and considers only the valence electrons. Comparison between theory and experiment shows that theory is able to describe low-energy electron scattering from these polyatomic targets quite well.

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.

INFERENCE FOR EIGENVALUES AND EIGENVECTORS OF GAUSSIAN SYMMETRIC MATRICES

This article presents maximum likelihood estimators (MLEs) and log-likelihood ratio (LLR) tests for the eigenvalues and eigenvectors of Gaussian random symmetric matrices of arbitrary dimension, where the observations are independent repeated samples from one or two populations. These inference problems are relevant in the analysis of diffusion tensor imaging data and polarized cosmic background radiation data, where the observations are, respectively, 3 x 3 and 2 x 2 symmetric positive definite matrices. The parameter sets involved in the inference problems for eigenvalues and eigenvectors are subsets of Euclidean space that are either affine subspaces, embedded submanifolds that are invariant under orthogonal transformations or polyhedral convex cones. We show that for a class of sets that includes the ones considered in this paper, the MLEs of the mean parameter do not depend on the covariance parameters if and only if the covariance structure is orthogonally invariant. Closed-form expressions for the MLEs and the associated LLRs are derived for this covariance structure.

TESTING STATISTICAL HYPOTHESIS ON RANDOM TREES AND APPLICATIONS TO THE PROTEIN CLASSIFICATION PROBLEM

Efficient automatic protein classification is of central importance in genomic annotation. As an independent way to check the reliability of the classification, we propose a statistical approach to test if two sets of protein domain sequences coming from two families of the Pfam database are significantly different. We model protein sequences as realizations of Variable Length Markov Chains (VLMC) and we use the context trees as a signature of each protein family. Our approach is based on a Kolmogorov-Smirnov-type goodness-of-fit test proposed by Balding et at. [Limit theorems for sequences of random trees (2008), DOI: 10.1007/s11749-008-0092-z]. The test statistic is a supremum over the space of trees of a function of the two samples; its computation grows, in principle, exponentially fast with the maximal number of nodes of the potential trees. We show how to transform this problem into a max-flow over a related graph which can be solved using a Ford-Fulkerson algorithm in polynomial time on that number. We apply the test to 10 randomly chosen protein domain families from the seed of Pfam-A database (high quality, manually curated families). The test shows that the distributions of context trees coming from different families are significantly different. We emphasize that this is a novel mathematical approach to validate the automatic clustering of sequences in any context. We also study the performance of the test via simulations on Galton-Watson related processes.

Heritability of cardiovascular risk factors in a Brazilian population: Baependi Heart Study

Background: The heritability of cardiovascular risk factors is expected to differ between populations because of the different distribution of environmental risk factors, as well as the genetic make-up of different human populations. Methods: The purpose of this analysis was to evaluate genetic and environmental influences on cardiovascular risk factor traits, using a variance component approach, by estimating the heritability of these traits in a sample of 1,666 individuals in 81 families ascertained randomly from a highly admixed population of a city in a rural area in Brazil. Results: Before adjustment for sex, age, age(2), and age x sex interaction, polygenic heritability of systolic (SBP) and diastolic (DBP) blood pressure were 15.0% and 16.4%, waist circumference 26.1%, triglycerides 25.7%, fasting glucose 32.8%, HDL-c 31.2%, total cholesterol 28.6%, LDL-c 26.3%, BMI 39.1%. Adjustment for covariates increased polygenic heritability estimates for all traits mainly systolic and diastolic blood pressure (25.9 and 26.2%, respectively), waist circumference (40.1%), and BMI (51.0%). Conclusion: Heritability estimates for cardiovascular traits in the Brazilian population are high and not significantly different from other studied worldwide populations. Mapping efforts to identify genetic loci associated with variability of these traits are warranted.

ON GENERIC ROTATIONLESS DIFFEOMORPHISMS OF THE ANNULUS

Let f be a C(r)-diffeomorphism of the closed annulus A that preserves the orientation, the boundary components and the Lebesgue measure. Suppose that f has a lift (f) over tilde to the infinite strip (A) over tilde which has zero Lebesgue measure rotation number. If the rotation number of f restricted to both boundary components of (f) over tilde is positive, then for such a generic f (r >= 16), zero is an interior point of its rotation set. This is a partial solution to a conjecture of P. Boyland.