Solution phase treatments of Sb2Se3 heterojunction photocathodes for improved water splitting performance.

Antimony selenide (Sb2Se3) is an auspicious material for solar energy conversion that has seen rapid improvement over the past ten years, but the photovoltage deficit remains a challenge. Here, simple and low-temperature treatments of the p-n heterojunction interface of Sb2Se3/TiO2-based photocathodes for photoelectrochemical water splitting were explored to address this challenge. The FTO/Ti/Au/Sb2Se3 (substrate configuration) stack was treated with (NH4)2S as an etching solution, followed by CuCl2 treatment prior to deposition of the TiO2 by atomic layer deposition. The different treatments show different mechanisms of action compared to similar reported treatments of the back Au/Sb2Se3 interface in superstrate configuration solar cells. These treatments collectively increased the onset potential from 0.14 V to 0.28 V vs. reversible hydrogen electrode (RHE) and the photocurrent from 13 mA cm-2 to 18 mA cm-2 at 0 V vs. RHE as compared to the untreated Sb2Se3 films. From SEM and XPS studies, it is clear that the etching treatment induces a morphological change and removes the surface Sb2O3 layer, which eliminates the Fermi-level pinning that the oxide layer generates. CuCl2 further enhances the performance due to the passivation of the surface defects, as supported by density functional theory molecular dynamics (DFT-MD) calculations, improving charge separation at the interface. The simple and low-cost semiconductor synthesis method combined with these facile, low-temperature treatments further increases the practical potential of Sb2Se3 for large-scale water splitting.

Two Novel Dinuclear Cobalt Polypyridyl Complexes in Electro- and Photocatalysis for Hydrogen Production: Cooperativity Increases Performance.

Syntheses and mechanisms of two dinuclear Co-polypyridyl catalysts for the H2 evolution reaction (HER) were reported and compared to their mononuclear analogue (R1). In both catalysts, two di-(2,2'-bipyridin-6-yl)-methanone units were linked by either 2,2'-bipyridin-6,6'-yl or pyrazin-2,5-yl. Complexation with CoII gave dinuclear compounds bridged by pyrazine (C2) or bipyridine (C1). Photocatalytic HER gave turnover numbers (TONs) of up to 20000 (C2) and 7000 (C1) in water. Electrochemically, C1 was similar to the R1, whereas C2 showed electronic coupling between the two Co centers. The E(CoII/I ) split by 360 mV into two separate waves. Proton reduction in DMF was investigated for R1 with [HNEt3 ](BF4 ) by simulation, foot of the wave analysis, and linear sweep voltammetry (LSV) with in-line detection of H2 . All methods agreed well with an (E)ECEC mechanism and the first protonation being rate limiting (≈104  m-1 s-1 ). The second reduction was more anodic than the first one. pKa values of around 10 and 7.5 were found for the two protonations. LSV analysis with H2 detection for all catalysts and acids with different pKa values [HBF4 , pKa (DMF)≈3.4], intermediate {[HNEt3 ](BF4 ), pKa (DMF)≈9.2} to weak [AcOH, pKa (DMF)≈13.5] confirmed electrochemical H2 production, distinctly dependent on the pKa values. Only HBF4 protonated CoI intermediates. The two metals in the dualcore C2 cooperated with an increase in rate to a competitive 105  m-1 s-1 with [HNEt3 ](BF4 ). The overpotential decreased compared to R1 by 100 mV. Chronoamperometry established high stabilities for all catalysts with TONlim of 100 for R1 and 320 for C1 and C2.

Closer Look at Inverse Electron Demand Diels-Alder and Nucleophilic Addition Reactions on s-Tetrazines Using Enhanced Sampling Methods.

Inverse electron demand [4+2] Diels-Alder (iEDDA) reactions as well as unprecedented nucleophilic (azaphilic) additions of R-substituted silyl-enol ethers (where R is Phenyl, Methyl, and Hydrogen) to 1,2,4,5-tetrazine (s-tetrazine) catalyzed by BF 3 have recently been discovered (Simon et al. in Org Lett 23(7):2426-2430, 2021), where static calculations were employed for calculation of activation energies. In order to have a more realistic dynamic description of these reactions in explicit solution at ambient conditions, in this work we use a semiempirical tight-binding method combined with enhanced sampling techniques to calculate free energy surfaces (FESs) of the iEDDA and azaphilic addition reactions. Relevant products of not only s-tetrazine but also its derivatives such as BF 3 -mediated s-tetrazine adducts are investigated. We reconstruct the FESs of the iEDDA and azaphilic addition reactions using metadynamics and blue moon ensemble, and compare the ability of different collective variables (CVs) including bond distances, Social PeRmutation INvarianT (SPRINT) coordinates, and path-CV to describe the reaction pathway. We find that when a bulky Phenyl is used as a substituent at the dienophile the azaphilic addition is preferred over the iEDDA reaction. In addition, we also investigate the effect of BF 3 in the diene and steric hindrance in the dienophile on the competition between the iEDDA and azaphilic addition reactions, providing chemical insight for reaction design. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s11244-021-01516-y.

Machine Learning-Assisted Discovery of Hidden States in Expanded Free Energy Space.

Collective variables (CVs) are crucial parameters in enhanced sampling calculations and strongly impact the quality of the obtained free energy surface. However, many existing CVs are unique to and dependent on the system they are constructed with, making the developed CV non-transferable to other systems. Herein, we develop a non-instructor-led deep autoencoder neural network (DAENN) for discovering general-purpose CVs. The DAENN is used to train a model by learning molecular representations upon unbiased trajectories that contain only the reactant conformers. The prior knowledge of nonconstraint reactants coupled with the here-introduced topology variable and loss-like penalty function are only required to make the biasing method able to expand its configurational (phase) space to unexplored energy basins. Our developed autoencoder is efficient and relatively inexpensive to use in terms of a priori knowledge, enabling one to automatically search for hidden CVs of the reaction of interest.

Water-Assisted Chemical Route Towards the Oxygen Evolution Reaction at the Hydrated (110) Ruthenium Oxide Surface: Heterogeneous Catalysis via DFT-MD and Metadynamics Simulations.

Notwithstanding that RuO2 is a promising catalyst for the oxygen evolution reaction (OER), a plethora of fundamental details on its catalytic properties are still elusive, severely limiting its large-scale deployment. It is also established experimentally that corrosion and wettability of metal oxides can, in fact, enhance the catalytic activity for OER owing to the formation of a hydrated surface layer. However, the mechanistic interplay between surface wettability, interfacial water dynamics and OER across RuO2 , and what degree these processes are correlated are still debated. Herein, spin-polarized Density Functional Theory Molecular Dynamics (DFT-MD) simulations, coupled with advanced enhanced sampling methods in the well-tempered metadynamics framework, are applied to gain a global understanding of RuO2 aqueous interface (explicit water solvent) in catalyzing the OER, and hence possibly help in the design of novel catalysts in the context of photochemical water oxidation. The present study quantitatively assesses the free-energy barriers behind the OER at the (110)-RuO2 catalyst surface revealing plausible pathways composing the reaction network of the O2 evolution. In particular, OER is investigated at room temperature when such a surface is exposed to both gas-phase and liquid-phase water. Albeit a unique efficient pathway has been identified in the gas-phase OER, a surprisingly lowest-free-energy-requiring reaction route is possible when (110)-RuO2 is in contact with explicit liquid water. By estimating the free-energy surfaces associated to these processes, we reveal a noticeable water-assisted OER mechanism which involves a crucial proton-transfer-step assisted by the local water environment. These findings pave the way toward the systematic usage of DFT-MD coupled with metadynamics techniques for the fine assessment of the activity of catalysts, considering finite-temperature and explicit-solvent effects.

Ions Tune Interfacial Water Structure and Modulate Hydrophobic Interactions at Silica Surfaces.

The structure and ultrafast dynamics of the electric double layer (EDL) are central to chemical reactivity and physical properties at solid/aqueous interfaces. While the Gouy-Chapman-Stern model is widely used to describe EDLs, it is solely based on the macroscopic electrostatic attraction of electrolytes for the charged surfaces. Structure and dynamics in the Stern layer are, however, more complex because of competing effects due to the localized surface charge distribution, surface-solvent-ion correlations, and the interfacial hydrogen bonding environment. Here, we report combined time-resolved vibrational sum frequency generation (TR-vSFG) spectroscopy with ab initio DFT-based molecular dynamics simulations (AIMD/DFT-MD) to get direct access to the molecular-level understanding of how ions change the structure and dynamics of the EDL. We show that innersphere adsorbed ions tune the hydrophobicity of the silica-aqueous interface by shifting the structural makeup in the Stern layer from dominant water-surface interactions to water-water interactions. This drives an initially inhomogeneous interfacial water coordination landscape observed at the neat interface toward a homogeneous, highly interconnected in-plane 2D hydrogen bonding (2D-HB) network at the ionic interface, reminiscent of the canonical, hydrophobic air-water interface. This ion-induced transformation results in a characteristic decrease of the vibrational lifetime (T1) of excited interfacial O-H stretching modes from T1 ∼ 600 fs to T1 ∼ 250 fs. Hence, we propose that the T1 determined by TR-vSFG in combination with DFT-MD simulations can be widely used for a quantitative spectroscopic probe of the ion kosmotropic/chaotropic effect at aqueous interfaces as well as of the ion-induced surface hydrophobicity.

Enhanced conductivity of water at the electrified air-water interface: a DFT-MD characterization.

DFT-based molecular dynamics simulations of the electrified air-liquid water interface are presented, where a homogeneous field is applied parallel to the surface plane. We unveil the field intensity for the onset of proton transfer and molecular dissociation; the protonic current/proton conductivity is measured as a function of the field intensity/voltage. The air-water interface is shown to exhibit a proton conductivity twice the one in the liquid water for field intensities below 0.40 V Å-1. We show that this difference arises from the very specific organization of water in the binding interfacial layer (BIL, i.e. the air-water interface region) into a 2D-HBond-network that is maintained and enforced at the electrified interface. Beyond fields of 0.40 V Å-1, water in the BIL and in the bulk liquid are aligned in the same way by the rather intense fields, hence leading to the same proton conductivity in both BIL and bulk water.

DFT-MD of the (110)-Co3O4 cobalt oxide semiconductor in contact with liquid water, preliminary chemical and physical insights into the electrochemical environment.

Within the general context of the electrochemical oxygen evolution reaction of the water oxidation/electrolysis, we focus on one essential aspect of electrochemical interfaces, i.e., the comprehension of the interaction and organisation of liquid water at the (semiconductor) (110)-Co3O4 surface using density functional theory-molecular dynamics simulations. A detailed characterization of the chemical and physical properties of the aqueous interface is provided in terms of structure, dynamics, electric field, work function, and spectroscopy, as a preliminary step into the modelling of the (110)-Co3O4 aqueous surface in more relevant electrochemical conditions. The water at the aqueous B-termination is, in particular, shown more dynamical than that at the A-termination and more "undisciplined": the water is indeed mostly an HB-acceptor with the solid, with an orientation of their dipole moments found opposite the field generated by the negative surface charge. At both aqueous interfaces, the work function is twice lower than that at the bare (non-hydroxylated) surfaces. The SFG (Sum Frequency Generation) spectroscopy is shown dominated by the water in the diffuse layer, while the SFG signal from the binding interfacial layer reflects the single orientation of water at the aqueous A-termination and the two orientations of water at the aqueous B-termination.

Ionic diffusion and proton transfer in aqueous solutions of alkali metal salts.

We report on a series of ab initio molecular dynamics investigations on LiCl, NaCl, and KCl aqueous solutions under the effect of static electric fields. We have found that although in low-to-moderate field intensity regimes the well-known sequence of cationic mobilities µ(K+) > µ(Na+) > µ(Li+) (i.e., the bigger the cation the higher the mobility) is recovered, from intense field strengths this intuitive rule is no longer verified. In fact, field-induced water molecular dissociations lead to more complex phenomena regulating the standard migration properties of the simplest monovalent cations. The water dissociation threshold is lowered from 0.35 V Å-1 to 0.25 V Å-1 by the presence of charged species in all samples. However, notwithstanding a one-stage process of water ionization and proton conduction takes place at 0.25 V Å-1 in the electrolyte solutions where "structure maker" cations are present (i.e., LiCl and NaCl), the KCl aqueous solution shows some hindrance in establishing a proton conductive regime, which is characterized by the same proton conduction threshold of neat water (i.e., 0.35 V Å-1). In addition, it turns out that protons flow easily in the LiCl (σp = 3.0 S cm-1) solution and then - in descending order - in the NaCl (σp = 2.5 S cm-1) and KCl (σp = 2.3 S cm-1) electrolyte solutions. The protonic conduction efficiency is thus inversely proportional to the ionic radii of the cations present in the samples. Moreover, Cl- anions act as a sort of "protonic well" for high field intensities, further lowering the overall proton transfer efficiency of the aqueous solutions. As a consequence, all the recorded protonic conductivities are lower than that for neat water (σp = 7.8 S cm-1), which strongly indicates that devices exploiting the proton transfer ability should be designed so as to minimize the presence of ionic impurities.

Ab initio molecular dynamics study of an aqueous NaCl solution under an electric field.

We report on an ab initio molecular dynamics study of an aqueous NaCl solution under the effect of static electric fields. We found that at low-to-moderate field intensity regimes chlorine ions have a greater mobility than sodium ions which, being a sort of "structure makers", are able to drag their own coordination shells. However, for field strengths exceeding 0.15 V Å(-1) the mobility of sodium ions overcomes that of chlorine ions as both types of ions do actually escape from their respective hydration cages. The presence of charged particles lowers the water dissociation threshold (i.e., the minimum field strength which induces a transfer of protons) from 0.35 V Å(-1) to 0.25 V Å(-1); moreover, a protonic current was also recorded at the estimated dissociation threshold of the solution. The behaviour of the current-voltage diagram of the protonic response to the external electric field is Ohmic as in pure water, with a resulting protonic conductivity of about 2.5 S cm(-1). This value is approximately one third of that estimated in pure water (7.8 S cm(-1)), which shows that the partial breaking of hydrogen bonds induced by the solvated ions hinders the migration of protonic defects. Finally, the conductivity of Na(+) and Cl(-) ions (0.2 S cm(-1)) is in fair agreement with the available experimental data for a solution molarity of 1.7 M.