1
|
Zhou Z, Della Sala F, Parker SM. Minimal Auxiliary Basis Set Approach for the Electronic Excitation Spectra of Organic Molecules. J Phys Chem Lett 2023; 14:1968-1976. [PMID: 36787711 DOI: 10.1021/acs.jpclett.2c03698] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
We report a minimal auxiliary basis model for time-dependent density functional theory (TDDFT) with hybrid density functionals that can accurately reproduce excitation energies and absorption spectra from TDDFT while reducing cost by about 2 orders of magnitude. Our method, dubbed TDDFT-ris, employs the resolution-of-the-identity technique with just one s-type auxiliary basis function per atom for the linear response operator, where the Gaussian exponents are parametrized across the periodic table using tabulated atomic radii with a single global scaling factor. By tuning on a small test set, we determine a single functional-independent scale factor that balances errors in excitation energies and absorption spectra. Benchmarked on organic molecules and compared to standard TDDFT, TDDFT-ris has an average energy error of only 0.06 eV and yields absorption spectra in close agreement with TDDFT. Thus, TDDFT-ris enables simulation of realistic absorption spectra in large molecules that would be inaccessible from standard TDDFT.
Collapse
Affiliation(s)
- Zehao Zhou
- Department of Chemistry, Case Western Reserve University, 10900 Euclid Avenue, Cleveland, Ohio 44106, United States
| | - Fabio Della Sala
- Center for Biomolecular Nanotechnologies, Istituto Italiano di Tecnologia, Via Barsanti 14, 73010 Arnesano (LE), Italy
- Institute for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, 73100 Lecce, Italy
| | - Shane M Parker
- Department of Chemistry, Case Western Reserve University, 10900 Euclid Avenue, Cleveland, Ohio 44106, United States
| |
Collapse
|
2
|
Pei Z, Mao Y, Shao Y, Liang W. Analytic high-order energy derivatives for metal nanoparticle-mediated infrared and Raman scattering spectra within the framework of quantum mechanics/molecular mechanics model with induced charges and dipoles. J Chem Phys 2022; 157:164110. [PMID: 36319412 PMCID: PMC9616608 DOI: 10.1063/5.0118205] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2022] [Accepted: 09/30/2022] [Indexed: 11/14/2022] Open
Abstract
This work is devoted to deriving and implementing analytic second- and third-order energy derivatives with respect to the nuclear coordinates and external electric field within the framework of the hybrid quantum mechanics/molecular mechanics method with induced charges and dipoles (QM/DIM). Using these analytic energy derivatives, one can efficiently compute the harmonic vibrational frequencies, infrared (IR) and Raman scattering (RS) spectra of the molecule in the proximity of noble metal clusters/nanoparticles. The validity and accuracy of these analytic implementations are demonstrated by the comparison of results obtained by the finite-difference method and the analytic approaches and by the full QM and QM/DIM calculations. The complexes formed by pyridine and two sizes of gold clusters (Au18 and Au32) at varying intersystem distances of 3, 4, and 5 Å are used as the test systems, and Raman spectra of 4,4'-bipyridine in the proximity of Au2057 and Ag2057 metal nanoparticles (MNP) are calculated by the QM/DIM method and compared with experimental results as well. We find that the QM/DIM model can well reproduce the IR spectra obtained from full QM calculations for all the configurations, while although it properly enhances some of the vibrational modes, it artificially overestimates RS spectral intensities of several modes for the systems with very short intersystem distance. We show that this could be improved, however, by incorporating the hyperpolarizability of the gold metal cluster in the evaluation of RS intensities. Additionally, we address the potential impact of charge migration between the adsorbate and MNPs.
Collapse
Affiliation(s)
- Zheng Pei
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Collaborative Innovation Center of Chemistry for Energy Materials, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Department of Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People’s Republic of China
| | - Yuezhi Mao
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, USA
| | - WanZhen Liang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Collaborative Innovation Center of Chemistry for Energy Materials, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Department of Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People’s Republic of China
| |
Collapse
|
3
|
Cheng Y, Verstraelen T. A new framework for frequency-dependent polarizable force fields. J Chem Phys 2022; 157:124106. [PMID: 36182425 DOI: 10.1063/5.0115151] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A frequency-dependent extension of the polarizable force field "Atom-Condensed Kohn-Sham density functional theory approximated to the second-order" (ACKS2) [Verstraelen et al., J. Chem. Phys. 141, 194114 (2014)] is proposed, referred to as ACKS2ω. The method enables theoretical predictions of dynamical response properties of finite systems after partitioning of the frequency-dependent molecular response function. Parameters in this model are computed simply as expectation values of an electronic wavefunction, and the hardness matrix is entirely reused from ACKS2 as an adiabatic approximation is used. A numerical validation shows that accurate models can already be obtained with atomic monopoles and dipoles. Absorption spectra of 42 organic and inorganic molecular monomers are evaluated using ACKS2ω, and our results agree well with the time-dependent DFT calculations. Also for the calculation of C6 dispersion coefficients, ACKS2ω closely reproduces its TDDFT reference. When parameters for ACKS2ω are derived from a PBE/aug-cc-pVDZ ground state, it reproduces experimental values for 903 organic and inorganic intermolecular pairs with an MAPE of 3.84%. Our results confirm that ACKS2ω offers a solid connection between the quantum-mechanical description of frequency-dependent response and computationally efficient force-field models.
Collapse
Affiliation(s)
- YingXing Cheng
- Center for Molecular Modeling (CMM), Ghent University, Technologiepark-Zwijnaarde 46, B-9052 Gent, Belgium
| | - Toon Verstraelen
- Center for Molecular Modeling (CMM), Ghent University, Technologiepark-Zwijnaarde 46, B-9052 Gent, Belgium
| |
Collapse
|
4
|
Liang W, Pei Z, Mao Y, Shao Y. Evaluation of molecular photophysical and photochemical properties using linear response time-dependent density functional theory with classical embedding: Successes and challenges. J Chem Phys 2022; 156:210901. [DOI: 10.1063/5.0088271] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Abstract
Time-dependent density functional theory (TDDFT) based approaches have been developed in recent years to model the excited-state properties and transition processes of the molecules in the gas-phase and in a condensed medium, such as in a solution and protein microenvironment or near semiconductor and metal surfaces. In the latter case, usually, classical embedding models have been adopted to account for the molecular environmental effects, leading to the multi-scale approaches of TDDFT/polarizable continuum model (PCM) and TDDFT/molecular mechanics (MM), where a molecular system of interest is designated as the quantum mechanical region and treated with TDDFT, while the environment is usually described using either a PCM or (non-polarizable or polarizable) MM force fields. In this Perspective, we briefly review these TDDFT-related multi-scale models with a specific emphasis on the implementation of analytical energy derivatives, such as the energy gradient and Hessian, the nonadiabatic coupling, the spin–orbit coupling, and the transition dipole moment as well as their nuclear derivatives for various radiative and radiativeless transition processes among electronic states. Three variations of the TDDFT method, the Tamm–Dancoff approximation to TDDFT, spin–flip DFT, and spin-adiabatic TDDFT, are discussed. Moreover, using a model system (pyridine–Ag20 complex), we emphasize that caution is needed to properly account for system–environment interactions within the TDDFT/MM models. Specifically, one should appropriately damp the electrostatic embedding potential from MM atoms and carefully tune the van der Waals interaction potential between the system and the environment. We also highlight the lack of proper treatment of charge transfer between the quantum mechanics and MM regions as well as the need for accelerated TDDFT modelings and interpretability, which calls for new method developments.
Collapse
Affiliation(s)
- WanZhen Liang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Collaborative Innovation Center of Chemistry for Energy Materials, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Department of Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People’s Republic of China
| | - Zheng Pei
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Collaborative Innovation Center of Chemistry for Energy Materials, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Department of Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People’s Republic of China
| | - Yuezhi Mao
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, USA
| |
Collapse
|
5
|
Makarov DE. Perspective: Mechanochemistry of biological and synthetic molecules. J Chem Phys 2016; 144:030901. [PMID: 26801011 DOI: 10.1063/1.4939791] [Citation(s) in RCA: 59] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Abstract
Coupling of mechanical forces and chemical transformations is central to the biophysics of molecular machines, polymer chemistry, fracture mechanics, tribology, and other disciplines. As a consequence, the same physical principles and theoretical models should be applicable in all of those fields; in fact, similar models have been invoked (and often repeatedly reinvented) to describe, for example, cell adhesion, dry and wet friction, propagation of cracks, and action of molecular motors. This perspective offers a unified view of these phenomena, described in terms of chemical kinetics with rates of elementary steps that are force dependent. The central question is then to describe how the rate of a chemical transformation (and its other measurable properties such as the transition path) depends on the applied force. I will describe physical models used to answer this question and compare them with experimental measurements, which employ single-molecule force spectroscopy and which become increasingly common. Multidimensionality of the underlying molecular energy landscapes and the ensuing frequent misalignment between chemical and mechanical coordinates result in a number of distinct scenarios, each showing a nontrivial force dependence of the reaction rate. I will discuss these scenarios, their commonness (or its lack), and the prospects for their experimental validation. Finally, I will discuss open issues in the field.
Collapse
Affiliation(s)
- Dmitrii E Makarov
- Department of Chemistry and Institute for Computational Engineering and Sciences, University of Texas at Austin, Austin, Texas 78712, USA
| |
Collapse
|
6
|
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
|
7
|
Davari N, Haghdani S, Åstrand PO, Schatz GC. Local electric field factors by a combined charge-transfer and point–dipole interaction model. RSC Adv 2015. [DOI: 10.1039/c5ra04183j] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
A model for the local electric field as a linear response to a frequency-dependent external electric field is presented based on a combined charge-transfer and point–dipole interaction force-field model.
Collapse
Affiliation(s)
- Nazanin Davari
- Department of Chemistry
- Norwegian University of Science and Technology (NTNU)
- Trondheim
- Norway
| | - Shokouh Haghdani
- Department of Chemistry
- Norwegian University of Science and Technology (NTNU)
- Trondheim
- Norway
| | - Per-Olof Åstrand
- Department of Chemistry
- Norwegian University of Science and Technology (NTNU)
- Trondheim
- Norway
| | | |
Collapse
|
8
|
Haghdani S, Davari N, Sandnes R, Åstrand PO. Complex Frequency-Dependent Polarizability through the π → π* Excitation Energy of Azobenzene Molecules by a Combined Charge-Transfer and Point-Dipole Interaction Model. J Phys Chem A 2014; 118:11282-92. [DOI: 10.1021/jp507639z] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Shokouh Haghdani
- Department of Chemistry, Norwegian University of Science and Technology (NTNU), N-7491 Trondheim, Norway
| | - Nazanin Davari
- Department of Chemistry, Norwegian University of Science and Technology (NTNU), N-7491 Trondheim, Norway
| | - Runar Sandnes
- Department of Chemistry, Norwegian University of Science and Technology (NTNU), N-7491 Trondheim, Norway
| | - Per-Olof Åstrand
- Department of Chemistry, Norwegian University of Science and Technology (NTNU), N-7491 Trondheim, Norway
| |
Collapse
|