Water Structure in the First Layers on TiO2: A Key Factor for Boosting Solar-Driven Water-Splitting Performances.

The water hydrogen-bonded network is strongly perturbed in the first layers in contact with the semiconductor surface. Even though this aspect influences the outer-sphere electron transfer, it was not recognized that it is a crucial factor impacting the solar-driven water-splitting performances. To fill this gap, we have selected two TiO2 anatase samples (with and without B-doping), and by extensive experimental and computational investigations, we have demonstrated that the remarkable 5-fold increase in water-splitting photoactivity of the B-doped sample cannot be ascribed to effects typically associated to enhanced photocatalytic properties, such as band gap, heterojunctions, crystal facets, and other aspects. Studying these samples by combining FTIR measurements under controlled humidity with first-principles simulations sheds light on the role and nature of the first-layer water structure in contact with the photocatalyst surfaces. It turns out that the doping hampers the percolation of tetrahedrally coordinated water molecules while enhancing the population of topological H-bond defects forming approximately linear H-bonded chains. This work unveils how doping the semiconductor surface affects the local electric field, determining the water splitting rate by influencing the H-bond topologies in the first water layers. This evidence opens new prospects for designing efficient photocatalysts for water splitting.

Stability and degradation of (oxy)nitride photocatalysts for solar water splitting.

Advancing towards alternative technologies for the sustainable production of hydrogen is a necessity for the successful integration of this potentially green fuel in the future. Photocatalytic and photoelectrochemical water splitting are promising concepts in this context. Over the past decades, researchers have successfully explored several materials classes, such as oxides, nitrides, and oxynitrides, in their quest for suitable photocatalysts with a focus on reaching higher efficiencies. However, to pave the way towards practicability, understanding degradation processes and reaching stability is essential, a domain where research has been scarcer. This perspective aims at providing an overview on recent progress concerning stability and degradation with a focus on (oxy)nitride photocatalysts and at providing insights into the opportunities and challenges coming along with the investigation of degradation processes and the attempts to improve the stability of photocatalysts.

Propagated (fragment) Pipek-Mezey Wannier functions in real-time time-dependent density functional theory.

Localization procedures are an important tool for analysis of complex systems in quantum chemistry, since canonical molecular orbitals are delocalized and can, therefore, be difficult to align with chemical intuition and obscure information at the local level of the system. This especially applies to calculations obeying periodic boundary conditions. The most commonly used approach to localization is Foster-Boys Wannier functions, which use a unitary transformation to jointly minimize the second moment of the orbitals. This procedure has proven to be robust and fast but has a side effect of often mixing σ- and π-type orbitals. σ/π-separation is achieved by the Pipek-Mezey Wannier function (PMWF) approach [Lehtola and Jónsson, J. Chem. Theory Comput. 10, 642 (2014) and Jónsson et al., J. Chem. Theory Comput. 13, 460 (2017)], which defines the spread functional in terms of partial charges instead. We have implemented a PMWF algorithm in the CP2K software package using the Cardoso-Souloumiac algorithm to enable their application to real-time time-dependent density functional theory. The method is demonstrated on stacked CO2 molecules, linear acetylenic carbon, boron and nitrogen co-doped graphene, and nitrogen-vacancy doped diamond. Finally, we discuss its computational scaling and recent efforts to improve it with fragment approaches.

The Role of Aqueous Solvation on the Intersystem Crossing of Nitrophenols.

The photochemistry of nitrophenols is a source of smog as nitrous acid is formed from their photolysis. Nevertheless, computational studies of the photochemistry of these widespread toxic molecules are scarce. In this work, the initial photodeactivation of ortho-nitrophenol and para-nitrophenol is modeled, both in gas phase and in aqueous solution to simulate atmospheric and aerosol environments. A large number of excited states, six for ortho-nitrophenol and 11 for para-nitrophenol, have been included and were all populated during the decay. Moreover, periodic time-dependent density functional theory (TDDFT) is used for both the explicitly included solvent and the solute. A comparison to periodic QM/MM (TDDFT/MM), with electrostatic embedding, is made, showing notable differences between the decays of solvated nitrophenols simulated with QM/MM and full (TD)DFT. A reduced intersystem crossing in aqueous solution could be observed thanks to the surface hopping approach using explicit, periodic TDDFT solvation including spin-orbit couplings.

Chiral Spectroscopy of Bulk Systems with Propagated Localized Orbitals.

We present approaches for the simulation of electronic circular dichroism, Raman, and Raman optical activity (ROA) spectra for isolated and periodic systems as well as subsystem analysis thereof. The method is based on the use of time-dependent maximally localized Wannier functions in the CP2K package and accounts for origin dependencies inherent to the Gaussian and plane wave with pseudopotentials approach as well as the origin dependence of the magnetic dipole and electric quadrupole operators. Tests on the H-bonded enantiomers of alanine by harmonic normal-mode analysis and on an aqueous solution of l-alanine by ab initio molecular dynamics obeying periodic boundary conditions (PBCs) are presented as total and subsystem-resolved spectra. To our knowledge, this is the first instance of an ROA spectrum derived from real-time propagation obeying PBCs and the first ROA simulation considering off-, pre-, and on-resonance effects within PBCs.

A Local Diabatisation Method for Two-State Adiabatic Conical Intersections.

A methodology to locally characterize conical intersections (CIs) between two adiabatic electronic states for which no nonadiabatic coupling (NAC) vectors are available is presented. Based on the Hessian and gradient at the CI, the branching space coordinates are identified. The potential energy surface around the CI in the branching space is expressed in the diabatic representation, from which the NAC vectors can be calculated in a wave-function-free, energy-based approach. To demonstrate the universality of the developed methodology, the minimum-energy CI (MECI) between the first (S1) and second (S2) singlet excited states of formamide is investigated at the state-averaged complete active space self-consistent field (SA-CASSCF) and extended multistate complete active space second-order perturbation theory (XMS-CASPT2) levels of theory. In addition, the asymmetrical MECI between the ground state (S0) and S1 of cyclopropanone is evaluated using SA-CASSCF, as well as (ME)CIs between the S1 and S2 states of benzene using SA-CASSCF and time-dependent density functional theory (TDDFT). Finally, a CI between the S1 and S2 excited states of thiophene was analyzed using TDDFT.

Enabling Direct Photoelectrochemical H2 Production using Alternative Oxidation Reactions on WO3.

The efficient and inexpensive conversion of solar energy into chemical bonds, such as in H2 via the photoelectrochemical splitting of H2O, is a promising route to produce green industrial feedstocks and renewable fuels, which is a key goal of the NCCR Catalysis. However, the oxidation product of the water splitting reaction, O2, has little economic or industrial value. Thus, upgrading key chemical species using alternative oxidation reactions is an emerging trend. WO3 has been identified as a unique photoanode material for this purpose since it performs poorly in the oxygen evolution reaction in H2O. Herein we highlight a collaboration in the NCCR Catalysis that has gained insights at the atomic level of the WO3 surface with ab initio computational methods that help to explain its unique catalytic activity. These computational efforts give new context to experimental results employing WO3 photoanodes for the direct photoelectrochemical oxidation of biomass-derived 5-(hydroxymethyl) furfural. While yield for the desired product, 2,5-furandicarboxylic acid is low, insights into the reaction rate constants using kinetic modelling and an electrochemical technique called derivative voltammetry, give indications on how to improve the system.

The position operator problem in periodic calculations with an emphasis on theoretical spectroscopy.

In this article, we present the challenges that arise when carrying out spectroscopic simulations within periodic boundary conditions. We present approaches which were proposed in the literature for the calculation of the extension of the electric dipole moment to periodic systems. Further, we describe the challenges arising for the simulation of magnetic properties within periodic boundary conditions and for the simulation of nuclear magnetic resonance shielding tensors and related quantities. Furthermore, issues arising in periodic implementations of vibrational circular dichroism spectroscopy are described, especially for the case of atom-centered basis functions and nuclear velocity perturbation theory.

Vibrational circular dichroism spectra of natural products by means of the nuclear velocity perturbation theory.

We present the application of the recently implemented nuclear velocity perturbation theory, using the combined Gaussian and plane waves approach in CP2K, to the vibrational circular dichroism (VCD) spectra of a set of natural products. Even though the calculations were carried out for isolated molecules in the gas-phase limit, neglecting inter-molecular interactions and anharmonic effects, the match between simulated and experimental spectra is reasonable. We also study the influence of different density functionals on the conformational search and the resulting VCD spectra via group coupling matrices (GCMs). The GCM analysis reveals that the VCD signal can in some cases arise from moieties which are close to each other and in other cases from moieties far from each other. Differences in spectra obtained using different exchange-correlation density functionals can be attributed to interaction terms between different moieties in the molecules changing their sign.

Magnetic Interactions in a [Co(II)3Er(III)(OR)4] Model Cubane through Forefront Multiconfigurational Methods.

Strong electron correlation effects are one of the major challenges in modern quantum chemistry. Polynuclear transition metal clusters are peculiar examples of systems featuring such forms of electron correlation. Multireference strategies, often based on but not limited to the concept of complete active space, are adopted to accurately account for strong electron correlation and to resolve their complex electronic structures. However, transition metal clusters already containing four magnetic centers with multiple unpaired electrons make conventional active space based strategies prohibitively expensive, due to their unfavorable scaling with the size of the active space. In this work, forefront techniques, such as density matrix renormalization group (DMRG), full configuration interaction quantum Monte Carlo (FCIQMC), and multiconfiguration pair-density functional theory (MCPDFT), are employed to overcome the computational limitation of conventional multireference approaches and to accurately investigate the magnetic interactions taking place in a [Co(II)3Er(III)(OR)4] (chemical formula [Co(II)3Er(III)(hmp)4(µ2-OAc)2(OH)3(H2O)], hmp = 2-(hydroxymethyl)-pyridine) model cubane water oxidation catalyst. Complete active spaces with up to 56 electrons in 56 orbitals have been constructed for the seven energetically lowest different spin states. Relative energies, local spin, and spin-spin correlation values are reported and provide crucial insights on the spin interactions for this model system, pivotal in the rationalization of the catalytic activity of this system in the water-splitting reaction. A ferromagnetic ground state is found with a very small, ∼50 cm-1, highest-to-lowest spin gap. Moreover, for the energetically lowest states, S = 3-6, the three Co(II) sites exhibit parallel aligned spins, and for the lower states, S = 0-2, two Co(II) sites retain strong parallel spin alignment.

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.

BODIPY-Based Photothermal Agents with Excellent Phototoxic Indices for Cancer Treatment.

Here, we report six novel, easily accessible BODIPY-based agents for cancer treatment. In contrast to established photodynamic therapy (PDT) agents, these BODIPY-based compounds show additional photothermal activity and their cytotoxicity is not dependent on the generation of reactive oxygen species (ROS). The agents show high photocytotoxicity upon irradiation with light and low dark toxicity in different cancer cell lines in 2D culture as well as in 3D multicellular tumor spheroids (MCTSs). The ratio of dark to light toxicity (phototoxic index, PI) of these agents reaches striking values exceeding 830,000 after irradiation with energetically low doses of light at 630 nm. The oxygen-dependent mechanism of action (MOA) of established photosensitizers (PSs) hampers effective clinical deployment of these agents. Under hypoxic conditions (0.2% O2), which are known to limit the efficiency of conventional PSs in solid tumors, photocytotoxicity was induced at the same concentration levels, indicating an oxygen-independent photothermal MOA. With a PI exceeding 360,000 under hypoxic conditions, both PI values are the highest reported to date. We anticipate that small molecule agents with a photothermal MOA, such as the BODIPY-based compounds reported in this work, may overcome this barrier and provide a new avenue to cancer therapy.

Automatic purpose-driven basis set truncation for time-dependent Hartree-Fock and density-functional theory.

Real-time time-dependent density-functional theory (RT-TDDFT) and linear response time-dependent density-functional theory (LR-TDDFT) are two important approaches to simulate electronic spectra. However, the basis sets used in such calculations are usually the ones designed mainly for electronic ground state calculations. In this work, we propose a systematic and robust scheme to truncate the atomic orbital (AO) basis set employed in TDDFT and TD Hartree-Fock (TDHF) calculations. The truncated bases are tested for both LR- and RT-TDDFT as well as RT-TDHF approaches, and provide an acceleration up to an order of magnitude while the shifts of excitation energies of interest are generally within 0.2 eV. The procedure only requires one extra RT calculation with 1% of the total propagation time and a simple modification on basis set file, which allows an instant application in any quantum chemistry package supporting RT-/LR-TDDFT calculations. Aside from the reduced computational effort, this approach also offers valuable insight into the effect of different basis functions on computed electronic excitations and further ideas on the design of basis sets for special purposes.

On the vibrations of formic acid predicted from first principles.

In this article, we review recent first principles, anharmonic studies on the molecular vibrations of gaseous formic acid in its monomer form. Transitions identified as fundamentals for both cis- and trans form reported in these studies are collected and supported by results from high-resolution experiments. Attention is given to the effect of coordinate coupling on the convergence of the computed vibrational states.

DeepCV: A Deep Learning Framework for Blind Search of Collective Variables in Expanded Configurational Space.

We present Deep learning for Collective Variables (DeepCV), a computer code that provides an efficient and customizable implementation of the deep autoencoder neural network (DAENN) algorithm that has been developed in our group for computing collective variables (CVs) and can be used with enhanced sampling methods to reconstruct free energy surfaces of chemical reactions. DeepCV can be used to conveniently calculate molecular features, train models, generate CVs, validate rare events from sampling, and analyze a trajectory for chemical reactions of interest. We use DeepCV in an example study of the conformational transition of cyclohexene, where metadynamics simulations are performed using DAENN-generated CVs. The results show that the adopted CVs give free energies in line with those obtained by previously developed CVs and experimental results. DeepCV is open-source software written in Python/C++ object-oriented languages, based on the TensorFlow framework and distributed free of charge for noncommercial purposes, which can be incorporated into general molecular dynamics software. DeepCV also comes with several additional tools, i.e., an application program interface (API), documentation, and tutorials.

Comparison of Length, Velocity, and Symmetric Gauges for the Calculation of Absorption and Electric Circular Dichroism Spectra with Real-Time Time-Dependent Density Functional Theory.

A velocity and symmetric gauge implementation for real-time time-dependent density functional theory (RT-TDDFT) in the CP2K package using a Gaussian and plane wave approach is presented, including the explicit gauge-transformed contributions due to the nonlocal part of pseudopotentials. Absorption spectra of gas-phase α-pinene are calculated in length and velocity gauges in the long-wavelength approximation for the application of a δ pulse in linear and full order. The velocity gauge implementation is also applied to a solvated uracil molecule to showcase its use within periodic boundary conditions (PBC). For the calculation of the expectation value of the electric dipole moment in PBC, both the velocity representation and the modern theory of polarization give equivalent absorption spectra if a distributed reference point is used for the nonlocal term of the velocity operator. The discussion of linear response theory takes place in a unified framework in terms of linear response functions in propagator notation, distinguishing the parts of the linear response functions associated with perturbation and response. To further investigate gauge dependence, electric circular dichroism (ECD) spectra of α-pinene were calculated either as magnetic response to an electric field perturbation, in length or velocity gauge, or as electric response to a magnetic field perturbation in the symmetric gauge. Both approaches, electric and magnetic perturbations, have been found to yield equivalent ECD spectra.

Synthesis of Mesodiphenylhelianthrene from 1-Aminoanthraquinone and the Structural Elucidation of Its Endoperoxide Species after Irradiation.

A safe, five-step synthetic route to yield the reliable chemical actinometer, mesodiphenylhelianthrene (MDH), is reported from a commercially available compound. Full characterization of the intermediates of the synthetic route and the final product MDH are presented together with four crystal structures of intermediates and MDH. The usage of the actinometer is described, and finally, the structure of the endoperoxide species (MDHPO), which is formed after irradiation of MDH, has been elucidated experimentally and theoretically.

Spin-Orbit Couplings for Nonadiabatic Molecular Dynamics at the ΔSCF Level.

A procedure for the calculation of spin-orbit coupling (SOC) at the delta self-consistent field (ΔSCF) level of theory is presented. Singlet and triplet excited electronic states obtained with the ΔSCF method are expanded into a linear combination of singly excited Slater determinants composed of ground electronic state Kohn-Sham orbitals. This alleviates the nonorthogonality between excited and ground electronic states and introduces a framework, similar to the auxiliary wave function at the time-dependent density functional theory (TD-DFT) level, for the calculation of observables. The ΔSCF observables of the formaldehyde system were compared to reference TD-DFT values. Our procedure gives all components (energies, gradients, nonadiabatic couplings, and SOC terms) at the ΔSCF level of theory for conducting efficient, full-atomistic nonadiabatic molecular dynamics with intersystem crossing, particularly in condensed phase systems.