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.

K*(0) production in Cu plus Cu and Au plus Au collisions at root s(NN)=62.4 GeV and 200 GeV

We report on K*(0) production at midrapidity in Au + Au and Cu + Cu collisions at root s(NN) = 62.4 and 200 GeV collected by the Solenoid Tracker at the Relativistic Heavy Ion Collider detector. The K*(0) is reconstructed via the hadronic decays K*(0) -> K(+)pi(-) and (K*(0)) over bar -> K(+)pi(-). Transverse momentum, p(T), spectra are measured over a range of p(T) extending from 0.2 GeV/c up to 5 GeV/c. The center-of-mass energy and system size dependence of the rapidity density, dN/dy, and the average transverse momentum, < p(T)>, are presented. The measured N(K*(0))/N(K) and N(phi)/N(K*(0)) ratios favor the dominance of rescattering of decay daughters of K*(0) over the hadronic regeneration for the K*(0) production. In the intermediate p(T) region (2.0 < p(T) < 4.0 GeV/c), the elliptic flow parameter, v(2), and the nuclear modification factor, R(CP), agree with the expectations from the quark coalescence model of particle production.

Charged and strange hadron elliptic flow in Cu plus Cu collisions at root s(NN)=62.4 and 200 GeV

We present the results of an elliptic flow, v(2), analysis of Cu + Cu collisions recorded with the solenoidal tracker detector (STAR) at the BNL Relativistic Heavy Ion Collider at root s(NN) = 62.4 and 200 GeV. Elliptic flow as a function of transverse momentum, v(2)(p(T)), is reported for different collision centralities for charged hadrons h(+/-) and strangeness-ontaining hadrons K(S)(0), Lambda, Xi, and phi in the midrapidity region vertical bar eta vertical bar < 1.0. Significant reduction in systematic uncertainty of the measurement due to nonflow effects has been achieved by correlating particles at midrapidity, vertical bar eta vertical bar < 1.0, with those at forward rapidity, 2.5 < vertical bar eta vertical bar < 4.0. We also present azimuthal correlations in p + p collisions at root s = 200 GeV to help in estimating nonflow effects. To study the system-size dependence of elliptic flow, we present a detailed comparison with previously published results from Au + Au collisions at root s(NN) = 200 GeV. We observe that v(2)(p(T)) of strange hadrons has similar scaling properties as were first observed in Au + Au collisions, that is, (i) at low transverse momenta, p(T) < 2 GeV/c, v(2) scales with transverse kinetic energy, m(T) - m, and (ii) at intermediate p(T), 2 < p(T) < 4 GeV/c, it scales with the number of constituent quarks, n(q.) We have found that ideal hydrodynamic calculations fail to reproduce the centrality dependence of v(2)(p(T)) for K(S)(0) and Lambda. Eccentricity scaled v(2) values, v(2)/epsilon, are larger in more central collisions, suggesting stronger collective flow develops in more central collisions. The comparison with Au + Au collisions, which go further in density, shows that v(2)/epsilon depends on the system size, that is, the number of participants N(part). This indicates that the ideal hydrodynamic limit is not reached in Cu + Cu collisions, presumably because the assumption of thermalization is not attained.

Identified particle production, azimuthal anisotropy, and interferometry measurements in Au plus Au collisions at root s(NN)=9.2 GeV

We present the first measurements of identified hadron production, azimuthal anisotropy, and pion interferometry from Au + Au collisions below the nominal injection energy at the BNL Relativistic Heavy-Ion Collider (RHIC) facility. The data were collected using the large acceptance solenoidal tracker at RHIC (STAR) detector at root s(NN) = 9.2 GeV from a test run of the collider in the year 2008. Midrapidity results on multiplicity density dN/dy in rapidity y, average transverse momentum < p(T)>, particle ratios, elliptic flow, and Hanbury-Brown-Twiss (HBT) radii are consistent with the corresponding results at similar root s(NN) from fixed-target experiments. Directed flow measurements are presented for both midrapidity and forward-rapidity regions. Furthermore the collision centrality dependence of identified particle dN/dy, < p(T)>, and particle ratios are discussed. These results also demonstrate that the capabilities of the STAR detector, although optimized for root s(NN) = 200 GeV, are suitable for the proposed QCD critical-point search and exploration of the QCD phase diagram at RHIC.

Sensitivity to noise and ergodicity of an assembly line of cellular automata that classifies density

We investigate the sensitivity of the composite cellular automaton of H. Fuks [Phys. Rev. E 55, R2081 (1997)] to noise and assess the density classification performance of the resulting probabilistic cellular automaton (PCA) numerically. We conclude that the composite PCA performs the density classification task reliably only up to very small levels of noise. In particular, it cannot outperform the noisy Gacs-Kurdyumov-Levin automaton, an imperfect classifier, for any level of noise. While the original composite CA is nonergodic, analyses of relaxation times indicate that its noisy version is an ergodic automaton, with the relaxation times decaying algebraically over an extended range of parameters with an exponent very close (possibly equal) to the mean-field value.

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.

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.

Excitation energies from ground-state density-functionals by means of generator coordinates

The generator-coordinate method is a flexible and powerful reformulation of the variational principle. Here we show that by introducing a generator coordinate in the Kohn-Sham equation of density-functional theory, excitation energies can be obtained from ground-state density functionals. As a viability test, the method is applied to ground-state energies and various types of excited-state energies of atoms and ions from the He and the Li isoelectronic series. Results are compared to a variety of alternative DFT-based approaches to excited states, in particular time-dependent density-functional theory with exact and approximate potentials.

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.

Near infrared spectroscopy for soil bulk density assessment

P>Soil bulk density values are needed to convert organic carbon content to mass of organic carbon per unit area. However, field sampling and measurement of soil bulk density are labour-intensive, costly and tedious. Near-infrared reflectance spectroscopy (NIRS) is a physically non-destructive, rapid, reproducible and low-cost method that characterizes materials according to their reflectance in the near-infrared spectral region. The aim of this paper was to investigate the ability of NIRS to predict soil bulk density and to compare its performance with published pedotransfer functions. The study was carried out on a dataset of 1184 soil samples originating from a reforestation area in the Brazilian Amazon basin, and conventional soil bulk density values were obtained with metallic ""core cylinders"". The results indicate that the modified partial least squares regression used on spectral data is an alternative method for soil bulk density predictions to the published pedotransfer functions tested in this study. The NIRS method presented the closest-to-zero accuracy error (-0.002 g cm-3) and the lowest prediction error (0.13 g cm-3) and the coefficient of variation of the validation sets ranged from 8.1 to 8.9% of the mean reference values. Nevertheless, further research is required to assess the limits and specificities of the NIRS method, but it may have advantages for soil bulk density predictions, especially in environments such as the Amazon forest.

Evaluation of grinding methods for pellets preparation aiming at the analysis of plant materials by laser induced breakdown spectrometry

It has been demonstrated that laser induced breakdown spectrometry (LIBS) can be used as an alternative method for the determination of macro (P, K. Ca, Mg) and micronutrients (B, Fe, Cu, Mn, Zn) in pellets of plant materials. However, information is required regarding the sample preparation for plant analysis by LIBS. In this work, methods involving cryogenic grinding and planetary ball milling were evaluated for leaves comminution before pellets preparation. The particle sizes were associated to chemical sample properties such as fiber and cellulose contents, as well as to pellets porosity and density. The pellets were ablated at 30 different sites by applying 25 laser pulses per site (Nd:YAG@1064 nm, 5 ns, 10 Hz, 25J cm(-2)). The plasma emission collected by lenses was directed through an optical fiber towards a high resolution echelle spectrometer equipped with an ICCD. Delay time and integration time gate were fixed at 2.0 and 4.5 mu s, respectively. Experiments carried out with pellets of sugarcane, orange tree and soy leaves showed a significant effect of the plant species for choosing the most appropriate grinding conditions. By using ball milling with agate materials, 20 min grinding for orange tree and soy, and 60 min for sugarcane leaves led to particle size distributions generally lower than 75 mu m. Cryogenic grinding yielded similar particle size distributions after 10 min for orange tree, 20 min for soy and 30 min for sugarcane leaves. There was up to 50% emission signal enhancement on LIBS measurements for most elements by improving particle size distribution and consequently the pellet porosity. (C) 2011 Elsevier B.V. All rights reserved.