Search

International audience; We present an analytical thermodynamic model of a surface in contact with a gas phase which enables us to determine the surface coverage of the adsorbate depending on the temperature, pressure and chemical potential of the gas. This model is applied to both the W(110) and W(100) surface of tungsten in contact with hydrogen. The thermodynamic model is built upon data computed by density functional theory that provide the complete electronic and vibrational energetics of both surfaces. It is further compared to experimental measurements of hydrogen coverage during isobar exposure at various temperature acquired via low energy ion scattering and direct recoil spectroscopy techniques. On W(110), the agreement is quantitative provided that an additional degree of freedom is added to the model. This degree of freedom accounts for the translational motion of the adsorbate along 1D channel on the surface. On W(100), surface reconstruction makes the energetics of the system more complex; the full details of the experimental isobar are not well reproduced, although the overall consistency is obtained. We end-up with a thermodynamic model based on DFT data having accurate predictive capabilities to determine the hydrogen coverage of tungsten at any p and T.

International audience; Herein, we investigate the saturation limits of hydrogen on the (110) and (100) surfaces of tungsten via Density Functional Theory (DFT) and complement our findings with experimental measurements. We present a detailed study of the various stable configurations that hydrogen can adopt upon the surfaces at coverage ratios starting below 1.0, up to the point of their experimental coverage ratios, and beyond. Our findings allow us to estimate that the saturation limit on each surface exists with one monolayer of hydrogen atoms adsorbed. In the case of (110) this corresponds to a coverage ratio of one hydrogen atom per tungsten atom, while in the case of (100) a full monolayer is present at a coverage ratio of 2.0. Preliminary Low Energy Ion Scattering (LEIS) and Direct Recoil Spectroscopy (DRS) measurements complement these results and tend to confirm the findings obtained by DFT. In particular, the preferred adsorption sites on both surfaces at any coverage, the reconstruction of the (100) surface and the saturation limits agree well. We show that depending on the coverage, hydrogen surface binding energies can be of the same magnitude as binding energies to defects like vacancies. As a consequence, surface effects should be included in models aiming to simulate retention and desorption of hydrogen from the bulk.

International audience; Tungsten will be used as a plasma facing material in the next generation of fusion reactors. To aid in understanding atomic scale H-W interactions, we investigated hydrogen coverage on the tungsten (110) surface via periodic density-functional theory, providing the most stable configurations that hydrogen forms on the surface at coverage ratios of interval 0.25, step-wise, up to a full mono-layer of hydrogen. We then calculate the Gibbs free energy for the stable configurations in the presence of hydrogen gas at specified temperature and pressure. It follows that the configuration, and corresponding coverage ratio, which yields the lowest Gibbs free energy is used to estimate the macroscopic surface state. Our findings based on the model compare well to low energy electron diffraction (LEED) measurements, primarily the presence of well-ordered phases with coverage ratios 0.5, 0.75, and 1.0, respectively, and that no ordered phases are expected as temperature increases and surface depletion occurs. 2

International audience; In this work, a kinetic model is presented to describe hydrogen absorption and desorption from tungsten at different surface coverages. Activation energies for hydrogen absorption into the bulk and desorption from the surface of tungsten are modelled by functions that depend explicitly and continuously on the hydrogen surface coverage. A steady-state model is developed to derive these activation energies from experimental data. The newly developed coverage dependent activation energies are then implemented in the non steady-state rate-equation code MHIMS. Published experimental results on D uptake and retention of self-damaged tungsten exposed to 0.28 eV deuterium atoms at different temperatures ranging from 450 K to 1000 K can be successfully described with this approach. Finally, the steady-state model is applied to determine surface concentration, bulk concentration and migration depths of hydrogen isotopes in tungsten exposed to various atomic fluxes and temperatures ranging from milder conditions in laboratory experiments to divertor strike point conditions in tokamaks.

Herein we investigate absorption and desorption of hydrogen in the sub-surface of tungsten via Density Functional Theory. Both the near-surface diffusion and recombination of a bulk hydrogen atom with a hydrogen atom adsorbed upon the W(110) and W(100) surfaces are investigated at various surface adsorption coverage ratios. This study intends to model the desorption processes occurring during Thermal-Desorption Spectroscopy experiments and the absorption of hydrogen during gaseous or low energy atomic exposure. Since the diffusion and recombination processes are expected to change as the hydrogen coverage of the surface varies, different coverage ratios were investigated on both surfaces. We found that at saturation coverage of hydrogen on both surfaces, the activation barriers for the recombination of molecular hydrogen are below 0.8 eV. On the contrary, below saturation, the activation barriers for recombination rise to 1.35 eV and 1.51 eV depending on the coverage and on the orientation of the surface. Regarding the absorption of atomic hydrogen from the surface into the bulk, the activation barrier raises from less than 1.0 eV at saturation to around 1.7 eV below saturation on both surfaces. These results indicate that surface mechanisms certainly play a significant role in the kinetics of desorption of hydrogen from tungsten; it is also expected that surface mechanisms affect the total amount on hydrogen absorbed in tungsten during implantation.

International audience; In this work, we combine Density Functional Theory data with a Thermodynamic and a kinetic model to determine the total concentration of hydrogen implanted in the sub-surface of tungsten exposed to a hydrogen flux. The sub-surface hydrogen concentration is calculated given a flux of hydrogen, a temperature of implantation, and the energy of the incoming hydrogen ions as independent variables. This global model is built step by step; an equilibrium between atomic hydrogen within bulk tungsten and a molecular hydrogen gas phase is first considered, and the calculated solubility is compared with experimental results. Subsequently, a kinetic model is used to determine the chemical potential for hydrogen in the sub-surface of tungsten. Combining both these models, two regimes are established in which hydrogen is preferentially trapped at either interstitial sites or in vacancies. We deduce from our model that the existence of these two regimes is driven by the temperature of the implanted tungsten sample; above a threshold or transition temperature is the interstitial regime, below is the vacancy regime in which super-saturated layers form within tenths of angstrom below the surface. A simple analytical expression is derived for the coexistence of the two regimes depending on the implantation temperature, the incident energy and the flux of the hydrogen ions which we use to plot the corresponding phase diagram.

We report a joint photoelectron spectroscopy and theoretical investigation of the gaseous Au2I3− cluster, which is found to exhibit two types of isomers due to competition between Au–I covalent bonding and Au–Au aurophilic interactions. The covalent bonding favors a bent IAuIAuI− structure with an obtuse Au–I–Au angle (100.7°), while aurophilic interactions pull the two Au atoms much closer, leading to an acutely bent structure (72.0°) with an Au–Au distance of 3.08 A. The two isomers are separated by a small barrier and are nearly degenerate with the obtuse isomer being slightly more stable. At low temperature, only the obtuse isomer is observed; distinct experimental evidence is observed for the co-existence of a combination of isomers with both acute and obtuse bending angles at room temperature. The two bond-bending isomers of Au2I3− reveal a unique example of one molecule being able to oscillate between different structures as a result of two competing chemical forces.

The behavior of hydrogen in perfect wurtzite beryllium oxide is herein investigated by means of electronic structure calculations based on density functional theory. The formation energies of the following set of states of hydrogen (H0, H+, H-, H2, [Formula: see text], [Formula: see text]) are computed and their solubility is established as a function of temperature and pressure with emphasis given to conditions relevant for hydrogen-implanted materials. It is found that all magnetic states H0, [Formula: see text], [Formula: see text] are unstable, while the relative stability of the non-magnetic states depends on the thermodynamic conditions: H2 prevails above temperatures around 900 K at standard pressure, which is the lowest temperature in experiments measuring the diffusion coefficient of hydrogen in wurtzite beryllium oxide. Under hydrogen implantation, the total concentration of hydrogen is fixed by the implantation source; it is found that molecular hydrogen prevails starting from a very low total concentration in hydrogen (as low as 10-40 at.fr.). Finally, the diffusion coefficient of H2 in beryllium oxide is calculated and the results are compared with previous experimental data.

Boron is an interesting element with unusual polymorphism. While three-dimensional (3D) structural motifs are prevalent in bulk boron, atomic boron clusters are found to have planar or quasi-planar structures, stabilized by localized two-center-two-electron (2c-2e) σ bonds on the periphery and delocalized multicenter-two-electron (nc-2e) bonds in both σ and π frameworks. Electron delocalization is a result of boron's electron deficiency and leads to fluxional behavior, which has been observed in B13(+) and B19(-). A unique capability of the in-plane rotation of the inner atoms against the periphery of the cluster in a chosen direction by employing circularly polarized infrared radiation has been suggested. Such fluxional behaviors in boron clusters are interesting and have been proposed as molecular Wankel motors. The concepts of aromaticity and antiaromaticity have been extended beyond organic chemistry to planar boron clusters. The validity of these concepts in understanding the electronic structures of boron clusters is evident in the striking similarities of the π-systems of planar boron clusters to those of polycyclic aromatic hydrocarbons, such as benzene, naphthalene, coronene, anthracene, or phenanthrene. Chemical bonding models developed for boron clusters not only allowed the rationalization of the stability of boron clusters but also lead to the design of novel metal-centered boron wheels with a record-setting planar coordination number of 10. The unprecedented highly coordinated borometallic molecular wheels provide insights into the interactions between transition metals and boron and expand the frontier of boron chemistry. Another interesting feature discovered through cluster studies is boron transmutation. Even though it is well-known that B(-), formed by adding one electron to boron, is isoelectronic to carbon, cluster studies have considerably expanded the possibilities of new structures and new materials using the B(-)/C analogy. It is believed that the electronic transmutation concept will be effective and valuable in aiding the design of new boride materials with predictable properties. The study of boron clusters with intermediate properties between those of individual atoms and bulk solids has given rise to a unique opportunity to broaden the frontier of boron chemistry. Understanding boron clusters has spurred experimentalists and theoreticians to find new boron-based nanomaterials, such as boron fullerenes, nanotubes, two-dimensional boron, and new compounds containing boron clusters as building blocks. Here, a brief and timely overview is presented addressing the recent progress made on boron clusters and the approaches used in the authors' laboratories to determine the structure, stability, and chemical bonding of size-selected boron clusters by joint photoelectron spectroscopy and theoretical studies. Specifically, key findings on all-boron hydrocarbon analogues, metal-centered boron wheels, and electronic transmutation in boron clusters are summarized.

Small boron clusters are known to be planar, and may be used as ligands to form novel coordination complexes with transition metals. Here we report a combined photoelectron spectroscopy and ab initio study of CoB12(-) and RhB12(-). Photoelectron spectra of the two doped-B12 clusters show similar spectral patterns, suggesting they have similar structures. Global minimum searches reveal that both CoB12(-) and RhB12(-) possess half-sandwich-type structures with the quasi-planar B12 moiety coordinating to the metal atom. The B12 ligand is found to have similar structure as the bare B12 cluster with C3v symmetry. Structures with Co or Rh inserted into the quasi-planar boron framework are found to be much higher in energy. Chemical bonding analyses of the two B12 half sandwiches reveal two sets of σ bonds on the boron unit: nine classical two-center-two-electron (2c-2e) σ bonds on the periphery of the B12 unit and four 3c-2e σ bonds within the boron unit. Both σ and π bonds are found between the metal and the B12 ligand: three M-B single σ bonds and one delocalized 4c-2e π bond. The exposed metal sites in these complexes can be further coordinated by other ligands or become reaction centers as model catalysts.

We report the observation of two transition-metal-centered nine-atom boron rings, RhⓒB(9)(-) and IrⓒB(9)(-). These two doped-boron clusters are produced in a laser-vaporization supersonic molecular beam and characterized by photoelectron spectroscopy and ab initio calculations. Large HOMO-LUMO gaps are observed in the anion photoelectron spectra, suggesting that neutral RhⓒB(9) and IrⓒB(9) are highly stable, closed shell species. Theoretical calculations show that RhⓒB(9) and IrⓒB(9) are of D(9h) symmetry. Chemical bonding analyses reveal that these complexes are doubly aromatic, each with six completely delocalized π and σ electrons, which describe the bonding between the central metal atom and the boron ring. This work establishes firmly the metal-doped B rings as a new class of novel aromatic molecular wheels.

ConspectusBoron is an interesting element with unusual polymorphism. While three-dimensional (3D) structural motifs are prevalent in bulk boron, atomic boron clusters are found to have planar or quasi-planar structures, stabilized by localized two-center–two-electron (2c–2e) σ bonds on the periphery and delocalized multicenter–two-electron (nc–2e) bonds in both σ and π frameworks. Electron delocalization is a result of boron’s electron deficiency and leads to fluxional behavior, which has been observed in B13+ and B19–. A unique capability of the in-plane rotation of the inner atoms against the periphery of the cluster in a chosen direction by employing circularly polarized infrared radiation has been suggested. Such fluxional behaviors in boron clusters are interesting and have been proposed as molecular Wankel motors. The concepts of aromaticity and antiaromaticity have been extended beyond organic chemistry to planar boron clusters. The validity of these concepts in understanding the electronic structu...

After the discovery of fullerene-C60, it took almost two decades for the possibility of boron-based fullerene structures to be considered. So far, there has been no experimental evidence for these nanostructures, in spite of the progress made in theoretical investigations of their structure and bonding. Here we report the observation, by photoelectron spectroscopy, of an all-boron fullerene-like cage cluster at B40(-) with an extremely low electron-binding energy. Theoretical calculations show that this arises from a cage structure with a large energy gap, but that a quasi-planar isomer of B40(-) with two adjacent hexagonal holes is slightly more stable than the fullerene structure. In contrast, for neutral B40 the fullerene-like cage is calculated to be the most stable structure. The surface of the all-boron fullerene, bonded uniformly via delocalized σ and π bonds, is not perfectly smooth and exhibits unusual heptagonal faces, in contrast to C60 fullerene.

The structure and chemical bonding of the 24-atom boron cluster are investigated using photoelectron spectroscopy and ab initio calculations. The joint experimental and theoretical investigation shows that B24 − possesses a quasi-planar structure containing fifteen outer and nine inner atoms with six of the inner atoms forming a filled pentagonal moiety. The central atom of the pentagonal moiety is puckered out of plane by 0.9 A, reminiscent of the six-atom pentagonal caps of the wellknown B12 icosahedral unit. The next closest isomer at the ROCCSD(T) level of theory has a tubular double-ring structure. Comparison of the simulated spectra with the experimental data shows that the global minimum quasi-planar B24 − isomer is the major contributor to the observed photoelectron spectrum, while the tubular isomer has no contribution to the experiment. Chemical bonding analyses reveal that the periphery of the quasi-planar B24 constitutes 15 classical 2c-2e B-B σ-bonds, whereas delocalized σ- and π-bonds are found in the interior of the cluster with one unique 6c-2e π-bond responsible for bonding in the B-centered pentagon. The current work suggests that the 24atom boron cluster continues to be quasi-2D, albeit the tendency to form filled pentagonal units, characteristic of 3D cage-like structures of bulk boron, is observed. © 2013 AIP Publishing LLC.

Boron is carbon's neighbour in the periodic table and has similar valence orbitals. However, boron cannot form graphene-like structures with a honeycomb hexagonal framework because of its electron deficiency. Computational studies suggest that extended boron sheets with partially filled hexagonal holes are stable; however, there has been no experimental evidence for such atom-thin boron nanostructures. Here, we show experimentally and theoretically that B36 is a highly stable quasiplanar boron cluster with a central hexagonal hole, providing the first experimental evidence that single-atom layer boron sheets with hexagonal vacancies are potentially viable. Photoelectron spectroscopy of B36(-) reveals a relatively simple spectrum, suggesting a symmetric cluster. Global minimum searches for B36(-) lead to a quasiplanar structure with a central hexagonal hole. Neutral B36 is the smallest boron cluster to have sixfold symmetry and a perfect hexagonal vacancy, and it can be viewed as a potential basis for extended two-dimensional boron sheets.

Clusters of boron atoms exhibit intriguing size-dependent structures and chemical bonding that are different from bulk boron and may lead to new boron-based nanostructures. We report a combined photoelectron spectroscopic and ab initio study of the 22- and 23-atom boron clusters. The joint experimental and theoretical investigation shows that B(22)(-) and B(23)(-) possess quasi-planar and planar structures, respectively. The quasi-planar B(22)(-) consists of fourteen peripheral atoms and eight interior atoms in a slightly buckled triangular lattice. Chemical bonding analyses of the closed-shell B(22)(2-) species reveal seven delocalized π orbitals, which are similar to those in anthracene. B(23)(-) is a perfectly planar and heart-shaped cluster with a pentagonal cavity and a π-bonding pattern similar to that in phenanthrene. Thus, B(22)(-) and B(23)(-), the largest negatively charged boron clusters that have been characterized experimentally to date, can be viewed as all-boron analogues of anthracene and phenanthrene, respectively. The current work shows not only that boron clusters are planar at very large sizes but also that they continue to yield surprises and novel chemical bonding analogous to specific polycyclic aromatic hydrocarbons.

The transition-metal-centered boron clusters M©B-9 (M: Rh, Ir) are produced in a laser-vaporization supersonic molecular beam and characterized by PES and quantum chemical ab initio calculations at different levels of theory.

Photoelectron spectroscopy and ab initio calculations have been carried out to probe the structures and chemical bonding of the B27 (-) cluster. Comparison between the experimental spectrum and the theoretical results reveals a two-dimensional (2D) global minimum with a triangular lattice containing a tetragonal defect (I) and two low-lying 2D isomers (II and III), each with a hexagonal vacancy. All three 2D isomers have 16 peripheral boron atoms and 11 inner boron atoms. Isomer I is shown to be mainly responsible for the observed photoelectron spectrum with isomers II and III as minor contributors. Chemical bonding analyses of these three isomers show that they all feature 16 localized peripheral B-B σ-bonds. Additionally, isomer I possesses 16 delocalized σ bonds and nine delocalized π bonds, while isomers II and III each contain 17 delocalized σ bonds and eight delocalized π bonds. It is found that the hexagonal vacancy is associated generally with an increase of delocalized σ bonds at the expense of delocalized π bonds in 2D boron clusters. The hexagonal vacancy, characteristic of borophenes, is found to be a general structural feature for mid-sized boron clusters. The current study shows that B27 (-) is the first boron cluster, where a hexagonal vacancy appears among the low-lying isomers accessible experimentally.

Photoelectron spectroscopy and ab initio calculations are used to investigate the structures and chemical bonding of the B25(-) cluster. Global minimum searches reveal a dense potential energy landscape with 13 quasi-planar structures within 10 kcal/mol at the CCSD(T)/6-311+G(d) level of theory. Three quasi-planar isomers (I, II, and III) are lowest in energy and nearly degenerate at the CCSD(T) level of theory, with II and III being 0.8 and 0.9 kcal/mol higher, respectively, whereas at two density functional levels of theory isomer III is the lowest in energy (8.4 kcal/mol more stable than I at PBE0/6-311+G(2df) level). Comparison with experimental photoelectron spectroscopic data shows isomer II to be the major contributor while isomers I and III cannot be ruled out as minor contributors to the observed spectrum. Theoretical analyses reveal similar chemical bonding in I and II, both involving peripheral 2c-2e B-B σ-bonding and delocalized interior σ- and π-bonding. Isomer III has an interesting elongated ribbon-like structure with a π-bonding pattern analogous to those of dibenzopentalene. The high density of low-lying isomers indicates the complexity of the medium-sized boron clusters; the method dependency of predicting relative energies of the low-lying structures for B25(-) suggests the importance of comparison with experiment in determining the global minima of boron clusters at this size range. The appearance of many low-lying quasi-planar structures containing a hexagonal hole in B25(-) suggests the importance of this structural feature in maintaining planarity of larger boron clusters.

The structures and chemical bonding of the B(21)(-) cluster have been investigated by a combined photoelectron spectroscopy and ab initio study. The photoelectron spectrum at 193 nm revealed a very high adiabatic electron binding energy of 4.38 eV for B(21)(-) and a congested spectral pattern. Extensive global minimum searches were conducted using two different methods, followed by high-level calculations of the low-lying isomers. The global minimum of B(21)(-) was found to be a quasiplanar structure with the next low-lying planar isomer only 1.9 kcal/mol higher in energy at the CCSD(T)/6-311-G* level of theory. The calculated vertical detachment energies for the two isomers were found to be in good agreement with the experimental spectrum, suggesting that they were both present experimentally and contributed to the observed spectrum. Chemical bonding analyses showed that both isomers consist of a 14-atom periphery, which is bonded by classical two-center two-electron bonds, and seven interior atoms in the planar structures. A localized two-center two-electron bond is found in the interior of the two planar isomers, in addition to delocalized multi-center σ and π bonds. The structures and the delocalized bonding of the two lowest lying isomers of B(21)(-) were found to be similar to those in the two lowest energy isomers in B(19)(-).

We report an investigation of the electronic structure and chemical bonding of AuH2− using photoelectron spectroscopy and ab initio calculations. We obtained vibrationally resolved photoelectron spectra of AuH2− at several photon energies. Six electronic states of AuH2 were observed and assigned according to the theoretical calculations. The ground state of AuH2− is known to be linear, while that of neutral AuH2 is bent with a ∠H–Au–H equilibrium bond angle of 129°. This large geometry change results in a very broad bending vibrational progression in the photoelectron spectra for the ground-state transition. The electron affinity of AuH2 is measured to be 3.030 ± 0.020 eV. A short bending vibrational progression is also observed in the second photodetachment band, suggesting a slightly bent structure for the first excited state of AuH2. The linear geometry is a saddle point for the ground and first excited states of AuH2, resulting in double-well potentials for these states along the bending coordinate. Spectroscopic evidence is observed for the detachment transitions to the double-well potentials of the ground and first excited states of AuH2. Higher excited states of AuH2 due to detachment from the nonbonding Au 5d electrons are all linear, similar to the anion ground state. Kohn–Sham molecular orbital analyses reveal surprising participation of H 2p orbitals in the Au–H chemical bonding and an unprecedented weak Au 5dπ to H 2pπ back donation. The simplicity of the linear AuH2− anion and its novel spectroscopic features make it a textbook example for understanding the covalent bonding properties and relativistic effects of Au.

A class of transition-metal-centered aromatic boron wheels (D(nh)-M©B(n)(q-)) have been recently produced and characterized according to an electronic design principle. Here we investigate the interplay between electronic and geometric requirements for the molecular wheels using the case of VB(10)(-), which is isoelectronic to the decacoordinated molecular wheels, Ta©B(10)(-) and Nb©B(10)(-). Photoelectron spectra of VB(10)(-) are observed to be broad and complicated with relatively low electron binding energies, in contrast to the simple and high electron binding energies observed for the molecular wheels of its heavier congeners. An unbiased global minimum search found the most stable isomer of VB(10)(-) to be a singlet "boat"-like structure (C(2)), in which the V atom is coordinated to a quasi-planar B(10) unit. A similar triplet C(2v) boat-like isomer is found to be almost degenerate to the C(2) structure, whereas the beautiful molecular wheel structure, D(10h)-V©B(10)(-), is significantly higher in energy on the potential energy surface. Therefore, even though the VB(10)(-) system fulfills the electronic requirement to form a D(10h)-M©B(10)(-) aromatic molecular wheel, the V atom is too small to stabilize the ten-membered boron ring.

International audience; Fusion fuel retention (trapping) and release (desorption) from plasma-facing components are critical issues for ITER and for any future industrial demonstration reactors such as DEMO. Therefore, understanding the fundamental mechanisms behind the retention of hydrogen isotopes in first wall and divertor materials is necessary. We developed an approach that couples dedicated experimental studies with modelling at all relevant scales, from microscopic elementary steps to macroscopic observables, in order to build a reliable and predictive fusion reactor wall model. This integrated approach is applied to the ITER divertor material (tungsten), and advances in the development of the wall model are presented. An experimental dataset, including focused ion beam scanning electron microscopy, isothermal desorption, temperature programmed desorption, nuclear reaction analysis and Auger electron spectroscopy, is exploited to initialize a macroscopic rate equation wall model. This model includes all elementary steps of modelled experiments: implantation of fusion fuel, fuel diffusion in the bulk or towards the surface, fuel trapping on defects and release of trapped fuel during a thermal excursion of materials. We were able to show that a single-trap-type single-detrapping-energy model is not able to reproduce an extended parameter space study of a polycrystalline sample exhibiting a single desorption peak. It is therefore justified to use density functional theory to guide the initialization of a more complex model. This new model still contains a single type of trap, but includes the density functional theory findings that the detrapping energy varies as a function of the number of hydrogen isotopes bound to the trap. A better agreement of the model with experimental results is obtained when grain boundary defects are included, as is consistent with the polycrystalline nature of the studied sample. Refinement of this grain boundary model is discussed as well as the inclusion in the model of a thin defective oxide layer following the experimental observation of the presence of an oxygen layer on the surface even after annealing to 1300 K.