Surface Passivation and Detrimental Heat-Induced Diffusion Effects in RbF-Treated Cu(In,Ga)Se2 Solar Cell Absorbers.

Alkali postdeposition treatments of Cu(In,Ga)Se2 absorbers with KF, RbF, and CsF have led to remarkable efficiency improvements for chalcopyrite thin film solar cells. However, the effect of such treatments on the electronic properties and defect physics of the chalcopyrite absorber surfaces are not yet fully understood. In this work, we use scanning tunneling spectroscopy and X-ray photoelectron spectroscopy to compare the surface defect electronic properties and chemical composition of RbF-treated and nontreated absorbers. We find that the RbF treatment is effective in passivating electronic defect levels at the surface by preventing surface oxidation. Our X-ray photoelectron spectroscopy (XPS) data points to the presence of chemisorbed Rb on the surface with a bonding configuration similar to that of a RbInSe2 bulk compound. Yet, a quantitative analysis indicates Rb coverage in the submonolayer regime, which is likely causing the surface passivation. Furthermore, ab initio calculations confirm that RbF-treated surfaces are less prone to oxidation (in the form of Ga, In, and Se oxides) than bare chalcopyrite surfaces. In addition, elemental diffusion of Rb along with Na, Cu, and Ga is found to occur when the samples are annealed under ultrahigh vacuum conditions. Magnetic sector secondary ion mass spectrometry measurements indicate that there is a homogeneous spatial distribution of Rb on the surface both before and after annealing, albeit with an increased concentration at the surface after heat treatment. Depth-resolved magnetic sector secondary ion mass spectrometry measurements show that Rb diffusion within the bulk occurs predominantly along grain boundaries. Scanning tunneling and XPS measurements after subsequent annealing steps demonstrate that the Rb accumulation at the surface leads to the formation of metallic Rb phases, involving a significant increase of electronic defect levels and/or surface dipole formation. These results strongly suggest a deterioration of the absorber-window interface because of increased recombination losses after the heat-induced diffusion of Rb toward the interface.

Photocatalytic Water Splitting Reaction Catalyzed by Ion-Exchanged Salts of Potassium Poly(heptazine imide) 2D Materials.

Potassium poly (heptazine imide) (K-PHI), a crystalline two-dimensional carbon-nitride material, is an active photocatalyst for water splitting. The potassium ions in K-PHI can be exchanged with other ions to change the properties of the material and eventually to design the catalysts. We report here the electronic structures of several ion-exchanged salts of K-PHI (K, H, Au, Ru, and Mg) and their feasibility as water splitting photocatalysts, which were determined by density functional theory (DFT) calculations. The DFT results are complemented by experiments where the performances in the photocatalytic hydrogen evolution reaction (HER) were recorded. We show that due to its narrow band gap, Ru-PHI is not a suitable photocatalyst. The water oxidation potentials are straddled between the band edge potentials of H-PHI, Au-PHI, and Mg-PHI; thus, these are active photocatalysts for both the oxygen and hydrogen evolution reactions, whereas K-PHI is active only for the HER. The experimental data show that these are active HER photocatalysts, in agreement with the DFT results. Furthermore, Mg-PHI has shown remarkable performance in the HER, with a rate of 539 µmol/(h·g) and a quantum efficiency of 7.14% at 410 nm light irradiation, which could be due to activation of the water molecule upon adsorption, as predicted by our DFT calculations.

Electrochemical N2 Reduction to Ammonia Using Single Au/Fe Atoms Supported on Nitrogen-Doped Porous Carbon.

The electrochemical nitrogen reduction reaction (NRR) to ammonia (NH3) is a promising alternative route for an NH3 synthesis at ambient conditions to the conventional high temperature and pressure Haber-Bosch process without the need for hydrogen gas. Single metal ions or atoms are attractive candidates for the catalytic activation of non-reactive nitrogen (N2), and for future targeted improvement of NRR catalysts, it is of utmost importance to get detailed insights into structure-performance relationships and mechanisms of N2 activation in such structures. Here, we report density functional theory studies on the NRR catalyzed by single Au and Fe atoms supported in graphitic C2N materials. Our results show that the metal atoms present in the structure of C2N are the reactive sites, which catalyze the aforesaid reaction by strong adsorption and activation of N2. We further demonstrate that a lower onset electrode potential is required for Fe-C2N than for Au-C2N. Thus, Fe-C2N is theoretically predicted to be a potentially better NRR catalyst at ambient conditions than Au-C2N owing to the larger adsorption energy of N2 molecules. Furthermore, we have experimentally shown that single sites of Au and Fe supported on nitrogen-doped porous carbon are indeed active NRR catalysts. However, in contrast to our theoretical results, the Au-based catalyst performed slightly better with a Faradaic efficiency (FE) of 10.1% than the Fe-based catalyst with an FE of 8.4% at -0.2 V vs. RHE. The DFT calculations suggest that this difference is due to the competitive hydrogen evolution reaction and higher desorption energy of ammonia.

In silico investigation of Cu(In,Ga)Se2-based solar cells.

Photovoltaics is one of the most promising and fastest-growing renewable energy technologies. Although the price-performance ratio of solar cells has improved significantly over recent years, further systematic investigations are needed to achieve higher performance and lower cost for future solar cells. In conjunction with experiments, computer simulations are powerful tools to investigate the thermodynamics and kinetics of solar cells. Over the last few years, we have developed and employed advanced computational techniques to gain a better understanding of solar cells based on copper indium gallium selenide (Cu(In,Ga)Se2). Furthermore, we have utilized state-of-the-art data-driven science and machine learning for the development of photovoltaic materials. In this Perspective, we review our results along with a survey of the field.

On the Possibility of Helium Adsorption in Nitrogen Doped Graphitic Materials.

The potassium salt of polyheptazine imide (K-PHI) is a promising photocatalyst for various chemical reactions. From powder X-ray diffraction data an idealized structural model of K-PHI has been derived. Using atomic coordinates of this model we defined an energetically optimized K-PHI structure, in which the K ions are present in the pore and between the PHI-planes. The distance between the anion framework and K+ resembles a frustrated Lewis pair-like structure, which we denote as frustrated Coulomb pair that results in an interesting adsorption environment for otherwise non-adsorbing, non-polar gas molecules. We demonstrate that even helium (He) gas molecules, which are known to have the lowest boiling point and the lowest intermolecular interactions, can be adsorbed in this polarized environment with an adsorption energy of  - 4.6 kJ mol-1 per He atom. The interaction between He atoms and K-PHI is partially originating from charge transfer, as disclosed by our energy decomposition analysis based on absolutely localized molecular orbitals. Due to very small charge transfer interactions, He gas adsorption saturates at 8 at%, which however can be subject to further improvement by cation variation.

Interfacing the Core-Shell or the Drude Polarizable Force Field With Car-Parrinello Molecular Dynamics for QM/MM Simulations.

We report a quantum mechanics/polarizable-molecular mechanics (QM/p-MM) potential based molecular dynamics (MD) technique where the core-shell (or the Drude) type polarizable MM force field is interfaced with the plane-wave density functional theory based QM force field which allows Car-Parrinello MD for the QM subsystem. In the QM/p-MM Lagrangian proposed here, the shell (or the Drude) MM variables are treated as extended degrees of freedom along with the Kohn-Sham (KS) orbitals describing the QM wavefunction. The shell and the KS orbital degrees of freedom are then adiabatically decoupled from the nuclear degrees of freedom. In this respect, we also present here the Nosé-Hoover Chain thermostat implementation for the dynamical subsystems. Our approach is then used to investigate the effect of MM polarization on the QM/MM results. Especially, the consequence of MM polarization on reaction free energy barriers, defect formation energy, and structural and dynamical properties are investigated. A low point charge polarizable potential (p-MZHB) for pure siliceous systems is also reported here.

CPMD/GULP QM/MM interface for modeling periodic solids: Implementation and its application in the study of Y-zeolite supported Rhn clusters.

We report here the development of hybrid quantum mechanics/molecular mechanics (QM/MM) interface between the plane-wave density functional theory based CPMD code and the empirical force-field based GULP code for modeling periodic solids and surfaces. The hybrid QM/MM interface is based on the electrostatic coupling between QM and MM regions. The interface is designed for carrying out full relaxation of all the QM and MM atoms during geometry optimizations and molecular dynamics simulations, including the boundary atoms. Both Born-Oppenheimer and Car-Parrinello molecular dynamics schemes are enabled for the QM part during the QM/MM calculations. This interface has the advantage of parallelization of both the programs such that the QM and MM force evaluations can be carried out in parallel to model large systems. The interface program is first validated for total energy conservation and parallel scaling performance is benchmarked. Oxygen vacancy in α-cristobalite is then studied in detail and the results are compared with a fully QM calculation and experimental data. Subsequently, we use our implementation to investigate the structure of rhodium cluster (Rhn ; n = 2 to 6) formed from Rh(C2 H4 )2 complex adsorbed within a cavity of Y-zeolite in a reducible atmosphere of H2 gas. © 2016 Wiley Periodicals, Inc.

A potential with low point charges for pure siliceous zeolites.

A modified CHARMM force-field (ZHB potential) with low point charges for silica was previously proposed by Zimmerman et al. (J. Chem. Theory Comput. 2011, 7, 1695). The ZHB potential is advantageous for quantum mechanics/molecular mechanics simulations as it minimizes the electron spill-out problems. In the same spirit, here we propose a modified ZHB potential (MZHB) by reformulating its bonding potential, while retaining the nonbonding potential as in the ZHB force-field. We show that several structural and dynamic properties of silica, like the IR spectrum, distribution functions, mechanical properties, and negative thermal expansion computed using the MZHB potential agree well with experimental data. Further, transferability of MZHB is also tested for reproducing the crystallographic structures of several polymorphs of silica. © 2015 Wiley Periodicals, Inc.