Combining Monte Carlo simulation and density-functional theory to describe the spectral changes of Na(2) in liquid helium

Spectral changes of Na(2) in liquid helium were studied using the sequential Monte Carlo-quantum mechanics method. Configurations composed by Na(2) surrounded by explicit helium atoms sampled from the Monte Carlo simulation were submitted to time-dependent density-functional theory calculations of the electronic absorption spectrum using different functionals. Attention is given to both line shift and line broadening. The Perdew, Burke, and Ernzerhof (PBE1PBE, also known as PBE0) functional, with the PBE1PBE/6-311++G(2d,2p) basis set, gives the spectral shift, compared to gas phase, of 500 cm(-1) for the allowed X (1)Sigma(+)(g) -> B (1)Pi(u) transition, in very good agreement with the experimental value (700 cm(-1)). For comparison, cluster calculations were also performed and the first X (1)Sigma(+)(g) -> A (1)Sigma(+)(u) transition was also considered.

Mn-doped cubic BN as an atomiclike memory device: A density functional study

We investigate the electronic properties of Mn(B) substitutional doping in cubic boron nitride (BN), for different charge states, using density functional theory (DFT) calculations. We show that the neutral Mn has a nonmagnetic ground state (S=0). Upon charge injection, it is unambiguously shown that the Mn(B)(-) has a high-spin configuration with a strong, localized magnetic moment of 5 mu(Bohr). We developed a simple model, parameterized by the DFT results, that allows us to interpret the rules played by the crystal-field and exchange-correlation splitting in the magnetization process.

Edge effects in bilayer graphene nanoribbons: Ab initio total-energy density functional theory calculations

We show that the ground state of zigzag bilayer graphene nanoribbons is nonmagnetic. It also possesses a finite gap, which has a nonmonotonic dependence with the width as a consequence of the competition between bulk and strongly attractive edge interactions. All results were obtained using ab initio total-energy density functional theory calculations with the inclusion of parametrized van der Waals interactions.

Approximation to density functional theory for the calculation of band gaps of semiconductors

The local-density approximation (LDA) together with the half occupation (transitionstate) is notoriously successful in the calculation of atomic ionization potentials. When it comes to extended systems, such as a semiconductor infinite system, it has been very difficult to find a way to half ionize because the hole tends to be infinitely extended (a Bloch wave). The answer to this problem lies in the LDA formalism itself. One proves that the half occupation is equivalent to introducing the hole self-energy (electrostatic and exchange correlation) into the Schrodinger equation. The argument then becomes simple: The eigenvalue minus the self-energy has to be minimized because the atom has a minimal energy. Then one simply proves that the hole is localized, not infinitely extended, because it must have maximal self-energy. Then one also arrives at an equation similar to the self- interaction correction equation, but corrected for the removal of just 1/2 electron. Applied to the calculation of band gaps and effective masses, we use the self- energy calculated in atoms and attain a precision similar to that of GW, but with the great advantage that it requires no more computational effort than standard LDA.

Designing real nanotube-based gas sensors

Using a combination of density functional theory and recursive Green's functions techniques, we present a full description of a large scale sensor, accounting for disorder and different coverages. Here, we use this method to demonstrate the functionality of nitrogen-rich carbon nanotubes as ammonia sensors as an example. We show how the molecules one wishes to detect bind to the most relevant defects on the nanotube, describe how these interactions lead to changes in the electronic transport properties of each isolated defect, and demonstrate that there are significative resistance changes even in the presence of disorder, elucidating how a realistic nanosensor works.

Communication structure of cortical networks

Large-scale cortical networks exhibit characteristic topological properties that shape communication between brain regions and global cortical dynamics. Analysis of complex networks allows the description of connectedness, distance, clustering, and centrality that reveal different aspects of how the network's nodes communicate. Here, we focus on a novel analysis of complex walks in a series of mammalian cortical networks that model potential dynamics of information flow between individual brain regions. We introduce two new measures called absorption and driftness. Absorption is the average length of random walks between any two nodes, and takes into account all paths that may diffuse activity throughout the network. Driftness is the ratio between absorption and the corresponding shortest path length. For a given node of the network, we also define four related measurements, namely in-and out-absorption as well as in-and out-driftness, as the averages of the corresponding measures from all nodes to that node, and from that node to all nodes, respectively. We find that the cat thalamo-cortical system incorporates features of two classic network topologies, Erdos-Renyi graphs with respect to in-absorption and in-driftness, and configuration models with respect to out-absorption and out-driftness. Moreover, taken together these four measures separate the network nodes based on broad functional roles (visual, auditory, somatomotor, and frontolimbic).

Atomic scale characterization of Mn doped InAs/GaAs quantum dots

Several growth procedures for doping InAs/GaAs quantum dots (QDs) with manganese (Mn) have been investigated with cross-sectional scanning tunneling microscopy. It is found that expulsion of Mn out of the QDs and subsequent segregation makes it difficult to incorporate Mn in the QDs even at low growth temperatures of T=320 degrees C and high Mn fluxes. Mn atoms in and around QDs have been observed with strain and potential confinement changing the appearance of the Mn features.

Ground-state functional for quantum spin chains with bond defects

We have developed a nonlocal functional of the exchange interaction for the ground-state energy of quantum spin chains described by the Heisenberg Hamiltonian. An alternating chain is used to obtain the correlation energy and a local unit-cell approximation is defined in the context of the density-functional theory. The agreement with our exact numerical data, for small chains, is significantly better than a previous formulation, even for chains with several ferromagnetic or antiferromagnetic bond defects. The results can be particularly relevant in the study of finite spin-1/2 Heisenberg chains, with exchange couplings changing, magnitude, or even sign, from bond-to-bond.

Spin gaps and spin-flip energies in density-functional theory

Energy gaps are crucial aspects of the electronic structure of finite and extended systems. Whereas much is known about how to define and calculate charge gaps in density-functional theory (DFT), and about the relation between these gaps and derivative discontinuities of the exchange-correlation functional, much less is known about spin gaps. In this paper we give density-functional definitions of spin-conserving gaps, spin-flip gaps and the spin stiffness in terms of many-body energies and in terms of single-particle (Kohn-Sham) energies. Our definitions are as analogous as possible to those commonly made in the charge case, but important differences between spin and charge gaps emerge already on the single-particle level because unlike the fundamental charge gap spin gaps involve excited-state energies. Kohn-Sham and many-body spin gaps are predicted to differ, and the difference is related to derivative discontinuities that are similar to, but distinct from, those usually considered in the case of charge gaps. Both ensemble DFT and time-dependent DFT (TDDFT) can be used to calculate these spin discontinuities from a suitable functional. We illustrate our findings by evaluating our definitions for the Lithium atom, for which we calculate spin gaps and spin discontinuities by making use of near-exact Kohn-Sham eigenvalues and, independently, from the single-pole approximation to TDDFT. The many-body corrections to the Kohn-Sham spin gaps are found to be negative, i.e., single-particle calculations tend to overestimate spin gaps while they underestimate charge gaps.

Density-functional description of superconducting and magnetic proximity effects across a tunneling barrier

A density-functional formalism for superconductivity and magnetism is presented. The resulting relations unify previously derived Kohn-Sham equations for superconductors and for noncollinear magnetism. The formalism, which discriminates Cooper-pair singlets from triplets, is applied to two quantum liquids coupled by tunneling through a barrier. An exact expression is derived, relating the eigenstates and eigenvalues of the Kohn-Sham equations, unperturbed by tunneling, on one side of the barrier to the proximity-induced ordering potential on the other.

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.

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.

Determination of the molecular weight of proteins in solution from a single small-angle X-ray scattering measurement on a relative scale

This paper describes a new and simple method to determine the molecular weight of proteins in dilute solution, with an error smaller than similar to 10%, by using the experimental data of a single small-angle X-ray scattering (SAXS) curve measured on a relative scale. This procedure does not require the measurement of SAXS intensity on an absolute scale and does not involve a comparison with another SAXS curve determined from a known standard protein. The proposed procedure can be applied to monodisperse systems of proteins in dilute solution, either in monomeric or multimeric state, and it has been successfully tested on SAXS data experimentally determined for proteins with known molecular weights. It is shown here that the molecular weights determined by this procedure deviate from the known values by less than 10% in each case and the average error for the test set of 21 proteins was 5.3%. Importantly, this method allows for an unambiguous determination of the multimeric state of proteins with known molecular weights.

Regional scale analysis of landform configuration with base-level (isobase) maps

Base-level maps (or ""isobase maps"", as originally defined by Filosofov, 1960), express a relationship between valley order and topography. The base-level map can be seen as a ""simplified"" version of the original topographic surface, from which the ""noise"" of the low-order stream erosion was removed. This method is able to identify areas with possible tectonic influence even within lithologically uniform domains. Base-level maps have been recently applied in semi-detail scale (e.g., 1:50 000 or larger) morphotectonic analysis. In this paper, we present an evaluation of the method's applicability in regional-scale analysis (e.g., 1:250 000 or smaller). A test area was selected in northern Brazil, at the lower course of the Araguaia and Tocantins rivers. The drainage network extracted from SRTM30_PLUS DEMs with spatial resolution of approximately 900 m was visually compared with available topographic maps and considered to be compatible with a 1:1,000 000 scale. Regarding the interpretation of regional-scale morphostructures, the map constructed with 2nd and 3rd-order valleys was considered to present the best results. Some of the interpreted base-level anomalies correspond to important shear zones and geological contacts present in the 1:5 000 000 Geological Map of South America. Others have no correspondence with mapped Precambrian structures and are considered to represent younger, probably neotectonic, features. A strong E-W orientation of the base-level lines over the inflexion of the Araguaia and Tocantins rivers, suggest a major drainage capture. A N-S topographic swath profile over the Tocantins and Araguaia rivers reveals a topographic pattern which, allied with seismic data showing a roughly N-S direction of extension in the area, lead us to interpret this lineament as an E-W, southward-dipping normal fault. There is also a good visual correspondence between the base-level lineaments and geophysical anomalies. A NW-SE lineament in the southeast of the study area partially corresponds to the northern border of the Mosquito lava field, of Jurassic age, and a NW-SE lineament traced in the northeastern sector of the study area can be interpreted as the Picos-Santa Ines lineament, identifiable in geophysical maps but with little expression in hypsometric or topographic maps.

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.