Encapsulation of small magnetic clusters in fullerene cages: A density functional theory investigation within van der Waals corrections

The encapsulation of magnetic transition-metal (TM) clusters inside carbon cages (fullerenes, nanotubes) has been of great interest due to the wide range of applications, which spread from medical sensors in magnetic resonance imaging to photonic crystals. Several theoretical studies have been reported; however, our atomistic understanding of the physical properties of encapsulated magnetic TM 3d clusters is far from satisfactory. In this work, we will report general trends, derived from density functional theory within the generalized gradient approximation proposed by Perdew, Burke, and Ernzerhof (PBE), for the encapsulation properties of the TMm@C-n (TM = Fe, Co, Ni; m = 2-6, n = 60,70,80,90) systems. Furthermore, to understand the role of the van der Waals corrections to the physical properties, we employed the empirical Grimme's correction (PBE + D2). We found that both PBE and PBE + D2 functionals yield almost the same geometric parameters, magnetic and electronic properties, however, PBE + D2 strongly enhances the encapsulation energy. We found that the center of mass of the TMm clusters is displaced towards the inside C-n surfaces, except for large TMm clusters (m = 5 and 6). For few cases, e. g., Co-4 and Fe-4, the encapsulation changes the putative lowest-energy structure compared to the isolated TMm clusters. We identified few physical parameters that play an important role in the sign and magnitude of the encapsulation energy, namely, cluster size, fullerene equatorial diameter, shape, curvature of the inside C-n surface, number of TM atoms that bind directly to the inside C-n surface, and the van der Waals correction. The total magnetic moment of encapsulated TMm clusters decreases compared with the isolated TMm clusters, which is expected due to the hybridization of the d-p states, and strongly depends on the size and shape of the fullerene cages.

Ethanol and Water Adsorption on Close-Packed 3d, 4d, and 5d Transition-Metal Surfaces: A Density Functional Theory Investigation with van der Waals Correction

Nowadays, there is a great interest in the economic success of direct ethanol fuel cells; however, our atomistic understanding of the designing of stable and low-cost catalysts for the steam reforming of ethanol is still far from satisfactory, in particular due to the large number of undesirable intermediates. In this study, we will report a first-principles investigation of the adsorption properties of ethanol and water at low coverage on close-packed transition-metal (TM) surfaces, namely, Fe(110), Co(0001), Ni(111), Cu(111), Ru(0001), Rh(111), Pd(111), Ag(111), Os(0001), Ir(111), Pt(111), and Au(111), employing density functional theory (DFT) calculations. We employed the generalized gradient approximation with the formulation proposed by Perdew, Burke, and Erzenholf (PBE) to the exchange correlation functional and the empirical correction proposed by S. Grimme (DFT+D3) for the van der Waals correction. We found that both adsorbates binds preferentially near or on the on top sites of the TM surfaces through the 0 atoms. The PBE adsorption energies of ethanol and water decreases almost linearly with the increased occupation of the 4d and 5d d-band, while there is a deviation for the 3d systems. The van der Waals correction affects the linear behavior and increases the adsorption energy for both adsorbates, which is expected as the van der Waals energy due to the correlation effects is strongly underestimated by DFT-PBE for weak interacting systems. The geometric parameters for water/TM are not affected by the van der Waals correction, i.e., both DFT and DFT+D3 yield an almost parallel orientation for water on the TM surfaces; however, DFT+D3 changes drastically the ethanol orientation. For example, DFT yields an almost perpendicular orientation of the C-C bond to the TM surface, while the C-C bond is almost parallel to the surface using DFT +D3 for all systems, except for ethanol/Fe(110). Thus, the van der Waals correction decreases the distance of the C atoms to the TM surfaces, which might contribute to break the C-C bond. The work function decreases upon the adsorption of ethanol and water, and both follow the same trends, however, with different magnitude (larger for ethanol/TM) due to the weak binding of water to the surface. The electron density increases mainly in the region between the topmost layer and the adsorbates, which explains the reduction of the substrate work function.

Adsorption of NO on the Rh-13, Pd-13, Ir-13, and Pt-13 Clusters: A Density Functional Theory Investigation

The adsorption of NO on transition-metal (TM) surfaces has been widely studied by experimental and theoretical techniques; however, our atomistic understanding of the interaction of nitrogen monoxide (NO) with small TM clusters is far from satisfactory, which compromises a deep understanding of real catalyst devices. In this study, we report a density functional theory study of the adsorption properties of NO on the TM13 (TM = Rh, Pd, Ir, Pt) clusters employing the projected augmented wave method. We found that the interaction of NO with TM13 is much more complex than that for NO/TM(111). In particular, for low symmetry TM13 clusters, there is a strong rearrangement of the electronic charge density upon NO adsorption and, as a consequence, the adsorption energy shows a very complex dependence even for adsorption sites with the same local effective coordination. We found a strong enhancement of the binding energy of NO to the TM13 clusters compared with the TM(111) surfaces, as the antibonding NO states are not occupied for NO/TM13, and the general relationship based on the d-band model between adsorption energy and the center of gravity of the occupied d-states does not hold for the studied TM13 clusters, in particular, for clusters with low symmetry. In contrast with the adsorption energy trends, the geometric NO/TM13 parameters and the vibrational N-O frequencies for different coordination sites follow the same trend as for the respective TM(111) surfaces, while the changes in the frequencies between different surfaces and TM13 clusters reflect the strong NO-TM13 interaction.

Sulfur Radicals as Tethers for the Adsorption of Aromatic Molecules on Silicon Surface

In this work we employ the state of the art pseudopotential method, within a generalized gradient approximation to the density functional theory, to investigate the adsorption process of benzenethiol and diphenyl disulfide with the silicon (001) surface. A direct comparison of different adsorption structures with Fourier transform infrared spectroscopy (FTIR), X-ray photoelectron spectroscopy (XPS) allow us to identify that benzenethiol and diphenyl disulfide dissociatively adsorb on the silicon surface. In addition, theoretically obtained data suggests that the C6H5SH:Si(001) presents a higher Schottky barrier height contact when compared to other similar aromatic molecules.

The T-2 phase in the Nb-Si-B system studied by ab initio calculations and synchrotron X-ray diffraction

The solid solution based on Nb5Si3 (Cr5B3 structure type, D8(l), tl32, 14/mcm, No140, a=6.5767 angstrom, c=11.8967 angstrom) in the Nb-Si-B system was studied from the structural and thermodynamic point of view both experimentally and by ab initio calculations. Rietveld refinement of powder X-ray synchrotron data allowed to determine the boron to silicon substitution mechanism and the structural parameters. Ab initio calculations of different ordered compounds and selected disordered alloys allowed to obtain in addition to the enthalpy of formation of the solution, substitution mechanism and structural parameters which are in excellent agreement with the experimental data. The stability of the phase is discussed. (C) 2012 Elsevier Inc. All rights reserved.

Platinum-Based Nanoalloys PtnTM55-n (TM=Co, Rh, Au): A Density Functional Theory Investigation

Structural and electronic properties of the PtnTM55-n (TM = Co, Rh, Au) nanoalloys are investigated using density functional theory within the generalized gradient approximation and employing the all-electron projected augmented wave method. For TM = Co and Rh, the excess energy, which measures the relative energy stability of the nanoalloys, is negative for all Pt compositions. We found that the excess energy has similar values for a wide range of Pt compositions, i.e., n = 20-42 and n = 28-42 for Co and Rh, respectively, with the core shell icosahedron-like configuration (n = 42) being slightly more stable for both Co and Rh systems because of the larger release of the strain energy due to the smaller atomic size of the Co and Rh atoms. For TM = Au, the excess energy is positive for all compositions, except for n = 13, which is energetically favorable due to the formation of the core-shell structure (Pt in the core and Au atoms at the surface). Thus, our calculations confirm that the formation of core-shell structures plays an important role to increase the stability of nanoalloys. The center of gravity of the occupied d-states changes almost linearly as a function of the Pt composition, and hence, based on the d-band model, the magnitude of the adsorption energy of an adsorbate can be tuned by changing the Pt composition. The magnetic moments of PtnCo55-n decrease almost linearly as a function of the Pt composition; however, the same does not hold for PtRh and PtAu. We found an enhancement of the magnetic moments of PtRh by a few times by increasing Pt composition, which we explain by the compression effects induced by the large size of the Pt atoms compared with the Rh atoms.

Electronic structure and transport properties of atomic NiO spinvalves

Ab initio quantum transport calculations show that short NiO chains suspended in Ni nanocontacts present a very strong spin-polarization of the conductance.The generalized gradient approximation we use here predicts a similar polarization of the conductance as the one previously computed with non-local exchange, confirming the robustness of the result. Their use as nanoscopic spinvalves is proposed.

Derivation of the spin Hamiltonians for Fe in MgO

A method to calculate the effective spin Hamiltonian for a transition metal impurity in a non-magnetic insulating host is presented and applied to the paradigmatic case of Fe in MgO. In the first step we calculate the electronic structure employing standard density functional theory (DFT), based on generalized gradient approximation (GGA), using plane waves as a basis set. The corresponding basis of atomic-like maximally localized Wannier functions is derived and used to represent the DFT Hamiltonian, resulting in a tight-binding model for the atomic orbitals of the magnetic impurity. The third step is to solve, by exact numerical diagonalization, the N electron problem in the open shell of the magnetic atom, including both effects of spin–orbit and Coulomb repulsion. Finally, the low energy sector of this multi-electron Hamiltonian is mapped into effective spin models that, in addition to the spin matrices S, can also include the orbital angular momentum L when appropriate. We successfully apply the method to Fe in MgO, considering both the undistorted and Jahn–Teller (JT) distorted cases. Implications for the influence of Fe impurities on the performance of magnetic tunnel junctions based on MgO are discussed.

Strain-induced structural instability in FeRh

We perform density functional calculations to investigate the structure of the intermetallic alloy FeRh under epitaxial strain. Bulk FeRh exhibits a metamagnetic transition from a low-temperature antiferromagnetic (AFM) phase to a ferromagnetic phase at 350 K, and its strain dependence is of interest for tuning the transition temperature to the room-temperature operating conditions of typical memory devices. We find an unusually strong dependence of the structural energetics on the choice of exchange-correlation functional, with the usual local density approximation yielding the wrong ground-state structure, and generalized gradient (GGA) extensions being in better agreement with the bulk experimental structure. Using the GGA we show the existence of a metastable face-centered-cubic-like AFM structure that is reached from the ground-state body-centered-cubic-like AFM structure by following the epitaxial Bain path. We show that the behavior is well described using nonlinear elasticity theory, which captures the softening and eventual sign change of the orthorhombic shear modulus under compressive strain, consistent with this structural instability. Finally, we predict the existence of an additional unit-cell-doubling lattice instability, which should be observable at low temperature.

Estudo in silico das interações da protease NS3-NS2B de DENV-2 com o inibidor peptídico Bz-nKRR-H com finalidades terapêuticas

Dengue virus is an important patogen that causes Dengue desease in all world, and belongs to Flavivirus gender. The virus consists of enveloped RNA with a single strand positive sense, 11Kb genome. The RNA is translated into a polyprotein precursor, wich is cleaved into 3 structural proteins (C, prM e E) and 7 non-structural proteins (NS1, NS2A, NS2B, NS3, NS4A, NS4B e NS5). The NS3 is a multifunctional protein, that besides to promote the polyprotein precursor cleavage, also have NTPase, helicase and RTPase activity. The NS3 needs a hydrophilic segment of 40 residues from the transmembrane NS2B protein (who acts like cofator) to realize this functions. Actually, there's no vacines available on the market, and the treatment are just symptomatic. The tetrapeptide inhibitor Bz-Nle-Lys-Arg-Arg-H (Ki de 5,8-7,0 M) was showed as a potent inhibitor μ for NS3prot in Dengue virus. That is a inteligent alternative to treat the dengue desease. The present work aimed analyse the interactions of the ligand bounded to the activity site to provid a clear and depth vision of that interaction. For this purpouse, it was conducted an in silico study, by using quantum mechanical calculations based on Density Functional Theory (DFT), with Generalized Gradient approximation (GGA) to describe the effects of exchange and correlation. The interaction energy of each amino acid belonging to the binding site to the ligand was calculated the using the method of molecular fragmentation with conjugated caps (MFCC). Besides energy, we calculated the distances, types of molecular interactions and atomic groups involved. The theoretical models used were satisfactory and show a more accurate description when the dielectric constant = 20 ε and 80 was used. The results demonstrate that the interaction energy of the system reached convergence at 13.5 A. Within a radius of 13,5A the most important residues were identified. Met49, Met84 and Asp81 perform interactions of hydrogen with the ligant. The Asp79 and Asp75 residues present high energy of attraction. Arg54, Arg85 and Lys 131 perform hydrogen interactions with the ligand, however, appear in BIRD graph having high repulsion energy with the inhibitor. The data also emphasizes the importance of residue Tyr161 and the involvement of the catalytic triad composed by Asp75, His51 and Ser135

Low valence cation doping of bulk Cr2O3: Charge compensation and oxygen vacancy formation

The different oxidation states of chromium allow its bulk oxide form to be reducible, facilitating the oxygen vacancy formation process, which is a key property in applications such as catalysis. Similar to other useful oxides such as TiO2, and CeO2, the effect of substitutional metal dopants in bulk Cr2O3 and its effect on the electronic structure and oxygen vacancy formation are of interest, particularly in enhancing the latter. In this paper, density functional theory (DFT) calculations with a Hubbard + U correction (DFT+U) applied to the Cr 3d and O 2p states, are carried out on pure and metal-doped bulk Cr2O3 to examine the effect of doping on the electronic and geometric structure. The role of dopants in enhancing the reducibility of Cr2O3 is examined to promote oxygen vacancy formation. The dopants are Mg, Cu, Ni, and Zn, which have a formal +2 oxidation state in their bulk oxides. Given this difference in host and, dopant oxidation states, we show that to predict the correct ground state two metal dopants charge compensated with an oxygen vacancy are required. The second oxygen atom removed is termed "the active" oxygen vacancy and it is the energy required to remove this atom that is related to the reduction process. In all cases, we find that substitutional doping improves the oxygen vacancy formation of bulk Cr2O3 by lowering the energy cost.

First-Principles Calculations of the Electronic and Structural Properties of GaSb

In this paper, we carried out first-principles calculations in order to investigate the structural and electronic properties of the binary compound gallium antimonide (GaSb). This theoretical study was carried out using the Density Functional Theory within the plane-wave pseudopotential method. The effects ofexchange and correlation (XC) were treated using the functional Local Density Approximation (LDA), generalized gradient approximation (GGA): Perdew–Burke–Ernzerhof (PBE), Perdew-Burke-Ernzerhof revised for solids (PBEsol), Perdew-Wang91 (PW91), revised Perdew–Burke–Ernzerhof (rPBE), Armiento–Mattson 2005 (AM05) and meta-generalized gradient approximation (meta-GGA): Tao–Perdew– Staroverov–Scuseria (TPSS) and revised Tao–Perdew–Staroverov–Scuseria (RTPSS) and modified Becke-Johnson (MBJ). We calculated the densities of state (DOS) and band structure with different XC potentials identified and compared them with the theoretical and experimental results reported in the literature. It was discovered that functional: LDA, PBEsol, AM05 and RTPSS provide the best results to calculate the lattice parameters (a) and bulk modulus (B0); while for the cohesive energy (Ecoh), functional: AM05, RTPSS and PW91 are closer to the values obtained experimentally. The MBJ, Rtpss and AM05 values found for the band gap energy is slightly underestimated with those values reported experimentally.

The generalized Onsager model for a binary gas mixture

The Onsager model for the secondary flow field in a high-speed rotating cylinder is extended to incorporate the difference in mass of the two species in a binary gas mixture. The base flow is an isothermal solid-body rotation in which there is a balance between the radial pressure gradient and the centrifugal force density for each species. Explicit expressions for the radial variation of the pressure, mass/mole fractions, and from these the radial variation of the viscosity, thermal conductivity and diffusion coefficient, are derived, and these are used in the computation of the secondary flow. For the secondary flow, the mass, momentum and energy equations in axisymmetric coordinates are expanded in an asymptotic series in a parameter epsilon = (Delta m/m(av)), where Delta m is the difference in the molecular masses of the two species, and the average molecular mass m(av) is defined as m(av) = (rho(w1)m(1) + rho(w2)m(2))/rho(w), where rho(w1) and rho(w2) are the mass densities of the two species at the wall, and rho(w) = rho(w1) + rho(w2). The equation for the master potential and the boundary conditions are derived correct to O(epsilon(2)). The leading-order equation for the master potential contains a self-adjoint sixth-order operator in the radial direction, which is different from the generalized Onsager model (Pradhan & Kumaran, J. Fluid Mech., vol. 686, 2011, pp. 109-159), since the species mass difference is included in the computation of the density, viscosity and thermal conductivity in the base state. This is solved, subject to boundary conditions, to obtain the leading approximation for the secondary flow, followed by a solution of the diffusion equation for the leading correction to the species mole fractions. The O(epsilon) and O(epsilon(2)) equations contain inhomogeneous terms that depend on the lower-order solutions, and these are solved in a hierarchical manner to obtain the O(epsilon) and O(epsilon(2)) corrections to the master potential. A similar hierarchical procedure is used for the Carrier-Maslen model for the end-cap secondary flow. The results of the Onsager hierarchy, up to O(epsilon(2)), are compared with the results of direct simulation Monte Carlo simulations for a binary hard-sphere gas mixture for secondary flow due to a wall temperature gradient, inflow/outflow of gas along the axis, as well as mass and momentum sources in the flow. There is excellent agreement between the solutions for the secondary flow correct to O(epsilon(2)) and the simulations, to within 15 %, even at a Reynolds number as low as 100, and length/diameter ratio as low as 2, for a low stratification parameter A of 0.707, and when the secondary flow velocity is as high as 0.2 times the maximum base flow velocity, and the ratio 2 Delta m/(m(1) + m(2)) is as high as 0.5. Here, the Reynolds number Re = rho(w)Omega R-2/mu, the stratification parameter A = root m Omega R-2(2)/(2k(B)T), R and Omega are the cylinder radius and angular velocity, m is the molecular mass, rho(w) is the wall density, mu is the viscosity and T is the temperature. The leading-order solutions do capture the qualitative trends, but are not in quantitative agreement.

Generalized power method for sparse principal component analysis

In this paper we develop a new approach to sparse principal component analysis (sparse PCA). We propose two single-unit and two block optimization formulations of the sparse PCA problem, aimed at extracting a single sparse dominant principal component of a data matrix, or more components at once, respectively. While the initial formulations involve nonconvex functions, and are therefore computationally intractable, we rewrite them into the form of an optimization program involving maximization of a convex function on a compact set. The dimension of the search space is decreased enormously if the data matrix has many more columns (variables) than rows. We then propose and analyze a simple gradient method suited for the task. It appears that our algorithm has best convergence properties in the case when either the objective function or the feasible set are strongly convex, which is the case with our single-unit formulations and can be enforced in the block case. Finally, we demonstrate numerically on a set of random and gene expression test problems that our approach outperforms existing algorithms both in quality of the obtained solution and in computational speed. © 2010 Michel Journée, Yurii Nesterov, Peter Richtárik and Rodolphe Sepulchre.