Abstract: Carbon nitride nanospheres have been synthesized into copper by simultaneous high fluence (10(18) at. cm(-2)) implantations of C-12 and N-15 ions. The composition of the implanted region has been measured using C-12(d,p(0))C-13 and N-15(d,alpha(0))C-13 nuclear reactions induced by a 1.05 MeV deuteron beam. The C-12 and N-15 depth profiles are very close and the retained doses into copper are relatively high, which indicates that carbon and nitrogen diffusion processes are likely limited during implantation. High resolution transmission electron microscopy (HRTEM) observations and electron diffraction (ED) analyses have been carried out to determine the structure of the nanospheres formed during implantation. Some consist in small hollow amorphous nanocapsules with sizes ranging from 30 to 100 nm. Large gas bubbles with diameters up to 300 mn have also been observed in the copper matrix. Electron energy-loss spectroscopy (EELS) measurements performed on the small nanocapsules indicate that their shells are composed of carbon and nitrogen. (c) 2007 Elsevier B.V. All rights reserved.

Abstract: We investigate the role of transition metal atoms of group V-b (V, Nb, Ta) and VI-b (Cr, Mo, W) as n- or p-type dopants in anatase TiO$_2$ using thermodynamic
principles and density functional theory with the Heyd-Scuseria-Ernzerhof HSE06 hybrid functional. The HSE06 functional provides a realistic value for the band gap, which ensures a correct classification of dopants as shallow or deep donors or acceptors. Defect formation energies and thermodynamic transition levels are calculated taking into account the constraints imposed by the stability of TiO$_2$ and the solubility limit of the impurities.
Nb, Ta, W and Mo are identified as shallow donors. Although W provides two electrons, Nb and Ta show a considerable lower formation energy, in particular under O-poor conditions. Mo donates in principle one electron, but under specific conditions can turn into a double donor. V impurities are deep donors and Cr
shows up as an amphoteric defect, thereby acting as an electron trapping center in n-type TiO$_2$ especially under O-rich conditions. A comparison with the available experimental data yields excellent agreement.

Abstract: In this study, we investigated the diffusion of H-atoms to the subsurface and their further diffusion into the bulk of a Ni(111) crystal by means of density functional theory calculations in the context of thermal and plasma-assisted catalysis. The H-atoms at the surface can originate from the dissociative adsorption of H2 or CH4 molecules, determining the surface H-coverage. When a threshold H-coverage is passed, corresponding to 1.00 ML for the crystalline Ni(111) surface, the surface-bound H-atoms start to diffuse to the subsurface. A similar threshold coverage is observed for the interstitial H-coverage. Once the interstitial sites are filled up with a coverage above 1.00 ML of H, dissolution of interstitial H-atoms to the layer below the interstitial sites will be initiated. Hence, by applying a high pressure or inducing a reactive plasma and high temperature, increasing the H-flux to the surface, a large amount of hydrogen can diffuse in a crystalline metal like Ni and can be absorbed. The formation of metal hydride may modify the entire reaction kinetics of the system. Equivalently, the H-atoms in the bulk can easily go back to the surface and release a large amount of heat. In a plasma process, H-atoms are formed in the plasma, and therefore the energy barrier for dissociative adsorption is dismissed, thus allowing achievement of the threshold coverage without applying a high pressure as in a thermal process. As a result, depending on the crystal plane and type of metal, a large number of H-atoms can be dissolved (absorbed) in the metal catalyst, explaining the high efficiency of plasma-assisted catalytic reactions. Here, the mechanism of H-dissolution is established as a chemical identifier for the investigation of the reaction kinetics of a chemical process.

Abstract: Using fixed-node diffusion quantum Monte Carlo (FN-DMC) simulation we investigate the electron correlation in all-metal aromatic clusters MAl4- (with M = Li, Na, K, Rb, Cu, Ag and Au). The electron detachment energies and electron affinities of the clusters are obtained. The vertical electron detachment energies obtained from the FN-DMC calculations are in very good agreement with the available experimental results. Calculations are also performed within the Hartree-Fock approximation, density-functional theory (DFT), and the couple-cluster (CCSD(T)) method. From the obtained results, we analyse the impact of the electron correlation effects in these bimetallic clusters and find that the correlation of the valence electrons contributes significantly to the detachment energies and electron affinities, varying between 20% and 50% of their total values. Furthermore, we discuss the electron correlation effects on the stability of the clusters as well as the accuracy of the DFT and CCSD(T) calculations in the present systems.

Abstract: Chalcopyrite semiconductors are of considerable interest for application as absorber layers in thin-film photovoltaic cells. When growing films of these compounds, however, they are often found to contain CuAu-like domains, a metastable phase of chalcopyrite. It has been reported that for CuInS2, the presence of the CuAu-like phase improves the short circuit current of the chalcopyrite-based photovoltaic cell. We investigate the thermodynamic stability of both phases for a selected list of I-III-VI2 materials using a first-principles density functional theory approach. For the CuIn-VI2 compounds, the difference in formation energy between the chalcopyrite and CuAu-like phase is found to be close to 2 meV per atom, indicating a high likelihood of the presence of CuAu-like domains. Next, we calculate the spectroscopic limited maximum efficiency (SLME) of the CuAu-like phase and compare the results with those of the corresponding chalcopyrite phase. We identify several candidates with a high efficiency, such as CuAu-like CuInS2, for which we obtain an SLME of 29% at a thickness of 500 nm. We observe that the SLME can have values above the Shockley-Queisser (SQ) limit, and show that this can occur because the SQ limit assumes the absorptivity to be a step function, thus overestimating the radiative recombination in the detailed balance approach. This means that it is possible to find higher theoretical efficiencies within this framework simply by calculating the J-V characteristic with an absorption spectrum. Finally, we expand our SLME analysis to indirect band gap absorbers by studying silicon, and find that the SLME quickly overestimates the reverse saturation current of indirect band gap materials, drastically lowering their calculated efficiency.

Abstract: Mg-Ti nanostructured samples with different Ti contents were prepared via compaction of nanoparticles grown by inert gas condensation with independent Mg and Ti vapour sources. The growth set-up offered the option to perform in situ hydrogen absorption before compaction. Structural and morphological characterisation was carried out by X-ray diffraction, energy dispersive spectroscopy and electron microscopy. The formation of an extended metastable solid solution of Ti in hcp Mg was detected up to 15 at% Ti in the as-grown nanoparticles, while after in situ hydrogen absorption, phase separation between MgH2 and TiH2 was observed. At a Ti content of 22 at%, a metastable Mg-Ti-H fcc phase was observed after in situ hydrogen absorption. The co-evaporation of Mg and Ti inhibited nanoparticle coalescence and crystallite growth in comparison with the evaporation of Mg only. In situ hydrogen absorption was beneficial to subsequent hydrogen behaviour, studied by high pressure differential scanning calorimetry and isothermal kinetics. A transformed fraction of 90% was reached within 100 s at 300 degrees C during both hydrogen absorption and desorption. The enthalpy of hydride formation was not observed to differ from bulk MgH2.

Abstract: Gold nanocluster properties exhibit unique size-dependence. In this contribution, we employ reactive molecular dynamics simulations to calculate the size- and temperature-dependent surface energies, strain energies and atomic displacements for icosahedral, cuboctahedral, truncated octahedral and decahedral Au-nanoclusters. The calculations demonstrate that the surface energy decreases with increasing cluster size at 0 K but increases with size at higher temperatures. The calculated melting curves as a function of cluster size demonstrate the Gibbs-Thomson effect. Atomic displacements and strain are found to strongly depend on the cluster size and both are found to increase with increasing cluster size. These results are of importance for understanding the size-and temperature-dependent surface processes on gold nanoclusters.

Abstract: van der Waals (vdW) forces play an important role in the adsorption of molecules on the surface of solids. However, the choice of the most suitable vdW functional for different systems is an essential problem which must be addressed for different systems. The lack of a systematic study on the performance of the vdW functionals in the bulk and adsorption properties of metal-oxides motivated us to examine different vdW approaches and compute the bulk and molecular adsorption properties of alpha-Cr2O3, alpha-Fe2O3, and alpha-Al2O3. For the bulk properties, we compared our results for the heat of formation, cohesive energy, lattice parameters and bond distances between the different vdW functionals and available experimental data. Next we studied the adsorption of benzene and CH3 molecules on top of different oxide surfaces. We employed different approximations to exchange and correlation within DFT, namely, the Perdew-Burke-Ernzerhof (PBE) GGA, (PBE)+U, and vdW density functionals [ DFT(vdW-DF/DF2/optPBE/optB86b/optB88)+U] as well as DFT-D2/D3(+U) methods of Grimme for the bulk calculations and optB86b-vdW(+U) and DFT-D2(+U) for the adsorption energy calculations. Our results highlight the importance of vdW interactions not only in the adsorption of molecules, but importantly also for the bulk properties. Although the vdW contribution in the adsorption of CH3 (as a chemisorption interaction) is less important compared to the adsorption of benzene (as a physisorption interaction), this contribution is not negligible. Also adsorption of benzene on ferryl/chromyl terminated surfaces shows an important chemisorption contribution in which the vdW interactions become less significant.

Abstract: We calculate bubble nucleation rates in a Lennard-Jones fluid through explicit molecular dynamics simulations. Our approach-based on a recent free energy method (dubbed reweighted Jarzynski sampling), transition state theory, and a simple recrossing correction-allows us to probe a fairly wide range of rates in several superheated and cavitation regimes in a consistent manner. Rate predictions from this approach bridge disparate independent literature studies on the same model system. As such, we find that rate predictions based on classical nucleation theory, direct brute force molecular dynamics simulations, and seeding are consistent with our approach and one another. Published rates derived from forward flux sampling simulations are, however, found to be outliers. This study serves two purposes: First, we validate the reliability of common modeling techniques and extrapolation approaches on a paradigmatic problem in materials science and chemical physics. Second, we further test our highly generic recipe for rate calculations, and establish its applicability to nucleation processes.

Abstract: Two-dimensional (2D) materials have attracted growing interest due to their potential use in the next generation of nanoelectronic and optoelectronic applications. On the basis of first-principles calculations based on density functional theory, we first investigate the electronic and mechanical properties of single layer boron phosphide (h-BP). Our calculations show that h-BP is a mechanically stable 2D material with a direct band gap of 0.9 eV at the K-point, promising for both electronic and optoelectronic applications. We next investigate the electron transport properties of a p-n junction constructed from single layer boron phosphide (h-BP) using the non-equilibrium Green's function formalism. The n-and p-type doping of BP are achieved by substitutional doping of B with C and P with Si, respectively. C(Si) substitutional doping creates donor (acceptor) states close to the conduction (valence) band edge of BP, which are essential to construct an efficient p-n junction. By modifying the structure and doping concentration, it is possible to tune the electronic and transport properties of the p-n junction which exhibits not only diode characteristics with a large current rectification but also negative differential resistance (NDR). The degree of NDR can be easily tuned via device engineering.

Abstract: Using first principles calculations we investigate the structural and electronic properties of interfaces between fluorographane and MoS2. Unsymmetrical functionalization of graphene with H and F results in an intrinsic dipole moment perpendicular to the plane of the buckled graphene skeleton. Depending on the orientation of this dipole moment, the electronic properties of a physically absorbed MoS2 monolayer can be switched from n-to p-type or vice versa. We show that one can realize vanishing n-type/p-type Schottky barrier heights when contacting MoS2 to fluorographane. By applying a perpendicular electric field, the size of the Schottky barrier and the degree of doping can be tuned. Our calculations indicate that a fluorographane monolayer is a promising candidate for bipolar doping of MoS2, which is vital in the design of novel technological applications based on two-dimensional materials.

Abstract: In this study, it is predicted by density functional theory calculations that graphene-like novel ultra-thin phases of manganese fluoride crystals, that have nonlayered structures in their bulk form, can be stabilized by fluorination of manganese dichalcogenide crystals. First, it is shown that substitution of fluorine atoms with chalcogens in the manganese dichalcogenide host lattice is favorable. Among possible crystal formations, three stable ultra-thin structures of manganese fluoride, 1H-MnF2, 1T-MnF2 and MnF3, are found to be stable by total energy optimization calculations. In addition, phonon calculations and Raman activity analysis reveal that predicted novel single-layers are dynamically stable crystal structures displaying distinctive characteristic peaks in their vibrational spectrum enabling experimental determination of the corresponding phases. Differing from 1H-MnF2 antiferromagnetic (AFM) large gap semiconductor, 1T-MnF2 and MnF3 single-layers are semiconductors with ferromagnetic (FM) ground state.

Abstract: We study the structure and phase behavior of fullerene-cubane C60·C8H8 by Monte Carlo simulation. Using a simple potential model capturing the icosahedral and cubic symmetries of its molecular constituents, we reproduce the experimentally observed phase transition from a cubic to an orthorhombic crystal lattice and the accompanying rotational freezing of the C60 molecules. We elaborate a scheme to identify the low-temperature orientations of individual molecules and to detect a pattern of orientational ordering similar to the arrangement of C60 molecules in solid C60. Our configuration of orientations supports a doubled periodicity along one of the crystal axes.

Abstract: In order to increase the accuracy of the linearized augmented plane wave (LAPW) method, we present a new approach where the plane wave basis function is augmented by two different atomic radial components constructed at two different linearization energies corresponding to two different electron bands (or energy windows). We demonstrate that this case can be reduced to the standard treatment within the LAPW paradigm where the usual basis set is enriched by the basis functions of the tight binding type, which go to zero with zero derivative at the sphere boundary. We show that the task is closely related with the problem of extended core states which is currently solved by applying the LAPW method with local orbitals (LAPW+LO). In comparison with LAPW+LO, the number of supplemented basis functions in our approach is doubled, which opens up a new channel for the extension of the LAPW and LAPW+LO basis sets. The appearance of new supplemented basis functions absent in the LAPW+LO treatment is closely related with the existence of the ul-component in the canonical LAPW method. We discuss properties of additional tight binding basis functions and apply the extended basis set for computation of electron energy bands of lanthanum (face and body centered structures) and hexagonal close packed lattice of cadmium. We demonstrate that the new treatment gives lower total energies in comparison with both canonical LAPW and LAPW+LO, with the energy difference more pronounced for intermediate and poor LAPW basis sets.

Abstract: The evaluation of nucleation rates from molecular dynamics trajectories is hampered by the slow nucleation time scale and impact of finite size effects. Here, we show that accurate nucleation rates can be obtained in a very general fashion relying only on the free energy barrier, transition state theory, and a simple dynamical correction for diffusive recrossing. In this setup, the time scale problem is overcome by using enhanced sampling methods, in casu metadynamics, whereas the impact of finite size effects can be naturally circumvented by reconstructing the free energy surface from an appropriate ensemble. Approximations from classical nucleation theory are avoided. We demonstrate the accuracy of the approach by calculating macroscopic rates of droplet nucleation from argon vapor, spanning 16 orders of magnitude and in excellent agreement with literature results, all from simulations of very small (512 atom) systems.