Simulating the solvation structure of low- and high-spin [Fe(bpy)3]2+: long-range dispersion and many-body effects.

When characterizing transition metal complexes and their functionalities, the importance of including the solvent as an active participant is becoming more and more apparent. Whereas many studies have evaluated long-range dispersion effects inside organic molecules and organometallics, less is known about their role in solvation. Here, we have analysed the components within solute-solvent and solvent-solvent interactions of one of the most studied iron-based photoswitch model systems, in two spin states. We find that long-range dispersion effects modulate the coordination significantly, and that this is accurately captured by density functional theory models including dispersion corrections. We furthermore correlate gas-phase relaxed complex-water clusters to thermally averaged molecular densities. This shows how the gas-phase interactions translate to solution structure, quantified through 3D molecular densities, angular distributions, and radial distribution functions. We show that finite-size simulation cells can cause the radial distribution functions to have artificially enlarged amplitudes. Finally, we quantify the effects of many-body interactions within the solvent shells, and find that almost a fifth of the total interaction energy of the solute-shell system in the high-spin state comes from many-body contributions, which cannot be captured by by pair-wise additive force field methods.

Taming the Complexity of Donor-Acceptor Stenhouse Adducts: Infrared Motion Pictures of the Complete Switching Pathway.

Switches that can be actively steered by external stimuli along multiple pathways at the molecular level are the basis for next-generation responsive material systems. The operation of commonly employed molecular photoswitches revolves around one key structural coordinate. Photoswitches with functionalities that depend on and can be addressed along multiple coordinates would offer novel means to tailor and control their behavior and performance. The recently developed donor-acceptor Stenhouse adducts (DASAs) are versatile switches suitable for such applications. Their photochemistry is well understood, but is only responsible for part of their overall photoswitching mechanism. The remaining thermal switching pathways are to date unknown. Here, rapid-scan infrared absorption spectroscopy is used to obtain transient fingerprints of reactions occurring on the ground state potential energy surface after reaching structures generated through light absorption. The spectroscopic data are interpreted in terms of structural transformations using kinetic modeling and quantum chemical calculations. Through this combined experimental-theoretical approach, we are able to unravel the complexity of the multidimensional ground-state potential energy surface explored by the photoswitch and use this knowledge to predict, and subsequently confirm, how DASA switches can be guided along this potential energy surface. These results break new ground for developing user-geared DASA switches but also shed light on the development of novel photoswitches in general.

Photoprogramming Allostery in Human Serum Albumin.

Developing strategies to interfere with allosteric interactions in proteins not only promises to deepen our understanding of vital cellular processes but also allows their regulation using external triggers. Light is particularly attractive as a trigger being spatiotemporally selective and compatible with the physiological environment. Here, we engineered a hybrid protein in which irradiation with light opens a new allosteric communication route that is not inherent to the natural system. We select human serum albumin, a promiscuous protein responsible for transporting a variety of ligands in plasma, and show that by covalently incorporating a synthetic photoswitch to subdomain IA we achieve optical control of the ligand binding in subdomain IB. Molecular dynamics simulations confirm the allosteric nature of the interactions between IA and IB in the engineered protein. Specifically, upon illumination, photoconversion of the switch is found to correlate with a less-coordinated motion of the two subdomains and an increased flexibility of the binding pocket in subdomain IB, whose fluctuations are cooperatively enhanced by the presence of ligands, ultimately facilitating their release. Our combined experimental and computational work demonstrates how harnessing artificial molecular switches enables photoprogramming the allosteric regulation of binding activities in such a prominent protein.

Crystal Field in Rare-Earth Complexes: From Electrostatics to Bonding.

The flexibility of first-principles (ab initio) calculations with the SO-CASSCF (complete active space self-consistent field theory with a treatment of the spin-orbit (SO) coupling by state interaction) method is used to quantify the electrostatic and covalent contributions to crystal field parameters. Two types of systems are chosen for illustration: 1) The ionic and experimentally well-characterized PrCl3 crystal; this study permits a revisitation of the partition of contributions proposed in the early days of crystal field theory; and 2) a series of sandwich molecules [Ln(ηn -Cn Hn )2 ]q , with Ln=Dy, Ho, Er, and Tm and n=5, 6, and 8, in which the interaction between LnIII and the aromatic ligands is more difficult to describe within an electrostatic approach. It is shown that a model with three layers of charges reproduces the electrostatic field generated by the ligands and that the covalency plays a qualitative role. The one-electron character of crystal field theory is discussed and shown to be valuable, although it is not completely quantitative. This permits a reduction of the many-electron problem to a discussion of the energy of the seven 4f orbitals.

Redshifted and thermally bistable one-way quantitative hemithioindigo-derived photoswitches enabled by isomer-specific excited state intramolecular proton transfer.

We report bistable indole-containing hemithioindigos (HTIs) with one-way quantitative photoswitching properties. Supported by state-averaged CASPT2/CASSCF calculations, we propose a mechanism for the observed one-way photoswitching that involves an isomer-specific excited state intramolecular proton transfer (ESIPT). Additionally, we developed a thermally bistable oligomer-inspired bipyrrole-containing HTI, which displays large band separation and bidirectional near-quantitative photoisomerization in the near-infrared, bio-optical window.

Multiple-Resonance Local Wave Functions for Accurate Excited States in Quantum Monte Carlo.

We introduce a novel class of local multideterminant Jastrow-Slater wave functions for the efficient and accurate treatment of excited states in quantum Monte Carlo. The wave function is expanded as a linear combination of excitations built from multiple sets of localized orbitals that correspond to the bonding patterns of the different Lewis resonance structures of the molecule. We capitalize on the concept of orbital domains of local coupled-cluster methods, which is here applied to the active space to select the orbitals to correlate and construct the important transitions. The excitations are further grouped into classes, which are ordered in importance and can be systematically included in the Jastrow-Slater wave function to ensure a balanced description of all states of interest. We assess the performance of the proposed wave function in the calculation of vertical excitation energies and excited-state geometry optimization of retinal models whose π → π* state has a strong intramolecular charge-transfer character. We find that our multiresonance wave functions recover the reference values of the total energies of the ground and excited states with only a small number of excitations and that the same expansion can be flexibly used at very different geometries. Furthermore, significant computational saving can also be gained in the orbital optimization step by selectively mixing occupied and virtual orbitals based on spatial considerations without loss of accuracy on the excitation energy. Our multiresonance wave functions are therefore compact, accurate, and very promising for the calculation of multiple excited states of different character in large molecules.

Introducing QMC/MMpol: Quantum Monte Carlo in Polarizable Force Fields for Excited States.

We present for the first time a quantum mechanics/molecular mechanics scheme which combines quantum Monte Carlo with the reaction field of classical polarizable dipoles (QMC/MMpol). In our approach, the optimal dipoles are self-consistently generated at the variational Monte Carlo level and then used to include environmental effects in diffusion Monte Carlo. We investigate the performance of this hybrid model in describing the vertical excitation energies of prototypical small molecules solvated in water, namely, methylenecyclopropene and s-trans acrolein. Two polarization regimes are explored where either the dipoles are optimized with respect to the ground-state solute density (polGS) or different sets of dipoles are separately brought to equilibrium with the states involved in the electronic transition (polSS). By comparing with reference supermolecular calculations where both solute and solvent are treated quantum mechanically, we find that the inclusion of the response of the environment to the excitation of the solute leads to superior results than the use of a frozen environment (point charges or polGS), in particular, when the solute-solvent coupling is dominated by electrostatic effects which are well recovered in the polSS condition. QMC/MMpol represents therefore a robust scheme to treat important environmental effects beyond static point charges, combining the accuracy of QMC with the simplicity of a classical approach.