Improving intrinsic electrocatalytic activity of layered transition metal chalcogenides as electrocatalysts for water splitting

Electrochemical water splitting has been considered an important method for facilitating renewable and sustainable energy conversion. For the practical application of water electrocatalysis, it is important to develop a non-noble metal-based, earth-abundant, highly efficient, and stable electrocatalysts for water splitting. Among the various non-noble metal-based electrocatalysts, layered transition metal chalcogenides (TMCs) have emerged as fascinating materials for electrochemical water splitting. The unique structural and electronic properties of layered TMCs make them very attractive for understanding the fundamental principles of electrocatalysis, as well as for developing highly efficient and stable electrocatalysts for the practical application of water electrocatalysis. In this mini review, we present a comprehensive overview of recent developments to improve the intrinsic electrocatalytic activity of layered transition metal chalcogenide (TMC)-based electrocatalysts for practical applications in water splitting.     

Computational Assessment of MLCT versus MC Stabilities in First-to-Third-Row d6 Pseudo-Octahedral Transition Metal Complexes

An accurate modeling of metal-to-ligand-charge-transfer (MLCT) and metal-centered (MC) excited state energies is key to predict the photoinduced response in transition metal complexes (TMCs). Herein, the importance of the ground state and excited state reference geometries is addressed for three-prototype d⁶ pseudo-octahedral TMCs, each displaying a different potential energy landscape of MLCT versus MC relative stabilities. Several functionals are used within the time-dependent density functional theory (TDDFT), as well as multireference wave-function theory (MS-CASPT2), applied to [Mn(im)(CO)₃(phen)]⁺, [Ru(im)₂(bpy)₂]²⁺, and [Re(im)(CO)₃(phen)]⁺, (im: imidazole, bpy: bypiridine, phen: phenantroline). The results revel that TDDFT is robust except when using B3LYP functional for first-row d⁶ TMCs. In contrast, MS-CASPT2 calculations are strongly biased in those cases with competitive MLCT/MC states. The results reinforce the reliability of B3LYP to describe the excited states in d⁶ TMCs, but question the validity of assessing the density functional theory (DFT)/TDDFT performance via direct comparison with MS-CASPT2 performed at the same DFT reference geometry as a standard strategy. © 2019 Wiley Periodicals, Inc.     

Studies on the Mechanism for the Bimolecular Metathesis Reaction CH3 + HCI ⇔︁ CH4 + CI

The geometry optimization of the transition state, the precursor complex and the successor complex was performed at the 6–311G* basis set level. From the analysis of the vibrational frequency of the precursor complex, transition state, successor complex and the isolated state, the reaction mechanism was derived which was complicated with the bond‐rupture electron transfer theory. The atom H in molecule HCI attacks the atom C, forming a transition state via the precursor complex and the electron‐transfer happens in precursor complex. And the active energy, electronic coupling matrix element, the reorganization energy, and the reaction rate are obtained.     

Natural transition orbitals for complex two-component excited state calculations

While the natural transition orbital (NTO) method has allowed electronic excitations from time-dependent Hartree-Fock and density functional theory to be viewed in a traditional orbital picture, the extension to multicomponent molecular orbitals such as those used in relativistic two-component methods or generalized Hartree-Fock (GHF) or generalized Kohn-Sham (GKS) is less straightforward due to mixing of spin-components and the inherent inclusion of spin-flip transitions in time-dependent GHF/GKS. An extension of single-component NTOs to the two-component framework is presented, in addition to a brief discussion of the practical aspects of visualizing two-component complex orbitals. Unlike the single-component analog, the method explicitly describes the spin and frequently obtains solutions with several significant orbital pairs. The method is presented using calculations on a mercury atom and a CrO₂Cl₂ complex.     

Magnetism in Simple Metal and 4<Emphasis Type="Italic">d</Emphasis> Transition Metal Clusters

In this review article I discuss two aspects of magnetism in small metal clusters. The first question discussed is whether simple metal clusters, that obey electronic shell models and mimic properties of elemental atoms, also obey Hund’s rule of maximum spin multiplicity. The second question is whether small clusters of 4d transition metal atoms, that are non-magnetic in the bulk, have magnetic ground states. The question arises because calculations showed that small V clusters are magnetic although the bulk metal is not. We discuss known results on Rh clusters in detail to show that small clusters are generally magnetic, but it is difficult to unequivocally identify the ground state due to the presence of many isomers and spin states that are very close in energy.     

Accurate redox potentials of mononuclear iron, manganese, and nickel model complexes*

Density functional theory (DFT) was combined with solution of the Poisson equation for continuum dielectric media to compute accurate redox potentials for several mononuclear transition metal complexes (TMCs) involving iron, manganese, and nickel. Progress was achieved by altering the B3LYP DFT functional (B4(XQ3)LYP-approach) and supplementing it with an empirical correction term G(X) having three additional adjustable parameters, which is applied after the quantum-chemical DFT computations. This method was used to compute 58 redox potentials of 48 different TMCs involving different pairs of redox states solvated in both protic and aprotic solvents. For the 58 redox potentials the root mean square deviation (RMSD) from experimental values is 65 mV. The reliability of the present approach is also supported by the observation that the energetic order of the spin multiplicities of the electronic ground states is fulfilled for all studied TMCs, if the influence from the solvent is considered as well.     

Intermolecular coulomb couplings from ab initio electrostatic potentials: application to optical transitions of strongly coupled pigments in photosynthetic antennae and reaction centers

An accurate and numerically efficient method for the calculation of intermolecular Coulomb couplings between charge densities of electronic states and between transition densities of electronic excitations is presented. The coupling of transition densities yields the F?rster type excitation energy transfer coupling, and from the charge density coupling, a shift in molecular excitation energies results. Starting from an ab initio calculation of the charge and transition densities, atomic partial charges are determined such as to fit the resulting electrostatic potentials of the different states and the transition. The different intermolecular couplings are then obtained from the Coulomb couplings between the respective atomic partial charges. The excitation energy transfer couplings obtained in the present TrEsp (transition charge from electrostatic potential) method are compared with couplings obtained from the simple point-dipole and extended dipole approximations and with those from the ab initio transition density cube method of Krüger, Scholes, and Fleming. The present method is of the same accuracy as the latter but computationally more efficient. The method is applied to study strongly coupled pigments in the light-harvesting complexes of green sulfur bacteria (FMO), purple bacteria (LH2), and higher plants (LHC-II) and the "special pairs" of bacterial reaction centers and reaction centers of photosystems I and II. For the pigment dimers in the antennae, it is found that the mutual orientation of the pigments is optimized for maximum excitonic coupling. A driving force for this orientation is the Coulomb coupling between ground-state charge densities. In the case of excitonic couplings in the "special pairs", a breakdown of the point-dipole approximation is found for all three reaction centers, but the extended dipole approximation works surprisingly well, if the extent of the transition dipole is chosen larger than assumed previously. For the "special pairs", a large shift in local transition energies is found due to charge density coupling.     

Potential-energy surface, van der Waals energy spectrum, and vibronic transitions in s-tetrazine-argon complex

The van der Waals vibrational states and the structure of the vibronic spectrum of s-tetrazine-argon complex have been studied by the ab initio methods. The potential-energy surface of the ground S(0) electronic state of the complex has been constructed by fitting the analytical many-body expansion to a large set of the interaction energy values computed using the second-order M?ller-Plesset perturbation theory combined with the standard aug-cc-pVDZ basis set. The equilibrium structure of the complex found is that with argon located above the tetrazine ring at a distance of 3.394 A. The calculated dissociation energy of 354 cm(-1) is compatible with the experiment. The van der Waals energy spectrum calculated from the potential-energy surface is explained analyzing a correlation with a simpler energy spectrum of benzene-argon. A new assignment of the S(0)-S(1) vibronic spectrum is proposed on the basis of the rigorous selection rules, vibrational energy levels in S(0) and S(1) states and vibronic transition intensities calculated from the electronic transition dipole moment surfaces.     

Multireference space without first solving the configuration interaction problem

We further develop an idea to generate a compact multireference space without first solving the configuration interaction problem previously proposed for the ground state (GS) (Glushkov, Chem. Phys. Lett. 1995, 244, 1). In the present contribution, our attention is focused on low‐lying excited states (ESs) with the same symmetry as the GS which can be adequately described in terms of an high‐spin open‐shell formalism. Two references Møller–Plesset (MP) like perturbation theory for ESs is developed. It is based on: (1) a main reference configuration constructed from the parent molecular orbitals adjusted to a given ES and (2) secondary double excitation configuration built on the GS like orbitals determined by the Hartree–Fock equations subject to some orthogonality constraints. It is shown how to modify the MP zeroth‐order Hamiltonian so that the reference configurations and corresponding excitations are eigenfunctions of it and are compatible with orthogonality conditions for the GS and ES. Intruder states appearance is also discussed. The proposed scheme is applied to the GS, ES, and excitation energies of small molecules to illustrate and calibrate our calculations. © 2013 Wiley Periodicals, Inc.     

Rational design of small indolic squaraine dyes with large two-photon absorption cross section

Small organic dyes with large two-photon absorption (TPA) cross sections (δ) are more desirable in many applications compared with large molecules. Herein, we proposed a facile theoretical method for the fast screening of small organic molecules as potential TPA dyes. This method is based on a theoretical analysis to the natural transition orbitals (NTOs) directly associated with the TPA transition. Experimental results on the small indolic squaraine dyes (ISD) confirmed that their TPA cross sections is strongly correlated to the delocalization degree of the NTOs of the S₂ excited states. Aided by this simple and intuitive method, we have successfully designed and synthesized a small indolic squaraine dye (ISD) with a remarkable δ value above 8000 GM at 780 nm. The ISD dye also exhibits a high singlet oxygen generation quantum yield about 0.90. The rationally designed TPA dye was successfully applied in both two-photon excited fluorescence cell imaging and in vivo cerebrovascular blood fluid tracing.     

Ab initio calculation of the electronic structure of the Ni(CN)4 ion

The electronic structure of the ion Ni(CN)₄²⁻ in its ground and first excited states is described by an ab initio SCF MO calculation. Formation of the complex in its ground state implies a large charge transfer corresponding to σ-bonding, but only a very small charge transfer due to π back-donation. The transitions of lowest energy are found to be d-d transitions which would appear as rather unfavourable on the basis of orbital energy considerations.     

Fast hydrogen elimination from the [Ru(PH3)3(CO)(H)2] complex in the first singlet excited states. A quantum dynamics study

The photodissociation dynamics of [Ru(PH₃)₃(CO)(H)₂] complex in the lowest two singlet excited electronic states has been theoretically analyzed. Reduced two-dimensional potential energy surfaces (PES) are built up by combining a time-dependent method to calculate the excited states energy with DFT (B3LYP) electronic calculations of the ground state. By means of a Fast Fourier Transform (FFT) algorithm the time evolution of the wavefunction upon vertical transition from the minimum of the ground state to both diabatic states has been followed. The propagation in S₁, the lowest in energy at the vertical transition point and the one with a larger transition probability from the ground state, discloses that the system is not evolving from the initial position at least in the time spanned by the calculations. Conversely the H₂ elimination is very fast (about 37 fs) in the S₂ state. In this state the vertical transition puts the system in a purely dissociative zone of the PES. In that state FFT results indicate that the lengthening of the Ru–H₂ distance and the shortening of the H–H one are taking place almost simultaneously.     

Noncovalent interactions in the gas phase: the anisole-phenol complex

The present paper reports on an integrated spectroscopic study of the anisole-phenol complex in a molecular beam environment. Combining REMPI and HR-LIF spectroscopy experimental data with density functional computations (TD-M05-2X/M05-2X//N07D) and first principle spectra simulations, it was possible to locate the band origin of the S(1) ← S(0) electronic transition and determine the equilibrium structure of the complex, both in the S(0) and S(1) electronic states. Experimental and computational evidence indicates that the observed band origin is due to an electronic transition localized on the phenol frame, while it was not possible to localize experimentally another band origin due to the electronic transition localized on the anisole molecule. The observed structure of the complex is stabilized by a hydrogen bond between the phenol, acting as a proton donor, and the anisole molecule, acting as an acceptor through the lone pairs of the oxygen atom. A secondary interaction involving the hydrogen atoms of the anisole methyl group and the π electron system of the phenol molecule stabilizes the complex in a nonplanar configuration. Additional insights about the landscapes of the potential energy surfaces governing the ground and first excited electronic states of the anisole-phenol complex, with the issuing implications on the system photodynamic, can be extracted from the combined experimental and computational studies.     

A model electronic Hamiltonian to study low-lying electronic states of [Fe(bpy)3]2+ in aqueous solution

A simple model electronic Hamiltonian to describe the potential energy surfaces of several low-lying d-d states of the [Fe(bpy)(3)](2+) complex is developed for use in molecular dynamics (MD) simulation studies. On the basis of a method proposed previously for first-row transition metal ions in aqueous solution, the model Hamiltonian is constructed using density functional theory calculations for the lowest singlet and quintet states. MD simulations are then carried out for the two spin states in aqueous solution in order to examine the performance of the model Hamiltonian. The simulation results indicate that the present model electronic Hamiltonian reasonably describes the potential energy surfaces of the two spin states of the aqueous [Fe(bpy)(3)](2+) system, while retaining sufficient simplicity for application in simulation studies on excited state dynamics.     

Experimental and theoretical characterization of the aromatization, epimerization, and fragmentation reactions of bi-2H-azirin-2-yls prepared from 1,4-diazidobuta-1,3-dienes

1,4‐Diazidobuta‐1,3‐dienes (Z,Z)‐ 10 , 17 , and 21 were photolyzed and thermolyzed to yield the pyridazines 13 , 20 , and 23 , respectively. To explain these aromatic final products, the generation of highly strained bi‐2H‐azirin‐2‐yls 12 , 19 , and 22 and their valence isomerization were postulated. In the case of meso‐ and rac‐ 22 , nearly quantitative formation from diazide 21 , isolation as stable solids, and complete characterization were possible. On the thermolysis of 22 , aromatization to 23 was only a side reaction, whereas equilibration of meso‐ and rac‐ 22 and fragmentation, which led to alkyne 24 and acetonitrile, dominated. Prolonged irradiation of 22 gave mainly the pyrimidine 25 . The change of the configuration at C‐2 of the 2H‐azirine unit was observed not only in the case of bi‐2H‐azirin‐2‐yls 22 but also for simple spirocyclic 2H‐azirines 29 at a relatively low temperature (75 °C). The fragmentation of rac‐ 22 to give alkyne 24 and two molecules of acetonitrile was also studied by high‐level quantum chemical calculations. For a related model system 30 (methyl instead of phenyl groups), two transition states TS‐ 30 – 31 of comparable energy with multiconfigurational electronic states could be localized on the energy hypersurface for this one‐step conversion. The symmetrical transition state complies with the definition of a coarctate mechanism.     

State-selective optimization of local excited electronic states in extended systems

Standard implementations of time-dependent density-functional theory (TDDFT) for the calculation of excitation energies give access to a number of the lowest-lying electronic excitations of a molecule under study. For extended systems, this can become cumbersome if a particular excited state is sought-after because many electronic transitions may be present. This often means that even for systems of moderate size, a multitude of excited states needs to be calculated to cover a certain energy range. Here, we present an algorithm for the selective determination of predefined excited electronic states in an extended system. A guess transition density in terms of orbital transitions has to be provided for the excitation that shall be optimized. The approach employs root-homing techniques together with iterative subspace diagonalization methods to optimize the electronic transition. We illustrate the advantages of this method for solvated molecules, core-excitations of metal complexes, and adsorbates at cluster surfaces. In particular, we study the local π→π(?) excitation of a pyridine molecule adsorbed at a silver cluster. It is shown that the method works very efficiently even for high-lying excited states. We demonstrate that the assumption of a single, well-defined local excitation is, in general, not justified for extended systems, which can lead to root-switching during optimization. In those cases, the method can give important information about the spectral distribution of the orbital transition employed as a guess.