1
|
Petrusevich EF, Bousquet MHE, Ośmiałowski B, Jacquemin D, Luis JM, Zaleśny R. Cost-Effective Simulations of Vibrationally-Resolved Absorption Spectra of Fluorophores with Machine-Learning-Based Inhomogeneous Broadening. J Chem Theory Comput 2023; 19:2304-2315. [PMID: 37096370 PMCID: PMC10134414 DOI: 10.1021/acs.jctc.2c01285] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/08/2023]
Abstract
The results of electronic and vibrational structure simulations are an invaluable support for interpreting experimental absorption/emission spectra, which stimulates the development of reliable and cost-effective computational protocols. In this work, we contribute to these efforts and propose an efficient first-principle protocol for simulating vibrationally-resolved absorption spectra, including nonempirical estimations of the inhomogeneous broadening. To this end, we analyze three key aspects: (i) a metric-based selection of density functional approximation (DFA) so to benefit from the computational efficiency of time-dependent density function theory (TD-DFT) while safeguarding the accuracy of the vibrationally-resolved spectra, (ii) an assessment of two vibrational structure schemes (vertical gradient and adiabatic Hessian) to compute the Franck-Condon factors, and (iii) the use of machine learning to speed up nonempirical estimations of the inhomogeneous broadening. In more detail, we predict the absorption band shapes for a set of 20 medium-sized fluorescent dyes, focusing on the bright ππ★ S0 → S1 transition and using experimental results as references. We demonstrate that, for the studied 20-dye set which includes structures with large structural variability, the preselection of DFAs based on an easily accessible metric ensures accurate band shapes with respect to the reference approach and that range-separated functionals show the best performance when combined with the vertical gradient model. As far as band widths are concerned, we propose a new machine-learning-based approach for determining the inhomogeneous broadening induced by the solvent microenvironment. This approach is shown to be very robust offering inhomogeneous broadenings with errors as small as 2 cm-1 with respect to genuine electronic-structure calculations, with a total CPU time reduced by 98%.
Collapse
Affiliation(s)
- Elizaveta F. Petrusevich
- Faculty of Chemistry, Wrocław University of Science and Technology, Wyb. Wyspiańskiego 27, PL-50370 Wrocław, Poland
- Institute of Computational Chemistry and Catalysis and Department of Chemistry, University of Girona, Campus de Montilivi, 17003 Girona, Catalonia, Spain
| | | | - Borys Ośmiałowski
- Faculty of Chemistry, Nicolaus Copernicus University, Gagarina Street 7, PL-87-100 Toruń, Poland
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
- Institut Universitaire de France (IUF), F-75005 Paris, France
| | - Josep M. Luis
- Institute of Computational Chemistry and Catalysis and Department of Chemistry, University of Girona, Campus de Montilivi, 17003 Girona, Catalonia, Spain
| | - Robert Zaleśny
- Faculty of Chemistry, Wrocław University of Science and Technology, Wyb. Wyspiańskiego 27, PL-50370 Wrocław, Poland
| |
Collapse
|
2
|
Chołuj M, Alam MM, Beerepoot MTP, Sitkiewicz SP, Matito E, Ruud K, Zaleśny R. Choosing Bad versus Worse: Predictions of Two-Photon-Absorption Strengths Based on Popular Density Functional Approximations. J Chem Theory Comput 2022; 18:1046-1060. [PMID: 35080389 PMCID: PMC8830054 DOI: 10.1021/acs.jctc.1c01056] [Citation(s) in RCA: 26] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
We present a benchmark
study of density functional approximation
(DFA) performances in predicting the two-photon-absorption strengths
in π-conjugated molecules containing electron-donating/-accepting
moieties. A set of 48 organic molecules is chosen for this purpose,
for which the two-photon-absorption (2PA) parameters are evaluated
using different DFAs, including BLYP, PBE, B3LYP, PBE0, CAM-B3LYP,
LC-BLYP, and optimally tuned LC-BLYP. Minnesota functionals and ωB97X-D
are also used, applying the two-state approximation, for a subset
of molecules. The efficient resolution-of-identity implementation
of the coupled-cluster CC2 model (RI-CC2) is used as a reference for
the assessment of the DFAs. Two-state models within the framework
of both DFAs and RI-CC2 are used to gain a deeper insight into the
performance of different DFAs. Our results give a clear picture of
the performance of the density functionals in describing the two-photon
activity in dipolar π-conjugated systems. The results show that
global hybrids are best suited to reproduce the absolute values of
2PA strengths of donor–acceptor molecules. The range-separated
functionals CAM-B3LYP and optimally tuned LC-BLYP, however, show the
highest linear correlations with the reference RI-CC2 results. Hence,
we recommend the latter DFAs for structure–property studies
across large series of dipolar compounds.
Collapse
Affiliation(s)
- Marta Chołuj
- Faculty of Chemistry, Wrocław University of Science and Technology, Wyb. Wyspiańskiego 27, PL-50370 Wrocław, Poland
| | - Md Mehboob Alam
- Department of Chemistry, Indian Institute of Technology Bhilai, Sejbahar, Raipur, Chhattisgarh 492015, India
| | - Maarten T P Beerepoot
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Sebastian P Sitkiewicz
- Donostia International Physics Center (DIPC), Manuel Lardizabal Ibilbidea 4, 20018 Donostia, Euskadi, Spain.,Kimika Fakultatea, Euskal Herriko Unibertsitatea (UPV/EHU), 20080 Donostia, Euskadi, Spain
| | - Eduard Matito
- Donostia International Physics Center (DIPC), Manuel Lardizabal Ibilbidea 4, 20018 Donostia, Euskadi, Spain.,Ikerbasque Foundation for Science, Plaza Euskadi 5, 48009 Bilbao, Euskadi, Spain
| | - Kenneth Ruud
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Robert Zaleśny
- Faculty of Chemistry, Wrocław University of Science and Technology, Wyb. Wyspiańskiego 27, PL-50370 Wrocław, Poland
| |
Collapse
|
3
|
Murugan NA, Zaleśny R. Multiscale Modeling of Two-Photon Probes for Parkinson's Diagnostics Based on Monoamine Oxidase B Biomarker. J Chem Inf Model 2020; 60:3854-3863. [PMID: 32786232 PMCID: PMC7458361 DOI: 10.1021/acs.jcim.0c00423] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
Monoamine oxidase B (MAO-B) is a potential biomarker for Parkinson's disease (PD), a neurodegenerative disease associated with the loss of motor activities in human subjects. The disease state is associated with dopamine deprival, and so the inhibitors of MAO-B can serve as therapeutic drugs for PD. Since the expression level of MAO-B directly correlates to the disease progress, the distribution and population of this enzyme can be employed to monitor disease development. One of the approaches available for estimating the population is two-photon imaging. The ligands used for two-photon imaging should have high binding affinity and binding specificity toward MAO-B along with significant two-photon absorption cross sections when they are bound to the target. In this article, we study using a multiscale modeling approach, the binding affinity and spectroscopic properties (one- and two-photon absorption) of three (Flu1, Flu2, Flu3) of the currently available probes for monitoring the MAO-B level. We report that the binding affinity of the probes can be explained using the molecular size and binding cavity volume. The experimentally determined one-photon absorption spectrum is well reproduced by the employed QM/MM approaches, and the most accurate spectral shifts, on passing from one probe to another, are obtained at the coupled-cluster (CC2) level of theory. An important conclusion from this study is also the demonstration that intrinsic molecular two-photon absorption strengths (δ2PA) increase in the order δ2PA (Flu1) > δ2PA (Flu2) > δ2PA (Flu3). This is in contrast with experimental data, which predict similar values of two-photon absorption cross sections for Flu1 and Flu3. We demontrate, based on the results of electronic-structure calculations for Flu1 that this discrepancy cannot be explained by an explicit account for neighboring residues (which could lead to charge transfer between a probe and neighboring aromatic amino acids thus boosting δ2PA). In summary, we show that the employed multiscale approach not only can optimize two-photon absorption properties and verify binding affinity, but it can also help in detailed analyses of experimental data.
Collapse
Affiliation(s)
- N Arul Murugan
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, S-106 91 Stockholm, Sweden
| | - Robert Zaleśny
- Department of Physical and Quantum Chemistry, Faculty of Chemistry, Wrocław University of Science and Technology, Wyb. Wyspiańskiego 27, PL-50370 Wrocław, Poland
| |
Collapse
|
4
|
Bednarska J, Zaleśny R, Tian G, Murugan NA, Ågren H, Bartkowiak W. Nonempirical Simulations of Inhomogeneous Broadening of Electronic Transitions in Solution: Predicting Band Shapes in One- and Two-Photon Absorption Spectra of Chalcones. Molecules 2017; 22:E1643. [PMID: 28973973 PMCID: PMC6151831 DOI: 10.3390/molecules22101643] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2017] [Revised: 09/20/2017] [Accepted: 09/21/2017] [Indexed: 11/23/2022] Open
Abstract
We have examined several approaches relying on the Polarizable Embedding (PE) scheme to predict optical band shapes for two chalcone molecules in methanol solution. The PE-TDDFT and PERI-CC2 methods were combined with molecular dynamics simulations, where the solute geometry was kept either as rigid, flexible or partly-flexible (restrained) body. The first approach, termed RBMD-PE-TDDFT, was employed to estimate the inhomogeneous broadening for subsequent convolution with the vibrationally-resolved spectra of the molecule in solution determined quantum-mechanically (QM). As demonstrated, the RBMD-PE-TDDFT/QM-PCM approach delivers accurate band widths, also reproducing their correct asymmetric shapes. Further refinement can be obtained by the estimation of the inhomogeneous broadening using the RBMD-PERI-CC2 method. On the other hand, the remaining two approaches (FBMD-PE-TDDFT and ResBMD-PE-TDDFT), which lack quantum-mechanical treatment of molecular vibrations, lead to underestimated band widths. In this study, we also proposed a simple strategy regarding the rapid selection of the exchange-correlation functional for the simulations of vibrationally-resolved one- and two-photon absorption spectra based on two easy-to-compute metrics.
Collapse
Affiliation(s)
- Joanna Bednarska
- Department of Physical and Quantum Chemistry, Faculty of Chemistry, Wrocław University of Science and Technology, Wyb. Wyspiańskiego 27, PL-50370 Wrocław, Poland.
| | - Robert Zaleśny
- Department of Physical and Quantum Chemistry, Faculty of Chemistry, Wrocław University of Science and Technology, Wyb. Wyspiańskiego 27, PL-50370 Wrocław, Poland.
| | - Guangjun Tian
- Hebei Key Laboratory of Microstructural Material Physics, School of Science, Yanshan University, Qinhuangdao 066004, China.
| | - Natarajan Arul Murugan
- Division of Theoretical Chemistry and Biology, School of Biotechnology, Royal Institute of Technology, SE-10691 Stockholm, Sweden.
| | - Hans Ågren
- Division of Theoretical Chemistry and Biology, School of Biotechnology, Royal Institute of Technology, SE-10691 Stockholm, Sweden.
| | - Wojciech Bartkowiak
- Department of Physical and Quantum Chemistry, Faculty of Chemistry, Wrocław University of Science and Technology, Wyb. Wyspiańskiego 27, PL-50370 Wrocław, Poland.
| |
Collapse
|
5
|
Zhang C, Zou LY, Wang D, Zhang L, Ren AM. A theoretical investigation and design of fluorescent probes with large two-photon absorption cross-sections and dual signaling behavior for mercury ions in biological system. J Photochem Photobiol A Chem 2017. [DOI: 10.1016/j.jphotochem.2016.09.016] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
6
|
Bednarska J, Zaleśny R, Arul Murugan N, Bartkowiak W, Ågren H, Odelius M. Elucidating the Mechanism of Zn(2+) Sensing by a Bipyridine Probe Based on Two-Photon Absorption. J Phys Chem B 2016; 120:9067-75. [PMID: 27494451 DOI: 10.1021/acs.jpcb.6b04949] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
In this work, we examine, by means of computational methods, the mechanism of Zn(2+) sensing by a bipyridine-centered, D-π-A-π-D-type ratiometric molecular probe. According to recently published experimental data [Divya, K. P.; Sreejith, S.; Ashokkumar, P.; Yuzhan, K.; Peng, Q.; Maji, S. K.; Tong, Y.; Yu, H.; Zhao, Y.; Ramamurthy, P.; Ajayaghosh, A. A ratiometric fluorescent molecular probe with enhanced two-photon response upon Zn(2+) binding for in vitro and in vivo bioimaging. Chem. Sci. 2014, 5, 3469-3474], after coordination to zinc ions the probe exhibits a large enhancement of the two-photon absorption cross section. The goal of our investigation was to elucidate the mechanism behind this phenomenon. For this purpose, linear and nonlinear optical properties of the unbound (cation-free) and bound probe were calculated, including the influence of solute-solvent interactions, implicitly using a polarizable continuum model and explicitely employing the QM/MM approach. Because the results of the calculations indicate that many conformers of the probe are energetically accessible at room temperature in solution and hence contribute to the signal, structure-property relationships were also taken into account. Results of our simulations demonstrate that the one-photon absorption bands for both the unbound and bound forms correspond to the bright π → π* transition to the first excited state, which, on the other hand, exhibits negligible two-photon activity. On the basis of the results of the quadratic response calculations, we put forward a notion that it is the second excited state that gives the strong signal in the experimental nonlinear spectrum. To explain the differences in the two-photon absorption activity for the two lowest-lying excited states and nonlinear response enhancement upon binding, we employed the generalized few-state model including the ground, first, and second excited states. The analysis of the optical channel suggests that the large two-photon response is due to the coordination-induced increase of the transition moment from the first to the second excited state.
Collapse
Affiliation(s)
- Joanna Bednarska
- Department of Physical and Quantum Chemistry, Faculty of Chemistry, Wrocław University of Science and Technology , Wyb. Wyspiańskiego 27, PL-50370 Wrocław, Poland
| | - Robert Zaleśny
- Department of Physical and Quantum Chemistry, Faculty of Chemistry, Wrocław University of Science and Technology , Wyb. Wyspiańskiego 27, PL-50370 Wrocław, Poland
| | - N Arul Murugan
- Division of Theoretical Chemistry and Biology, School of Biotechnology, Royal Institute of Technology , SE-10691 Stockholm, Sweden
| | - Wojciech Bartkowiak
- Department of Physical and Quantum Chemistry, Faculty of Chemistry, Wrocław University of Science and Technology , Wyb. Wyspiańskiego 27, PL-50370 Wrocław, Poland
| | - Hans Ågren
- Division of Theoretical Chemistry and Biology, School of Biotechnology, Royal Institute of Technology , SE-10691 Stockholm, Sweden
| | - Michael Odelius
- Division of Chemical Physics, Department of Physics, Stockholm University , SE-106 91 Stockholm, Sweden
| |
Collapse
|
7
|
Zaleśny R, Murugan NA, Tian G, Medved’ M, Ågren H. First-Principles Simulations of One- and Two-Photon Absorption Band Shapes of the Bis(BF2) Core Complex. J Phys Chem B 2016; 120:2323-32. [DOI: 10.1021/acs.jpcb.5b09726] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- Robert Zaleśny
- Department of Physical and
Quantum Chemistry, Faculty of
Chemistry, Wrocław University of Technology, Wyb. Wyspiańskiego 27, PL-50370 Wrocław, Poland
| | - N. Arul Murugan
- Division of Theoretical Chemistry and Biology, School of Biotechnology, Royal Institute of Technology, SE-10691 Stockholm, Sweden
| | - Guangjun Tian
- Division of Theoretical Chemistry and Biology, School of Biotechnology, Royal Institute of Technology, SE-10691 Stockholm, Sweden
| | - Miroslav Medved’
- Department of Chemistry, Faculty of Natural Sciences, Matej Bel University, Tajovského 40, 974 01 Banská Bystrica, Slovak Republic
| | - Hans Ågren
- Division of Theoretical Chemistry and Biology, School of Biotechnology, Royal Institute of Technology, SE-10691 Stockholm, Sweden
| |
Collapse
|
8
|
Alam MM, Kundi V, Thankachan PP. Solvent effects on static polarizability, static first hyperpolarizability and one- and two-photon absorption properties of functionalized triply twisted Möbius annulenes: a DFT study. Phys Chem Chem Phys 2016; 18:21833-42. [DOI: 10.1039/c6cp02732f] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Solvent effects on the polarizability (α), static first hyperpolarizability (β) and one- and two-photon absorption (OPA and TPA) properties of triply twisted Möbius annulenes.
Collapse
Affiliation(s)
- Md Mehboob Alam
- Laboratoire de Chimie Quantique
- Institute de Chimie
- CNRS/Université de Strasbourg
- 67000 Strasbourg
- France
| | - Varun Kundi
- Theoretical Chemistry Lab
- Department of Chemistry
- Indian Institute of Technology Roorkee
- Roorkee-247667
- India
| | | |
Collapse
|
9
|
Zaleśny R, Tian G, Hättig C, Bartkowiak W, Ågren H. Toward assessment of density functionals for vibronic coupling in two-photon absorption: A case study of 4-nitroaniline. J Comput Chem 2015; 36:1124-31. [DOI: 10.1002/jcc.23903] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2014] [Revised: 02/24/2015] [Accepted: 02/28/2015] [Indexed: 12/13/2022]
Affiliation(s)
- Robert Zaleśny
- Department of Physical and Quantum Chemistry, Faculty of Chemistry; Wrocław University of Technology; Wyb. Wyspiańskiego 27 PL-50370 Wrocław Poland
- Division of Theoretical Chemistry and Biology; School of Biotechnology, Royal Institute of Technology; SE-10691 Stockholm Sweden
| | - Guangjun Tian
- Division of Theoretical Chemistry and Biology; School of Biotechnology, Royal Institute of Technology; SE-10691 Stockholm Sweden
| | - Christof Hättig
- Lehrstuhl fuer Theoretische Chemie; Ruhr-Universität Bochum, Universitätsstrasse 150; D-44780 Bochum Germany
| | - Wojciech Bartkowiak
- Department of Physical and Quantum Chemistry, Faculty of Chemistry; Wrocław University of Technology; Wyb. Wyspiańskiego 27 PL-50370 Wrocław Poland
| | - Hans Ågren
- Division of Theoretical Chemistry and Biology; School of Biotechnology, Royal Institute of Technology; SE-10691 Stockholm Sweden
| |
Collapse
|
10
|
Matczyszyn K, Olesiak-Banska J, Nakatani K, Yu P, Murugan NA, Zaleśny R, Roztoczyńska A, Bednarska J, Bartkowiak W, Kongsted J, Ågren H, Samoć M. One- and two-photon absorption of a spiropyran-merocyanine system: experimental and theoretical studies. J Phys Chem B 2015; 119:1515-22. [PMID: 25531561 DOI: 10.1021/jp5071715] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We report on the nonlinear optical properties measurements and quantum-chemical calculations of a well-known photochromic system consisting of spiropyran and the merocyanine photoproduct. The study of nonlinear absorption and refraction properties of the molecules dissolved in chloroform were performed with the Z-scan technique, using femtosecond pulses in a wide range of wavelengths. Maxima in the two-photon absorption spectrum at 700 and 1050 nm were found for the merocyanine form, and the corresponding two-photon absorption cross section is 80 GM and 20 GM, respectively. The latter feature does not vanish completely in the nonlinear spectrum of the spiropyran form, possibly because of the existence of some photoconversion caused by the laser beam during the measurements. A nonlinear absorption peak at 900 nm is found in the spiropyran form with an effective cross section of about 20 GM; it is likely due to three-photon absorption or to absorption by some intermediate species. The experimental data are supported by calculations performed with the use of a hybrid quantum mechanics-molecular mechanics approach.
Collapse
Affiliation(s)
- Katarzyna Matczyszyn
- Institute of Physical and Theoretical Chemistry, Wrocław University of Technology , Wyb. Wyspiańskiego 27, PL-50-370 Wrocław, Poland
| | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
11
|
Zaleśny R, Murugan NA, Gel’mukhanov F, Rinkevicius Z, Ośmiałowski B, Bartkowiak W, Ågren H. Toward Fully Nonempirical Simulations of Optical Band Shapes of Molecules in Solution: A Case Study of Heterocyclic Ketoimine Difluoroborates. J Phys Chem A 2014; 119:5145-52. [DOI: 10.1021/jp5094417] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
| | | | | | | | - Borys Ośmiałowski
- Faculty
of Chemical Technology and Engineering, UTP University of Science and Technology, Seminaryjna 3, PL-85326 Bydgoszcz, Poland
| | - Wojciech Bartkowiak
- Institute
of Physical and Theoretical Chemistry, Wroclaw University of Technology, Wyb. Wyspiańskiego 27, PL-50370 Wrocław, Poland
| | | |
Collapse
|