Density functional theory investigation of 3d, 4d, and 5d 13-atom metal clusters

The knowledge of the atomic structure of clusters composed by few atoms is a basic prerequisite to obtain insights into the mechanisms that determine their chemical and physical properties as a function of diameter, shape, surface termination, as well as to understand the mechanism of bulk formation. Due to the wide use of metal systems in our modern life, the accurate determination of the properties of 3d, 4d, and 5d metal clusters poses a huge problem for nanoscience. In this work, we report a density functional theory study of the atomic structure, binding energies, effective coordination numbers, average bond lengths, and magnetic properties of the 3d, 4d, and 5d metal (30 elements) clusters containing 13 atoms, M(13). First, a set of lowest-energy local minimum structures (as supported by vibrational analysis) were obtained by combining high-temperature first- principles molecular-dynamics simulation, structure crossover, and the selection of five well-known M(13) structures. Several new lower energy configurations were identified, e. g., Pd(13), W(13), Pt(13), etc., and previous known structures were confirmed by our calculations. Furthermore, the following trends were identified: (i) compact icosahedral-like forms at the beginning of each metal series, more opened structures such as hexagonal bilayerlike and double simple-cubic layers at the middle of each metal series, and structures with an increasing effective coordination number occur for large d states occupation. (ii) For Au(13), we found that spin-orbit coupling favors the three-dimensional (3D) structures, i.e., a 3D structure is about 0.10 eV lower in energy than the lowest energy known two-dimensional configuration. (iii) The magnetic exchange interactions play an important role for particular systems such as Fe, Cr, and Mn. (iv) The analysis of the binding energy and average bond lengths show a paraboliclike shape as a function of the occupation of the d states and hence, most of the properties can be explained by the chemistry picture of occupation of the bonding and antibonding states.

Antiferromagnetism with spin polarization of GaN-based diluted magnetic semiconductors

The diluted magnetic semiconductors are promising materials for spintronic applications. Usually one intents to find the ferromagnetic state but recently the antiferromagnetism (AFM) was proposed to have some advantages. In this work, we verify the possibility to obtain spin polarization with an AFM state. In particular, we studied GaN 5% double doped with two different transition metals atoms (Mn and Co or Cr and Ni), forming the Mn(x)Co(0.056-x)Ga(0.944)N and Cr(x)Ni(0.056-x)Ga(0.944)N quaternary alloys. In order to simulate these systems in a more realistic way, and take into account composition fluctuations, we adapted the generalized quasichemical approach to diluted alloys, which is used in combination with spin density-functional theory. We find that is possible to obtain an AFM ground state up to 70% spin polarization.

Density functional theory and ab initio molecular dynamics study of NO adsorption on Pd(111) and Pt(111) surfaces

The origin of the unique geometry for nitric oxide (NO) adsorption on Pd(111) and Pt(111) surfaces as well as the effect of temperature were studied by density functional theory calculations and ab initio molecular dynamics at finite temperature. We found that at low coverage, the adsorption geometry is determined by electronic interactions, depending sensitively on the adsorption sites and coverages, and the effect of temperature on geometries is significant. At coverage of 0.25 monolayer (ML), adsorbed NO at hollow sites prefer an upright configuration, while NO adsorbed at top sites prefer a tilting configuration. With increase in the coverage up to 0.50 ML, the enhanced steric repulsion lead to the tilting of hollow NO. We found that the tilting was enhanced by the thermal effects. At coverage of 0.75 ML with p(2 x 2)-3NO(fcc+hcp+top) structure, we found that there was no preferential orientation for tilted top NO. The interplay of the orbital hybridization, thermal effects, steric repulsion, and their effects on the adsorption geometries were highlighted at the end.

Theoretical investigation of atomic and electronic structures of Ga(2)O(3)(ZnO)(6)

Transparent conducting oxides (TCO) are widely used in technological applications ranging from photovoltaics to thin-film transparent field-effect transistors. In this work we report a first-principles investigation, based on density-functional theory, of the atomic and electronic properties of Ga(2)O(3)(ZnO)(6) (GZO(6)), which is a promising candidate to be used as host oxide for wide band gap TCO applications. We identify a low-energy configuration for the coherent distribution of the Ga and Zn atoms in the cation positions within the experimentally reported orthorhombic GZO(6) structure. Four Ga atoms are located in four-fold sites, while the remaining 12 Ga atoms in the unit cell form four shared Ga agglomerates (a motif of four atoms). The Zn atoms are distributed in the remaining cation sites with effective coordination numbers from 3.90 to 4.50. Furthermore, we identify the natural formation of twin-boundaries in GZO(6), which can explain the zigzag modulations observed experimentally by high-resolution transmission electron microscopy in GZO(n) (n=9). Due to the intrinsic twin-boundary formation, polarity inversion in the ZnO tetrahedrons is present which is facilitated by the formation of the Ga agglomerates. Our analysis shows that the formation of fourfold Ga sites and Ga agglomerates are stabilized by the electronic octet rule, while the distribution of Ga atoms and the formation of the twin-boundary help alleviate excess strain. Finally we identify that the electronic properties of GZO(6) are essentially determined by the electronic properties of ZnO, i.e., there are slight changes in the band gap and optical absorption properties.

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.

Feature selection environment for genomic applications

Background: Feature selection is a pattern recognition approach to choose important variables according to some criteria in order to distinguish or explain certain phenomena (i.e., for dimensionality reduction). There are many genomic and proteomic applications that rely on feature selection to answer questions such as selecting signature genes which are informative about some biological state, e. g., normal tissues and several types of cancer; or inferring a prediction network among elements such as genes, proteins and external stimuli. In these applications, a recurrent problem is the lack of samples to perform an adequate estimate of the joint probabilities between element states. A myriad of feature selection algorithms and criterion functions have been proposed, although it is difficult to point the best solution for each application. Results: The intent of this work is to provide an open-source multiplataform graphical environment for bioinformatics problems, which supports many feature selection algorithms, criterion functions and graphic visualization tools such as scatterplots, parallel coordinates and graphs. A feature selection approach for growing genetic networks from seed genes ( targets or predictors) is also implemented in the system. Conclusion: The proposed feature selection environment allows data analysis using several algorithms, criterion functions and graphic visualization tools. Our experiments have shown the software effectiveness in two distinct types of biological problems. Besides, the environment can be used in different pattern recognition applications, although the main concern regards bioinformatics tasks.

An algorithmic Friedman-Pippenger theorem on tree embeddings and applications

An (n, d)-expander is a graph G = (V, E) such that for every X subset of V with vertical bar X vertical bar <= 2n - 2 we have vertical bar Gamma(G)(X) vertical bar >= (d + 1) vertical bar X vertical bar. A tree T is small if it has at most n vertices and has maximum degree at most d. Friedman and Pippenger (1987) proved that any ( n; d)- expander contains every small tree. However, their elegant proof does not seem to yield an efficient algorithm for obtaining the tree. In this paper, we give an alternative result that does admit a polynomial time algorithm for finding the immersion of any small tree in subgraphs G of (N, D, lambda)-graphs Lambda, as long as G contains a positive fraction of the edges of Lambda and lambda/D is small enough. In several applications of the Friedman-Pippenger theorem, including the ones in the original paper of those authors, the (n, d)-expander G is a subgraph of an (N, D, lambda)-graph as above. Therefore, our result suffices to provide efficient algorithms for such previously non-constructive applications. As an example, we discuss a recent result of Alon, Krivelevich, and Sudakov (2007) concerning embedding nearly spanning bounded degree trees, the proof of which makes use of the Friedman-Pippenger theorem. We shall also show a construction inspired on Wigderson-Zuckerman expander graphs for which any sufficiently dense subgraph contains all trees of sizes and maximum degrees achieving essentially optimal parameters. Our algorithmic approach is based on a reduction of the tree embedding problem to a certain on-line matching problem for bipartite graphs, solved by Aggarwal et al. (1996).

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.

Functional microarray analysis suggests repressed cell-cell signaling and cell survival-related modules inhibit progression of head and neck squamous cell carcinoma

Background: Cancer shows a great diversity in its clinical behavior which cannot be easily predicted using the currently available clinical or pathological markers. The identification of pathways associated with lymph node metastasis (N+) and recurrent head and neck squamous cell carcinoma (HNSCC) may increase our understanding of the complex biology of this disease. Methods: Tumor samples were obtained from untreated HNSCC patients undergoing surgery. Patients were classified according to pathologic lymph node status (positive or negative) or tumor recurrence (recurrent or non-recurrent tumor) after treatment (surgery with neck dissection followed by radiotherapy). Using microarray gene expression, we screened tumor samples according to modules comprised by genes in the same pathway or functional category. Results: The most frequent alterations were the repression of modules in negative lymph node (N0) and in non-recurrent tumors rather than induction of modules in N+ or in recurrent tumors. N0 tumors showed repression of modules that contain cell survival genes and in non-recurrent tumors cell-cell signaling and extracellular region modules were repressed. Conclusions: The repression of modules that contain cell survival genes in N0 tumors reinforces the important role that apoptosis plays in the regulation of metastasis. In addition, because tumor samples used here were not microdissected, tumor gene expression data are represented together with the stroma, which may reveal signaling between the microenvironment and tumor cells. For instance, in non-recurrent tumors, extracellular region module was repressed, indicating that the stroma and tumor cells may have fewer interactions, which disable metastasis development. Finally, the genes highlighted in our analysis can be implicated in more than one pathway or characteristic, suggesting that therapeutic approaches to prevent tumor progression should target more than one gene or pathway, specially apoptosis and interactions between tumor cells and the stroma.

Spectroscopic, structural, and functional characterization of the alternative low-spin state of horse heart cytochrome c

The alternative low-spin states of Fe3+ and Fe2+ cytochrome c induced by SDS or AOT/hexane reverse micelles exhibited the heme group in a less rhombic symmetry and were characterized by electron paramagnetic resonance, UV-visible, CD, magnetic CD, fluorescence, and Raman resonance. Consistent with the replacement of Met 80 by another strong field ligand at the sixth heme iron coordination position, Fe3+ ALSScytc exhibited 1-nm Soret band blue shift and e enhancement accompanied by disappearance of the 695-nm charge transfer band. The Raman resonance, CD, and magnetic CD spectra of Fe3+ and Fe2+ ALSScytc exhibited significant changes suggestive of alterations in the heme iron microenvironment and conformation and should not be assigned to unfold because the Trp(59) fluorescence remained quenched by the neighboring heme group. ALSScytc was obtained with His(33) and His(26) carboxyethoxylated horse cytochrome c and with tuna cytochrome c (His(33) replaced by Asn) pointing out Lys(79) as the probable heme iron ligand. Fe3+ ALSScytc retained the capacity to cleave tert-butylhydroperoxide and to be reduced by dithiothreitol and diphenylacetaldehyde but not by ascorbate. Compatible with a more open heme crevice, ALSScytc exhibited a redox potential similar to 200 mV lower than the wild-type protein (1220 mV) and was more susceptible to the attack of free radicals.

Structure, stability, depolarized light scattering, and vibrational spectra of fullerenols from all-electron density-functional-theory calculations

We have investigated the stability, electronic properties, Rayleigh (elastic), and Raman (inelastic) depolarization ratios, infrared and Raman absorption vibrational spectra of fullerenols [C(60)(OH)(n)] with different degrees of hydroxylation by using all-electron density-functional-theory (DFT) methods. Stable arrangements of these molecules were found by means of full geometry optimizations using Becke's three-parameter exchange functional with the Lee, Yang, and Parr correlation functional. This DFT level has been combined with the 6-31G(d,p) Gaussian-type basis set, as a compromise between accuracy and capability to treat highly hydroxylated fullerenes, e.g., C(60)(OH)(36). Thus, the molecular properties of fullerenols were systematically analyzed for structures with n=1, 2, 3, 4, 8, 10, 16, 18, 24, 32, and 36. From the electronic structure analysis of these molecules, we have evidenced an important effect related to the weak chemical reactivity of a possible C(60)(OH)(24) isomer. To investigate Raman scattering and the vibrational spectra of the different fullerenols, frequency calculations are carried out within the harmonic approximation. In this case a systematic study is only performed for n=1-4, 8, 10, 16, 18, and 24. Our results give good agreements with the expected changes in the spectral absorptions due to the hydroxylation of fullerenes.

Bulk structures of PtO and PtO(2) from density functional calculations

Platinum plays an important role in catalysis and electrochemistry, and it is known that the direct interaction of oxygen with Pt surfaces can lead to the formation of platinum oxides (PtO(x)), which can affect the reactivity. To contribute to the atomistic understanding of the atomic structure of PtO(x), we report a density functional theory study of the atomic structure of bulk PtO(x) (1 <= x <= 2). From our calculations, we identified a lowest-energy structure (GeS type, space group Pnma) for PtO, which is 0.181 eV lower in energy than the structure suggested by W. J. Moore and L. Pauling [J. Am. Chem. Soc. 63, 1392 (1941)] (PtS type). Furthermore, two atomic structures were identified for PtO(2), which are almost degenerate in energy with the lowest-energy structure reported so far for PtO(2) (CaCl(2) type). Based on our results and analysis, we suggest that Pt and O atoms tend to form octahedron motifs in PtO(x) even at lower O composition by the formation of Pt-Pt bonds.

Functional specialization of Eucalyptus fine roots: contrasting potential uptake rates for nitrogen, potassium and calcium tracers at varying soil depths

1. Little is known about the role of deep roots in the nutrition of forest trees and their ability to provide a safety-net service taking up nutrients leached from the topsoil. 2. To address this issue, we studied the potential uptake of N, K and Ca by Eucalyptus grandis trees (6 years of age - 25 m mean height), in Brazil, as a function of soil depth, texture and water content. We injected NO(3)(-)- (15)N, Rb(+) (analogue of K(+)) and Sr(2+) (analogue of Ca(2+)) tracers simultaneously in a solution through plastic tubes at 10, 50, 150 and 300 cm in depth in a sandy and a clayey Ferralsol soil. A complete randomized design was set up with three replicates of paired trees per injection depth and soil type. Recently expanded leaves were sampled at various times after tracer injection in the summer, and the experiment was repeated in the winter. Soil water contents were continuously monitored at the different depths in the two soils. 3. Determination of foliar Rb and Sr concentrations and (15)N atom % made it possible to estimate the relative uptake potential (RUP) of tracer injections from the four soil depths and the specific RUP (SRUP), defined as RUP, per unit of fine root length density in the corresponding soil layer. 4. The highest tracer uptake rates were found in the topsoil, but contrasting RUP distributions were observed for the three tracers. Whilst the RUP was higher for NO(3)(-)- (15)N than for Rb(+) and Sr(2+) in the upper 50 cm of soil, the highest SRUP values for Sr(2+) and Rb(+) were found at a depth of 300 cm in the sandy soil, as well as in the clayey soil when gravitational solutions reached that depth. 5. Our results suggest that the fine roots of E. grandis trees exhibit contrasting potential uptake rates with depth depending on the nutrient. This functional specialization of roots might contribute to the high growth rates of E. grandis trees, efficiently providing the large amounts of nutrients required throughout the development of these fast-growing plantations.

Economic viability of doses and split-applications of nitrogen fertilization in corn crop in a eutrophic Red Latosol

Nitrogen is the nutrient that is most absorbed by the corn crop, with the most complex management, and has the highest share on the cost of corn production. The objective of this work was to evaluate the economic viability of different rates and split-applications of nitrogen fertilization, as such as urea, in the corn crop in a eutrophic Red Latosol (Oxisol). The study was carried out in the Experimental Station of the Regional Pole of the Sao Paulo Northwest Agribusiness Development (APTA), in Votuporanga, State of Sao Paulo, Brazil. The experimental design was randomized complete blocks with nine treatments and four replications, consisting of five N rates: 0, 55, 95, 135 and 175 kg ha(-1), 15 kg ha-l applied in the seeding and the remainder in top dressing: 40 and 80 kg ha(-1) N at forty days after seeding (DAS), or 1/2 + 1/2 at 20 and 40 DAS; 120 kg ha-1 N split in 1/2 + 1/2 or 1/3 + 1/3 + 1/3 at 20, 40 or 60 DAS; 160 kg ha(-1) N split in 1/4 + 3/8 + 3/8 or 114 + 1/4 + 1/4 + 1/4 at 20, 40, 60 and 80 DAS. The application of 135 kg ha-l of N split in three times provided the best benefit/cost ratio. The non-application of N provided the lowest economic return, proving to be unviable.

Fractionation of Zn, Cd and Pb in a Tropical Soil After Nine-Year Sewage Sludge Applications

A long-term field experiment was carried out in the experiment farm of the Sao Paulo State University, Brazil, to evaluate the phytoavailability of Zn, Cd and Pb in a Typic Eutrorthox soil treated with sewage sludge for nine consecutive years, using the sequential extraction and organic matter fractionation methods. During 2005-2006, maize (Zea mays L.) was used as test plants and the experimental design was in randomized complete blocks with four treatments and five replicates. The treatments consisted of four sewage sludge rates (in a dry basis): 0.0 (control, with mineral fertilization), 45.0, 90.0 and 127.5 t ha(-1), annually for nine years. Before maize sowing, the sewage sludge was manually applied to the soil and incorporated at 10 cm depth. Soil samples (0-20 cm layer) for Zn, Cd and Pb analysis were collected 60 days after sowing. The successive applications of sewage sludge to the soil did not affect heavy metal (Cd and Pb) fractions in the soil, with exception of Zn fractions. The Zn, Cd and Pb distributions in the soil were strongly associated with humin and residual fractions, which are characterized by stable chemical bonds. Zinc, Cd and Pb in the soil showed low phytoavailability after nine-year successive applications of sewage sludge to the soil.