Ab initio study of the EFG at the N sites in imidazole

We study the nuclear quadrupole interaction at the nitrogen sites in the molecular and crystalline phases of the imidazole compound. We use PAW which is a state-of-the-art method to calculate the electronic structure and electric field gradient at the nucleus in the framework of the density functional theory. The quadrupole frequencies at both imino and amino N sites are in excellent agreement with measurements. This is the first time that the electric field gradient at crystalline imidazole is correctly treated by an ab initio theoretical approach.

Driving forces for the adsorption of cyclopentene on InP(001)

In this work the interaction of cyclopentene with a set of InP(001) surfaces is investigated by means of the density functional theory. We propose a simple approach for evaluating the surface strain and based on it we have found a linear relation between bond and strain energies and the adsorption energy. Our results also indicate that the higher the bond energy, the more disperse the charge distribution is around the adsorption site associated to the high occupied state, a key feature that characterizes the adsorption process. Different adsorption coverages are used to evaluate the proposed equation. Our results suggest that the proposed approach might be extended to other systems where the interaction of the semiconductor surface and the molecule is restricted to first neighbor sites. (C) 2011 Elsevier B.V. All rights reserved.

Semiconductor nanoparticle modeling via density functional theory

In this work, a 2.0 nm nanoparticle (low limit synthesized system) is compared to possible simplified models: passivated clusters, small (1.3 nm) nanoparticles and sets of plane surfaces. Our density functional theory results suggest that even when geometric aspects are properly described by the simplifications considered, electronic properties might be very different, especially when edge atoms are not properly taken into account in the nanoparticle`s modeling. In addition, we propose a protocol that might help future theoretical descriptions of nanoparticles.

The role of carbon impurities on the Si(001)-c(4 x 4) surface reconstruction: Theoretical calculations

In this work we employ the state-of-the-art pseudopotential method, within a generalized gradient approximation to the density functional theory, combined with a recently developed method for the calculation of HREELS spectra to study a series of different proposed models for carbon incorporation on the silicon (001) surface. A fully discussion on the geometry, energetics and specially the comparison between experimental and theoretical STM images and electron energy loss spectra indicate that the Si(100)-c(4 x 4) is probably induced by Si-C surface dinners, in agreement with recent experimental findings. (C) 2009 Elsevier B.V. All rights reserved.

Furan interaction with the Si(001)-(2 x 2) surface: structural, energetics, and vibrational spectra from first-principles

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 furan on the silicon (001) surface. A direct comparison of different adsorption structures with x-ray photoelectron spectroscopy (XPS), ultra-violet photoelectron spectroscopy (UPS), high resolution electron energy loss spectroscopy (HREELS), near edge x-ray absorption fine structure (NEXAFS), and high resolution spectroscopy experimental data allows us to identify the [4 + 2] cycloaddition reaction as the most probable adsorbate. In addition, theoretical scanning tunnelling microscopy (STM) images are presented, with a view to contributing to further experimental investigations.

Proximal methods for nonlinear programming: double regularization and inexact subproblems

This paper describes the first phase of a project attempting to construct an efficient general-purpose nonlinear optimizer using an augmented Lagrangian outer loop with a relative error criterion, and an inner loop employing a state-of-the art conjugate gradient solver. The outer loop can also employ double regularized proximal kernels, a fairly recent theoretical development that leads to fully smooth subproblems. We first enhance the existing theory to show that our approach is globally convergent in both the primal and dual spaces when applied to convex problems. We then present an extensive computational evaluation using the CUTE test set, showing that some aspects of our approach are promising, but some are not. These conclusions in turn lead to additional computational experiments suggesting where to next focus our theoretical and computational efforts.

On Approximate KKT Condition and its Extension to Continuous Variational Inequalities

In this work, we introduce a necessary sequential Approximate-Karush-Kuhn-Tucker (AKKT) condition for a point to be a solution of a continuous variational inequality, and we prove its relation with the Approximate Gradient Projection condition (AGP) of Garciga-Otero and Svaiter. We also prove that a slight variation of the AKKT condition is sufficient for a convex problem, either for variational inequalities or optimization. Sequential necessary conditions are more suitable to iterative methods than usual punctual conditions relying on constraint qualifications. The AKKT property holds at a solution independently of the fulfillment of a constraint qualification, but when a weak one holds, we can guarantee the validity of the KKT conditions.

A relaxed constant positive linear dependence constraint qualification and applications

In this work we introduce a relaxed version of the constant positive linear dependence constraint qualification (CPLD) that we call RCPLD. This development is inspired by a recent generalization of the constant rank constraint qualification by Minchenko and Stakhovski that was called RCRCQ. We show that RCPLD is enough to ensure the convergence of an augmented Lagrangian algorithm and that it asserts the validity of an error bound. We also provide proofs and counter-examples that show the relations of RCRCQ and RCPLD with other known constraint qualifications. In particular, RCPLD is strictly weaker than CPLD and RCRCQ, while still stronger than Abadie's constraint qualification. We also verify that the second order necessary optimality condition holds under RCRCQ.

Hybrid density functional study of small Rh-n (n=2-15) clusters

The physical properties of small rhodium clusters, Rh-n, have been in debate due to the shortcomings of density functional theory (DFT). To help in the solution of those problems, we obtained a set of putative lowest energy structures for small Rh-n (n = 2-15) clusters employing hybrid-DFT and the generalized gradient approximation (GGA). For n = 2-6, both hybrid and GGA functionals yield similar ground-state structures (compact), however, hybrid favors compact structures for n = 7-15, while GGA favors open structures based on simple cubic motifs. Thus, experimental results are crucial to indicate the correct ground-state structures, however, we found that a unique set of structures (compact or open) is unable to explain all available experimental data. For example, the GGA structures (open) yield total magnetic moments in excellent agreement with experimental data, while hybrid structures (compact) have larger magnetic moments compared with experiments due to the increased localization of the 4d states. Thus, we would conclude that GGA provides a better description of the Rh-n clusters, however, a recent experimental-theoretical study [ Harding et al., J. Chem. Phys. 133, 214304 (2010)] found that only compact structures are able to explain experimental vibrational data, while open structures cannot. Therefore, it indicates that the study of Rh-n clusters is a challenging problem and further experimental studies are required to help in the solution of this conundrum, as well as a better description of the exchange and correlation effects on the Rh n clusters using theoretical methods such as the quantum Monte Carlo method.

Site localization of Cd impurities in sapphire

By combining first-principles electronic structure calculations and existing time-differential gamma-gamma perturbed-angular-correlation experiments we studied the site localization, the local environment, and the electronic structure of Cd impurities in sapphire (alpha-Al2O3) single crystals in different charged states. The ab initio calculations were performed with the full-potential augmented plane wave plus local orbitals method and the projector augmented wave method. Comparing the calculated electric-field-gradient tensor at the Cd nuclei in the alpha-Al2O3 host lattice and the corresponding available experimental values, we have seen that it is equally possible for Cd to replace an Al atom (in a negative charge state) or to be placed in an interstitial site (in a neutral charge state). To finally address the issue of the Cd impurity localization, we performed formation energy calculations. These results have shown that Cd placed in the substitutional Al site, in the negatively charged state, is the most probable configuration.

Evaluating bound-constrained minimization software

Bound-constrained minimization is a subject of active research. To assess the performance of existent solvers, numerical evaluations and comparisons are carried on. Arbitrary decisions that may have a crucial effect on the conclusions of numerical experiments are highlighted in the present work. As a result, a detailed evaluation based on performance profiles is applied to the comparison of bound-constrained minimization solvers. Extensive numerical results are presented and analyzed.

Quantum Monte Carlo study of small aluminum clusters Al-n (n=2-13)

Using fixed node diffusion quantum Monte Carlo (FN-DMC) simulations and density functional theory (DFT) within the generalized gradient approximations, we calculate the total energies of the relaxed and unrelaxed neutral, cationic, and anionic aluminum clusters, Al-n (n = 1-13). From the obtained total energies, we extract the ionization potential and electron detachment energy and compare with previous theoretical and experimental results. Our results for the electronic properties from both the FN-DMC and DFT calculations are in reasonably good agreement with the available experimental data. A comparison between the FN-DMC and DFT results reveals that their differences are a few tenths of electron volt for both the ionization potential and the electron detachment energy. We also observe two distinct behaviors in the electron correlation contribution to the total energies from smaller to larger clusters, which could be assigned to the structural transition of the clusters from planar to three-dimensional occurring at n = 4 to 5.

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.