Zhi-Pan Liu

Observation of an all-boron fullerene

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.

Mechanism and Tafel Lines of Electro-Oxidation of Water to Oxygen on RuO2(110)

How to efficiently oxidize H(2)O to O(2) (H(2)O → 1/2O(2) + 2H(+) + 2e(-)) is a great challenge for electrochemical/photo water splitting owing to the high overpotential and catalyst corrosion. Here extensive periodic first-principles calculations integrated with modified-Poisson-Boltzmann electrostatics are utilized to reveal the physical origin of the high overpotential of the electrocatalytic oxygen evolution reaction (OER) on RuO(2)(110). By determining the surface phase diagram, exploring the possible reaction channels, and computing the Tafel lines, we are able to elucidate some long-standing puzzles on the OER kinetics from the atomic level. We show that OER occurs directly on an O-terminated surface phase above 1.58 V vs NHE, but indirectly on a OH/O mixed phase below 1.58 V by converting first the OH/O mixed phase to the O-terminated phase locally. The rate-determining step of OER involves an unusual water oxidation reaction following a Eley-Rideal-like mechanism, where a water molecule from solution breaks its OH bond over surface Os with concurrent new O-OH bond formation. The free energy barrier is 0.74 eV at 1.58 V, and it decreases linearly with the increase of potential above 1.58 V (a slope of 0.56). In contrast, the traditionally regarded surface oxygen coupling reaction with a Langmuir-Hinshelwood mechanism is energetically less favored and its barrier is weakly affected by the potential. Fundamentally, we show that the empirical linear barrier~potential relation is caused by the linear structural response of the solvated transition state to the change of potential. Finally, the general strategy for finding better OER anode is also presented.

Mechanism and Activity of Photocatalytic Oxygen Evolution on Titania Anatase in Aqueous Surroundings

Due to its high overpotential and low efficiency, the conversion of water to O(2) using solar energy remains a bottleneck for photocatalytic water splitting. Here the microscopic mechanisms of the oxygen evolution reaction (OER) on differently structured anatase surfaces in aqueous surroundings, namely, (101), (001), and (102), are determined and compared systematically by combining first-principles density functional theory calculations and a parallel periodic continuum solvation model. We show that OER involves the sequential removal of protons from surface oxidative species, forming surface peroxo and superoxo intermediates. The initiating step, the first proton removal, dictates the high overpotential. Only at an overpotential of 0.7 V (1.93 V vs SHE) does this rate-controlling step become surmountable at room temperature: the free energy change of the step is 0.69, 0.63, and 0.61 eV for (101), (102), and (001) surfaces, respectively. We therefore conclude that (i) OER is not sensitive to the local surface structure of anatase and (ii) visible light (<∼590 nm) is, in principle, capable of driving the photocatatlytic OER on anatase kinetically. By co-doping high-valent elements into the anatase subsurface, we demonstrate that the high overpotential of the OER can be significantly reduced, with extra occupied levels above the valence band.

Aggregation-induced phosphorescent emission (AIPE) of iridium(III) complexes

A novel aggregation-induced phosphorescent emission (AIPE) was observed for iridium(III) complexes. This interesting phenomenon was attributed to the intermolecular packing, resulting in a switch from the non-emissive 3LX excited state to the emissive 3MLLCT transition, which is confirmed by X-ray diffraction studies as well as theoretical calculations.

Particle size, shape and activity for photocatalysis on titania anatase nanoparticles in aqueous surroundings.

TiO(2) nanoparticles have been widely utilized in photocatalysis, but the atomic level understanding on their working mechanism falls much short of expectations. In particular, the correlation between the particle structure and the photocatalytic activity is not established yet, although it was observed that the activity is sensitive to the particle size and shape. This work, by investigating a series of TiO(2) anatase nanoparticles with different size and shape as the photocatalyst for water oxidation, correlates quantitatively the particle size and shape with the photocatalytic activity of the oxygen evolution reaction (OER). Extensive density functional theory (DFT) calculations combined with the periodic continuum solvation model have been utilized to compute the electronic structure of nanoparticles in aqueous solution and provide the reaction energetics for the key elementary reaction. We demonstrate that the equilibrium shape of nanoparticle is sensitive to its size from 1 to 30 nm, and the sharp crystals possess much higher activity than the flat crystals in OER, which in combination lead to the morphology dependence of photocatalytic activity. The conventionally regarded quantum size effect is excluded as the major cause. The physical origin for the shape-activity relationship is identified to be the unique spatial separation/localization of the frontier orbitals in the sharp nanoparticles, which benefits the adsorption of the key reaction intermediate (i.e., OH) in OER on the exposed five-coordinated Ti of {101} facet. The theoretical results here provide a firm basis for maximizing photocatalytic activity via nanostructure engineering and are also of significance for understanding photocatalysis on nanomaterials in general.

Origin and Activity of Gold Nanoparticles as Aerobic Oxidation Catalysts in Aqueous Solution

Whether gold is catalytically active on its own has been hotly debated since the discovery of gold-based catalysis in the 1980s. One of the central controversies is on the O(2) activation mechanism. This work, by investigating aerobic phenylethanol oxidation on gold nanoparticles in aqueous solution, demonstrates that gold nanoparticles are capable to activate O(2) at the solid-liquid interface. Extensive density functional theory (DFT) calculations combined with the periodic continuum solvation model have been utilized to provide a complete reaction network of aerobic alcohol oxidation. We show that the adsorption of O(2) is very sensitive to the environment: the presence of water can double the O(2) adsorption energy to ~0.4 eV at commonly available edge sites of nanoparticles (~4 nm) because of its strongly polarized nature in adsorption. In alcohol oxidation, the hydroxyl bond of alcohol can break only with the help of an external base at ambient conditions, while the consequent α-C-H bond breaking occurs on pure Au, both on edges and terraces, with a reaction barrier of 0.7 eV, which is the rate-determining step. The surface H from the α-C-H bond cleavage can be easily removed by O(2) and OOH via a H(2)O(2) pathway without involving atomic O. We find that Au particles become negatively charged at the steady state because of a facile proton-shift equilibrium on surface, OOH + OH ↔ O(2) + H(2)O. The theoretical results are utilized to rationalize experimental findings and provide a firm basis for utilizing nanoparticle gold as aerobic oxidation catalysts in aqueous surroundings.

General trends in the barriers of catalytic reactions on transition metal surfaces

A catalyst preparation by design is one of the ultimate goals in chemistry. The first step towards this goal is to understand the origin of reaction barriers. In this study, we have investigated several catalytic reactions on some transition metal surfaces, using density functional theory. All the reaction barriers have been determined. By detailed analyses we obtain some insight into the reaction barrier. Each barrier is related to (i) the potential energy surface of reactants on the surface, (ii) the total chemisorption energy of reactants, and (iii) the metal d orbital occupancy and the reactant valency.

Is the uniform electron gas limit important for small Ag clusters? Assessment of different density functionals for Agn (n⩽4)

Twenty-three density functional theory (DFT) methods, including the second- and the third-generation functionals, are tested in conjunction with two basis sets (LANL2DZ and SDD) for studying the properties of neutral and ionic silver clusters. We find that DFT methods incorporating the uniform electron gas limit in the correlation functional, namely, those with Perdews correlation functionals (PW91, PBE, P86, and TPSS), Beckes B95, and the Van Voorhis-Scuseria functional VSXC, generally perform better than the other group of functionals, e.g., those incorporating the LYP correlation functional and variations of the B97 functional. Strikingly, these two groups of functionals can produce qualitatively different results for the Ag3 and Ag4 clusters. The energetic properties and vibrational frequencies of Ag(n) are also evaluated by the different functionals. The present study shows that the choice of DFT methods for heavy metals may be critical. It is found that the exact-exchange-incorporated PBE functional (PBE1PBE) is among the best for predicting the range of properties.

Stochastic Surface Walking Method for Structure Prediction and Pathway Searching

We propose an unbiased general-purpose potential energy surface (PES) searching method for both the structure and the pathway prediction of a complex system. The method is based on the idea of bias-potential-driven dynamics and Metropolis Monte Carlo. A central feature of the method is able to perturb smoothly a structural configuration toward a new configuration and simultaneously has the ability to surmount the high barrier in the path. We apply the method for locating the global minimum (GM) of short-ranged Morse clusters up to 103 atoms starting from a random structure without using extra information from the system. In addition to GM searching, the method can identify the pathways for chemical reactions with large dimensionality, as demonstrated in a nanohelix transformation containing 222 degrees of freedoms.

Origin of Selectivity Switch in Fischer−Tropsch Synthesis over Ru and Rh from First-Principles Statistical Mechanics Studies

For its unique position in chemical industry, Fischer-Tropsch (FT) synthesis has been a hot subject in heterogeneous catalysis. Due to its great complexity in product distribution, it remains unclear how to maximally convert syngas to long-chain hydrocarbons. By combining extensive DFT calculations with grand canonical Monte Carlo simulations, this work examines the key elementary steps in FT synthesis over Ru and Rh surfaces, including CO dissociation, C/C coupling, and hydrogenations. The origin of the relationship between activity and selectivity of catalysts is revealed based on the calculated reaction rate at working temperatures, in which the catalytic role of surface steps as the center of accumulating surface CH(x) species is highlighted. This theoretical work demonstrates that the ability to dissociate CO under carbon-rich conditions is the key requirement for a good FT catalyst. The RC + C (R = alkyl or H) pathway occurring at surface steps may be a general mechanism for FT chain propagation on transition metals.