Enrico Berardo

Modeling Excited States in TiO2 Nanoparticles: On the Accuracy of a TD-DFT Based Description.

We have investigated the suitability of Time-Dependent Density Functional Theory (TD-DFT) to describe vertical low-energy excitations in naked and hydrated titanium dioxide nanoparticles. Specifically, we compared TD-DFT results obtained using different exchange-correlation (XC) potentials with those calculated using Equation-of-Motion Coupled Cluster (EOM-CC) quantum chemistry methods. We demonstrate that TD-DFT calculations with commonly used XC potentials (e.g., B3LYP) and EOM-CC methods give qualitatively similar results for most TiO2 nanoparticles investigated. More importantly, however, we also show that, for a significant subset of structures, TD-DFT gives qualitatively different results depending upon the XC potential used and that only TD-CAM-B3LYP and TD-BHLYP calculations yield results that are consistent with those obtained using EOM-CC theory. Moreover, we demonstrate that the discrepancies for such structures originate from a particular combination of defects that give rise to charge-transfer excitations, which are poorly described by XC potentials that do not contain sufficient Hartree–Fock like exchange. Finally, we consider that such defects are readily healed in the presence of ubiquitously present water and that, as a result, the description of vertical low-energy excitations for hydrated TiO2 nanoparticles is nonproblematic.

Amine Molecular Cages as Supramolecular Fluorescent Explosive Sensors: A Computational Perspective

We investigate using a computational approach the physical and chemical processes underlying the application of organic (macro)molecules as fluorescence quenching sensors for explosives sensing. We concentrate on the use of amine molecular cages to sense nitroaromatic analytes, such as picric acid and 2,4-dinitrophenol, through fluorescence quenching. Our observations for this model system hold for many related systems. We consider the different possible mechanisms of fluorescence quenching: Förster resonance energy transfer, Dexter energy transfer and photoinduced electron transfer, and show that in the case of our model system, the fluorescence quenching is driven by the latter and involves stable supramolecular sensor-analyte host-guest complexes. Furthermore, we demonstrate that the experimentally observed selectivity of amine molecular cages for different explosives can be explained by the stability of these host-guest complexes and discuss how this is related to the geometry of the binding site in the sensor. Finally, we discuss what our observations mean for explosive sensing by fluorescence quenching in general and how this can help in future rational design of new supramolecular detection systems.

DFT modeling of 45S5 and 77S soda-lime phospho-silicate glass surfaces: clues on different bioactivity mechanism.

The reactivity of bioglasses, which is related to the dissolution of cations and orthosilicate groups in the physiological fluid, strongly depends on the key structural features present at the glass surfaces. On the basis of the composition and the synthetic routes employed to make the glass, surfaces with very different characteristics and thus presenting different mechanisms of dissolution can be observed. In this paper, the surface structures of two very different bioglass compositions, namely 45S5 (46.1 SiO2, 24.4 Na2O, 26.9 CaO, and 2.6 P2O5 mol %) and 77S (80.0 SiO2, 16.0 CaO, and 4.0 P2O5 mol %), have been investigated by means of periodic DFT calculations based on a PBE functional and localized Gaussian basis set as encoded in the CRYSTAL code. Our calculations show that the two glass surfaces differ by the relative amount of key structural sites such as NBOs, exposed ions, orthosilicate units, and small rings. We have demonstrated how the number of these sites affects the surface stability and reactivity (bioactivity).

Modelling materials for solar fuel synthesis by artificial photosynthesis; predicting the optical, electronic and redox properties of photocatalysts

In this mini-review, we discuss what insight computational modelling can provide into the working of photocatalysts for solar fuel synthesis and how calculations can be used to screen for new promising materials for photocatalytic water splitting and carbon dioxide reduction. We will extensively discuss the different relevant (material) properties and the computational approaches (DFT, TD-DFT, GW/BSE) available to model them. We illustrate this with examples from the literature, focussing on polymeric and nanoparticle photocatalysts. We finish with a perspective on the outstanding conceptual and computational challenges.

Probing the fate of interstitial water in bulk bioactive glass by ab initio simulations

As also observed for conventional silicate glasses, water can be incorporated in the bulk interstitial regions of a bioactive glass (BG) matrix during the glass preparation and/or upon exposure to an aqueous environment. However, in the case of BGs, very little is known about the effect of hydration on the bulk structure, and then on key properties of these materials, such as biodegradation and bioactivity, that depend on the bulk structure itself. Here we employ a combination of atomistic simulation techniques to explore the nature and effects of water–BG interactions in the bulk of a bioactive glass. The fate of water inserted in the bulk interstitial region of 45S5 bioglass has been studied by ab initio geometry relaxations and Molecular Dynamics (AIMD) simulations. We probed the interaction of a water molecule with silica rings and cages of different size, as well as the stability of potentially relevant configurations involving manually dissociated water and opened rings. The local stability of selected configurations was further assessed by subjecting them to AIMD runs, in order to overcome possible kinetic barriers for water diffusion and dissociation. Small rings do not appear as favourable absorption sites in the bulk of a bioactive glass as they are for bioinert glasses. Moreover, water dissociation through rupture of Si–O bonds of silica rings formed in the bulk was thermodynamically unfavourable. However, a high-temperature AIMD run led to a dissociated state involving no broken Si–O bonds and a free hydroxyl: because re-optimization of this state produced the most favourable hydration energy identified in this study, dissociative absorption through this mechanism appears a likely outcome of the water–45S5 interaction at low water content. We discuss the structural and dynamical basis for the stability of this and other water–glass adducts identified, and the potential consequences of these interactions for the behaviour of the glass in a biological context.

stk: A python toolkit for supramolecular assembly

A tool for the automated assembly, molecular optimization and property calculation of supramolecular materials is presented. stk is a modular, extensible and open‐source Python library that provides a simple Python API and integration with third party computational codes. stk currently supports the construction of linear polymers, small linear oligomers, organic cages in multiple topologies and covalent organic frameworks (COFs) in multiple framework topologies, but is designed to be easy to extend to new, unrelated, supramolecules or new topologies. Extension to metal–organic frameworks (MOFs), metallocycles or supramolecules, such as catenanes, would be straightforward. Through integration with third party codes, stk offers the user the opportunity to explore the potential energy landscape of the assembled supramolecule and then calculate the supramolecules structural features and properties. stk provides support for high‐throughput screening of large batches of supramolecules at a time. The source code of the program can be found at https://github.com/supramolecular-toolkit/stk.

Maximising the hydrogen evolution activity in organic photocatalysts by co-polymerisation

The hydrogen evolution activity of a polymeric photocatalyst was maximised by co-polymerisation, using both experimental and computational screening, for a family of 1,4-phenylene/2,5-thiophene co-polymers. The photocatalytic activity is the product of multiple material properties that are affected in different ways by the polymer composition and microstructure. For the first time, the photocatalytic activity was shown to be a function of the arrangement of the building blocks in the polymer chain as well as the overall composition. The maximum in hydrogen evolution for the co-polymer series appears to result from a trade-off between the fraction of light absorbed and the thermodynamic driving force for proton reduction and sacrificial electron donor oxidation, with the co-polymer of p-terphenyl and 2,5-thiophene showing the highest activity.

High-Throughput Screening Approach for the Optoelectronic Properties of Conjugated Polymers

We propose a general high-throughput virtual screening approach for the optical and electronic properties of conjugated polymers. This approach makes use of the recently developed xTB family of low-computational-cost density functional tight-binding methods from Grimme and co-workers, calibrated here to (Time-Dependent) Density Functional Theory ((TD)DFT) data computed for a representative diverse set of (co)polymers. Parameters drawn from the resulting calibration using a linear model can then be applied to the xTB derived results for new polymers, thus generating near DFT-quality data with orders of magnitude reduction in computational cost. As a result, after an initial computational investment for calibration, this approach can be used to quickly and accurately screen on the order of thousands of polymers for target applications. We also demonstrate that the (opto)electronic properties of the conjugated polymers show only a very minor variation when considering different conformers and that the results of high-throughput screening are therefore expected to be relatively insensitive with respect to the conformer search methodology applied.