Grand Canonical Ensemble Approaches in CP2K for Modeling Electrochemistry at Constant Electrode Potentials.

In electrochemical experiments, the number of electrons of the electrode immersed in the electrolyte is usually variable. Additionally, the numbers of adsorbed substances on the surface of the electrode, the solvent molecules, and counter charge ions in the near-surface region can also vary. Treating electrochemical solid-liquid interfaces with the typical fixed electron number density functional theory (DFT) approach tends to be a challenge. This can be addressed by using grand canonical ensemble approaches. We present the implementation of two grand canonical ensemble approaches in the open-source computational chemistry software CP2K that go beyond the existing canonical ensemble paradigm. The first approach is based on implicit solvent models and explicit atomistic solute (electrode with/without adsorbed species) models, and includes two recent developments: (a) grand canonical self-consistent field (GC-SCF) method (J. Chem. Phys. 2017, 146, 114104) allowing the electron number of the system to fluctuate naturally and accordingly with the experimental electrode potential, (b) planar counter charge (J. Chem. Phys. 2019, 150, 041722, Phys. Rev. B 2003, 68, 245416) salt model completely screening the net charge of the electrode model. In contrast with previous studies, in our implementation, the work function (WF) (absolute electrode potential if the potential drop at the electrolyte-vacuum interface is omitted) is the constrained quantity during an SCF optimization instead of the Fermi energy. The chemical potential of electrons (negative WF) is a natural variable of the grand potential in the GC ensemble of electronic states, and this method can easily achieve stable SCF convergence and obtain an electronic structure that precisely corresponds to a user-specified WF. The second approach referred to as the GC DFT molecular dynamics (DFT-MD) simulation scheme (Phys. Rev. Lett. 2002, 88, 213002, J. Chem. Phys. 2005, 122, 234505, J. Am. Chem. Soc. 2004, 126 (12), 3928-3938) is based on fully explicit modeling the solvent molecules and the ions and is used to calculate the electron chemical potential corresponding to an equilibrium electrochemical half-reaction (M(n+m)+ + ne- ⇌ Mm+) which involves DFT-MD, by allowing the number of electrons to vary during the DFT-MD simulation process. This opens the way for forefront electrochemical calculations in CP2K for a broad range of systems.

Machine Learning Interatomic Potentials for Catalysis.

Atomistic modeling can provide insights into the design of novel catalysts in modern industries of chemistry, materials science, and biology. Classical force fields and ab initio calculations have been widely adopted in molecular simulations. How- ever, these methods suffer from the drawbacks of either low accuracy or high cost. Recently, the development of machine learning interatomic potentials (MLIPs) has become more and more popular as they can tackle the problems in question and can deliver rather accurate results at significantly lower computational cost. In this review, the atomistic modeling of catalytic systems with the aid of MLIPs is discussed, showcasing recently developed MLIP models and selected applications for the modeling of catalytic systems. We also highlight the best practices, and challenges for MLIPs and give an outlook for future works on MLIPs in the field of catalysis.

Water-air interface revisited by means of path-integral ab initio molecular dynamics.

Although nuclear quantum effects (NQEs) have been considered on bulk liquid water, the impact of these latter on the air-water interface has not yet been reported. Herein, by performing and comparing ab initio molecular dynamics (AIMD) and path integral AIMD (PI-AIMD) simulations, we reveal the impact of NQEs on structural, dynamical and electronic properties as well as IR spectra of the air-water interface at room temperature. NQEs, being able to describe a more accurate proton delocalization in H-bonded system than AIMD, reveal a different structural arrangement and dynamical behaviour of both bulk and interfacial water molecules in comparison to AIMD results. A more de-structured and de-bound water arrangement and coordination are identified when the quantum nature of nuclei are considered for both bulk and interfacial water molecules. Structural properties, such as inter-/intra-molecular bond lengths, coordination numbers and H-bonding angles of bulk and interfacial water molecules here calculated, are affected by NQEs mitigating the overstructured description given by AIMD. Further evidences of an AIMD overstructured description of bulk water are in the computed IR spectra, where an increased absorption peak intensity and an increased strength of the hydrogen-bond network are alleviated by NQEs. In addition, NQEs show a valuable impact on the electronic structure of the air-water interface, reducing the total valence bandwidth and the electronic energy band-gap when passing from bulk to interfacial water. This work proves how NQEs significantly affect properties and features of the air-water interface, that are essential to accurately describe H-bonded interfacial systems.

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.

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.

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.

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.

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.

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.

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.

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.