1
|
Hrivnák T, Medveď M, Bartkowiak W, Zaleśny R. Hyperpolarizabilities of Push-Pull Chromophores in Solution: Interplay between Electronic and Vibrational Contributions. Molecules 2022; 27:8738. [PMID: 36557870 PMCID: PMC9783928 DOI: 10.3390/molecules27248738] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2022] [Revised: 11/25/2022] [Accepted: 11/29/2022] [Indexed: 12/14/2022] Open
Abstract
Contemporary design of new organic non-linear optical (NLO) materials relies to a large extent on the understanding of molecular and electronic structure-property relationships revealed during the years by available computational approaches. The progress in theory-hand-in-hand with experiment-has enabled us to identify and analyze various physical aspects affecting the NLO responses, such as the environmental effects, molecular vibrations, frequency dispersion, and system dynamics. Although it is nowadays possible to reliably address these effects separately, the studies analyzing their mutual interplay are still very limited. Here, we employ density functional theory (DFT) methods in combination with an implicit solvent model to examine the solvent effects on the electronic and harmonic as well as anharmonic vibrational contributions to the static first hyperpolarizability of a series of push-pull α,ω-diphenylpolyene oligomers, which were experimentally shown to exhibit notable second-order NLO responses. We demonstrate that the magnitudes of both vibrational and electronic contributions being comparable in the gas phase significantly increase in solvents, and the enhancement can be, in some cases, as large as three- or even four-fold. The electrical and mechanical anharmonic contributions are not negligible but cancel each other out to a large extent. The computed dynamic solute NLO properties of the studied systems are shown to be in a fair agreement with those derived from experimentally measured electric-field-induced second-harmonic generation (EFISHG) signals. Our results substantiate the necessity to consider concomitantly both solvation and vibrational effects in modeling static NLO properties of solvated systems.
Collapse
Affiliation(s)
- Tomáš Hrivnák
- Department of Molecular Simulations of Polymers, Polymer Institute, Slovak Academy of Sciences, SK-845 41 Bratislava, Slovakia
| | - Miroslav Medveď
- Department of Chemistry, Faculty of Natural Sciences, Matej Bel University, SK-974 00 Banská Bystrica, Slovakia
- Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacký University Olomouc, 783 71 Olomouc, Czech Republic
| | - 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
| | - Robert Zaleśny
- Faculty of Chemistry, Wrocław University of Science and Technology, Wyb. Wyspiańskiego 27, PL-50370 Wrocław, Poland
| |
Collapse
|
2
|
Yokogawa D, Suda K. Analytical second derivatives of the free energy in solution by the reference interaction site model self-consistent field explicitly including constrained spatial electron density distribution. J Chem Phys 2021; 155:204102. [PMID: 34852465 DOI: 10.1063/5.0067248] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The application of analytical derivative methods to solution systems is important because several chemical reactions occur in solution. The reference interaction site model (RISM) is one of the solvation theories used to study solution systems and has shown good performance, especially in the polar solvent systems. Although the analytical first derivative based on the RISM coupled with quantum methods (RISM-SCF) has already been derived, the analytical second derivative has not been proposed yet. Therefore, in this study, the analytical second derivative was derived using RISM-SCF explicitly including constrained spatial electron density distribution (RISM-SCF-cSED). The performance of this method was validated with the Hessian calculations of formaldehyde and para-nitroaniline in solution, and the results demonstrated that the method accurately calculated frequency values at a small computational cost.
Collapse
Affiliation(s)
- Daisuke Yokogawa
- Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro-ku, Tokyo 153-8902, Japan
| | - Kayo Suda
- Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro-ku, Tokyo 153-8902, Japan
| |
Collapse
|
3
|
Marrazzini G, Giovannini T, Egidi F, Cappelli C. Calculation of Linear and Non-linear Electric Response Properties of Systems in Aqueous Solution: A Polarizable Quantum/Classical Approach with Quantum Repulsion Effects. J Chem Theory Comput 2020; 16:6993-7004. [PMID: 33058671 PMCID: PMC8015238 DOI: 10.1021/acs.jctc.0c00674] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2020] [Indexed: 01/10/2023]
Abstract
We present a computational study of polarizabilities and hyperpolarizabilities of organic molecules in aqueous solutions, focusing on solute-water interactions and the way they affect a molecule's linear and non-linear electric response properties. We employ a polarizable quantum mechanics/molecular mechanics (QM/MM) computational model that treats the solute at the QM level while the solvent is treated classically using a force field that includes polarizable charges and dipoles, which dynamically respond to the solute's quantum-mechanical electron density. Quantum confinement effects are also treated by means of a recently implemented method that endows solvent molecules with a parametric electron density, which exerts Pauli repulsion forces upon the solute. By applying the method to a set of aromatic molecules in solution we show that, for both polarizabilities and first hyperpolarizabilities, observed solution values are the result of a delicate balance between electrostatics, hydrogen-bonding, and non-electrostatic solute solvent interactions.
Collapse
Affiliation(s)
- Gioia Marrazzini
- Scuola
Normale Superiore, Piazza dei Cavalieri 7, Pisa 56126, Italy
| | - Tommaso Giovannini
- Department
of Chemistry, Norwegian University of Science
and Technology, Trondheim 7491, Norway
| | - Franco Egidi
- Scuola
Normale Superiore, Piazza dei Cavalieri 7, Pisa 56126, Italy
| | - Chiara Cappelli
- Scuola
Normale Superiore, Piazza dei Cavalieri 7, Pisa 56126, Italy
| |
Collapse
|
4
|
Hrivnák T, Reis H, Neogrády P, Zaleśny R, Medved’ M. Accurate Nonlinear Optical Properties of Solvated para-Nitroaniline Predicted by an Electrostatic Discrete Local Field Approach. J Phys Chem B 2020; 124:10195-10209. [DOI: 10.1021/acs.jpcb.0c06046] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Tomáš Hrivnák
- Department of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University, Mlynská dolina, Ilkovičova 6, SK-842 15 Bratislava, Slovakia
| | - Heribert Reis
- Institute of Chemical Biology, National Hellenic Research Foundation, Vasileos Constantinou 48, GR-11635 Athens, Greece
| | - Pavel Neogrády
- Department of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University, Mlynská dolina, Ilkovičova 6, SK-842 15 Bratislava, Slovakia
| | - Robert Zaleśny
- Department of Physical and Quantum Chemistry, Faculty of Chemistry, Wroclaw University of Science and Technology, Wyb. Wyspiańskiego 27, PL-50370 Wrocław, Poland
| | - Miroslav Medved’
- Department of Chemistry, Faculty of Natural Sciences, Matej Bel University, Tajovského 40, SK-97400 Banská Bystrica, Slovakia
- Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacký University Olomouc, Šlechtitelů 27, 783 71 Olomouc, Czech Republic
| |
Collapse
|
5
|
Giovannini T, Riso RR, Ambrosetti M, Puglisi A, Cappelli C. Electronic transitions for a fully polarizable QM/MM approach based on fluctuating charges and fluctuating dipoles: Linear and corrected linear response regimes. J Chem Phys 2019; 151:174104. [PMID: 31703497 DOI: 10.1063/1.5121396] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Abstract
The fully polarizable Quantum Mechanics/Molecular Mechanics (QM/MM) approach based on fluctuating charges and fluctuating dipoles, named QM/FQFμ [T. Giovannini et al., J. Chem. Theory Comput. 15, 2233 (2019)], is extended to the calculation of vertical excitation energies of solvated molecular systems. Excitation energies are defined within two different solvation regimes, i.e., linear response (LR), where the response of the MM portion is adjusted to the QM transition density, and corrected-Linear Response (cLR) in which the MM response is adjusted to the relaxed QM density, thus being able to account for charge equilibration in the excited state. The model, which is specified in terms of three physical parameters (electronegativity, chemical hardness, and polarizability) is applied to vacuo-to-water solvatochromic shifts of aqueous solutions of para-nitroaniline, pyridine, and pyrimidine. The results show a good agreement with their experimental counterparts, thus highlighting the potentialities of this approach.
Collapse
Affiliation(s)
- Tommaso Giovannini
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| | | | | | | | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
6
|
Sharipov AS, Loukhovitski BI. Small atomic clusters: quantum chemical research of isomeric composition and physical properties. Struct Chem 2019. [DOI: 10.1007/s11224-019-01417-7] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
7
|
Giovannini T, Grazioli L, Ambrosetti M, Cappelli C. Calculation of IR Spectra with a Fully Polarizable QM/MM Approach Based on Fluctuating Charges and Fluctuating Dipoles. J Chem Theory Comput 2019; 15:5495-5507. [DOI: 10.1021/acs.jctc.9b00574] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Tommaso Giovannini
- Department of Chemistry, Norwegian University of Science and Technology, Trondheim 7491, Norway
| | - Laura Grazioli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, Pisa 56126, Italy
| | | | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, Pisa 56126, Italy
| |
Collapse
|
8
|
Di Remigio R, Giovannini T, Ambrosetti M, Cappelli C, Frediani L. Fully Polarizable QM/Fluctuating Charge Approach to Two-Photon Absorption of Aqueous Solutions. J Chem Theory Comput 2019; 15:4056-4068. [DOI: 10.1021/acs.jctc.9b00305] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Roberto Di Remigio
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Tromsø - The Arctic University of Norway, N-9037 Tromsø, Norway
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Tommaso Giovannini
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| | | | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Luca Frediani
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Tromsø - The Arctic University of Norway, N-9037 Tromsø, Norway
| |
Collapse
|
9
|
Giovannini T, Lafiosca P, Chandramouli B, Barone V, Cappelli C. Effective yet reliable computation of hyperfine coupling constants in solution by a QM/MM approach: Interplay between electrostatics and non-electrostatic effects. J Chem Phys 2019; 150:124102. [DOI: 10.1063/1.5080810] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | - Piero Lafiosca
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Balasubramanian Chandramouli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
- Compunet, Istituto Italiano di Tecnologia (IIT), Via Morego 30, 16163 Genova, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
10
|
Lao KU, Jia J, Maitra R, DiStasio RA. On the geometric dependence of the molecular dipole polarizability in water: A benchmark study of higher-order electron correlation, basis set incompleteness error, core electron effects, and zero-point vibrational contributions. J Chem Phys 2018; 149:204303. [PMID: 30501247 DOI: 10.1063/1.5051458] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, we investigate how geometric changes influence the static dipole polarizability (α) of a water molecule by explicitly computing the corresponding dipole polarizability surface (DPS) across 3125 total (1625 symmetry-unique) geometries using linear response coupled cluster theory including single, double, and triple excitations (LR-CCSDT) and the doubly augmented triple-ζ basis set (d-aug-cc-pVTZ). Analytical formulae based on power series expansions of this ab initio surface are generated using linear least-squares analysis and provide highly accurate estimates of this quantity as a function of molecular geometry (i.e., bond and angle variations) in a computationally tractable manner. An additional database, which consists of 25 representative molecular geometries and incorporates a more thorough treatment of both basis sets and core electron effects, is provided as a current benchmark for this quantity and the corresponding leading-order C 6 dispersion coefficient. This database has been utilized to assess the importance of these effects as well as the relative accuracy that can be obtained using several quantum chemical methods and a library of density functional approximations. In addition to high-level electron correlation methods (like CCSD) and our analytical least-squares formulae, we find that the SCAN0, PBE0, MN15, and B97-2 hybrid functionals yield the most accurate descriptions of the molecular polarizability tensor in H2O. Using first-order perturbation theory, we compute the zero-point vibrational correction to α at the CCSDT/d-aug-cc-pVTZ level and find that this correction contributes approximately 3% to the isotropic (α iso) and nearly 50% to the anisotropic (α aniso) polarizability values. In doing so, we find that α iso = 9.8307 bohr3, which is in excellent agreement with the experimental value of 9.83 ± 0.02 bohr3 provided by Russell and Spackman. The DPS reported herein provides a benchmark-quality quantum mechanical estimate of this fundamental quantity of interest and should find extensive use in the development (and assessment) of next-generation force fields and machine-learning based approaches for modeling water in complex condensed-phase environments.
Collapse
Affiliation(s)
- Ka Un Lao
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Junteng Jia
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Rahul Maitra
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Robert A DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
11
|
Battisti A, Ambrosetti M, Ruggeri G, Cappelli C, Pucci A. A 4,4'-bis(2-benzoxazolyl)stilbene luminescent probe: assessment of aggregate formation through photophysics experiments and quantum-chemical calculations. Phys Chem Chem Phys 2018; 20:26249-26258. [PMID: 30324198 DOI: 10.1039/c8cp04450c] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
A combination of experimental and quantum mechanical investigations is applied to the study of the optical features of 4,4'-bis(2-benzoxazolyl)stilbene (BBS) dissolved in solution or in a poly(l-lactic acid) (PLA) thermoplastic matrix at different concentrations. The experimental analyses allow the characterization of BBS solutions and dispersions in terms of absorption and emission features, along with the collection of some key parameters such as fluorescence quantum yield, anisotropy and lifetime, while the computational approach gives a detailed description of the photophysical behavior of BBS in the different environments. For the 10-5 M BBS solution, the fluorescence spectra show the expected peaks at 425 and 455 nm of the non-interacting BBS molecules with a single fluorescence lifetime of 0.85 ns without revealing any aggregation phenomena, prevented by the short lifetime and fast diffusion rate of the monomer. Moreover, the calculated spectra are in excellent agreement with the experiments, thus showing the reliability of the computational approach. In time-resolved emission experiments (TRES) on more concentrated solutions (10-4 M) and on BBS crystals, the presence of an excimer is revealed by the appearance of a broad peak around 540 nm, followed by the disappearance of the two main peaks at 460 nm on a time scale of about 10 ns. The computational analysis attributes this behavior to the formation of aggregates of different geometries. The BBS dispersions in PLA reveal the presence of different BBS architectures depending on the fluorophore content. Even at low concentrations, BBS is mainly dispersed as a monomer in the matrix, spheroid aggregates of about 800-900 nm in diameter are also present and the relevant fluorescence spectra arise from the combination of monomer and aggregate contributions. At higher concentrations, BBS starts forming crystals of a peculiar helicoidal shape, with a diameter of about 2 μm, variable length up to several hundreds of μm and emission spectra similar to those of isolated BBS crystals.
Collapse
Affiliation(s)
- Antonella Battisti
- NEST - Scuola Normale Superiore and Istituto Nanoscienze - CNR (CNR-NANO), Piazza San Silvestro 12, I-56127 Pisa, Italy. and M3-Village, Consorzio Interuniversitario Nazionale per la Scienza e Tecnologia dei Materiali - INSTM, p.zza San Silvestro 12, I-56127 Pisa, Italy
| | - Matteo Ambrosetti
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy.
| | - Giacomo Ruggeri
- Dipartimento di Chimica e Chimica Industriale, Universitá di Pisa, Via Moruzzi 13, I-56124 Pisa, Italy.
| | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy.
| | - Andrea Pucci
- Dipartimento di Chimica e Chimica Industriale, Universitá di Pisa, Via Moruzzi 13, I-56124 Pisa, Italy.
| |
Collapse
|
12
|
|
13
|
A polarizable embedding approach to second harmonic generation (SHG) of molecular systems in aqueous solutions. Theor Chem Acc 2018. [DOI: 10.1007/s00214-018-2247-7] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
14
|
Katariya SB, Patil D, Rhyman L, Alswaidan IA, Ramasami P, Sekar N. Triphenylamine-based fluorescent NLO phores with ICT characteristics: Solvatochromic and theoretical study. J Mol Struct 2017. [DOI: 10.1016/j.molstruc.2017.08.084] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
15
|
Rajeshirke M, Katariya SB, Sekar N. Carbazole-Based Mono and Bis-styryl NLOphores: Structure Property Correlations. J SOLUTION CHEM 2017. [DOI: 10.1007/s10953-017-0691-y] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
16
|
Giovannini T, Olszówka M, Egidi F, Cheeseman JR, Scalmani G, Cappelli C. Polarizable Embedding Approach for the Analytical Calculation of Raman and Raman Optical Activity Spectra of Solvated Systems. J Chem Theory Comput 2017; 13:4421-4435. [DOI: 10.1021/acs.jctc.7b00628] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
| | - Marta Olszówka
- Scuola Normale
Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Franco Egidi
- Scuola Normale
Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - James R. Cheeseman
- Gaussian
Inc., 340 Quinnipiac Street Building
40, Wallingford, Connecticut 06492, United States
| | - Giovanni Scalmani
- Gaussian
Inc., 340 Quinnipiac Street Building
40, Wallingford, Connecticut 06492, United States
| | - Chiara Cappelli
- Scuola Normale
Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
17
|
Loukhovitski BI, Sharipov AS, Starik AM. Quantum chemical study of small Al n B m clusters: Structure and physical properties. Chem Phys 2017. [DOI: 10.1016/j.chemphys.2017.06.006] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
18
|
Djorović A, Meyer M, Darby BL, Le Ru EC. Accurate Modeling of the Polarizability of Dyes for Electromagnetic Calculations. ACS OMEGA 2017; 2:1804-1811. [PMID: 31457544 PMCID: PMC6640948 DOI: 10.1021/acsomega.7b00171] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/13/2017] [Accepted: 04/25/2017] [Indexed: 06/10/2023]
Abstract
The wavelength-dependent complex linear polarizability of a dye is a crucial input for the modeling of the optical properties of dye-containing systems. We here propose and discuss methods to obtain an accurate polarizability model by combining absorption spectrum measurements, Kramers-Kronig (KK) tranformations, and density functional theory (DFT) calculations. We focus, in particular, on the real part of the polarizability and its link with static polarizability. In addition, we introduce simple KK-consistent analytic functions based on the theory of critical points as a much more accurate approach to model dye polarizabilities compared with existing models based on Lorentz oscillators. Accurate polarizability models based on critical points and DFT calculations of the static polarizability are derived for five commonly used dyes: Rhodamine 6G, Rhodamine 700, Crystal Violet, Nile Blue A, and Methylene Blue. Finally, we demonstrate explicitly, using examples of Mie Theory calculations of nanoparticle-dye interactions, how an inaccurate polarizability model can result in fundamentally different predictions, further emphasizing the importance of accurate models, such as the one proposed here.
Collapse
|
19
|
Egidi F, Williams-Young DB, Baiardi A, Bloino J, Scalmani G, Frisch MJ, Li X, Barone V. Effective Inclusion of Mechanical and Electrical Anharmonicity in Excited Electronic States: VPT2-TDDFT Route. J Chem Theory Comput 2017; 13:2789-2803. [PMID: 28453287 DOI: 10.1021/acs.jctc.7b00218] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present a reliable and cost-effective procedure for the inclusion of anharmonic effects in excited-state energies and spectroscopic intensities by means of second-order vibrational perturbation theory. This development is made possible thanks to a recent efficient implementation of excited-state analytic Hessians and properties within the time-dependent density functional theory framework. As illustrated in this work, by taking advantage of such algorithmic developments, it is possible to perform calculations of excited-state infrared spectra of medium-large isolated molecular systems, with anharmonicity effects included in both the energy and property surfaces. We also explore the use of this procedure for the inclusion of anharmonic effects in the simulation of vibronic bandshapes of electronic spectra and compare the results with previous, more approximate models.
Collapse
Affiliation(s)
- Franco Egidi
- Scuola Normale Superiore , Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - David B Williams-Young
- Department of Chemistry, University of Washington , Seattle, Washington 98195, United States
| | - Alberto Baiardi
- Scuola Normale Superiore , Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Julien Bloino
- Consiglio Nazionale delle Ricerche, Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR) , UOS di Pisa, Area della Ricerca CNR, Via G. Moruzzi 1, Pisa 56124, Italy
| | - Giovanni Scalmani
- Gaussian, Inc. , 340 Quinnipiac St., Bldg. 40, Wallingford, Connecticut 06492, United States
| | - Michael J Frisch
- Gaussian, Inc. , 340 Quinnipiac St., Bldg. 40, Wallingford, Connecticut 06492, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington , Seattle, Washington 98195, United States
| | - Vincenzo Barone
- Scuola Normale Superiore , Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
20
|
Cannelli O, Giovannini T, Baiardi A, Carlotti B, Elisei F, Cappelli C. Understanding the interplay between the solvent and nuclear rearrangements in the negative solvatochromism of a push–pull flexible quinolinium cation. Phys Chem Chem Phys 2017; 19:32544-32555. [DOI: 10.1039/c7cp07104c] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The main effects (solvation, vibronic progression) affecting the band position and shape of a push–pull flexible quinolinium cation OPA are highlighted.
Collapse
Affiliation(s)
- Oliviero Cannelli
- Scuola Normale Superiore
- Pisa
- Italy
- Laboratoire de Spectroscopie Ultrarapide
- ISIC-FSB
| | | | | | - Benedetta Carlotti
- Department of Chemistry
- Biology and Biotechnology and Centro di Eccellenza sui Materiali Innovativi Nanostrutturati (CEMIN)
- University of Perugia
- via Elce di Sotto 8
- 06123 Perugia
| | - Fausto Elisei
- Department of Chemistry
- Biology and Biotechnology and Centro di Eccellenza sui Materiali Innovativi Nanostrutturati (CEMIN)
- University of Perugia
- via Elce di Sotto 8
- 06123 Perugia
| | | |
Collapse
|
21
|
Giovannini T, Olszòwka M, Cappelli C. Effective Fully Polarizable QM/MM Approach To Model Vibrational Circular Dichroism Spectra of Systems in Aqueous Solution. J Chem Theory Comput 2016; 12:5483-5492. [PMID: 27704812 DOI: 10.1021/acs.jctc.6b00768] [Citation(s) in RCA: 44] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We propose a methodology, based on the combination of classical Molecular Dynamics (MD) simulations with a fully polarizable Quantum Mechanical (QM)/Molecular Mechanics (MM)/Polarizable Continuum Model (PCM) Hamiltonian, to calculate Vibrational Circular Dichroism (VCD) spectra of chiral systems in aqueous solution. Polarization effects are included in the MM force field by exploiting an approach based on Fluctuating Charges (FQ). By performing the MD, the description of the solvating environment is enriched by taking into account the dynamical aspects of the solute-solvent interactions. On the other hand, the QM/FQ/PCM calculation of the VCD spectrum ensures an accurate description of the electronic density of the solute and a proper account for the specific interactions in solution. The application of our approach to (R)-methyloxirane and (l)-alanine in aqueous solution gives calculated spectra in remarkable agreement with their experimental counterparts and a substantial improvement with respect to the same spectra calculated with the PCM.
Collapse
Affiliation(s)
| | - Marta Olszòwka
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa , Via Moruzzi 13, 56124 Pisa, Italy
| | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
22
|
Zaleśny R, Baranowska-Łączkowska A, Medveď M, Luis JM. Comparison of Property-Oriented Basis Sets for the Computation of Electronic and Nuclear Relaxation Hyperpolarizabilities. J Chem Theory Comput 2015; 11:4119-28. [PMID: 26575907 DOI: 10.1021/acs.jctc.5b00434] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
In the present work, we perform an assessment of several property-oriented atomic basis sets in computing (hyper)polarizabilities with a focus on the vibrational contributions. Our analysis encompasses the Pol and LPol-ds basis sets of Sadlej and co-workers, the def2-SVPD and def2-TZVPD basis sets of Rappoport and Furche, and the ORP basis set of Baranowska-Łączkowska and Łączkowski. Additionally, we use the d-aug-cc-pVQZ and aug-cc-pVTZ basis sets of Dunning and co-workers to determine the reference estimates of the investigated electric properties for small- and medium-sized molecules, respectively. We combine these basis sets with ab initio post-Hartree-Fock quantum-chemistry approaches (including the coupled cluster method) to calculate electronic and nuclear relaxation (hyper)polarizabilities of carbon dioxide, formaldehyde, cis-diazene, and a medium-sized Schiff base. The primary finding of our study is that, among all studied property-oriented basis sets, only the def2-TZVPD and ORP basis sets yield nuclear relaxation (hyper)polarizabilities of small molecules with average absolute errors less than 5.5%. A similar accuracy for the nuclear relaxation (hyper)polarizabilites of the studied systems can also be reached using the aug-cc-pVDZ basis set (5.3%), although for more accurate calculations of vibrational contributions, i.e., average absolute errors less than 1%, the aug-cc-pVTZ basis set is recommended. It was also demonstrated that anharmonic contributions to first and second hyperpolarizabilities of a medium-sized Schiff base are particularly difficult to accurately predict at the correlated level using property-oriented basis sets. For instance, the value of the nuclear relaxation first hyperpolarizability computed at the MP2/def2-TZVPD level of theory is roughly 3 times larger than that determined using the aug-cc-pVTZ basis set. We link the failure of the def2-TZVPD basis set with the difficulties in predicting the first-order field-induced coordinates. On the other hand, the aug-cc-pVDZ and ORP basis sets overestimate the property in question only by roughly 30%. In this study, we also propose a low-cost composite treatment of anharmonicity that relies on the combination of two basis sets, i.e., a large-sized basis set is employed to determine lowest-order derivatives with respect to the field-induced coordinates, and a medium-sized basis set is used to compute the higher-order derivatives. The results of calculations performed at the MP2 level of theory demonstrate that this approximate scheme is very successful at predicting nuclear relaxation hyperpolarizabilities.
Collapse
Affiliation(s)
- Robert Zaleśny
- Department of Chemistry, Faculty of Natural Sciences, Matej Bel University , Tajovského 40, 974 01 Banská Bystrica, Slovak Republic.,Department of Physical and Quantum Chemistry, Faculty of Chemistry, Wrocław University of Technology , Wyb. Wyspiańskiego 27, PL-50370 Wrocław, Poland
| | | | - Miroslav Medveď
- Department of Chemistry, Faculty of Natural Sciences, Matej Bel University , Tajovského 40, 974 01 Banská Bystrica, Slovak Republic
| | - Josep M Luis
- Institute of Computational Chemistry and Catalysis and Department of Chemistry, University of Girona , Campus de Montilivi, 17071 Girona, Catalonia, Spain
| |
Collapse
|
23
|
Davari N, Daub CD, Åstrand PO, Unge M. Local Field Factors and Dielectric Properties of Liquid Benzene. J Phys Chem B 2015; 119:11839-45. [DOI: 10.1021/acs.jpcb.5b07043] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Nazanin Davari
- Department
of Chemistry, Norwegian University of Science and Technology (NTNU), NO-7491 Trondheim, Norway
| | - Christopher D. Daub
- Department
of Chemistry, Norwegian University of Science and Technology (NTNU), NO-7491 Trondheim, Norway
| | - Per-Olof Åstrand
- Department
of Chemistry, Norwegian University of Science and Technology (NTNU), NO-7491 Trondheim, Norway
| | - Mikael Unge
- ABB Corporate Research, SE-72178 Västerås, Sweden
| |
Collapse
|
24
|
Loukhovitski BI, Sharipov AS, Starik AM. Physical and Thermodynamic Properties of AlnCm Clusters: Quantum-Chemical Study. J Phys Chem A 2015; 119:1369-80. [DOI: 10.1021/jp5108087] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
- Boris I. Loukhovitski
- Central Institute of Aviation Motors, Moscow, Russia Scientific Educational
Centre “Physical-Chemical Kinetics and Combustion”, Aviamotornaya 2, Moscow 111116, Russia
| | - Alexander S. Sharipov
- Central Institute of Aviation Motors, Moscow, Russia Scientific Educational
Centre “Physical-Chemical Kinetics and Combustion”, Aviamotornaya 2, Moscow 111116, Russia
| | - Alexander M. Starik
- Central Institute of Aviation Motors, Moscow, Russia Scientific Educational
Centre “Physical-Chemical Kinetics and Combustion”, Aviamotornaya 2, Moscow 111116, Russia
| |
Collapse
|
25
|
Benassi E, Egidi F, Barone V. General Strategy for Computing Nonlinear Optical Properties of Large Neutral and Cationic Organic Chromophores in Solution. J Phys Chem B 2015; 119:3155-73. [DOI: 10.1021/jp512342y] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Affiliation(s)
- Enrico Benassi
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Franco Egidi
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|