1
|
Imamura K, Yokogawa D, Sato H. Recent developments and applications of reference interaction site model self-consistent field with constrained spatial electron density (RISM-SCF-cSED): A hybrid model of quantum chemistry and integral equation theory of molecular liquids. J Chem Phys 2024; 160:050901. [PMID: 38341702 DOI: 10.1063/5.0190116] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2023] [Accepted: 01/04/2024] [Indexed: 02/13/2024] Open
Abstract
The significance of solvent effects in electronic structure calculations has long been noted, and various methods have been developed to consider this effect. The reference interaction site model self-consistent field with constrained spatial electron density (RISM-SCF-cSED) is a hybrid model that combines the integral equation theory of molecular liquids with quantum chemistry. This method can consider the statistically convergent solvent distribution at a significantly lower cost than molecular dynamics simulations. Because the RISM theory explicitly considers the solvent structure, it performs well for systems where hydrogen bonds are formed between the solute and solvent molecules, which is a challenge for continuum solvent models. Taking advantage of being founded on the variational principle, theoretical developments have been made in calculating various properties and incorporating electron correlation effects. In this review, we organize the theoretical aspects of RISM-SCF-cSED and its distinctions from other hybrid methods involving integral equation theories. Furthermore, we carefully present its progress in terms of theoretical developments and recent applications.
Collapse
Affiliation(s)
- Kosuke Imamura
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan
| | - Daisuke Yokogawa
- Graduate School of Arts and Science, The University of Tokyo, Komaba, Meguro-ku, Tokyo 153-8902, Japan
| | - Hirofumi Sato
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan
- Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8103, Japan
| |
Collapse
|
2
|
Reimann M, Kaupp M. Evaluation of an Efficient 3D-RISM-SCF Implementation as a Tool for Computational Spectroscopy in Solution. J Phys Chem A 2020; 124:7439-7452. [PMID: 32838530 DOI: 10.1021/acs.jpca.0c06322] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
The 3D-RISM-SCF solvent-model implementation of Gusarov et al. [ J. Phys. Chem. A 2006, 110, 6083-6090] in the Amsterdam density functional program has been improved and extended. In particular, an accurate yet efficient representation of the solute electrostatic potential is provided. The Coulomb-potential fitting of many DFT codes can be used advantageously in this context. The extra effort compared to a point-charge representation is small for a given SCF cycle and compensated by faster SCF convergence. This allows applications to large solutes, as demonstrated by evaluation of the solvatochromism of Reichardt's dye. In general, TDDFT applications to excitation energies in solution stand out and are highlighted. Applications to the 17O NMR chemical shifts of N-methylformamide in different solvents also demonstrate the distinct advantages of 3D-RISM over continuum solvents. Limitations are observed in this case for water solvent, where the solvent shielding is overestimated. This shortcoming applies also to the 17O gas-to-liquid shift of water, where we used localized molecular orbital analyses for a deeper understanding. For such cases of extremely strong solute-solvent interactions, couplings between solute and solvent orbitals induced by the magnetic perturbation are relevant. These clearly require a quantum-mechanical treatment of the most closely bound solvent molecules. Except for such extreme cases, 3D-RISM-SCF is very well suited to treat solvent effects on NMR parameters. More serious limitations pertain to the treatment of vibrational spectra, where the absence of the coupling between solute and solvent vibrational modes limits the accuracy of applications of 3D-RISM-SCF. The reported extended, efficient, and numerically accurate 3D-RISM-SCF implementation should provide a useful tool to study chemical and spectroscopic properties of molecules of appreciable size in a realistic solvent environment.
Collapse
Affiliation(s)
- Marc Reimann
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, D-10623, Berlin, Germany
| | - Martin Kaupp
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, D-10623, Berlin, Germany
| |
Collapse
|
3
|
Peng Y, Ye Y, Xiu X, Sun S. Mechanism of Excited-State Intramolecular Proton Transfer for 1,2-Dihydroxyanthraquinone: Effect of Water on the ESIPT. J Phys Chem A 2017; 121:5625-5634. [DOI: 10.1021/acs.jpca.7b03877] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Yajing Peng
- Department
of Physics, Bohai University, JinZhou 121013, China
- State
Key Laboratory of Molecular Reaction Dynamics Dalian Institute of
Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
| | - Yuqing Ye
- Department
of Physics, Bohai University, JinZhou 121013, China
| | - Xianming Xiu
- Department
of Physics, Bohai University, JinZhou 121013, China
| | - Shuang Sun
- Department
of Physics, Bohai University, JinZhou 121013, China
| |
Collapse
|
4
|
Yokogawa D. Time-dependent density functional theory (TD-DFT) coupled with reference interaction site model self-consistent field explicitly including spatial electron density distribution (RISM-SCF-SEDD). J Chem Phys 2016; 145:094101. [DOI: 10.1063/1.4962062] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- D. Yokogawa
- Department of Chemistry, Graduate School of Science, Nagoya University, Chikusa, Nagoya 464-8602, Japan
- Institute of Transformative Bio-Molecules (WPI-ITbM), Nagoya University, Chikusa, Nagoya 464-8602, Japan
| |
Collapse
|
5
|
Aono S, Mori T, Sakaki S. 3D-RISM-MP2 Approach to Hydration Structure of Pt(II) and Pd(II) Complexes: Unusual H-Ahead Mode vs Usual O-Ahead One. J Chem Theory Comput 2016; 12:1189-206. [PMID: 26863511 DOI: 10.1021/acs.jctc.5b01137] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Solvation of transition metal complexes with water has been one of the fundamental topics in physical and coordination chemistry. In particular, Pt(II) complexes have recently attracted considerable interest for their relation to anticancer activity in cisplatin and its analogues, yet the interaction of the water molecule and the metal center has been obscured. The challenge from a theoretical perspective remains that both the microscopic solvation effect and the dynamical electron correlation (DEC) effect have to be treated simultaneously in a reasonable manner. In this work we derive the analytical gradient for the three-dimensional reference interaction site model Møller-Plesset second order (3D-RISM-MP2) free energy. On the basis of the three-regions 3D-RISM self-consistent field (SCF) method recently proposed by us, we apply a new layer of the Z-vector method to the CP-RISM equation as well as point-charge approximation to the derivatives with respect to the density matrix elements in the RISM-CPHF equation to remarkably reduce the computational cost. This method is applied to study the interaction of H2O with the d(8) square planar transition metal complexes in aqueous solution, trans-[Pt(II)Cl2(NH3)(glycine)] (1a), [Pt(II)(NH3)4](2+) (1b), [Pt(II)(CN)4](2-) (1c), and their Pd(II) analogues 2a, 2b, and 2c, respectively, to elucidate whether the usual H2O interaction through O atom (O-ahead mode) or unusual one through H atom (H-ahead mode) is stable in these complexes. We find that the interaction energy of the coordinating water and the transition metal complex changes little when switching from gas to aqueous phase, but the solvation free energy differs remarkably between the two interaction modes, thereby affecting the relative stability of the H-ahead and O-ahead modes. Particularly, in contrast to the expectation that the O-ahead mode is preferred due to the presence of positive charges in 1b, the H-ahead mode is also found to be more stable. The O-ahead mode is found to be more stable than the H-ahead one only in 2b. The energy decomposition analysis (EDA) at the 3D-RISM-MP2 level revealed that the O-ahead mode is stabilized by the electrostatic (ES) interaction, whereas the H-ahead one is mainly stabilized by the DEC effect. The ES interaction is also responsible for the difference between the Pd(II) and Pt(II) complexes; because the electrostatic potential is more negative along the z-axis in the Pt(II) complex than in the Pd(II) one, the O-ahead mode prefers the Pd(II) complexes, whereas the H-ahead becomes predominant in the Pt(II) complexes.
Collapse
Affiliation(s)
- Shinji Aono
- Fukui Institute for Fundamental Chemistry, Kyoto University , Nishihiraki-cho, Takano, Sakyo-ku, Kyoto 606-8103, Japan
| | - Toshifumi Mori
- Institute for Molecular Science , Okazaki, Aichi 444-8585, Japan.,School of Physical Sciences, The Graduate University for Advanced Studies , Okazaki, Aichi 444-8585, Japan
| | - Shigeyoshi Sakaki
- Fukui Institute for Fundamental Chemistry, Kyoto University , Nishihiraki-cho, Takano, Sakyo-ku, Kyoto 606-8103, Japan
| |
Collapse
|
6
|
Azarias C, Budzák Š, Laurent AD, Ulrich G, Jacquemin D. Tuning ESIPT fluorophores into dual emitters. Chem Sci 2016; 7:3763-3774. [PMID: 29997864 PMCID: PMC6008603 DOI: 10.1039/c5sc04826e] [Citation(s) in RCA: 118] [Impact Index Per Article: 14.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2015] [Accepted: 02/23/2016] [Indexed: 12/24/2022] Open
Abstract
Using first-principle approaches, we show how ESIPT can be controlled by fine-tuning of substituents, hence leading to new dual emitters.
Dyes undergoing excited-state intramolecular proton transfer (ESIPT) are known to present large Stokes shifts as a result of the important geometrical reorganisation following photon absorption. When the ESIPT process is not quantitative, one can obtain dual emitters characterised by two distinct fluorescence bands, observed due to emissions from both the canonical and ESIPT isomers. However, dual emission generally requires to maintain a very specific balance, as the relative excited-state free energies of the two tautomers have to be within a narrow window to observe the phenomenon. Consequently, simple chemical intuition is insufficient to optimise dual emission. In the present contribution, we investigate, with the help of quantum-mechanical tools and more precisely, time-dependent density functional theory (TD-DFT) and algebraic diagrammatic construction (ADC), a wide panel of possible ESIPT/dual emitters with various substituents. The selected protocol is first shown to be very robust on a series of structures with known experimental behaviour, and next is applied to novel derivatives with various substituents located at different positions. This work encompasses the largest chemical library of potential ESIPT compounds studied to date. We pinpoint the most promising combinations for building dual emitters, highlight unexpected combination effects and rationalise the impact of the different auxochromes.
Collapse
Affiliation(s)
- Cloé Azarias
- CEISAM , UMR CNRS 6230 , Université de Nantes , BP 92208, 2, Rue de la Houssinière , 44322 Nantes, Cedex 3 , France . ; Tel: +33-2-51-12-55-64
| | - Šimon Budzák
- CEISAM , UMR CNRS 6230 , Université de Nantes , BP 92208, 2, Rue de la Houssinière , 44322 Nantes, Cedex 3 , France . ; Tel: +33-2-51-12-55-64
| | - Adèle D Laurent
- CEISAM , UMR CNRS 6230 , Université de Nantes , BP 92208, 2, Rue de la Houssinière , 44322 Nantes, Cedex 3 , France . ; Tel: +33-2-51-12-55-64
| | - Gilles Ulrich
- Institut de Chimie et Procédés pour l'Energie , l'Environnement et la Santé (ICPEES) UMR7515 CNRS Université de Strasbourg , 25 rue Becquerel , 67087 Strasbourg Cedex 02 , France
| | - Denis Jacquemin
- CEISAM , UMR CNRS 6230 , Université de Nantes , BP 92208, 2, Rue de la Houssinière , 44322 Nantes, Cedex 3 , France . ; Tel: +33-2-51-12-55-64.,Institut Universitaire de France , 1, rue Descartes , F-75005 Paris Cedex 05 , France
| |
Collapse
|
7
|
Laurent AD, Adamo C, Jacquemin D. Dye chemistry with time-dependent density functional theory. Phys Chem Chem Phys 2015; 16:14334-56. [PMID: 24548975 DOI: 10.1039/c3cp55336a] [Citation(s) in RCA: 233] [Impact Index Per Article: 25.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
In this perspective, we present an overview of the determination of excited-state properties of "real-life" dyes, and notably of their optical absorption and emission spectra, performed during the last decade with time-dependent density functional theory (TD-DFT). We discuss the results obtained with both vertical and adiabatic (vibronic) approximations, choosing relevant examples for several series of dyes. These examples include reproducing absorption wavelengths of numerous families of coloured molecules, understanding the specific band shape of amino-anthraquinones, optimising the properties of dyes used in solar cells, mimicking the fluorescence wavelengths of fluorescent brighteners and BODIPY dyes, studying optically active biomolecules and photo-induced proton transfer, as well as improving the properties of photochromes.
Collapse
Affiliation(s)
- Adèle D Laurent
- Chimie Et Interdisciplinarité, Synthèse, Analyse, Modélisation (CEISAM), UMR CNRS no. 6230, BP 92208, Université de Nantes, 2, Rue de la Houssinière, 44322 Nantes, Cedex 3, France.
| | | | | |
Collapse
|
8
|
Houari Y, Chibani S, Jacquemin D, Laurent AD. TD-DFT Assessment of the Excited State Intramolecular Proton Transfer in Hydroxyphenylbenzimidazole (HBI) Dyes. J Phys Chem B 2014; 119:2180-92. [DOI: 10.1021/jp505036d] [Citation(s) in RCA: 47] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Ymène Houari
- Laboratoire
CEISAM - UMR CNR 6230, Université de Nantes, 2 Rue de la
Houssinière, BP 92208, 44322 Nantes Cedex 3, France
| | - Siwar Chibani
- Laboratoire
CEISAM - UMR CNR 6230, Université de Nantes, 2 Rue de la
Houssinière, BP 92208, 44322 Nantes Cedex 3, France
| | - Denis Jacquemin
- Laboratoire
CEISAM - UMR CNR 6230, Université de Nantes, 2 Rue de la
Houssinière, BP 92208, 44322 Nantes Cedex 3, France
- Institut Universitaire
de France, 103, bd Saint-Michel, F-75005 Paris Cedex 05, France
| | - Adèle D. Laurent
- Laboratoire
CEISAM - UMR CNR 6230, Université de Nantes, 2 Rue de la
Houssinière, BP 92208, 44322 Nantes Cedex 3, France
| |
Collapse
|
9
|
State-specific solvation effect on the intramolecular charge transfer reaction in solution: A linear-response free energy TDDFT method. Chem Phys Lett 2014. [DOI: 10.1016/j.cplett.2014.05.104] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|