The Conformational Dynamics of the Ligands Determines the Electronic Circular Dichroism of the Chiral Au38(SC2H4Ph)24 Cluster.

Effects of the conformational dynamics of 2-PET protective ligands on the electronic circular dichroism (ECD) of the chiral Au38(SC2H4Ph)24 cluster are investigated. We adopt a computational protocol in which ECD spectra are calculated via the first principle polTDDFT approach on a series of conformations extracted from MD simulations by using Essential Dynamics (ED) analysis, and then properly weighted to predict the final spectrum. We find that the experimental spectral features are well reproduced, whereas significant discrepancies arise when the spectrum is calculated using the experimental X-ray structure. This result unambiguously demonstrates the need to account for the conformational effects in the ECD modeling of chiral nanoclusters. The present procedure proved to be able of capturing the essential conformational features of the dynamic Au38(SC2H4Ph)24 system, opening the possibility to model the ECD of soluble chiral nanoclusters in a realistic way.

Time-Resolved Excited-State Analysis of Molecular Electron Dynamics by TDDFT and Bethe-Salpeter Equation Formalisms.

In this work, a theoretical and computational set of tools to study and analyze time-resolved electron dynamics in molecules, under the influence of one or more external pulses, is presented. By coupling electronic-structure methods with the resolution of the time-dependent Schrödinger equation, we developed and implemented the time-resolved induced density of the electronic wavepacket, the time-resolved formulation of the differential projection density of states (ΔPDOS), and of transition contribution map (TCM) to look at the single-electron orbital occupation and localization change in time. Moreover, to further quantify the possible charge transfer, we also defined the energy-integrated ΔPDOS and the fragment-projected TCM. We have used time-dependent density-functional theory (TDDFT), as implemented in ADF software, and the Bethe-Salpeter equation, as provided by MolGW package, for the description of the electronic excited states. This suite of postprocessing tools also provides the time evolution of the electronic states of the system of interest. To illustrate the usefulness of these postprocessing tools, excited-state populations have been computed for HBDI (the chromophore of GFP) and DNQDI molecules interacting with a sequence of two pulses. Time-resolved descriptors have been applied to study the time-resolved electron dynamics of HBDI, DNQDI, LiCN (being a model system for dipole switching upon highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) electronic excitation), and Ag22. The computational analysis tools presented in this article can be employed to help the interpretation of fast and ultrafast spectroscopies on molecular, supramolecular, and composite systems.

Hybrid theoretical models for molecular nanoplasmonics.

The multidisciplinary nature of the research in molecular nanoplasmonics, i.e., the use of plasmonic nanostructures to enhance, control, or suppress properties of molecules interacting with light, led to contributions from different theory communities over the years, with the aim of understanding, interpreting, and predicting the physical and chemical phenomena occurring at molecular- and nano-scale in the presence of light. Multiscale hybrid techniques, using a different level of description for the molecule and the plasmonic nanosystems, permit a reliable representation of the atomistic details and of collective features, such as plasmons, in such complex systems. Here, we focus on a selected set of topics of current interest in molecular plasmonics (control of electronic excitations in light-harvesting systems, polaritonic chemistry, hot-carrier generation, and plasmon-enhanced catalysis). We discuss how their description may benefit from a hybrid modeling approach and what are the main challenges for the application of such models. In doing so, we also provide an introduction to such models and to the selected topics, as well as general discussions on their theoretical descriptions.

Manipulating azobenzene photoisomerization through strong light-molecule coupling.

The formation of hybrid light-molecule states (polaritons) offers a new strategy to manipulate the photochemistry of molecules. To fully exploit its potential, one needs to build a toolbox of polaritonic phenomenologies that supplement those of standard photochemistry. By means of a state-of-the-art computational photochemistry approach extended to the strong-coupling regime, here we disclose various mechanisms peculiar of polaritonic chemistry: coherent population oscillations between polaritons, quenching by trapping in dead-end polaritonic states and the alteration of the photochemical reaction pathway and quantum yields. We focus on azobenzene photoisomerization, that encompasses the essential features of complex photochemical reactions such as the presence of conical intersections and reaction coordinates involving multiple internal modes. In the strong coupling regime, a polaritonic conical intersection arises and we characterize its role in the photochemical process. Our chemically detailed simulations provide a framework to rationalize how the strong coupling impacts the photochemistry of realistic molecules.

Excitation energy-transfer in functionalized nanoparticles: Going beyond the Förster approach.

We develop a novel approach to treat excitation energy transfer in hybrid nanosystems composed by an organic molecule attached to a semiconductor nanoparticle. Our approach extends the customary Förster theory by considering interaction between transition multipole moments of the nanoparticle at all orders and a point-like transition dipole moment representing the molecule. Optical excitations of the nanoparticle are described through an envelope-function configuration interaction method for a single electron-hole pair. We applied the method to the prototypical case of a core/shell CdSe/ZnS semiconductor quantum dot which shows a complete suppression of the energy transfer for specific transitions which could not be captured by Förster theory.

Self-assembly of mono- and bidentate oligoarylene thiols onto polycrystalline Au.

Four thiolated oligoarylene molecules (i) 4-methoxy-terphenyl-4″-methanethiol (MTM), (ii) 4-methoxy-terphenyl-3″,5″-dimethanethiol (MTD), (iii) 4-nitro-terphenyl-4″-methanethiol (NTM), and (iv) 4-nitro-terphenyl-3″,5″-dimethanethiol (NTD) were synthesized and self-assembled as monolayers (SAMs) on polycrystalline Au electrodes of organic field-effect transistors (OFETs). SAMs were characterized by contact angle and AC/DC electrochemical measurements, whereas atomic force microscopy was used for imaging the pentacene films grown on the coated electrodes. The electrical properties of functionalized OFETs, the electrochemical SAMs features and the morphology of pentacene films were correlated to the molecular organization of the thiolated oligoarylenes on Au, as calculated by means of the density functional theory. This multi-methodological approach allows us to associate the systematic replacement of the SAM anchoring head group (viz. methanethiol and dimethanethiol) and/or terminal tail group (viz. nitro-, -NO2, and methoxy, -OCH3) with the change of the electrical features. The dimethanethiol head group endows SAMs with higher resistive features along with higher surface tensions compared with methanethiol. Furthermore, the different number of thiolated heads affects the kinetics of Au passivation as well as the pentacene morphology. On the other hand, the nitro group confers further distinctive properties, such as the positive shift of both threshold and critical voltages of OFETs with respect to the methoxy one. The latter experimental evidence arise from its electron-withdrawing capability, which has been verified by both DFT calculations and DC electrochemical measurements.

Load-induced confinement activates diamond lubrication by water.

Tribochemical reactions are chemical processes, usually involving lubricant or environment molecules, activated at the interface between two solids in relative motion. They are difficult to be monitored in situ, which leaves a gap in the atomistic understanding required for their control. Here we report the real-time atomistic description of the tribochemical reactions occurring at the interface between two diamond films in relative motion, by means of large scale ab initio molecular dynamics. We show that the load-induced confinement is able to catalyze diamond passivation by water dissociative adsorption. Such passivation decreases the energy of the contacting surfaces and increases their electronic repulsion. At sufficiently high coverages, the latter prevents surface sealing, thus lowering friction. Our findings elucidate effects of the nanoscale confinement on reaction kinetics and surface thermodynamics, which are important for the design of new lubricants.

GolP: an atomistic force-field to describe the interaction of proteins with Au(111) surfaces in water.

A classical atomistic force field to describe the interaction of proteins with gold (111) surfaces in explicit water has been devised. The force field is specifically designed to be easily usable in most common bio-oriented molecular dynamics codes, such as GROMACS and NAMD. Its parametrization is based on quantum mechanical (density functional theory [DFT] and second order Möller-Plesset perturbation theory [MP2]) calculations and experimental data on the adsorption of small molecules on gold. In particular, a systematic DFT survey of the interaction between Au(111) and the natural amino acid side chains has been performed to single out chemisorption effects. Van der Waals parameters have been instead fitted to experimental desorption energy data of linear alkanes and were also studied via MP2 calculations. Finally, gold polarization (image charge effects) is taken into account by a recently proposed procedure (Iori, F.; Corni, S. J Comp Chem 2008, 29, 1656). Preliminary validation results of GolP on an independent test set of small molecules show the good performances of the force field.

Including image charge effects in the molecular dynamics simulations of molecules on metal surfaces.

Combinatorial bio-techniques have demonstrated that proteins can be good and even selective binders for several inorganic surfaces, including metals. However, the understanding of the basic physical mechanisms that govern such interactions did not keep up with the success in these experiments. The comprehension of such mechanisms would greatly benefit from the computational investigation of the problem. Because of the complexity of the system, classical molecular dynamics simulations based on an atomistic description appear to be the best compromise between reliability and feasibility. For proteins interacting with metal surfaces, however, methodological improvements with respect to standard Molecular Dynamics (MD) of proteins are needed, since the polarization of the metal induced by the protein (and the surrounding water) is not generally negligible. In this article, we present a simple approach to introduce metal polarization effects (often termed image effects) in MD simulations by exploiting standard features of bio-oriented MD codes such as the widely used GROMACS and NAMD. Tests to show the reliability of the proposed methods are presented, and the results for a model application showing the importance of image effects are also discussed.

First-principles density-functional theory calculations of electron-transfer rates in azurin dimers.

We have conceived and implemented a new method to calculate transfer integrals between molecular sites, which exploits few quantities derived from density-functional theory electronic structure computations and does not require the knowledge of the exact transition state coordinate. The method uses a complete multielectron scheme, thus including electronic relaxation effects. Moreover, it makes no use of empirical parameters. The computed electronic couplings can then be combined with estimates of the reorganization energy to evaluate electron-transfer rates that are measured in kinetic experiments: the latter are the basis to interpret electron-transfer mechanisms. We have applied our approach to the study of the electron self-exchange reaction of azurin, an electron-transfer protein belonging to the family of cupredoxins. The transfer integral estimates provided by the proposed method have been compared with those resulting from other computational techniques, from empirical models, and with available experimental data.

Electronic excitation energies of molecules in solution within continuum solvation models: investigating the discrepancy between state-specific and linear-response methods.

In a recent article [R. Cammi, S. Corni, B. Mennucci, and J. Tomasi, J. Chem. Phys. 122, 104513 (2005)], we demonstrated that the state-specific (SS) and the linear-response (LR) approaches, two different ways to calculate solute excitation energies in the framework of quantum-mechanical continuum models of solvation, give different excitation energy expressions. In particular, they differ in the terms related to the electronic response of the solvent. In the present work, we further investigate this difference by comparing the excitation energy expressions of SS and LR with those obtained through a simple model for solute-solvent systems that bypasses one of the basic assumptions of continuum solvation models, i.e., the use of a single Hartree product of a solute and a solvent wave function to describe the total solute-solvent wave function. In particular, we consider the total solute-solvent wave function as a linear combination of the four products of two solute states and two solvent electronic states. To maximize the comparability with quantum-mechanical continuum model the resulting excitation energy expression is recast in terms of response functions of the solvent and quantities proper for the solvated molecule. The comparison of the presented expressions with the LR and SS ones enlightens the physical meaning of the terms included or neglected by these approaches and shows that SS agrees with the results of the four-level model, while LR includes a term classified as dispersion in previous treatments and neglects another related to electrostatic. A discussion on the possible origin of the LR flaw is finally given.

Retention of nativelike conformation by proteins embedded in high external electric fields.

In this Communication, we show that proteins embedded in high external electric fields are capable of retaining a nativelike fold pattern. We have tested the metalloprotein azurin, immobilized onto SiO2 substrates in air with proper electrode configuration, by applying static fields up to 10(6)-10(7) Vm. The effects on the conformational properties of protein molecules have been determined by means of intrinsic fluorescence measurements. Experimental results indicate that no significant field-induced conformational alteration occurs. Such results are also discussed and supported by theoretical predictions of the inner protein fields.

Electronic excitation energies of molecules in solution: state specific and linear response methods for nonequilibrium continuum solvation models.

We present a formal comparison between the two different approaches to the calculation of electronic excitation energies of molecules in solution within the continuum solvation model framework, taking also into account nonequilibrium effects. These two approaches, one based on the explicit evaluation of the excited state wave function of the solute and the other based on the linear response theory, are here proven to give formally different expressions for the excitation energies even when exact eigenstates are considered. Calculations performed for some illustrative examples show that this formal difference has sensible effects on absolute solvatochromic shifts (i.e., with respect to gas phase) while it has small effects on relative (i.e., nonpolar to polar solvent) solvatochromic shifts.