Toward a Unified Modeling of Environment and Bridge-Mediated Contributions to Electronic Energy Transfer: A Fully Polarizable QM/MM/PCM Approach

A Vision for the Future of Multiscale Modeling

The rise of modern computer science enabled physical chemistry to make enormous progresses in understanding and harnessing natural and artificial phenomena. Nevertheless, despite the advances achieved over past decades, computational resources are still insufficient to thoroughly simulate extended systems from first principles. Indeed, countless biological, catalytic and photophysical processes require ab initio treatments to be properly described, but the breadth of length and time scales involved makes it practically unfeasible. A way to address these issues is to couple theories and algorithms working at different scales by dividing the system into domains treated at different levels of approximation, ranging from quantum mechanics to classical molecular dynamics, even including continuum electrodynamics. This approach is known as multiscale modeling and its use over the past 60 years has led to remarkable results. Considering the rapid advances in theory, algorithm design, and computing power, we believe multiscale modeling will massively grow into a dominant research methodology in the forthcoming years. Hereby we describe the main approaches developed within its realm, highlighting their achievements and current drawbacks, eventually proposing a plausible direction for future developments considering also the emergence of new computational techniques such as machine learning and quantum computing. We then discuss how advanced multiscale modeling methods could be exploited to address critical scientific challenges, focusing on the simulation of complex light-harvesting processes, such as natural photosynthesis. While doing so, we suggest a cutting-edge computational paradigm consisting in performing simultaneous multiscale calculations on a system allowing the various domains, treated with appropriate accuracy, to move and extend while they properly interact with each other. Although this vision is very ambitious, we believe the quick development of computer science will lead to both massive improvements and widespread use of these techniques, resulting in enormous progresses in physical chemistry and, eventually, in our society.

Polarizable Embedding without Artificial Boundary Polarization

We present a fully self-consistent polarizable embedding (PE) model that does not suffer from unphysical boundary polarization. This is achieved through the use of the minimum-image convention (MIC) in the induced electrostatics. It is a simple yet effective approach that includes a more physically accurate description of the polarization throughout the molecular system. Using PE with MIC (PE-MIC), we shed new light on the limitations of commonly employed cutoff models, such as the droplet model, when used in PE calculations. Specifically, we investigate the effects of the unphysical polarization at the outer boundary by comparing induced dipoles and the associated electrostatic potentials, as well as some optical properties of solute-solvent and biomolecular systems. We show that the magnitude of the inaccuracies caused by the unphysical polarization depends on multiple parameters: the nature of the quantum subsystem and of the environment, the cutoff model and distance, and the calculated property.

The Influence of Electronic Polarization on Nonlinear Optical Spectroscopy

The environment surrounding a chromophore can dramatically affect the energy absorption and relaxation process, as manifested in optical spectra. Simulations of nonlinear optical spectroscopy, such as two-dimensional electronic spectroscopy (2DES) and transient absorption (TA), will be influenced by the computational model of the environment. We here compare a fixed point charge molecular mechanics model and a quantum mechanical (QM) model of the environment in computed 2DES and TA spectra of Nile red in water and the chromophore of photoactive yellow protein (PYP) in water and protein environments. In addition to simulating these nonlinear optical spectra, we directly juxtapose the computed excitation energy correlation function to the dynamic Stokes shift function often used to analyze environment dynamics. Overall, we find that for the three systems studied here the mutual electronic polarization provided by the QM environment manifests in broader 2DES signals, as well as a larger reorganization energy and a larger static Stokes shift due to stronger coupling between the chromophore and the environment.

Energy, Structures, and Response Properties with a Fully Coupled QM/AMOEBA/ddCOSMO Implementation

We present the implementation of a fully coupled polarizable QM/MM/continuum model based on the AMOEBA polarizable force field and the domain decomposition implementation of the conductor-like screening model. Energies, response properties, and analytical gradients with respect to both QM and MM nuclear positions are available, and a generic, atomistic cavity can be employed. The model is linear scaling in memory requirements and computational cost with respect to the number of classical atoms and is therefore suited to model large, complex systems. Using three variants of the green-fluorescent protein, we investigate the overall computational cost of such calculations and the effect of the continuum model on the convergence of the computed properties with respect to the size of the embedding. We also demonstrate the fundamental role of polarization effects by comparing polarizable and nonpolarizable embeddings to fully QM ones.

Vibronic and Environmental Effects in Simulations of Optical Spectroscopy

Including both environmental and vibronic effects is important for accurate simulation of optical spectra, but combining these effects remains computationally challenging. We outline two approaches that consider both the explicit atomistic environment and the vibronic transitions. Both phenomena are responsible for spectral shapes in linear spectroscopy and the electronic evolution measured in nonlinear spectroscopy. The first approach utilizes snapshots of chromophore-environment configurations for which chromophore normal modes are determined. We outline various approximations for this static approach that assumes harmonic potentials and ignores dynamic system-environment coupling. The second approach obtains excitation energies for a series of time-correlated snapshots. This dynamic approach relies on the accurate truncation of the cumulant expansion but treats the dynamics of the chromophore and the environment on equal footing. Both approaches show significant potential for making strides toward more accurate optical spectroscopy simulations of complex condensed phase systems.

Resonance energy transfer mediated by a chiral molecule

The problem of resonant energy transfer (RET) between an electric dipole donor, D, and an electric dipole acceptor, A, mediated by a passive, chiral third-body, T, is considered within the framework of molecular quantum electrodynamics theory. To account for the optical activity of the mediator, magnetic dipole and electric quadrupole coupling terms are included in addition to the leading electric dipole interaction term. Fourth-order diagrammatic time-dependent perturbation theory is used to obtain the matrix element. It is found that the Fermi golden rule rate depends on pure multipole moment polarizabilities and susceptibilities of T, as well as on various mixed electric and magnetic multipole moment response functions. The handedness of T manifests through mixed electric-magnetic dipole and mixed electric dipole-quadrupole polarizabilities, which affect the rate and, respectively, require the use of fourth-rank and sixth-rank Cartesian tensor averages over T, yielding non-vanishing isotropic rate formulae applicable to a chiral fluid medium. Terms of a similar order of magnitude proportional to the product of electric dipole polarizability and either magnetic dipole susceptibility or electric quadrupole polarizability of T are also computed for oriented and freely tumbling molecules. Migration rates dependent upon the product of the pure electric dipole or magnetic dipole polarizability with the mixed electric-magnetic or electric dipole-quadrupole analogs, which require fourth- and fifth-rank Cartesian tensor averaging, vanish for randomly oriented systems. Asymptotically limiting rate expressions are also evaluated. Insight is gained into RET occurring in complex media.

Bridge-Mediated RET between Two Chiral Molecules

Molecular quantum electrodynamics (QED) theory is employed to calculate the rate of resonance energy transfer (RET) between a donor, D, described by an electric dipole and quadrupole, and magnetic dipole coupling, and an identical acceptor molecule, A, that is mediated by a third body, T, which is otherwise inert. A single virtual photon propagates between D and T, and between T and A. Time-dependent perturbation theory is used to compute the matrix element, from which the transfer rate is evaluated using the Fermi golden rule. This extends previous studies that were limited to the electric dipole approximation only and admits the possibility of the exchange of excitation between a chiral emitter and absorber. Rate terms are computed for specific pure and mixed multipole-dependent contributions of D and A for both an oriented arrangement of the three particles and for the freely tumbling situation. Mixed multipole moment contributions, such as those involving electric–magnetic dipole or electric dipole–quadrupole coupling at one center, do not survive random orientational averaging. Interestingly, the mixed electric–magnetic dipole D and A rate term is non-vanishing and discriminatory, exhibiting a dependence on the chirality of the emitter and absorber, and is entirely retarded. It vanishes, however, if D and A are oriented perpendicularly to one another. Near- and far-zone asymptotes of isotropic contributions to the rate are also evaluated, demonstrating radiationless short-range transfer and inverse-square radiative exchange at very large separations.

Theoretical investigation of non-Förster exciton transfer mechanisms in perylene diimide donor, phenylene bridge, and terrylene diimide acceptor systems

The rates of exciton transfer within dyads of perylene diimide and terrylene diimide connected by oligophenylene bridge units have been shown to deviate significantly from those of Förster's resonance energy transfer theory, according to single molecule spectroscopy experiments. The present work provides a detailed computational and theoretical study investigating the source of such a discrepancy. Electronic spectroscopy data are calculated by time-dependent density functional theory and then compared with experimental results. Electronic couplings between the exciton donor and the acceptor are estimated based on both the transition density cube method and transition dipole approximation. These results confirm that the delocalization of the exciton to the bridge parts contributes to significant enhancement of donor-acceptor electronic coupling. Mechanistic details of exciton transfer are examined by estimating the contributions of the bridge electronic states, vibrational modes of the dyads commonly coupled to both donor and acceptor, inelastic resonance energy transfer mechanism, and dark exciton states. These analyses suggest that the contribution of common vibrational modes serves as the main source of deviation from Förster's spectral overlap expression.

Ultrafast Spectroscopy of Photoactive Molecular Systems from First Principles: Where We Stand Today and Where We Are Going

Computational spectroscopy is becoming a mandatory tool for the interpretation of the complex, and often congested, spectral maps delivered by modern non-linear multi-pulse techniques. The fields of Electronic Structure Methods, Non-Adiabatic Molecular Dynamics, and Theoretical Spectroscopy represent the three pillars of the virtual ultrafast optical spectrometer, able to deliver transient spectra in silico from first principles. A successful simulation strategy requires a synergistic approach that balances between the three fields, each one having its very own challenges and bottlenecks. The aim of this Perspective is to demonstrate that, despite these challenges, an impressive agreement between theory and experiment is achievable now regarding the modeling of ultrafast photoinduced processes in complex molecular architectures. Beyond that, some key recent developments in the three fields are presented that we believe will have major impacts on spectroscopic simulations in the very near future. Potential directions of development, pending challenges, and rising opportunities are illustrated.

Polariton mediated resonance energy transfer in a fluid

The focus of this work is on a microscopic quantum electrodynamical understanding of cumulative quantum effects in resonance energy transfer occurring in an isotropic and disordered medium. In particular, we consider quantum coherence, defined in terms of interferences between Feynman pathways, and analyze pure-amplitude and phase cross terms that appear in the Fermi golden rule rate equation that results from squaring the matrix element for mediated energy transfer. It is shown that pure-amplitude terms dominate in the near-zone when chromophores are close in proximity to one another (within a few nanometers), and phase cross terms dominate toward the far-zone when phase differences between different Feynman pathways begin to emerge. This can be understood in terms of physical attributes of the mediating photon, whose character becomes more real at long distances, coinciding with vanishing longitudinal components of the field, as transverse components begin to dominate.

Polarizable embedding QM/MM: the future gold standard for complex (bio)systems?

We provide a perspective of the induced dipole formulation of polarizable QM/MM, showing how efficient implementations will enable their application to the modeling of dynamics, spectroscopy, and reactivity in complex biosystems.

A Quantum Electrodynamics Description of Quantum Coherence and Damping in Condensed-Phase Energy Transfer

Quantum coherence in condensed-phase electronic resonance energy transfer (RET) is described within the context of quantum electrodynamics (QED) theory. Mediating dressed virtual photons (polaritons) are explicitly incorporated into the treatment, and coherence is understood within the context of interfering Feynman pathways connecting the initial and final states for the RET process. The model investigated is that of an oriented three-body donor, acceptor, and mediator RET system embedded within a dispersive and absorbing polarizable medium. We show how quantum coherence can significantly enhance the rate of RET and give a rigorous picture for subsequent decoherence that is driven by both phase and amplitude damping. Energy-conserving phase damping occurs as a result of geometric and dispersive effects and is associated with destructive interference between Feynman pathways. Dissipative amplitude damping, on the other hand, is attributed to vibronic relaxation and absorptivity of the medium and can be understood as virtual photons (polaritons) leaking into the environment. This model offers insights into the emergence of coherence and subsequent decoherence for energy transfer in photosynthetic systems.

Quantum Chemical Modeling of the Photoinduced Activity of Multichromophoric Biosystems

Multichromophoric biosystems represent a broad family with very diverse members, ranging from light-harvesting pigment-protein complexes to nucleic acids. The former are designed to capture, harvest, efficiently transport, and transform energy from sunlight for photosynthesis, while the latter should dissipate the absorbed radiation as quickly as possible to prevent photodamages and corruption of the carried genetic information. Because of the unique electronic and structural characteristics, the modeling of their photoinduced activity is a real challenge. Numerous approaches have been devised building on the theoretical development achieved for single chromophores and on model Hamiltonians that capture the essential features of the system. Still, a question remains: is a general strategy for the accurate modeling of multichromophoric systems possible? By using a quantum chemical point of view, here we review the advancements developed so far highlighting differences and similarities with the single chromophore treatment. Finally, we outline the important limitations and challenges that still need to be tackled to reach a complete and accurate picture of their photoinduced properties and dynamics.

Coupled Cluster Theory with Induced Dipole Polarizable Embedding for Ground and Excited States

In this work, we present the theory and implementation of the coupled cluster single and double excitations (CCSD) method combined with a classical polarizable molecular mechanics force field (MMPol) based on the induced dipole model. The method is developed to compute electronic excitation energies within the state specific (SS) and linear response (LR) formalisms for the interaction of the quantum mechanical and classical regions. Furthermore, we consider an approximate expression of the correlation energy, originally developed for CCSD with implicit solvation models, where the interaction term is linear in the coupled cluster density. This approximation allows us to include the explicit contribution of the environment to the CC equations without increasing the computational effort. The test calculations on microsolvated systems, where the CCSD/MMPol method is compared to full CCSD calculations, demonstrates the reliability of this computational protocol for all interaction schemes (errors < 2%). We also show that it is important to include induced dipoles on all atom centers of the classical region and that too diffuse functions in the basis set may be problematic due to too strong interaction with the environment.

Towards large scale hybrid QM/MM dynamics of complex systems with advanced point dipole polarizable embeddings

Hybrid DFT(Gaussian)/AMOEBA(Tinker/Tinker-HP) polarizable molecular dynamics including the QM/MM mutual polarization on large complex systems. Example of the thiazole orange dye buried in a DNA double helix, embedded in a sphere of water (16 500 atoms).

In this work, we present a general route to hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) Molecular Dynamics for complex systems using a polarizable embedding. We extend the capabilities of our hybrid framework, combining the Gaussian and Tinker/Tinker-HP packages in the context of the AMOEBA polarizable force field to treat large (bio)systems where the QM and the MM subsystems are covalently bound, adopting pseudopotentials at the boundaries between the two regions. We discuss in detail the implementation and demonstrate the global energy conservation of our QM/MM Born–Oppenheimer molecular dynamics approach using Density Functional Theory. Finally, the approach is assessed on the electronic absorption properties of a 16 500 atom complex encompassing an organic dye embedded in a DNA matrix in solution, extending the hybrid method to a time-dependent Density Functional Theory approach. The results obtained comparing different partitions between the quantum and the classical subsystems also suggest that large QM portions are not necessary if accurate polarizable force fields are used in a variational formulation of the embedding, properly including the QM/MM mutual polarization.

Mutually polarizable QM/MM model with in situ optimized localized basis functions

We extend our recently developed quantum-mechanical/molecular mechanics (QM/MM) approach [Dziedzic et al., J. Chem. Phys. 145, 124106 (2016)] to enable in situ optimization of the localized orbitals. The quantum subsystem is described with onetep linear-scaling density functional theory and the classical subsystem - with the AMOEBA polarizable force field. The two subsystems interact via multipolar electrostatics and are fully mutually polarizable. A total energy minimization scheme is employed for the Hamiltonian of the coupled QM/MM system. We demonstrate that, compared to simpler models using fixed basis sets, the additional flexibility offered by in situ optimized basis functions improves the accuracy of the QM/MM interface, but also poses new challenges, making the QM subsystem more prone to overpolarization and unphysical charge transfer due to increased charge penetration. We show how these issues can be efficiently solved by replacing the classical repulsive van der Waals term for QM/MM interactions with an interaction of the electronic density with a fixed, repulsive MM potential that mimics Pauli repulsion, together with a modest increase in the damping of QM/MM polarization. We validate our method, with particular attention paid to the hydrogen bond, in tests on water-ion pairs, the water dimer, first solvation shells of neutral and charged species, and solute-solvent interaction energies. As a proof of principle, we determine suitable repulsive potential parameters for water, K⁺, and Cl^-. The mechanisms we employed to counteract the unphysical overpolarization of the QM subsystem are demonstrated to be adequate, and our approach is robust. We find that the inclusion of explicit polarization in the MM part of QM/MM improves agreement with fully QM calculations. Our model permits the use of minimal size QM regions and, remarkably, yields good energetics across the well-balanced QM/MM interface.

QM/MM methods for free energies and photochemistry

Hybrid computational methods describing a small region of a biomolecular system with quantum mechanics and the bulk with molecular mechanics, referred to as QM/MM methods, are now a central part of computational biochemistry. This review considers developments in the QM/MM approach that make it easier to calculate free energies using accurate QM-based potential energy expressions. We also describe techniques to treat electronic coupling between the core region and the MM environment. Polarizability of the protein matrix is important but so is electronic coupling. Applications of these new methods, especially to photochemistry, are discussed.

Is energy transfer limiting multiphotochromism? answers from ab initio quantifications

Dithienylethenes (DTEs) can be assembled to form supramolecular multiphotochromic systems that are highly functional molecular architectures of potential interest for building complex optoelectronic devices. Yet even simple DTE dimers relying on an organic linker may suffer from a partial photoactivity, i.e., only one of the two switches does isomerise. Experimentally, this limited photochromism has been attributed to an excited state energy transfer (EET) between the two DTEs of the multimer; this EET taking place instead of the desired photoinduced cyclisation of the DTE. However, no clear evidences of this phenomenon have been provided so far. In this work we propose the first rationalisation of this potential parasite photoinduced event using a computational approach based on Time-Dependent Density Functional Theory (TD-DFT) for the calculation of the electronic coupling in DTE dimers. Besides quantifying EET in several systems, we dissect the role of through-bond and through-space mechanisms on this process and clarify their dependence on both the nature and length of the bridge separating the two photochromes. The theoretical data obtained in this framework are in full agreement with the experimental outcomes and pave the way toward a molecular design of coupled, yet fully functionals, DTE-based multiswitches.

Open-ended formulation of self-consistent field response theory with the polarizable continuum model for solvation

The study of high-order absorption properties of molecules is a field of growing importance. Quantum-chemical studies can help design chromophores with desirable characteristics. Given that most experiments are performed in solution, it is important to devise a cost-effective strategy to include solvation effects in quantum-chemical studies of these properties. We here present an open-ended formulation of self-consistent field (SCF) response theory for a molecular solute coupled to a polarizable continuum model (PCM) description of the solvent. Our formulation relies on the open-ended, density matrix-based quasienergy formulation of SCF response theory of Thorvaldsen, et al., [J. Chem. Phys., 2008, 129, 214108] and the variational formulation of the PCM, as presented by Lipparini et al., [J. Chem. Phys., 2010, 133, 014106]. Within the PCM approach to solvation, the mutual solute-solvent polarization is represented by means of an apparent surface charge (ASC) spread over the molecular cavity defining the solute-solvent boundary. In the variational formulation, the ASC is an independent, variational degree of freedom. This allows us to formulate response theory for molecular solutes in the fixed-cavity approximation up to arbitrary order and with arbitrary perturbation operators. For electric dipole perturbations, pole and residue analyses of the response functions naturally lead to the identification of excitation energies and transition moments. We document the implementation of this approach in the Dalton program package using a recently developed open-ended response code and the PCMSolver libraries and present results for one-, two-, three-, four- and five-photon absorption processes of three small molecules in solution.

Modelling excitation energy transfer in covalently linked molecular dyads containing a BODIPY unit and a macrocycle

We model the singlet–singlet Excitation Energy Transfer (EET) process in a panel of large BODIPY–macrocycle dyads, including some azacalixphyrin derivatives.

Electron and excitation energy transfers in covalently linked donor–acceptor dyads: mechanisms and dynamics revealed using quantum chemistry

A quantum chemical strategy is developed to study photoinduced electron and energy transfer processes in covalently linked dyads. The strategy is applied to two zinc and free-based porphyrin–naphthalenediimide dyads in solution.

Performance of the AMOEBA Water Model in the Vicinity of QM Solutes: A Diagnosis Using Energy Decomposition Analysis

The importance of incorporating solvent polarization effects into the modeling of solvation processes has been well-recognized, and therefore a new generation of hybrid quantum mechanics/molecular mechanics (QM/MM) approaches that accounts for this effect is desirable. We present a fully self-consistent, mutually polarizable QM/MM scheme using the AMOEBA force field, in which the total energy of the system is variationally minimized with respect to both the QM electronic density and the MM induced dipoles. This QM/AMOEBA model is implemented through the Q-Chem/LibEFP code interface and then applied to the evaluation of solute-solvent interaction energies for various systems ranging from the water dimer to neutral and ionic solutes (NH₃, NH₄⁺, CN^-) surrounded by increasing numbers of water molecules (up to 100). In order to analyze the resulting interaction energies, we also utilize an energy decomposition analysis (EDA) scheme which identifies contributions from permanent electrostatics, polarization, and van der Waals (vdW) interaction for the interaction between the QM solute and the solvent molecules described by AMOEBA. This facilitates a component-wise comparison against full QM calculations where the corresponding energy components are obtained via a modified version of the absolutely localized molecular orbitals (ALMO)-EDA. The results show that the present QM/AMOEBA model can yield reasonable solute-solvent interaction energies for neutral and cationic species, while further scrutiny reveals that this accuracy highly relies on the delicate balance between insufficiently favorable permanent electrostatics and softened vdW interaction. For anionic solutes where the charge penetration effect becomes more pronounced, the QM/MM interface turns out to be unbalanced. These results are consistent with and further elucidate our findings in a previous study using a slightly different QM/AMOEBA model ( Dziedzic et al. J. Chem. Phys. 2016 , 145 , 124106 ). The implications of these results for further refinement of this model are also discussed.

Modeling excitation energy transfer in multi-BODIPY architectures

The excitation energy transfer (EET) allowing the concentration of the energy has been investigated in several multi-BODIPY architectures with the help of an approach coupling time dependent density functional theory to an implicit solvation scheme, the polarizable continuum model. We have first considered several strategies to compute the electronic coupling in a dyad varying the size of the donor/acceptor units, the bridge, the geometries and conformations. We have next studied the electronic coupling in three different architectures for which the EET rate constants have been experimentally measured both from luminescence and transient absorption data and from Förster theory. A good agreement with experimental values was obtained. Finally, in an effort to further improve these systems, we have designed several series of BODIPY triads, investigating the effect of acidochromism, core modifications, the position of the linkage and chemical substitutions on the EET coupling and rate constant. We show that several architectures allow us to increase the EET rate by one order of magnitude compared to the original compound.

Theoretical description of protein field effects on electronic excitations of biological chromophores

Photoinitiated phenomena play a crucial role in many living organisms. Plants, algae, and bacteria absorb sunlight to perform photosynthesis, and convert water and carbon dioxide into molecular oxygen and carbohydrates, thus forming the basis for life on Earth. The vision of vertebrates is accomplished in the eye by a protein called rhodopsin, which upon photon absorption performs an ultrafast isomerisation of the retinal chromophore, triggering the signal cascade. Many other biological functions start with the photoexcitation of a protein-embedded pigment, followed by complex processes comprising, for example, electron or excitation energy transfer in photosynthetic complexes. The optical properties of chromophores in living systems are strongly dependent on the interaction with the surrounding environment (nearby protein residues, membrane, water), and the complexity of such interplay is, in most cases, at the origin of the functional diversity of the photoactive proteins. The specific interactions with the environment often lead to a significant shift of the chromophore excitation energies, compared with their absorption in solution or gas phase. The investigation of the optical response of chromophores is generally not straightforward, from both experimental and theoretical standpoints; this is due to the difficulty in understanding diverse behaviours and effects, occurring at different scales, with a single technique. In particular, the role played by ab initio calculations in assisting and guiding experiments, as well as in understanding the physics of photoactive proteins, is fundamental. At the same time, owing to the large size of the systems, more approximate strategies which take into account the environmental effects on the absorption spectra are also of paramount importance. Here we review the recent advances in the first-principle description of electronic and optical properties of biological chromophores embedded in a protein environment. We show their applications on paradigmatic systems, such as the light-harvesting complexes, rhodopsin and green fluorescent protein, emphasising the theoretical frameworks which are of common use in solid state physics, and emerging as promising tools for biomolecular systems.

Excitation Dynamics in Hetero-bichromophoric Calixarene Systems

In this work, the dynamics of electronic energy transfer (EET) in bichromophoric donor-acceptor systems, obtained by functionalizing a calix[4]arene scaffold with two dyes, was experimentally and theoretically characterized. The investigated compounds are highly versatile, due to the possibility of linking the dye molecules to the cone or partial cone structure of the calix[4]arene, which directs the two active units to the same or opposite side of the scaffold, respectively. The dynamics and efficiency of the EET process between the donor and acceptor units was investigated and discussed through a combined experimental and theoretical approach, involving ultrafast pump-probe spectroscopy and density functional theory based characterization of the energetic and spectroscopic properties of the system. Our results suggest that the external medium strongly determines the particular conformation adopted by the bichromophores, with a direct effect on the extent of excitonic coupling between the dyes and hence on the dynamics of the EET process itself.

LICHEM: A QM/MM program for simulations with multipolar and polarizable force fields

We introduce an initial implementation of the LICHEM software package. LICHEM can interface with Gaussian, PSI4, NWChem, TINKER, and TINKER-HP to enable QM/MM calculations using multipolar/polarizable force fields. LICHEM extracts forces and energies from unmodified QM and MM software packages to perform geometry optimizations, single-point energy calculations, or Monte Carlo simulations. When the QM and MM regions are connected by covalent bonds, the pseudo-bond approach is employed to smoothly transition between the QM region and the polarizable force field. A series of water clusters and small peptides have been employed to test our initial implementation. The results obtained from these test systems show the capabilities of the new software and highlight the importance of including explicit polarization. © 2016 Wiley Periodicals, Inc.

Quantum Chemical Studies of Light Harvesting

The design of optimal light-harvesting (supra)molecular systems and materials is one of the most challenging frontiers of science. Theoretical methods and computational models play a fundamental role in this difficult task, as they allow the establishment of structural blueprints inspired by natural photosynthetic organisms that can be applied to the design of novel artificial light-harvesting devices. Among theoretical strategies, the application of quantum chemical tools represents an important reality that has already reached an evident degree of maturity, although it still has to show its real potentials. This Review presents an overview of the state of the art of this strategy, showing the actual fields of applicability but also indicating its current limitations, which need to be solved in future developments.

The fragment molecular orbital method combined with density-functional tight-binding and the polarizable continuum model

The electronic gap in proteins is analyzed in detail, and it is shown that FMO-DFTB/PCM is efficient and accurate in describing the molecular structure of proteins in solution.

A curve-crossing model to rationalize and optimize diarylethene dyads

Extra crossing points play a key role in the photochemistry of diarylethene dyads.

Going from photochromic compounds presenting a single switchable function to multi-addressable photochromic multimers remains an extremely difficult task notably because the interactions of several photochromic units through a linker generally result in a substantial loss of photoactivity. Due to their size and the intrinsic complexity of their electronic structure, coupled photochromes also constitute a fundamental challenge for theoretical chemistry. We present here an effective curve-crossing model that, used in connection with easily accessible ab initio data, allows a first understanding of the difficulty to obtain efficient multiphotochromes. Indeed, we demonstrate that extra crossing points, specific to multiphotochromes, have to be passed to ensure reactivity. In addition, the proposed approach allows the definition of an intuitive tilt criterion that can be used to screen a large number of substitution patterns and hence help in the design of new compounds, an aspect that is also developed here. The compatibility of this tilt criterion with previously proposed static Franck–Condon parameters is discussed as well.

Quantum, classical, and hybrid QM/MM calculations in solution: general implementation of the ddCOSMO linear scaling strategy

We present the general theory and implementation of the Conductor-like Screening Model according to the recently developed ddCOSMO paradigm. The various quantities needed to apply ddCOSMO at different levels of theory, including quantum mechanical descriptions, are discussed in detail, with a particular focus on how to compute the integrals needed to evaluate the ddCOSMO solvation energy and its derivatives. The overall computational cost of a ddCOSMO computation is then analyzed and decomposed in the various steps: the different relative weights of such contributions are then discussed for both ddCOSMO and the fastest available alternative discretization to the COSMO equations. Finally, the scaling of the cost of the various steps with respect to the size of the solute is analyzed and discussed, showing how ddCOSMO opens significantly new possibilities when cheap or hybrid molecular mechanics/quantum mechanics methods are used to describe the solute.