Excited-State Properties for Extended Systems: Efficient Hybrid Density Functional MethodsClick to copy article linkArticle link copied!
- Anna-Sophia Hehn*Anna-Sophia Hehn*A.-S. Hehn. Email: [email protected]Department of Chemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, SwitzerlandMore by Anna-Sophia Hehn
- Beliz SertcanBeliz SertcanDepartment of Chemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, SwitzerlandMore by Beliz Sertcan
- Fabian BelleflammeFabian BelleflammeDepartment of Chemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, SwitzerlandMore by Fabian Belleflamme
- Sergey K. ChulkovSergey K. ChulkovSchool of Mathematics and Physics, University of Lincoln, Brayford Pool, Lincoln LN67TS, United KingdomMore by Sergey K. Chulkov
- Matthew B. WatkinsMatthew B. WatkinsSchool of Mathematics and Physics, University of Lincoln, Brayford Pool, Lincoln LN67TS, United KingdomMore by Matthew B. Watkins
- Jürg HutterJürg HutterDepartment of Chemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, SwitzerlandMore by Jürg Hutter
Abstract
Time-dependent density functional theory has become state-of-the-art for describing photophysical and photochemical processes in extended materials because of its affordable cost. The inclusion of exact exchange was shown to be essential for the correct description of the long-range asymptotics of electronic interactions and thus a well-balanced description of valence, Rydberg, and charge-transfer excitations. Several approaches for an efficient treatment of exact exchange have been established for the ground state, while implementations for excited-state properties are rare. Furthermore, the high computational costs required for excited-state properties in comparison to ground-state computations often hinder large-scale applications on periodic systems with hybrid functional accuracy. We therefore propose two approximate schemes for improving computational efficiency for the treatment of exact exchange. Within the auxiliary density matrix method (ADMM), exact exchange is estimated using a relatively small auxiliary basis and the introduced basis set incompleteness error is compensated by an exchange density functional correction term. Benchmark results for a test set of 35 molecules demonstrate that the mean absolute error introduced by ADMM is smaller than 0.3 pm for excited-state bond lengths and in the range of 0.02–0.04 eV for vertical excitation, adiabatic excitation, and fluorescence energies. Computational timings for a series of covalent-organic frameworks demonstrate that a speed-up of at least 1 order of magnitude can be achieved for excited-state geometry optimizations in comparison to conventional hybrid functionals. The second method is to use a semiempirical tight binding approximation for both Coulomb and exchange contributions to the excited-state kernel. This simplified Tamm–Dancoff approximation (sTDA) achieves an accuracy comparable to approximated hybrid density functional theory when referring to highly accurate coupled-cluster reference data. We find that excited-state bond lengths deviate by 1.1 pm on average and mean absolute errors in vertical excitation, adiabatic excitation, and fluorescence energies are in the range of 0.2–0.5 eV. In comparison to ADMM-approximated hybrid functional theory, sTDA accelerates the computation of broad-band excitation spectra by 1 order of magnitude, suggesting its potential use for large-scale screening purposes.
This publication is licensed under
License Summary*
You are free to share(copy and redistribute) this article in any medium or format within the parameters below:
Creative Commons (CC): This is a Creative Commons license.
Attribution (BY): Credit must be given to the creator.
Non-Commercial (NC): Only non-commercial uses of the work are permitted.
No Derivatives (ND): Derivative works may be created for non-commercial purposes, but sharing is prohibited.
*Disclaimer
This summary highlights only some of the key features and terms of the actual license. It is not a license and has no legal value. Carefully review the actual license before using these materials.
License Summary*
You are free to share(copy and redistribute) this article in any medium or format within the parameters below:
Creative Commons (CC): This is a Creative Commons license.
Attribution (BY): Credit must be given to the creator.
Non-Commercial (NC): Only non-commercial uses of the work are permitted.
No Derivatives (ND): Derivative works may be created for non-commercial purposes, but sharing is prohibited.
*Disclaimer
This summary highlights only some of the key features and terms of the actual license. It is not a license and has no legal value. Carefully review the actual license before using these materials.
License Summary*
You are free to share(copy and redistribute) this article in any medium or format within the parameters below:
Creative Commons (CC): This is a Creative Commons license.
Attribution (BY): Credit must be given to the creator.
Non-Commercial (NC): Only non-commercial uses of the work are permitted.
No Derivatives (ND): Derivative works may be created for non-commercial purposes, but sharing is prohibited.
*Disclaimer
This summary highlights only some of the key features and terms of the actual license. It is not a license and has no legal value. Carefully review the actual license before using these materials.
License Summary*
You are free to share(copy and redistribute) this article in any medium or format within the parameters below:
Creative Commons (CC): This is a Creative Commons license.
Attribution (BY): Credit must be given to the creator.
Non-Commercial (NC): Only non-commercial uses of the work are permitted.
No Derivatives (ND): Derivative works may be created for non-commercial purposes, but sharing is prohibited.
*Disclaimer
This summary highlights only some of the key features and terms of the actual license. It is not a license and has no legal value. Carefully review the actual license before using these materials.
License Summary*
You are free to share(copy and redistribute) this article in any medium or format within the parameters below:
Creative Commons (CC): This is a Creative Commons license.
Attribution (BY): Credit must be given to the creator.
Non-Commercial (NC): Only non-commercial uses of the work are permitted.
No Derivatives (ND): Derivative works may be created for non-commercial purposes, but sharing is prohibited.
*Disclaimer
This summary highlights only some of the key features and terms of the actual license. It is not a license and has no legal value. Carefully review the actual license before using these materials.
1. Introduction
2. Theoretical Background
2.1. The Tamm–Dancoff Approximation
2.2. Exact Exchange Using the Auxiliary Density Matrix Method (ADMM)
2.3. Semiempirical Coulomb and Exchange Contributions within the Simplified Tamm–Dancoff Approximation (sTDA)
2.4. Periodic Boundary Conditions
2.5. The Excited-State Lagrangian in the Tamm–Dancoff Approximation
3. Results and Discussion
3.1. Tests for Accuracy: Benchmark Results for 35 Main-Group Molecules
3.1.1. Excited-State Geometries: Impact of the Auxiliary Basis Set Size for ADMM Kernels
Figure 1
Figure 1. Error (in pm) for selected bond lengths for (a) optimized ground- and (b) first singlet excited-state geometries comparing conventional PBE0 and ADMM-PBE0 computations for different auxiliary (ABS) and primary basis sets (PBS) indicated as ABS/PBS.
Figure 2
Figure 2. Error (in deg) for selected angles for optimized first singlet excited-state geometries comparing conventional PBE0 and ADMM-PBE0 computations.
Bonds (pm) | Angles (deg) | Dihedrals (deg) | |||||||
---|---|---|---|---|---|---|---|---|---|
ABS/PBS | ME | MAE | STD | ME | MAE | STD | ME | MAE | STD |
dzp/DZVP | 1.31 | 1.31 | 0.55 | –0.74 | 0.85 | 1.19 | 2.07 | 2.07 | 1.20 |
dzp/TZVP | 1.37 | 1.37 | 0.57 | –0.75 | 0.85 | 1.13 | 1.94 | 1.94 | 1.11 |
dzp/TZV2P | 1.38 | 1.38 | 0.56 | –0.73 | 0.83 | 1.18 | 2.05 | 2.05 | 1.23 |
tzp/TZVP | 0.21 | 0.23 | 0.25 | –0.11 | 0.21 | 0.35 | 0.07 | 0.34 | 0.54 |
tzp/TZV2P | 0.25 | 0.26 | 0.24 | –0.13 | 0.20 | 0.32 | 0.30 | 0.36 | 0.49 |
tz2p/TZV2P | 0.21 | 0.22 | 0.26 | –0.10 | 0.18 | 0.28 | 0.02 | 0.13 | 0.16 |
Evert (eV) | Ead (eV) | Efl (eV) | |||||||
---|---|---|---|---|---|---|---|---|---|
ABS/PBS | ME | MAE | STD | ME | MAE | STD | ME | MAE | STD |
dzp/DZVP | –0.064 | 0.095 | 0.117 | –0.023 | 0.048 | 0.054 | –0.055 | 0.080 | 0.074 |
dzp/TZVP | –0.058 | 0.091 | 0.109 | –0.025 | 0.055 | 0.061 | –0.057 | 0.083 | 0.075 |
dzp/TZV2P | –0.062 | 0.093 | 0.109 | –0.028 | 0.054 | 0.058 | –0.058 | 0.081 | 0.073 |
tzp/TZVP | –0.012 | 0.036 | 0.057 | 0.004 | 0.022 | 0.031 | –0.002 | 0.022 | 0.030 |
tzp/TZV2P | –0.015 | 0.035 | 0.048 | 0.004 | 0.022 | 0.032 | 0.003 | 0.022 | 0.029 |
tz2p/TZV2P | –0.021 | 0.033 | 0.044 | –0.004 | 0.016 | 0.022 | –0.005 | 0.019 | 0.025 |
3.1.2. Excited-State Geometries: ADMM and sTDA Kernels in Comparison to EOM-CCSD References
Figure 3
Figure 3. Normal distributions based on the mean errors (ME) and standard deviations (STD) wrt EOM-CCSD reference geometries in (pm)/(deg)/(deg) for selected (a) bond lengths/(b) angles/(c) dihedral angles for optimized first singlet excited-state geometries comparing conventional PBE0, ADMM-PBE0, and sTDA computations using triple-ζ basis sets.
3.1.3. Vertical Excitation Energies
Bonds (pm) | Angles (deg) | Dihedrals (deg) | ||||||||
---|---|---|---|---|---|---|---|---|---|---|
Kernel | PBS+ABS | ME | MAE | STD | ME | MAE | STD | ME | MAE | STD |
PBE0 | TZVP | –0.15 | 0.93 | 1.18 | 0.32 | 1.24 | 2.20 | 0.80 | 2.34 | 3.21 |
ADMM-PBE0 | TZVP+tzp | 0.06 | 0.92 | 1.18 | 0.20 | 1.18 | 2.09 | 0.86 | 2.31 | 3.46 |
sTDA@PBE0 | TZVP | 0.10 | 1.06 | 1.50 | –0.03 | 1.60 | 2.57 | 3.98 | 4.06 | 3.95 |
sTDA@ADMM-PBE0 | TZVP+tzp | 0.28 | 1.03 | 1.35 | –0.05 | 1.51 | 2.34 | 4.03 | 4.14 | 4.06 |
Evert (eV) | Ead (eV) | Efl (eV) | ||||||||
---|---|---|---|---|---|---|---|---|---|---|
ME | MAE | STD | ME | MAE | STD | ME | MAE | STD | ||
PBE0 | DZVP | –0.02 | 0.05 | 0.06 | –0.16 | 0.30 | 0.42 | –0.05 | 0.06 | 0.05 |
PBE0 | TZVP | –0.01 | 0.03 | 0.03 | –0.15 | 0.29 | 0.41 | –0.04 | 0.05 | 0.04 |
PBE0 | TZV2P | 0.00 | 0.02 | 0.03 | –0.13 | 0.28 | 0.40 | –0.02 | 0.02 | 0.02 |
ADMM-PBE0 | TZVP+tzp | –0.01 | 0.03 | 0.04 | –0.14 | 0.29 | 0.41 | –0.05 | 0.07 | 0.06 |
sTDA@PBE0 | TZVP | 0.45 | 0.45 | 0.18 | 0.23 | 0.34 | 0.39 | 0.27 | 0.30 | 0.21 |
sTDA@ADMM-PBE0 | TZVP+tzp | 0.34 | 0.34 | 0.17 | 0.10 | 0.33 | 0.44 | 0.14 | 0.18 | 0.16 |
Errors in the geometrical data and Ead are calculated with respect to the EOM-CCSD reference data of ref (51); errors in Evert and Efl are referring to PBE0/def2-QZVPP/TURBOMOLE reference computations.
Figure 4
Figure 4. Normal distributions based on the MEs and STDs for vertical excitation energies (in eV) depicting (a) the error introduced by ADMM-PBE0 in comparison to conventional PBE0 computations (upper left plot), (b) the performance of PBE0, ADMM-PBE0, and sTDA kernel computations using double- and triple-ζ basis sets in comparison to a PBE0/TZV2P/CP2K reference (upper right plot) as well as (c) an analogous assessment of the kernels with respect to a PBE0/def2-QZVPP/Turbomole reference (lower plot).
3.1.4. Adiabatic Excitation Energies and Fluorescence
Figure 5
Figure 5. ADMM error distribution in (a) adiabatic excitation and (b) fluorescence energies comparing approximated ADMM-PBE0 and conventional PBE0 results for different auxiliary basis set sizes (in eV).
Figure 6
Figure 6. Comparison of (a) adiabatic excitation energies with respect to EOM-CCSD reference data and (b) fluorescence energies with respect to PBE0/def2-QZVPP/Turbomole reference data for PBE0, ADMM-PBE0, and sTDA kernels (in eV).
3.2. Increasing Computational Efficiency: Treating Extended Systems Using ADMM-PBE0 and sTDA Kernels
PBE0 | ADMM-PBE0 | sTDA@ADMM-PBE0 | |
---|---|---|---|
GS SCF | 301 | 137 | 134 |
ES energy and gradient | 1919 | 403 | 148 |
ERI for ES energy | 687 | 86 | 33 |
DERI for ES gradient | 1265 | 11 | 7 |
Total computation time | 2220 | 540 | 282 |
Computations were performed using a ccGRB-T primary basis (2076 basis functions) and a tzp auxiliary basis (1164 basis functions). Timings are reported for a Intel Xeon E5-2670 processor using in total 4608 cores, analysing the cost for converging the GS SCF computation, for calculating the ES energy and gradient as well as for the computation of the ERIs needed for the ES energy calculation and the DERIs needed for the ES gradient.
Absorption spectra | ES optimization | |||||||
---|---|---|---|---|---|---|---|---|
ADMM | sTDA | ADMM | sTDA | |||||
# Atoms | # Basis | Time | Time | Time | Time/Step | Time | Time/Step | |
05000N2 | 84 | 2076 | 15758 | 2157 | 9590 | 461 | 8937 | 329 |
05001N2 | 192 | 4848 | 33678 | 2472 | 77332 | 1556 | 33837 | 637 |
15100N2 | 300 | 7440 | 36907 | 2623 | 155696 | 1729 | 37804 | 433 |
15101N2 | 192 | 5028 | 30931 | 1932 | 67962 | 1680 | 29258 | 417 |
20610N2 | 240 | 5700 | 30014 | 2119 | 59603 | 1086 | 16627 | 266 |
20611N2 | 262 | 6188 | 27965 | 1993 | 136122 | 1818 | 35955 | 362 |
Timings are reported for a Intel Xeon E5-2670 processor using in total 2304 cores.
Evert | Ead | Efl | λ | |||||
---|---|---|---|---|---|---|---|---|
ADMM | sTDA | ADMM | sTDA | ADMM | sTDA | ADMM | sTDA | |
05000N2 | 4.71 | 4.95 | 4.40 | 3.87 | 4.49 | 4.16 | 0.22 | 0.79 |
05001N2 | 3.77 | 3.44 | 2.60 | 2.14 | 2.97 | 2.56 | 0.80 | 0.88 |
15100N2 | 3.54 | 3.08 | 2.27 | 1.85 | 2.43 | 1.98 | 1.11 | 1.10 |
15101N2 | 3.10 | 2.67 | 2.04 | 1.46 | 2.18 | 1.68 | 0.92 | 0.98 |
20610N2 | 3.97 | 3.67 | 2.82 | 2.45 | 3.08 | 2.74 | 0.89 | 0.93 |
20611N2 | 3.21 | 3.10 | 2.47 | 1.96 | 2.65 | 2.25 | 0.56 | 0.85 |
MAE | 0.31 | 0.42 | 0.40 | 0.17 |
Figure 7
Figure 7. Comparison of ADMM-PBE0 and sTDA absorption spectra for a series of fluorescent COFs.
4. Conclusions
Supporting Information
The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jctc.2c00144.
Treatment of the virtual space, structural information for molecular benchmark set, geometrical data on optimized molecular geometries, statistical analysis for ADMM and sTDA vertical excitation energies, adiabatic excitation energies and statistical analysis, fluorescence energies and statistical analysis, and structural information for covalent organic frameworks (PDF)
Terms & Conditions
Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.
Acknowledgments
We thank Dr. Augustin Bussy (University of Zurich) for his kind support and the valuable discussions. This work was supported in parts by the MARVEL National Centre for Competence in Research funded by the Swiss National Science Foundation (grant agreement ID 51NF40-182892). We thank the Swiss National Supercomputing Center (CSCS) for providing computational resources (Grant No. S1109). Furthermore, this work has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement no. 798196. M.W. and S.C. were assisted by membership of the UK’s HEC Materials Chemistry Consortium, which is funded by EPSRC (EP/R029431).
References
This article references 79 other publications.
- 1González, L.; Escudero, D.; Serrano-Andrés, L. Progress and Challenges in the Calculation of Electronic Excited States. ChemPhysChem 2012, 13, 28– 51, DOI: 10.1002/cphc.201100200Google Scholar1https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC38XjvVGntw%253D%253D&md5=710b04d12866b33b409e5377a8d3be72Progress and Challenges in the Calculation of Electronic Excited StatesGonzalez, Leticia; Escudero, Daniel; Serrano-Andres, LuisChemPhysChem (2012), 13 (1), 28-51CODEN: CPCHFT; ISSN:1439-4235. (Wiley-VCH Verlag GmbH & Co. KGaA)A review. A detailed understanding of the properties of electronic excited states and the reaction mechanisms that mols. undergo after light irradn. is a fundamental ingredient for following light-driven natural processes and for designing novel photonic materials. The aim of this review is to present an overview of the ab initio quantum chem. and time-dependent d. functional theory methods that can be used to model spectroscopy and photochem. in mol. systems. The applicability and limitations of the different methods as well as the main frontiers are discussed. To illustrate the progress achieved by excited-state chem. in the recent years as well as the main challenges facing computational chem., three main applications that reflect the authors' experience are addressed: the UV/Vis spectroscopy of org. mols., the assignment of absorption and emission bands of organometallic complexes, and finally, the obtainment of non-adiabatic photoinduced pathways mediated by conical intersections. In the latter case, special emphasis is put on the photochem. of DNA. These applications show that the description of electronically excited states is a rewarding but challenging area of research.
- 2Dreuw, A.; Head-Gordon, M. Single-Reference ab Initio Methods for the Calculation of Excited States of Large Molecules. Chem. Rev. 2005, 105, 4009– 4037, DOI: 10.1021/cr0505627Google Scholar2https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2MXhtVKmtLrK&md5=31bf0c10a174b6d5f2aa914bb0118957Single-Reference ab Initio Methods for the Calculation of Excited States of Large MoleculesDreuw, Andreas; Head-Gordon, MartinChemical Reviews (Washington, DC, United States) (2005), 105 (11), 4009-4037CODEN: CHREAY; ISSN:0009-2665. (American Chemical Society)A review of single-ref. ab initio excited state methods, which are applicable to large mols. and do not explicitly include correlation through the ground-state wave function. CIS, TDHF, and TDDFT are rigorously introduced by outlining their derivations and theor. footing, where special emphasis is put on their relations to each other. Different methods for the anal. of complicated electronically excited states are reviewed and comparatively discussed. The applicability of the presented methods is outlined, their limitations are show, and out their strengths and weaknesses are pointed out.
- 3Laurent, A. D.; Jacquemin, D. TD-DFT benchmarks: A review. Int. J. Quantum Chem. 2013, 113, 2019– 2039, DOI: 10.1002/qua.24438Google Scholar3https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3sXlsVGgt74%253D&md5=2d81e67d518499c1a10f72ec369ec6aeTD-DFT benchmarks: A reviewLaurent, Adele D.; Jacquemin, DenisInternational Journal of Quantum Chemistry (2013), 113 (17), 2019-2039CODEN: IJQCB2; ISSN:0020-7608. (John Wiley & Sons, Inc.)A review. Time-Dependent D. Functional Theory (TD-DFT) has become the most widely-used theor. approach to simulate the optical properties of both org. and inorg. mols. In this contribution, we review TD-DFT benchmarks that have been performed during the last decade. The aim is often to pinpoint the most accurate or adequate exchange-correlation functional(s). We present both the different strategies used to assess the functionals and the main results obtained in terms of accuracy. In particular, we discuss both vertical and adiabatic benchmarks and comparisons with both exptl. and theor. ref. transition energies. More specific benchmarks (oscillator strengths, excited-state geometries, dipole moments, vibronic shapes, etc.) are summarized as well. © 2013 Wiley Periodicals, Inc.
- 4de Wergifosse, M.; Grimme, S. Perspective on Simplified Quantum Chemistry Methods for Excited States and Response Properties. J. Phys. Chem. A 2021, 125, 3841– 3851, DOI: 10.1021/acs.jpca.1c02362Google Scholar4https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3MXpsl2ltLo%253D&md5=279421848df210a3de326d8bd50e783cPerspective on Simplified Quantum Chemistry Methods for Excited States and Response PropertiesDe wergifosse, Marc; Grimme, StefanJournal of Physical Chemistry A (2021), 125 (18), 3841-3851CODEN: JPCAFH; ISSN:1089-5639. (American Chemical Society)A review. We review recent developments in the framework of simplified quantum chem. for excited state and optical response properties (sTD-DFT) and present future challenges for new method developments to improve accuracy and extend the range of application. In recent years, the scope of sTD-DFT was extended to mol. response calcns. of the polarizability, optical rotation, first hyperpolarizability, two-photon absorption (2PA), and excited-state absorption for large systems with hundreds to thousands of atoms. The recently introduced spin-flip simplified time-dependent d. functional theory (SF-sTD-DFT) variant enables an ultrafast treatment for diradicals and related strongly correlated systems. A few drawbacks were also identified, specifically for the computation of 2PA cross sections. We propose solns. to this problem and how to generally improve the accuracy of simplified schemes. New possible simplified schemes are also introduced for strongly correlated systems, e.g., with a second-order perturbative correlation correction. Interpretation tools that can ext. chem. structure-property relationships from excited state or response calcns. are also discussed. In particular, the recently introduced method-agnostic RespA approach based on natural response orbitals (NROs) as the key concept is employed.
- 5Grimme, S.; Schreiner, P. R. Computational Chemistry: The Fate of Current Methods and Future Challenges. Angew. Chem., Int. Ed. 2018, 57, 4170– 4176, DOI: 10.1002/anie.201709943Google Scholar5https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhvFGitrzO&md5=4890cef4031e2238fd8db2f01049cbedComputational Chemistry: The Fate of Current Methods and Future ChallengesGrimme, Stefan; Schreiner, Peter R.Angewandte Chemie, International Edition (2018), 57 (16), 4170-4176CODEN: ACIEF5; ISSN:1433-7851. (Wiley-VCH Verlag GmbH & Co. KGaA)In this essay, we attempt to make predictions about the fate and development of the computational mol. sciences. Of course, it is not the first time that the future challenges for computational org. chem. and biochem. are considered; these were outlined in complementary contexts recently. In this article, the authors take a somewhat different perspective and emphasize the changes expected for chem. that are triggered by the rapid developments and increasingly stronger influences from theory, algorithms, and data-driven technologies. The authors of this essay are about the same age and have a general overview of a period of about 25 years during which they have actively contributed to the field of computational chem. Hence, it appears sensible to make predictions extending 25 years into the future, roughly to the year 2043 (when both authors will long be retired).
- 6Dierksen, M.; Grimme, S. The Vibronic Structure of Electronic Absorption Spectra of Large Molecules: A Time-Dependent Density Functional Study on the Influence of Exact Hartree-Fock Exchange. J. Phys. Chem. A 2004, 108, 10225– 10237, DOI: 10.1021/jp047289hGoogle Scholar6https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2cXotlKktLk%253D&md5=b85a6bb8357480b91217f853af23c6a3The Vibronic Structure of Electronic Absorption Spectra of Large Molecules: A Time-Dependent Density Functional Study on the Influence of "Exact" Hartree-Fock ExchangeDierksen, Marc; Grimme, StefanJournal of Physical Chemistry A (2004), 108 (46), 10225-10237CODEN: JPCAFH; ISSN:1089-5639. (American Chemical Society)The functional dependence of excited-state geometries and normal modes calcd. with time-dependent d. functional theory (TDDFT) is studied from vibronic structure calcns. of the absorption spectra of large mols. For a set of mols. covering a wide range of different structures including org. dyes, biol. chromophores, and mols. of importance in material science, quantum mech. simulations of the vibronic structure are performed. In total over 40 singlet-singlet transitions of neutral closed-shell compds. and doublet-doublet transitions of neutral radicals, radical cations, and anions are considered. Calcns. with different std. d. functionals show that the predicted vibronic structure critically depends on the fraction of the exact Hartree-Fock exchange (EEX) included in hybrid functionals. The effect can been traced back to a large influence of EEX on the geometrical displacement upon excitation. On the contrary, the dependence of the results on the choice of the local exchange-correlation functional is rather small. From detailed comparisons with exptl. spectra conclusions are drawn concerning the optimum amt. of EEX mixing for a proper description of the excited-state properties. The relation of the quality of the simulated spectra with the errors for 0-0 transition energies is discussed. For the studied singlet-singlet π → π* transitions and the 1st strongly dipole-allowed transitions of PAH radical cations some rules of thumb concerning the optimum portion of EEX are derived. However, in general no universal amt. of EEX seems to exist that gives a uniformly good description for all systems and states. Nevertheless an inclusion of ∼30-40% of EEX in the functional is found empirically to yield in most cases simulated spectra that compare very well with those from expt. and thus seems to be necessary for an accurate description of the excited-state geometry. Pure d. functionals that are computationally more efficient provide less accurate spectra in most cases and their application is recommended solely for comparison purposes to obtain ests. for the reliability of the theor. predictions.
- 7Merlot, P.; Izsák, R.; Borgoo, A.; Kjærgaard, T.; Helgaker, T.; Reine, S. Charge-constrained auxiliary-density-matrix methods for the Hartree-Fock exchange contribution. J. Chem. Phys. 2014, 141, 094104, DOI: 10.1063/1.4894267Google Scholar7https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXhsFSjtrzP&md5=0d138a12fe1e543edb43a3e92af05f10Charge-constrained auxiliary-density-matrix methods for the Hartree-Fock exchange contributionMerlot, Patrick; Izsak, Robert; Borgoo, Alex; Kjaergaard, Thomas; Helgaker, Trygve; Reine, SimenJournal of Chemical Physics (2014), 141 (9), 094104/1-094104/11CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Three new variants of the auxiliary-d.-matrix method (ADMM) of Guidon, Hutter, and VandeVondele [J. Chem. Theory Comput. 6, 2348 (2010)] are presented with the common feature that they have a simplified constraint compared with the full orthonormality requirement of the earlier ADMM1 method. All ADMM variants are tested for accuracy and performance in all-electron B3LYP calcns. with several commonly used basis sets. The effect of the choice of the exchange functional for the ADMM exchange-correction term is also investigated. (c) 2014 American Institute of Physics.
- 8Guidon, M.; Hutter, J.; VandeVondele, J. Auxiliary Density Matrix Methods for Hartree-Fock Exchange Calculations. J. Chem. Theory Comput. 2010, 6, 2348– 2364, DOI: 10.1021/ct1002225Google Scholar8https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXot1als7Y%253D&md5=7a84c086d49f17bab239fd3f8d004c6eAuxiliary Density Matrix Methods for Hartree-Fock Exchange CalculationsGuidon, Manuel; Hutter, Jurg; Vande Vondele, JoostJournal of Chemical Theory and Computation (2010), 6 (8), 2348-2364CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)The calcn. of Hartree-Fock exchange (HFX) is computationally demanding for large systems described with high-quality basis sets. In this work, we show that excellent performance and good accuracy can nevertheless be obtained if an auxiliary d. matrix is employed for the HFX calcn. Several schemes to derive an auxiliary d. matrix from a high-quality d. matrix are discussed. Key to the accuracy of the auxiliary d. matrix methods (ADMM) is the use of a correction based on std. generalized gradient approxns. for HFX. ADMM integrates seamlessly in existing HFX codes and, in particular, can be employed in linear scaling implementations. Demonstrating the performance of the method, the effect of HFX on the structure of liq. water is investigated in detail using Born-Oppenheimer mol. dynamics simulations (300 ps) of a system of 64 mols. Representative for large systems are calcns. on a solvated protein (Rubredoxin), for which ADMM outperforms the corresponding std. HFX implementation by approx. a factor 20.
- 9Kumar, C.; Fliegl, H.; Jensen, F.; Teale, A. M.; Reine, S.; Kjærgaard, T. Accelerating Kohn-Sham response theory using density fitting and the auxiliary-density-matrix method. Int. J. Quantum Chem. 2018, 118, e25639 DOI: 10.1002/qua.25639Google ScholarThere is no corresponding record for this reference.
- 10Rebolini, E.; Izsák, R.; Reine, S. S.; Helgaker, T.; Pedersen, T. B. Comparison of Three Efficient Approximate Exact-Exchange Algorithms: The Chain-of-Spheres Algorithm, Pair-Atomic Resolution-of-the-Identity Method, and Auxiliary Density Matrix Method. J. Chem. Theory Comput. 2016, 12, 3514– 3522, DOI: 10.1021/acs.jctc.6b00074Google Scholar10https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28Xos1Sku7s%253D&md5=dccdeec570871ad4ed29846f8f6d70a4Comparison of Three Efficient Approximate Exact-Exchange Algorithms: The Chain-of-Spheres Algorithm, Pair-Atomic Resolution-of-the-Identity Method, and Auxiliary Density Matrix MethodRebolini, Elisa; Izsak, Robert; Reine, Simen Sommerfelt; Helgaker, Trygve; Pedersen, Thomas BondoJournal of Chemical Theory and Computation (2016), 12 (8), 3514-3522CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)We compare the performance of three approx. methods for speeding up evaluation of the exchange contribution in Hartree-Fock and hybrid Kohn-Sham calcns.: the chain-of-spheres algorithm, and the auxiliary d. matrix method. Both the efficiency relative to that of a conventional linear-scaling algorithm and the accuracy of total, atomization, and orbital energies are compared for a subset contg. 25 of the 200 mols. in the Rx200 set using double-, triple-, and quadruple-ζ basis sets. The accuracy of relative energies is further compared for small alkane conformers (ACONF test set) and Diels-Alder reactions (DARC test set). Overall, we find that the COSX method provides good accuracy for orbital energies as well as total and relative energies, and the method delivers a satisfactory speedup. The PARI-K and in particular ADMM algorithms require further development and optimization to fully exploit their indisputable potential.
- 11Merlot, P.; Kjærgaard, T.; Helgaker, T.; Lindh, R.; Aquilante, F.; Reine, S.; Pedersen, T. B. Attractive electron–electron interactions within robust local fitting approximations. J. Comput. Chem. 2013, 34, 1486– 1496, DOI: 10.1002/jcc.23284Google Scholar11https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3sXltlOntrg%253D&md5=d51f5c5ca0807fb8dd24244873b2a0d2Attractive electron-electron interactions within robust local fitting approximationsMerlot, Patrick; Kjaergaard, Thomas; Helgaker, Trygve; Lindh, Roland; Aquilante, Francesco; Reine, Simen; Pedersen, Thomas BondoJournal of Computational Chemistry (2013), 34 (17), 1486-1496CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)An anal. of Dunlap's robust fitting approach reveals that the resulting two-electron integral matrix is not manifestly pos. semidefinite when local fitting domains or non-Coulomb fitting metrics are used. We present a highly local approx. method for evaluating four-center two-electron integrals based on the resoln.-of-the-identity (RI) approxn. and apply it to the construction of the Coulomb and exchange contributions to the Fock matrix. In this pair-at. resoln.-of-the-identity (PARI) approach, at.-orbital (AO) products are expanded in auxiliary functions centered on the two atoms assocd. with each product. Numerical tests indicate that in 1% or less of all Hartree-Fock and Kohn-Sham calcns., the indefinite integral matrix causes nonconvergence in the self-consistent-field iterations. In these cases, the two-electron contribution to the total energy becomes neg., meaning that the electronic interaction is effectively attractive, and the total energy is dramatically lower than that obtained with exact integrals. In the vast majority of our test cases, however, the indefiniteness does not interfere with convergence. The total energy accuracy is comparable to that of the std. Coulomb-metric RI method. The speed-up compared with conventional algorithms is similar to the RI method for Coulomb contributions; exchange contributions are accelerated by a factor of up to eight with a triple-zeta quality basis set. A pos. semidefinite integral matrix is recovered within PARI by introducing local auxiliary basis functions spanning the full AO product space, as may be achieved by using Cholesky-decompn. techniques. Local completion, however, slows down the algorithm to a level comparable with or below conventional calcns. © 2013 Wiley Periodicals, Inc.
- 12Manzer, S. F.; Epifanovsky, E.; Head-Gordon, M. Efficient Implementation of the Pair Atomic Resolution of the Identity Approximation for Exact Exchange for Hybrid and Range-Separated Density Functionals. J. Chem. Theory Comput. 2015, 11, 518– 527, DOI: 10.1021/ct5008586Google Scholar12https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXitFehtrrN&md5=bc28289b8fa1caf00b97bdda95c12c88Efficient Implementation of the Pair Atomic Resolution of the Identity Approximation for Exact Exchange for Hybrid and Range-Separated Density FunctionalsManzer, Samuel F.; Epifanovsky, Evgeny; Head-Gordon, MartinJournal of Chemical Theory and Computation (2015), 11 (2), 518-527CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)An efficient new MO basis algorithm is reported implementing the pair at. resoln. of the identity approxn. (PARI) to evaluate the exact exchange contribution (K) to SCF methods, such as hybrid and range-sepd. hybrid d. functionals. The PARI approxn., in which AO basis function pairs are expanded using auxiliary basis functions centered only on their two resp. atoms, was recently investigated by Merlot et al. Our algorithm is significantly faster than quartic scaling RI-K, with an asymptotic exchange speedup for hybrid functionals of (1 + X/N), where N and X are the AO and auxiliary basis dimensions. The asymptotic speedup is 2 + 2X/N for range sepd. hybrids such as CAM-B3LYP, ωB97X-D, and ωB97X-V which include short- and long-range exact exchange. The obsd. speedup for exchange in ωB97X-V for a C68 graphene fragment in the cc-pVTZ basis is 3.4 relative to RI-K. Like conventional RI-K, our method greatly outperforms conventional integral evaluation in large basis sets; a speedup of 19 is obtained in the cc-pVQZ basis on a C54 graphene fragment. Negligible loss of accuracy relative to exact integral evaluation is demonstrated on databases of bonded and nonbonded interactions. We also demonstrate both anal. and numerically that the PARI-K approxn. is variationally stable.
- 13Neese, F.; Wennmohs, F.; Hansen, A.; Becker, U. Efficient, approximate and parallel Hartree–Fock and hybrid DFT calculations. A ‘chain-of-spheres’ algorithm for the Hartree–Fock exchange. Chem. Phys. 2009, 356, 98– 109, DOI: 10.1016/j.chemphys.2008.10.036Google Scholar13https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1MXhs1eru7g%253D&md5=0a7b02cc3eeb4cad044b9b726b8c64ddEfficient, approximate and parallel Hartree-Fock and hybrid DFT calculations. A 'chain-of-spheres' algorithm for the Hartree-Fock exchangeNeese, Frank; Wennmohs, Frank; Hansen, Andreas; Becker, UteChemical Physics (2009), 356 (1-3), 98-109CODEN: CMPHC2; ISSN:0301-0104. (Elsevier B.V.)In this paper, the possibility is explored to speed up Hartree-Fock and hybrid d. functional calcns. by forming the Coulomb and exchange parts of the Fock matrix by different approxns. For the Coulomb part the previously introduced Split-RI-J variant of the well-known d. fitting' approxn. is used. The exchange part is formed by semi-numerical integration techniques that are closely related to Friesner's pioneering pseudo-spectral approach. Our potentially linear scaling realization of this algorithm is called the 'chain-of-spheres exchange' (COSX). A combination of semi-numerical integration and d. fitting is also proposed. Both Split-RI-J and COSX scale very well with the highest angular momentum in the basis sets. It is shown that for extended basis sets speed-ups of up to two orders of magnitude compared to traditional implementations can be obtained in this way. Total energies are reproduced with an av. error of <0.3 kcal/mol as detd. from extended test calcns. with various basis sets on a set of 26 mols. with 20-200 atoms and up to 2000 basis functions. Reaction energies agree to within 0.2 kcal/mol (Hartree-Fock) or 0.05 kcal/mol (hybrid DFT) with the canonical values. The COSX algorithm parallelizes with a speedup of 8.6 obsd. for 10 processes. Min. energy geometries differ by less than 0.3 pm in the bond distances and 0.5° in the bond angles from their canonical values. These developments enable highly efficient and accurate SCF calcns. including nonlocal Hartree-Fock exchange for large mols. In combination with the RI-MP2 method and large basis sets, second-order many body perturbation energies can be obtained for medium sized mols. with unprecedented efficiency. The algorithms are implemented into the ORCA electronic structure system.
- 14Izsák, R.; Neese, F. An overlap fitted chain of spheres exchange method. J. Chem. Phys. 2011, 135, 144105, DOI: 10.1063/1.3646921Google Scholar14https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXht12hsr3P&md5=9c8a746be8f60e71deff75d4029477f9An overlap fitted chain of spheres exchange methodIzsak, Robert; Neese, FrankJournal of Chemical Physics (2011), 135 (14), 144105/1-144105/11CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The "chain of spheres" (COS) algorithm, as part of the RIJCOSX SCF procedure, approximates the exchange term by performing analytic integration with respect to the coordinates of only one of the two electrons, whereas for the remaining coordinates, integration is carried out numerically. In the present work, we attempt to enhance the efficiency of the method by minimizing numerical errors in the COS procedure. The main idea is based on the work of Friesner and consists of finding a fitting matrix, Q, which leads the numerical and anal. evaluated overlap matrixes to coincide. Using Q, the evaluation of exchange integrals can indeed be improved. Improved results and timings are obtained with the present default grid setup for both single point calcns. and geometry optimizations. The fitting procedure results in a redn. of grid sizes necessary for achieving chem. accuracy. We demonstrate this by testing a no. of grids and comparing results to the fully analytic and the earlier COS approxns. This turns out to be favorable for total and reaction energies, for which chem. accuracy can now be reached with a corresponding ∼30% speedup over the original RIJCOSX procedure for single point energies. Results are slightly less favorable for the accuracy of geometry optimizations, but the procedure is still shown to yield geometries with errors well below the method inherent errors of the employed theor. framework. (c) 2011 American Institute of Physics.
- 15Laqua, H.; Thompson, T. H.; Kussmann, J.; Ochsenfeld, C. Highly Efficient, Linear-Scaling Seminumerical Exact-Exchange Method for Graphic Processing Units. J. Chem. Theory Comput. 2020, 16, 1456– 1468, DOI: 10.1021/acs.jctc.9b00860Google Scholar15https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXivVejs7Y%253D&md5=606be25e20a70cc545ff71ebabf28152Highly Efficient, Linear-Scaling Seminumerical Exact-Exchange Method for Graphic Processing UnitsLaqua, Henryk; Thompson, Travis H.; Kussmann, Joerg; Ochsenfeld, ChristianJournal of Chemical Theory and Computation (2020), 16 (3), 1456-1468CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)We present a highly efficient and asymptotically linear-scaling graphic processing unit accelerated seminumerical exact-exchange method (sn-LinK). We go beyond our previous central processing unit-based method (H. Laqua et al., 2018) by employing our recently developed integral bounds (T.H. Thomson and C. Ochsenfeld, 2019) and high-accuracy numerical integration grid (H. Laqua et al., 2018). The accuracy is assessed for several established test sets, providing errors significantly below 1mEh for the smallest grid. Moreover, a comprehensive performance anal. for large mols. between 62 and 1347 atoms is provided, revealing the outstanding performance of our method, in particular, for large basis sets such as the polarized quadruple-zeta level with diffuse functions.
- 16Holzer, C. An improved seminumerical Coulomb and exchange algorithm for properties and excited states in modern density functional theory. J. Chem. Phys. 2020, 153, 184115, DOI: 10.1063/5.0022755Google Scholar16https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXitlelur7N&md5=cad44a4ffbef63e25cf320bbcd316688An improved seminumerical Coulomb and exchange algorithm for properties and excited states in modern density functional theoryHolzer, ChristofJournal of Chemical Physics (2020), 153 (18), 184115CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)A seminumerical algorithm capable of performing large-scale (time-dependent) d. functional theory (TD-DFT) calcns. to ext. excitation energies and other ground-state and excited-state properties is outlined. The algorithm uses seminumerical integral techniques for evaluating Coulomb and exchange parts for a set of d. matrixes as occurring in std. TD-DFT or similar methods for the evaluation of vibrational frequencies. A suitable optimized de-aliasing procedure is introduced. The latter does not depend on further auxiliary quantities and retains the symmetry of a given d. matrix. The algorithm is self-contained and applicable to any orbital basis set available without the need for further auxiliary basis sets or optimized de-aliasing grids. Relativistic two-component excited-state TD-DFT calcns. are reported for the first time using the developed seminumerical algorithm for std. and local hybrid d. functional approxns. Errors are compared with the widely used "resoln. of the identity" (RI) approxns. for Coulomb (RI-J) and exchange integrals (RI-K). The fully seminumerical algorithm does not exhibit an enlarged error for std. DFT functionals compared to the RI approxn. For the more involved local hybrid functionals and within strong external fields, accuracy is even considerably improved. (c) 2020 American Institute of Physics.
- 17Ohno, K. Some remarks on the Pariser-Parr-Pople method. Theor. Chim. Acta 1964, 2, 219– 227, DOI: 10.1007/BF00528281Google Scholar17https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaF2cXkt12htLk%253D&md5=b652ff406687369fbbdb1d53b19870e4The Pariser-Parr-Pople methodOhno, KimioTheoretica Chimica Acta (1964), 2 (3), 219-27CODEN: TCHAAM; ISSN:0040-5744.Basic assumptions which characterize the Pariser-Parr (CA 49, 10725f)-Pople (CA 48, 8639g) method of computing mol. electronic wave functions are examd. crit. By restricted variational calcn. of the valence state of C and N atoms and ions, it is demonstrated that the usual methods of evaluation of center Coulomb integrals and at. core energies are rather good. A semitheoretical means of estg. the core resonance integral is proposed and shown to give fair agreement with the empirical values for C-C, O-O, C-N, and C-O bonds.
- 18Klopman, G. A Semiempirical Treatment of molecular Structures. II. Molecular Terms and Application to diatomic Molecules. J. Am. Chem. Soc. 1964, 86, 4550– 4557, DOI: 10.1021/ja01075a008Google Scholar18https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaF2cXkvVygsr4%253D&md5=be4431018f6ae6937415a7bdf8370bc6A semiempirical treatment of molecular structures. II. Molecular terms and application to diatomic moleculesKlopman, G.Journal of the American Chemical Society (1964), 86 (21), 4550-7CODEN: JACSAT; ISSN:0002-7863.cf. CA 60, 12679e. A self-consistent semiempirical method which is designed for the calcn. of heats of formation and charge distribution of nonconjugated mols. is outlined. The method is based on an antisymmetrized product of mol. orbitals, simplified in such a way as to make the values of all involved integrals directly available from at. spectra (loc. cit.) and mol. bond distances. In a preliminary study, this method was used to calc. satisfactory values of bond energies and reasonable values of charge distributions in 80 diat. mols. (σ-bonded).
- 19Nishimoto, K.; Mataga, N. Z. Electronic Structure and Spectra of Some Nitrogen Heterocycles. Z. Phys. Chem. 1957, 12, 335, DOI: 10.1524/zpch.1957.12.5_6.335Google Scholar19https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaG2sXps1ymsg%253D%253D&md5=95bf1fdd3345cb66b0f8af10bf781f0bElectronic structure and spectra of some nitrogen heretocyclesNishimoto, Kitisuke; Mataga, NoboruZeitschrift fuer Physikalische Chemie (Muenchen, Germany) (1957), 12 (), 335-8CODEN: ZPCFAX; ISSN:0044-3336.The energies of lower excited states of pyridine, pyrazine, and sym-triazine and the oscillator strengths of the transitions to those excited states have been calcd. The results are listed and discussed.
- 20Bannwarth, C.; Grimme, S. A simplified time-dependent density functional theory approach for electronic ultraviolet and circular dichroism spectra of very large molecules. Comput. Theor. Chem. 2014, 1040–1041, 45– 53, DOI: 10.1016/j.comptc.2014.02.023Google Scholar20https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXksFylt7o%253D&md5=a51265489cf75ee7284bf9cd31ece615A simplified time-dependent density functional theory approach for electronic ultraviolet and circular dichroism spectra of very large moleculesBannwarth, Christoph; Grimme, StefanComputational & Theoretical Chemistry (2014), 1040-1041 (), 45-53CODEN: CTCOA5; ISSN:2210-271X. (Elsevier B.V.)We present a simplified time-dependent d. functional theory approach (sTD-DFT) that allows fast computation of electronic UV or CD spectra of mols. with 500-1000 atoms. The matrix elements are treated in the same way as in the recently proposed simplified Tamm-Dancoff approach but instead of applying the Tamm-Dancoff approxn., the std. linear-response d. functional theory problem is solved. Compared to sTDA, the method leads to an increase in computation time (typically a factor of 2-5 compared to the corresponding sTDA) which is justified since the resulting transition dipole moments are in general of higher quality. This becomes important if spectral intensities (e.g. single-photon oscillator and rotatory transition strengths) are of interest. Comparison of UV and CD spectra obtained from sTD-DFT and sTDA for some typical systems employing std. hybrid functionals shows that both yield very similar excitation energies but the advantage of using the former approach for transition moments. In order to show the applicability of sTD-DFT to systems which are far beyond the scope of conventional TD-DFT, we present the CD spectrum of a substituted, chiral fullerene over a range of almost 1200 excited states. We propose this method as a more reliable alternative for the prediction esp. of the more challenging CD spectra.
- 21Grimme, S. A simplified Tamm-Dancoff density functional approach for the electronic excitation spectra of very large molecules. J. Chem. Phys. 2013, 138, 244104, DOI: 10.1063/1.4811331Google Scholar21https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3sXpvFCit7s%253D&md5=1c591b50b78ab4df676b394c5754f9d3A simplified Tamm-Dancoff density functional approach for the electronic excitation spectra of very large moleculesGrimme, StefanJournal of Chemical Physics (2013), 138 (24), 244104/1-244104/14CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Two approxns. in the Tamm-Dancoff d. functional theory approach (TDA-DFT) to electronically excited states are proposed which allow routine computations for electronic UV- or CD spectra of mols. with 500-1000 atoms. Speed-ups compared to conventional time-dependent DFT (TD-DFT) treatments of about two to three orders of magnitude in the excited state part at only minor loss of accuracy are obtained. The method termed sTDA ("s" for simplified) employs atom-centered Loewdin-monopole based two-electron repulsion integrals with the asymptotically correct 1/R behavior and perturbative single excitation configuration selection. It is formulated generally for any std. global hybrid d. functional with given Fock-exchange mixing parameter ax. The method performs well for two std. benchmark sets of vertical singlet-singlet excitations for values of ax in the range 0.2-0.6. The mean abs. deviations from ref. data are only 0.2-0.3 eV and similar to those from std. TD-DFT. In three cases (two dyes and one polypeptide), good mutual agreement between the electronic spectra (up to 10-11 eV excitation energy) from the sTDA method and those from TD(A)-DFT is obtained. The computed UV- and CD-spectra of a few typical systems (e.g., C60, two transition metal complexes, [7]helicene, polyalanine, a supramol. aggregate with 483 atoms and about 7000 basis functions) compare well with corresponding exptl. data. The method is proposed together with medium-sized double- or triple-zeta type at.-orbital basis sets as a quantum chem. tool to investigate the spectra of huge mol. systems at a reliable DFT level. (c) 2013 American Institute of Physics.
- 22Grimme, S.; Bannwarth, C. Ultra-fast computation of electronic spectra for large systems by tight-binding based simplified Tamm-Dancoff approximation (sTDA-xTB). J. Chem. Phys. 2016, 145, 054103, DOI: 10.1063/1.4959605Google Scholar22https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28Xht1Klt7nK&md5=2ff48457d9d88c019084d3fdf6fee0e9Ultra-fast computation of electronic spectra for large systems by tight-binding based simplified Tamm-Dancoff approximation (sTDA-xTB)Grimme, Stefan; Bannwarth, ChristophJournal of Chemical Physics (2016), 145 (5), 054103/1-054103/20CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The computational bottleneck of the extremely fast simplified Tamm-Dancoff approximated (sTDA) time-dependent d. functional theory procedure [S. Grimme, J. Chem. Phys. 138, 244104 (2013)] for the computation of electronic spectra for large systems is the detn. of the ground state Kohn-Sham orbitals and eigenvalues. This limits such treatments to single structures with a few hundred atoms and hence, e.g., sampling along mol. dynamics trajectories for flexible systems or the calcn. of chromophore aggregates is often not possible. The aim of this work is to solve this problem by a specifically designed semiempirical tight binding (TB) procedure similar to the well established self-consistent-charge d. functional TB scheme. The new special purpose method provides orbitals and orbital energies of hybrid d. functional character for a subsequent and basically unmodified sTDA procedure. Compared to many previous semiempirical excited state methods, an advantage of the ansatz is that a general eigenvalue problem in a non-orthogonal, extended AO basis is solved and therefore correct occupied/virtual orbital energy splittings as well as Rydberg levels are obtained. A key idea for the success of the new model is that the detn. of at. charges (describing an effective electron-electron interaction) and the one-particle spectrum is decoupled and treated by two differently parametrized Hamiltonians/basis sets. The three-diagonalization-step composite procedure can routinely compute broad range electronic spectra (0-8 eV) within minutes of computation time for systems composed of 500-1000 atoms with an accuracy typical of std. time-dependent d. functional theory (0.3-0.5 eV av. error). An easily extendable parametrization based on coupled-cluster and d. functional computed ref. data for the elements H-Zn including transition metals is described. The accuracy of the method termed sTDA-xTB is first benchmarked for vertical excitation energies of open- and closed-shell systems in comparison to other semiempirical methods and applied to exemplary problems in electronic spectroscopy. As side products of the development, a robust and efficient valence electron TB method for the accurate detn. of at. charges as well as a more accurate calcn. scheme of dipole rotatory strengths within the Tamm-Dancoff approxn. is proposed. (c) 2016 American Institute of Physics.
- 23de Wergifosse, M.; Bannwarth, C.; Grimme, S. A Simplified Spin-Flip Time-Dependent Density Functional Theory Approach for the Electronic Excitation Spectra of Very Large Diradicals. J. Phys. Chem. A 2019, 123, 5815– 5825, DOI: 10.1021/acs.jpca.9b03176Google Scholar23https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXhtFyjsbzE&md5=e9769d3847d1c74864f2d04259fd5a9dA Simplified Spin-Flip Time-Dependent Density Functional Theory Approach for the Electronic Excitation Spectra of Very Large Diradicalsde Wergifosse, Marc; Bannwarth, Christoph; Grimme, StefanJournal of Physical Chemistry A (2019), 123 (27), 5815-5825CODEN: JPCAFH; ISSN:1089-5639. (American Chemical Society)Experimentalists working with diradicals are often facing the question of what kind of species among singlet or triplet diradicals or closed-shell mols. are obsd. To treat large diradicals with a high d. of electronic states, we propose a simplified version of the spin-flip time-dependent d. functional theory (SF-TD-DFT) method for a fast computation of their state energies and absorption spectra with an accuracy similar to the nonsimplified scheme. An ultrafast tight-binding variant called SF-sTD-DFT-xTB is also developed to treat even larger systems. For a benchmark set of nine diradicals, good agreement between simplified and conventional SF excitation energies for std. functionals is found. This shows that the proposed parameterization is robust for a wide range of Fock exchange mixing values. With the asymptotically correct response integrals used in SF-sTD-DFT and a correction factor of √(2) for the transition moments, the SF-sTD-DFT/B5050LYP/cc-pVDZ method even outperforms the nonsimplified scheme at drastically reduced computational effort when comparing to the exptl. absorption spectra for this set of diradicals. To showcase the actual performance of the method, absorption spectra of two μ-hydroxo-bridged dimers of corrole tape Ga(III) complex derivs. were computed and compared to the expt., providing good qual. agreement. Finally, a comparison with the high-spin triplet spectrum of a perylene bisimide biradical and the one detd. at the SF-sTD-DFT level showed that at room temp., mostly triplet diradicals are present in soln.
- 24De Wergifosse, M.; Grimme, S. Nonlinear-response properties in a simplified time-dependent density functional theory (sTD-DFT) framework: Evaluation of excited-state absorption spectra. J. Chem. Phys. 2019, 150, 094112, DOI: 10.1063/1.5080199Google Scholar24https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXktlersLY%253D&md5=8f0cea7b5ce0b88cc73e98a90cce4240Nonlinear-response properties in a simplified time-dependent density functional theory (sTD-DFT) framework: Evaluation of excited-state absorption spectrade Wergifosse, Marc; Grimme, StefanJournal of Chemical Physics (2019), 150 (9), 094112/1-094112/9CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The energy conversion efficiency of org. solar cells seems crucial for a clean future. The design of new light-harvesting devices needs an in-depth understanding of their optical properties, including the excited-state absorption (ESA). In biol., the optical characterization of photochem./phys. processes happening in photosynthetic pigments and proteins can be difficult to interpret due to their structural complexities. Exptl., an ultrafast transient absorption expt. can probe the excited state interaction with light. Quantum chem. could play an important role to model the transient absorption spectrum of excited states. However, systems that need to be investigated can be way too large for existent software implementations. In this contribution, we present the first sTDA/sTD-DFT (simplified time-dependent d. functional theory with and without Tamm Dancoff approxn.) implementation to evaluate the ESA of mols. The ultrafast ESA evaluation presents a negligible extra cost with respect to sTDA/sTD-DFT original schemes for std. ground state absorption. The sTD-DFT method shows ability to assign ESA spectra to the correct excited state. We showed that in the literature, wrong assignments were proposed as for the L34/L44 mixt. and N-methylfulleropyrrolidine. In addn., sTDA/sTD-DFT-xTB tight-binding variants are also available, allowing the evaluation of ESA for systems of a few thousands of atoms, e.g., the spectrum of the photoactive yellow protein composed of 1931 atoms. (c) 2019 American Institute of Physics.
- 25De Wergifosse, M.; Grimme, S. Nonlinear-response properties in a simplified time-dependent density functional theory (sTD-DFT) framework: Evaluation of the first hyperpolarizability. J. Chem. Phys. 2018, 149, 024108, DOI: 10.1063/1.5037665Google Scholar25https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXhtlWqurzI&md5=8a2bb42755257c0761582ab77e2efdc5Nonlinear-response properties in a simplified time-dependent density functional theory (sTD-DFT) framework: Evaluation of the first hyperpolarizabilityde Wergifosse, Marc; Grimme, StefanJournal of Chemical Physics (2018), 149 (2), 024108/1-024108/13CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Recent developments in nonlinear imaging microscopy show the need to implement new theor. tools, which are able to characterize nonlinear optical properties in an efficient way. For second-harmonic imaging microscopy (SHIM), quantum chem. could play an important role to design new exogenous dyes with enhanced first hyperpolarizabilities or to characterize the response origin in large endogenous biol. systems. Such methods should be able to screen a large no. of compds. while reproducing their trends and to treat large systems in reasonable computation times. To fulfill these requirements, we present a new simplified time-dependent d. functional theory (sTD-DFT) implementation to evaluate the first hyperpolarizability where the Coulomb and exchange integrals are approximated by short-range damped Coulomb interactions of transition d. monopoles. For an ultra-fast computation of the first hyperpolarizability, a tight-binding version (sTD-DFT-xTB) is also proposed. In our implementation, a sTD-DFT calcn. is more than 600 time faster with respect to a regular TD-DFT treatment, while the xTB version speeds up the entire calcn. further by at least two orders of magnitude. We challenge our implementation on three test cases: typical push-pull π-conjugated compds., fluorescent proteins, and a collagen model, which were selected to model requirements for SHIM applications. (c) 2018 American Institute of Physics.
- 26Bannwarth, C.; Seibert, J.; Grimme, S. Electronic Circular Dichroism of [16]Helicene With Simplified TD-DFT: Beyond the Single Structure Approach. Chirality 2016, 28, 365– 369, DOI: 10.1002/chir.22594Google Scholar26https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XlvV2mtrc%253D&md5=40b2daee5ca75233db6ffaa58bccbb58Electronic Circular Dichroism of [16]Helicene With Simplified TD-DFT: Beyond the Single Structure ApproachBannwarth, Christoph; Seibert, Jakob; Grimme, StefanChirality (2016), 28 (5), 365-369CODEN: CHRLEP; ISSN:0899-0042. (Wiley-Liss, Inc.)The electronic CD (ECD) spectrum of the recently synthesized [16]helicene and a deriv. comprising two triisopropylsilyloxy protection groups was computed by means of the very efficient simplified time-dependent d. functional theory (sTD-DFT) approach. Different from many previous ECD studies of helicenes, nonequil. structure effects were accounted for by computing ECD spectra on "snapshots" obtained from a mol. dynamics (MD) simulation including solvent mols. The trajectories are based on a mol. specific classical potential as obtained from the recently developed quantum chem. derived force field (QMDFF) scheme. The reduced computational cost in the MD simulation due to the use of the QMDFF (compared to ab-initio MD) as well as the sTD-DFT approach make realistic spectral simulations feasible for these compds. that comprise more than 100 atoms. While the ECD spectra of [16]helicene and its deriv. computed vertically on the resp. gas phase, equil. geometries show noticeable differences, these are "washed" out when nonequil. structures are taken into account. The computed spectra with two recommended d. functionals (ωB97X and BHLYP) and extended basis sets compare very well with the exptl. one. In addn. we provide an est. for the missing abs. intensities of the latter. The approach presented here could also be used in future studies to capture nonequil. effects, but also to systematically av. ECD spectra over different conformations in more flexible mols. Chirality 00:000-000, 2016. 2016 Wiley Periodicals, Inc.
- 27Risthaus, T.; Hansen, A.; Grimme, S. Excited states using the simplified Tamm-Dancoff-Approach for range-separated hybrid density functionals: development and application. Phys. Chem. Chem. Phys. 2014, 16, 14408– 14419, DOI: 10.1039/C3CP54517BGoogle Scholar27https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXhtVKgs7bM&md5=978384902d7b53a6ef2eefa57ef407feExcited states using the simplified Tamm-Dancoff-Approach for range-separated hybrid density functionals: development and applicationRisthaus, Tobias; Hansen, Andreas; Grimme, StefanPhysical Chemistry Chemical Physics (2014), 16 (28), 14408-14419CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)The recently introduced sTDA methodol. [S. Grimme, J. Chem. Phys., 2013, 138, 244104] to compute excitation spectra of huge mol. systems is extended to range-sepd. hybrid (RSH) d. functionals. The three empirical parameters of the method which describe a screened two-electron interaction are obtained for some common RSH functionals (ωB97 family, CAM-B3LYP, LC-BLYP) from a fit to theor. SCS-CC2 ref. vertical excitation energies for a set of small to medium-sized chromophores. The method is cross-validated on a set of inter- and intramol. charge transfer states and a set composed of typical valence transitions. Overall small deviations from ref. data of only about 0.2-0.4 eV are found with best performance for CAM-B3LYP and ωB97X-D3. To demonstrate versatility and robustness of the new methodol., applications (the UV/Vis spectrum of the pyridine polymer and the ECD spectrum of (P)-[11]helicene) and frequently used charge transfer examples are discussed. In one case, 11 000 + excited electronic states of a system contg. 330 atoms were calcd. We show that the asymptotically correct sTDA-RSH combination yields results often superior to those based on global hybrids and that it opens up new possibilities for the computation of excited states in materials science and bio-mol. systems.
- 28Elstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, T.; Suhai, S.; Seifert, G. Self-consistent-charge density-functional tight-binding method for simulations of complex materials properties. Phys. Rev. B 1998, 58, 7260– 7268, DOI: 10.1103/PhysRevB.58.7260Google Scholar28https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK1cXmtVamsb0%253D&md5=e3b3143c475c70203badf7db6ef668c3Self-consistent-charge density-functional tight-binding method for simulations of complex materials propertiesElstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, Th.; Suhai, S.; Seifert, G.Physical Review B: Condensed Matter and Materials Physics (1998), 58 (11), 7260-7268CODEN: PRBMDO; ISSN:0163-1829. (American Physical Society)We outline details about an extension of the tight-binding (TB) approach to improve total energies, forces, and transferability. The method is based on a second-order expansion of the Kohn-Sham total energy in d.-functional theory (DFT) with respect to charge-d. fluctuations. The zeroth-order approach is equiv. to a common std. non-self-consistent (TB) scheme, while at second-order a transparent, parameter-free, and readily calculable expression for generalized Hamiltonian matrix elements may be derived. These are modified by a self-consistent redistribution of Mulliken charges (SCC). Besides the usual "band structure" and short-range repulsive terms the final approx. Kohn-Sham energy addnl. includes a Coulomb interaction between charge fluctuations. At large distances this accounts for long-range electrostatic forces between two point charges and approx. includes self-interaction contributions of a given atom if the charges are located at one and the same atom. We apply the new SCC scheme to problems where deficiencies within the non-SCC std. TB approach become obvious. We thus considerably improve transferability.
- 29Seibert, J.; Pisarek, J.; Schmitz, S.; Bannwarth, C.; Grimme, S. Extension of the element parameter set for ultra-fast excitation spectra calculation (sTDA-xTB). Mol. Phys. 2019, 117, 1104– 1116, DOI: 10.1080/00268976.2018.1510141Google Scholar29https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXhsFCjs7bK&md5=57617dc9f5fcc3316f1edabcb4581b83Extension of the element parameter set for ultra-fast excitation spectra calculation (sTDA-xTB)Seibert, Jakob; Pisarek, Jana; Schmitz, Sarah; Bannwarth, Christoph; Grimme, StefanMolecular Physics (2019), 117 (9-12), 1104-1116CODEN: MOPHAM; ISSN:0026-8976. (Taylor & Francis Ltd.)The extension of the parameter set for an ultra-fast electronic excitation spectra calcn. is presented. The semiempirical theory based on a tight-binding approach, called extended tight-binding (xTB) in combination with the simplified Tamm-Dancoff approxn. (sTDA) shows remarkable accuracy at very low computational cost for the calcn. of vertical excitation energies of mols. It enables the possibility for computing even large systems up to thousands of atoms or sampling along mol. dynamic (MD) trajectories. The original publication of the sTDA-xTB method included parameters for the most important elements (H-Zn,Br,I). In this work, element parameters for 4d and 5d metals, and the missing ones in 4p, 5p and 6p element blocks are presented and analyzed for their quality. Comparisons to theory and expt. show that sTDA-xTB provides similar good results as for the elements in the original publication with an av. deviation of excitation energies of 0.3-0.5 eV.
- 30Cho, Y.; Bintrim, S. J.; Berkelbach, T. C. A simplified GW/BSE approach for charged and neutral excitation energies of large molecules and nanomaterials ; 2021, arXiv:2109.04421. DOI: 10.48550/arXiv.2109.04421 .Google ScholarThere is no corresponding record for this reference.
- 31Rüger, R.; van Lenthe, E.; Heine, T.; Visscher, L. Tight-binding approximations to time-dependent density functional theory ─ A fast approach for the calculation of electronically excited states. J. Chem. Phys. 2016, 144, 184103, DOI: 10.1063/1.4948647Google Scholar31https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XnslWru7g%253D&md5=0efdefe5a46f1be35cbf902b0cfc3726Tight-binding approximations to time-dependent density functional theory - A fast approach for the calculation of electronically excited statesRueger, Robert; van Lenthe, Erik; Heine, Thomas; Visscher, LucasJournal of Chemical Physics (2016), 144 (18), 184103/1-184103/12CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We propose a new method of calcg. electronically excited states that combines a d. functional theory based ground state calcn. with a linear response treatment that employs approxns. used in the time-dependent d. functional based tight binding (TD-DFTB) approach. The new method termed time-dependent d. functional theory TD-DFT+TB does not rely on the DFTB parametrization and is therefore applicable to systems involving all combinations of elements. We show that the new method yields UV/Vis absorption spectra that are in excellent agreement with computationally much more expensive TD-DFT calcns. Errors in vertical excitation energies are reduced by a factor of two compared to TD-DFTB. (c) 2016 American Institute of Physics.
- 32Asadi-Aghbolaghi, N.; Pototschnig, J.; Jamshidi, Z.; Visscher, L. Effects of ligands on (de-)enhancement of plasmonic excitations of silver, gold and bimetallic nanoclusters: TD-DFT+TB calculations. Phys. Chem. Chem. Phys. 2021, 23, 17929– 17938, DOI: 10.1039/D1CP03220HGoogle Scholar32https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3MXhslGms7%252FN&md5=00c70d22cdec2f0da8c7a1a1de065a96Effects of ligands on (de-)enhancement of plasmonic excitations of silver, gold and bimetallic nanoclusters: TD-DFT+TB calculationsAsadi-Aghbolaghi, Narges; Pototschnig, Johann; Jamshidi, Zahra; Visscher, LucasPhysical Chemistry Chemical Physics (2021), 23 (33), 17929-17938CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)Metal nanoclusters can be synthesized in various sizes and shapes and are typically protected with ligands to stabilize them. These ligands can also be used to tune the plasmonic properties of the clusters as the absorption spectrum of a protected cluster can be significantly altered compared to the bare cluster. In this paper, we computationally investigate the influence of thiolate ligands on the plasmonic intensity for silver, gold and alloy clusters. Using time-dependent d. functional theory with tight-binding approxns., TD-DFT+TB, we show that this level of theory can reproduce the broad exptl. spectra of Au144(SR)60 and Ag53Au91(SR)60 (R = CH3) compds. with satisfactory agreement. As TD-DFT+TB does not depend on atom-type parameters we were able to apply this approach on large ligand-protected clusters with various compns. With these calcns. we predict that the effect of ligands on the absorption can be a quenching as well as an enhancement. We furthermore show that it is possible to unambiguously identify the plasmonic peaks by the scaled Coulomb kernel technique and explain the influence of ligands on the intensity (de-)enhancement by analyzing the plasmonic excitations in terms of the dominant orbital contributions.
- 33Van Caillie, C.; Amos, R. D. Geometric derivatives of excitation energies using SCF and DFT. Chem. Phys. Lett. 1999, 308, 249– 255, DOI: 10.1016/S0009-2614(99)00646-6Google Scholar33https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK1MXksFWku70%253D&md5=b0a42882739b7c4d495c632cc8205f5aGeometric derivatives of excitation energies using SCF and DFTVan Caillie, Carole; Amos, Roger D.Chemical Physics Letters (1999), 308 (3,4), 249-255CODEN: CHPLBC; ISSN:0009-2614. (Elsevier Science B.V.)There is increasing interest in using the methods of time-dependent d. functional theory to calc. electronic excitation energies. We have implemented an analytic gradient method to find the geometric derivs. of the excitation energies. When added to the gradient for the ground state, this yields the excited-state energy derivs. This enables the efficient generation and searching of excited-state potential energy surfaces to obtain excited-state geometries and other properties. The initial implementation is for SCF methods and for the local d. approxn. Some examples of excited-state geometry optimizations are given.
- 34Van Caillie, C.; Amos, R. D. Geometric derivatives of density functional theory excitation energies using gradient-corrected functionals. Chem. Phys. Lett. 2000, 317, 159– 164, DOI: 10.1016/S0009-2614(99)01346-9Google Scholar34https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3cXotlyksw%253D%253D&md5=fb54eacfa34d1288afe2075376db3e3cGeometric derivatives of density functional theory excitation energies using gradient-corrected functionalsVan Caillie, C.; Amos, R. D.Chemical Physics Letters (2000), 317 (1,2), 159-164CODEN: CHPLBC; ISSN:0009-2614. (Elsevier Science B.V.)D. functional theory (DFT) is having increasing success in predicting excitation energies using the methods of time-dependent DFT. As a result, it should be possible to generate potential energy surfaces for excited states by adding the excitation energy, as a function of geometry, to the ground-state energy. It is easier to find stationary points such as min. and transition states if the gradient of the energy is known. The present Letter extends earlier work on the gradients on excited-state surfaces using SCF and LDA (local d. approxn.) methods, to use gradient-cor. and hybrid functionals. Some examples of geometry optimizations are given.
- 35Furche, F.; Ahlrichs, R. Adiabatic time-dependent density functional methods for excited state properties. J. Chem. Phys. 2002, 117, 7433– 7447, DOI: 10.1063/1.1508368Google Scholar35https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD38XnvVWrurY%253D&md5=061f0264a5993d772854715400d3d189Adiabatic time-dependent density functional methods for excited state propertiesFurche, Filipp; Ahlrichs, ReinhartJournal of Chemical Physics (2002), 117 (16), 7433-7447CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)This work presents theory, implementation, and validation of excited state properties obtained from time-dependent d. functional theory (TDDFT). Based on a fully variational expression for the excited state energy, a compact derivation of first order properties is given. We report an implementation of analytic excited state gradients and charge moments for local, gradient cor., and hybrid functionals, as well as for the CI singles (CIS) and time-dependent Hartree-Fock (TDHF) methods. By exploiting analogies to ground state energy and gradient calcns., efficient techniques can be transferred to excited state methods. Benchmark results demonstrate that, for low-lying excited states, geometry optimizations are not substantially more expensive than for the ground state, independent of the mol. size. We assess the quality of calcd. adiabatic excitation energies, structures, dipole moments, and vibrational frequencies by comparison with accurate exptl. data for a variety of excited states and mols. Similar trends are obsd. for adiabatic excitation energies as for vertical ones. TDDFT is more robust than CIS and TDHF, in particular, for geometries differing significantly from the ground state min. The TDDFT excited state structures, dipole moments, and vibrational frequencies are of a remarkably high quality, which is comparable to that obtained in ground state d. functional calcns. Thus, yielding considerably more accurate results at similar computational cost, TDDFT rivals CIS as a std. method for calcg. excited state properties in larger mols.
- 36Furche, F.; Ahlrichs, R. Erratum: ”Time-dependent density functional methods for excited state properties” [J. Chem. Phys. 117, 7433 (2002)]. J. Chem. Phys. 2004, 121, 12772– 12773, DOI: 10.1063/1.1824903Google Scholar36https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2cXhtFSnsrnE&md5=c2047481400c80542918b0de9189c0baAdiabatic time-dependent density functional methods for excited state properties. [Erratum to document cited in CA138:044999]Furche, Filipp; Ahlrichs, ReinhartJournal of Chemical Physics (2004), 121 (24), 12772-12773CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)A factor 1/2 is missing in the expression for the diagonal elements of the energy weighted difference d. matrix W throughout the paper. The cor. versions of Equations (24), (A4), (A10), and (A13) and Tables V, VI, and VIII are given. The implementation and results presented in the paper were based on the correct expressions. Some of the calcd. spectroscopic consts. of the 1 1.sum.u- state of N2, of the 1 1.sum.u+ state of Mg2, and of the 2 1.sum.+ state of CuH in Tables V, VI, and VII are incorrect and should be replaced by the values in the cor. equations. The conclusions are not affected by these changes.
- 37Hutter, J. Excited state nuclear forces from the Tamm-Dancoff approximation to time-dependent density functional theory within the plane wave basis set framework. J. Chem. Phys. 2003, 118, 3928– 3934, DOI: 10.1063/1.1540109Google Scholar37https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3sXht1CltLc%253D&md5=4099596d811bae697319ed324e79ce8cExcited state nuclear forces from the Tamm-Dancoff approximation to time-dependent density functional theory within the plane wave basis set frameworkHutter, JurgJournal of Chemical Physics (2003), 118 (9), 3928-3934CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)An efficient formulation of time-dependent linear response d. functional theory for the use within the plane wave basis set framework is presented. The method avoids the transformation of the Kohn-Sham matrix into the canonical basis and refs. virtual orbitals only through a projection operator. Using a Lagrangian formulation nuclear derivs. of excited state energies within the Tamm-Dancoff approxn. are derived. The algorithms were implemented into a pseudo potential/plane wave code and applied to the calcn. of adiabatic excitation energies, optimized geometries and vibrational frequencies of three low lying states of formaldehyde. An overall good agreement with other time-dependent d. functional calcns., multireference CI calcns. and exptl. data was found.
- 38Petrenko, T.; Kossmann, S.; Neese, F. Efficient time-dependent density functional theory approximations for hybrid density functionals: Analytical gradients and parallelization. J. Chem. Phys. 2011, 134, 054116, DOI: 10.1063/1.3533441Google Scholar38https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXhsFSqtrs%253D&md5=1eae10d380aea0e3920f1870fc2a181fEfficient time-dependent density functional theory approximations for hybrid density functionals: Analytical gradients and parallelizationPetrenko, Taras; Kossmann, Simone; Neese, FrankJournal of Chemical Physics (2011), 134 (5), 054116/1-054116/14CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We present the implementation of efficient approxns. to time-dependent d. functional theory (TDDFT) within the Tamm-Dancoff approxn. (TDA) for hybrid d. functionals. For the calcn. of the TDDFT/TDA excitation energies and anal. gradients, we combine the resoln. of identity (RI-J) algorithm for the computation of the Coulomb terms and the recently introduced "chain of spheres exchange" (COSX) algorithm for the calcn. of the exchange terms. It is shown that for extended basis sets, the RIJCOSX approxn. leads to speedups of up to 2 orders of magnitude compared to traditional methods, as demonstrated for hydrocarbon chains. The accuracy of the adiabatic transition energies, excited state structures, and vibrational frequencies is assessed on a set of 27 excited states for 25 mols. with the CI singles and hybrid TDDFT/TDA methods using various basis sets. Compared to the canonical values, the typical error in transition energies is of the order of 0.01 eV. Similar to the ground-state results, excited state equil. geometries differ by less than 0.3 pm in the bond distances and 0.5° in the bond angles from the canonical values. The typical error in the calcd. excited state normal coordinate displacements is of the order of 0.01, and relative error in the calcd. excited state vibrational frequencies is less than 1%. The errors introduced by the RIJCOSX approxn. are, thus, insignificant compared to the errors related to the approx. nature of the TDDFT methods and basis set truncation. For TDDFT/TDA energy and gradient calcns. on Ag-TB2-helicate (156 atoms, 2732 basis functions), it is demonstrated that the COSX algorithm parallelizes almost perfectly (speedup ∼26-29 for 30 processors). The exchange-correlation terms also parallelize well (speedup ∼27-29 for 30 processors). The soln. of the Z-vector equations shows a speedup of ∼24 on 30 processors. The parallelization efficiency for the Coulomb terms can be somewhat smaller (speedup ∼15-25 for 30 processors), but their contribution to the total calcn. time is small. Thus, the parallel program completes a Becke3-Lee-Yang-Parr energy and gradient calcn. on the Ag-TB2-helicate in less than 4 h on 30 processors. We also present the necessary extension of the Lagrangian formalism, which enables the calcn. of the TDDFT excited state properties in the frozen-core approxn. The algorithms described in this work are implemented into the ORCA electronic structure system. (c) 2011 American Institute of Physics.
- 39Grotjahn, R.; Furche, F.; Kaupp, M. Development and Implementation of Excited-State Gradients for Local Hybrid Functionals. J. Chem. Theory Comput. 2019, 15, 5508– 5522, DOI: 10.1021/acs.jctc.9b00659Google Scholar39https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXhs12jsbbE&md5=5943c43f7445f4cf7546917346d65f36Development and Implementation of Excited-State Gradients for Local Hybrid FunctionalsGrotjahn, Robin; Furche, Filipp; Kaupp, MartinJournal of Chemical Theory and Computation (2019), 15 (10), 5508-5522CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)Local hybrid functionals are a relatively recent class of exchange-correlation functionals that use a real-space dependent admixt. of exact exchange. Here, we present the first implementation of time-dependent d. functional theory excited-state gradients for these functionals. Based on the ansatz of a fully variational auxiliary Lagrangian of the excitation energy, the working equations for the case of a local hybrid functional are deduced. For the implementation, we derive the third-order functional derivs. used in the hyper-kernel and kernel-gradients following a seminumerical integration scheme. The first assessment for a test set of small mols. reveals competitive performance for excited-state structural parameters with typical mean abs. errors (MAEs) of 1.2 pm (PBE0: 1.4 pm) for bond lengths as well as for vibrational frequencies with typical MAEs of 81 cm-1 (PBE0: 76 cm-1). Excellent performance was found for adiabatic triplet excitation energies with typical MAEs of 0.08 eV (PBE0: 0.32 eV). In a detailed case anal. of the first singlet and triplet excited states of formaldehyde the conceptional (dis-)advantages of the local hybrid scheme for excited-state gradients are exposed.
- 40Kretz, B.; Egger, D. A. Accurate Molecular Geometries in Complex Excited-State Potential Energy Surfaces from Time-Dependent Density Functional Theory. J. Chem. Theory Comput. 2021, 17, 357– 366, DOI: 10.1021/acs.jctc.0c00858Google Scholar40https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXisVyrurjN&md5=74c6717ff8768412998953095ed42cf5Accurate Molecular Geometries in Complex Excited-State Potential Energy Surfaces from Time-Dependent Density Functional TheoryKretz, Bernhard; Egger, David A.Journal of Chemical Theory and Computation (2021), 17 (1), 357-366CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)The interplay of electronic excitations and structural changes in mols. impacts nonradiative decay and charge transfer in the excited state, thus influencing excited-state lifetimes and photocatalytic reaction rates in optoelectronic and energy devices. To capture such effects requires computational methods providing an accurate description of excited-state potential energy surfaces and geometries. We suggest time-dependent d. functional theory using optimally tuned range-sepd. hybrid (OT-RSH) functionals as an accurate approach to obtain excited-state mol. geometries. We show that OT-RSH provides accurate mol. geometries in excited-state potential energy surfaces that are complex and involve an interplay of local and charge-transfer excitations, for which conventional semilocal and hybrid functionals fail. At the same time, the nonempirical OT-RSH approach maintains the high accuracy of parametrized functionals (e.g., B3LYP) for predicting excited-state geometries of small org. mols. showing valence excited states.
- 41Sokolov, M.; Bold, B. M.; Kranz, J. J.; Höfener, S.; Niehaus, T. A.; Elstner, M. Analytical Time-Dependent Long-Range Corrected Density Functional Tight Binding (TD-LC-DFTB) Gradients in DFTB+: Implementation and Benchmark for Excited-State Geometries and Transition Energies. J. Chem. Theory Comput. 2021, 17, 2266– 2282, DOI: 10.1021/acs.jctc.1c00095Google Scholar41https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3MXlvFWit7Y%253D&md5=433eaf26623e19cbabe5f82254b17193Analytical Time-Dependent Long-Range Corrected Density Functional Tight Binding (TD-LC-DFTB) Gradients in DFTB+: Implementation and Benchmark for Excited-State Geometries and Transition EnergiesSokolov, Monja; Bold, Beatrix M.; Kranz, Julian J.; Hoefener, Sebastian; Niehaus, Thomas A.; Elstner, MarcusJournal of Chemical Theory and Computation (2021), 17 (4), 2266-2282CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)The absorption and emission of light is a ubiquitous process in chem. and biol. processes, making a theor. description inevitable for understanding and predicting such properties. Although ab initio and DFT methods are capable of describing excited states with good accuracy in many cases, the investigation of dynamical processes and the need to sample the phase space in complex systems often requires methods with reduced computational costs but still sufficient accuracy. In the present work, we report the derivation and implementation of anal. nuclear gradients for time-dependent long-range cor. d. functional tight binding (TD-LC-DFTB) in the DFTB+ program. The accuracy of the TD-LC-DFTB potential-energy surfaces is benchmarked for excited-state geometries and adiabatic as well as vertical transition energies. The benchmark set consists of more than 100 org. mols. taken as subsets from available benchmark sets. The reported method yields a mean deviation of 0.31 eV for adiabatic excitation energies with respect to CC2. In order to study more subtle effects, seminumerical second derivs. based on the anal. gradients are employed to simulate vibrationally resolved UV/vis spectra. This extensive test exhibits few problematic cases, which can be traced back to the parametrization of the repulsive potential.
- 42Kühne, T. D.; Iannuzzi, M.; Del Ben, M.; Rybkin, V. V.; Seewald, P.; Stein, F.; Laino, T.; Khaliullin, R. Z.; Schütt, O.; Schiffmann, F.; Golze, D.; Wilhelm, J.; Chulkov, S.; Bani-Hashemian, M. H.; Weber, V.; Borštnik, U.; Taillefumier, M.; Jakobovits, A. S.; Lazzaro, A.; Pabst, H.; Müller, T.; Schade, R.; Guidon, M.; Andermatt, S.; Holmberg, N.; Schenter, G. K.; Hehn, A.; Bussy, A.; Belleflamme, F.; Tabacchi, G.; Glöß, A.; Lass, M.; Bethune, I.; Mundy, C. J.; Plessl, C.; Watkins, M.; VandeVondele, J.; Krack, M.; Hutter, J. CP2K: An electronic structure and molecular dynamics software package - Quickstep: Efficient and accurate electronic structure calculations. J. Chem. Phys. 2020, 152, 194103, DOI: 10.1063/5.0007045Google Scholar42https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXhtVSgtrzF&md5=b9e5975bc402f0d53e2d99da998adf5fCP2K: An electronic structure and molecular dynamics software package - Quickstep: Efficient and accurate electronic structure calculationsKuehne, Thomas D.; Iannuzzi, Marcella; Del Ben, Mauro; Rybkin, Vladimir V.; Seewald, Patrick; Stein, Frederick; Laino, Teodoro; Khaliullin, Rustam Z.; Schuett, Ole; Schiffmann, Florian; Golze, Dorothea; Wilhelm, Jan; Chulkov, Sergey; Bani-Hashemian, Mohammad Hossein; Weber, Valery; Borstnik, Urban; Taillefumier, Mathieu; Jakobovits, Alice Shoshana; Lazzaro, Alfio; Pabst, Hans; Mueller, Tiziano; Schade, Robert; Guidon, Manuel; Andermatt, Samuel; Holmberg, Nico; Schenter, Gregory K.; Hehn, Anna; Bussy, Augustin; Belleflamme, Fabian; Tabacchi, Gloria; Gloess, Andreas; Lass, Michael; Bethune, Iain; Mundy, Christopher J.; Plessl, Christian; Watkins, Matt; VandeVondele, Joost; Krack, Matthias; Hutter, JuergJournal of Chemical Physics (2020), 152 (19), 194103CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)A review. CP2K is an open source electronic structure and mol. dynamics software package to perform atomistic simulations of solid-state, liq., mol., and biol. systems. It is esp. aimed at massively parallel and linear-scaling electronic structure methods and state-of-the-art ab initio mol. dynamics simulations. Excellent performance for electronic structure calcns. is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2K to perform efficient and accurate electronic structure simulations. The emphasis is put on d. functional theory and multiple post-Hartree-Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension. (c) 2020 American Institute of Physics.
- 43VandeVondele, J.; Krack, M.; Mohamed, F.; Parrinello, M.; Chassaing, T.; Hutter, J. Quickstep: Fast and accurate density functional calculations using a mixed Gaussian and plane waves approach. Comput. Phys. Commun. 2005, 167, 103– 128, DOI: 10.1016/j.cpc.2004.12.014Google Scholar43https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2MXjt1aitb4%253D&md5=8c5393031c9dbd341e0e73fcdacad486QUICKSTEP: fast and accurate density functional calculations using a mixed Gaussian and plane waves approachVandeVondele, Joost; Krack, Matthias; Mohamed, Fawzi; Parrinello, Michele; Chassaing, Thomas; Hutter, JuergComputer Physics Communications (2005), 167 (2), 103-128CODEN: CPHCBZ; ISSN:0010-4655. (Elsevier B.V.)We present the Gaussian and plane waves (GPW) method and its implementation in which is part of the freely available program package CP2K. The GPW method allows for accurate d. functional calcns. in gas and condensed phases and can be effectively used for mol. dynamics simulations. We show how derivs. of the GPW energy functional, namely ionic forces and the Kohn-Sham matrix, can be computed in a consistent way. The computational cost of computing the total energy and the Kohn-Sham matrix is scaling linearly with the system size, even for condensed phase systems of just a few tens of atoms. The efficiency of the method allows for the use of large Gaussian basis sets for systems up to 3000 atoms, and we illustrate the accuracy of the method for various basis sets in gas and condensed phases. Agreement with basis set free calcns. for single mols. and plane wave based calcns. in the condensed phase is excellent. Wave function optimization with the orbital transformation technique leads to good parallel performance, and outperforms traditional diagonalisation methods. Energy conserving Born-Oppenheimer dynamics can be performed, and a highly efficient scheme is obtained using an extrapolation of the d. matrix. We illustrate these findings with calcns. using commodity PCs as well as supercomputers.
- 44Lippert, G.; Hutter, J.; Parrinello, M. A hybrid Gaussian and plane wave density functional scheme. Mol. Phys. 1997, 92, 477– 488, DOI: 10.1080/00268979709482119Google Scholar44https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK2sXntVSrs7s%253D&md5=ece29370fc0b1b66ddd0b82b937bfd6bA hybrid Gaussian and plane wave density functional schemeLippert, Gerald; Hutter, Juerg; Parrinello, MicheleMolecular Physics (1997), 92 (3), 477-487CODEN: MOPHAM; ISSN:0026-8976. (Taylor & Francis)A d.-functional theory-based algorithm for periodic and nonperiodic ab initio calcns. is presented. This scheme uses pseudopotentials in order to integrate out the core electrons from the problem. The valence pseudo wave functions are expanded in Gaussian-type orbitals and the d. is represented in a plane wave auxiliary basis. The Gaussian basis functions make it possible to use the efficient anal. integration schemes and screening algorithms of quantum chem. Novel recursion relations are developed for the calcn. of the matrix elements of the d.-dependent Kohn-Sham self-consistent potential. At the same time the use of a plane wave basis for the electron d. permits efficient calcn. of the Hartree energy using fast Fourier transforms, thus circumventing one of the major bottlenecks of std. Gaussian based calcns. Furthermore, this algorithm avoids the fitting procedures that go along with intermediate basis sets for the charge d. The performance and accuracy of this new scheme are discussed and selected examples are given.
- 45Iannuzzi, M.; Chassaing, T.; Wallman, T.; Hutter, J. Ground and Excited State Density Functional Calculations with the Gaussian and Augmented-Plane-Wave Method. CHIMIA Int. J. Chem. 2005, 59, 499– 503, DOI: 10.2533/000942905777676164Google Scholar45https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2MXpvF2msL4%253D&md5=2ef7b146fbe6660372f872678e1bdc8eGround and excited state density functional calculations with the Gaussian and Augmented-Plane-Wave methodIannuzzi, Marcella; Chassaing, Thomas; Wallman, Thomas; Hutter, JurgChimia (2005), 59 (7-8), 499-503CODEN: CHIMAD; ISSN:0009-4293. (Swiss Chemical Society)The calcn. of the electronic structure of large systems by methods based on d. functional theory has recently gained a central role in mol. simulations. However, the extensive study of quantities like excited states and related properties is still out of reach due to high computational costs. We present a new implementation of a hybrid method, the Gaussian and APW (GAPW) method, where the electronic d. is partitioned in hard and soft contributions. The former are local terms naturally expanded in a Gaussian basis, whereas the soft contributions are expanded in plane-waves by using a low energy cutoff, without loss in accuracy, even for all-electron calcns. For the calcn. of excitation energies a recently developed, time-dependent d. functional response theory (TD-DFRT) technique is joined with the GAPW procedure. We demonstrate the accuracy of the method by comparison with std. quantum chem. calcns. for a set of small mols. To highlight the performance and efficiency of GAPW we show calcns. on systems with several thousands of basis functions.
- 46Strand, J.; Chulkov, S. K.; Watkins, M. B.; Shluger, A. L. First principles calculations of optical properties for oxygen vacancies in binary metal oxides. J. Chem. Phys. 2019, 150, 044702, DOI: 10.1063/1.5078682Google Scholar46https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXhvFOhtro%253D&md5=97b0d171cbcf065900b510bd2222a25eFirst principles calculations of optical properties for oxygen vacancies in binary metal oxidesStrand, Jack; Chulkov, Sergey K.; Watkins, Matthew B.; Shluger, Alexander L.Journal of Chemical Physics (2019), 150 (4), 044702/1-044702/10CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Using an advanced computational methodol. implemented in CP2K, a non-local PBE0-TC-LRC d. functional and the recently implemented linear response formulation of the Time-dependent D. Functional Theory equations, we test the interpretation of the optical absorption and photoluminescence signatures attributed by previous exptl. and theor. studies to O-vacancies in two widely used oxides-cubic MgO and monoclinic (m)-HfO2. The results obtained in large periodic cells including up to 1000 atoms emphasize the importance of accurate predictions of defect-induced lattice distortions. They confirm that optical transitions of O-vacancies in 0, +1, and +2 charge states in MgO all have energies close to 5 eV. We test the models of photoluminescence of O-vacancies proposed in the literature. The photoluminescence of V+2O centers in m-HfO2 is predicted to peak at 3.7 eV and originate from radiative tunneling transition between a V+1O center and a self-trapped hole created by the 5.2 eV excitation. (c) 2019 American Institute of Physics.
- 47Poli, E.; Elliott, J. D.; Chulkov, S. K.; Watkins, M. B.; Teobaldi, G. The Role of Cation-Vacancies for the Electronic and Optical Properties of Aluminosilicate Imogolite Nanotubes: A Non-local, Linear-Response TDDFT Study. Front. Chem. 2019, 7, 210, DOI: 10.3389/fchem.2019.00210Google Scholar47https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXhtlOnu7jK&md5=139ba03de31f5c581350d6fd5a16473dCation-vacancies for the electronic and optical properties of aluminosilicate imogolite nanotubes and a non-local, linear-response TDDFT studyPoli, Emiliano; Elliott, Joshua D.; Chulkov, Sergey K.; Watkins, Matthew B.; Teobaldi, GilbertoFrontiers in Chemistry (Lausanne, Switzerland) (2019), 7 (), 210CODEN: FCLSAA; ISSN:2296-2646. (Frontiers Media S.A.)We report a combined non-local (PBE-TC-LRC) D. Functional Theory (DFT) and linear-response time-dependent DFT (LR-TDDFT) study of the structural, electronic, and optical properties of the cation-vacancy based defects in aluminosilicate (AlSi) imogolite nanotubes (Imo-NTs) that have been recently proposed on the basis of NMR (NMR) expts. Following numerical detn. of the smallest AlSi Imo-NT model capable of accommodating the defect-induced relaxation with negligible finite-size errors, we analyze the defect-induced structural deformations in the NTs and ensuing changes in the NTs' electronic structure. The NMR-derived defects are found to introduce both shallow and deep occupied states in the pristine NTs' band gap (BG). These BG states are found to be highly localized at the defect site. No empty defect-state is modeled for any of the considered systems. LR-TDDFT simulation of the defects reveal increased low-energy optical absorbance for all but one defects, with the appearance of optically active excitations at energies lower than for the defect-free NT. These results enable interpretation of the low-energy tail in the exptl. UV-vis spectra for AlSi NTs as being due to the defects. Finally, the PBE-TC-LRC-approximated exciton binding energy for the defects' optical transitions is found to be substantially lower (up to 0.8 eV) than for the pristine defect-free NT's excitations (1.1 eV).
- 48Sternheimer, R. On Nuclear Quadrupole Moments. Phys. Rev. 1951, 84, 244– 253, DOI: 10.1103/PhysRev.84.244Google Scholar48https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaG38XhtlOltw%253D%253D&md5=dccfdcd6ac94c5ac879db40925286c28Nuclear quadrupole momentsSternheimer, R.Physical Review (1951), 84 (), 244-53CODEN: PHRVAO; ISSN:0031-899X.cf. C.A. 45, 443i. The correction to nuclear quadrupole moments on account of the quadrupole moment induced in the electron shells was obtained by solving the Schroedinger equation for the perturbed core wave functions for Li, Al, and Cl. The correction factor by which the average 〈l/r3〉 over the valence electron function in the equation for the quadrupole coupling should be multiplied to take account of the induced effect is 1.11, 0.83, and 0.68, resp. The previously described Thomas-Fermi calcn. of this effect was carried out for F, Cl, Cu, Br, Y, Ag, I, La, Lu, Pt, Tl, At, and U.
- 49Baroni, S.; de Gironcoli, S.; Dal Corso, A.; Giannozzi, P. Phonons and related crystal properties from density-functional perturbation theory. Rev. Mod. Phys. 2001, 73, 515– 562, DOI: 10.1103/RevModPhys.73.515Google Scholar49https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3MXlvFKrtLc%253D&md5=20ea8e1535ceb775168384a30fa2846dPhonons and related crystal properties from density-functional perturbation theoryBaroni, Stefano; De Gironcoli, Stefano; Dal Corso, Andrea; Giannozzi, PaoloReviews of Modern Physics (2001), 73 (2), 515-562CODEN: RMPHAT; ISSN:0034-6861. (American Physical Society)This article reviews with many refs. the current status of lattice-dynamical calcns. in crystals, using d.-functional perturbation theory, with emphasis on the plane-wave pseudopotential method. Several specialized topics are treated, including the implementation for metals, the calcn. of the response to macroscopic elec. fields and their relevance to long-wavelength vibrations in polar materials, the response to strain deformations, and higher-order responses. The success of this methodol. is demonstrated with a no. of applications existing in the literature.
- 50Sternheimer, R. M. Electronic Polarizabilities of Ions from the Hartree-Fock Wave Functions. Phys. Rev. 1954, 96, 951– 968, DOI: 10.1103/PhysRev.96.951Google Scholar50https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaG2MXhsVOktw%253D%253D&md5=629629b68b30bee1c7fe6de9cca98a06Electron polarizabilities of ions from the Hartree-Fock wave functionsSternheimer, R. M.Physical Review (1954), 96 (), 951-68CODEN: PHRVAO; ISSN:0031-899X.cf. ibid. 95, 655. The electronic polarizability α was calcd. for several ions by obtaining the perturbation of the wave functions by an external field from a numerical solution of the differential equation satisfied by the perturbation. For the He-like ions an analytic solution was obtained by using the wave functions of L.ovrddot.owdin (C.A. 47, 9135c). The calcd. values of α are, in general, 1 to 1.5 times the observed values. For several ions values were calcd. for the quadrupole polarizability, which measures the quadrupole moment induced in the ion by an external charge. The effect of the dipole moment induced in the ion on the elec. field at the nucleus is discussed.
- 51Budzak, S.; Scalmani, G.; Jacquemin, D. Accurate Excited-State Geometries: A CASPT2 and Coupled-Cluster Reference Database for Small Molecules. J. Chem. Theory Comput. 2017, 13, 6237– 6252, DOI: 10.1021/acs.jctc.7b00921Google Scholar51https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhvVWlsb3I&md5=9b7945dbdfd31d2b2055a7efbe1eacfbAccurate Excited-State Geometries: A CASPT2 and Coupled-Cluster Reference Database for Small MoleculesBudzak, Simon; Scalmani, Giovanni; Jacquemin, DenisJournal of Chemical Theory and Computation (2017), 13 (12), 6237-6252CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)We present an investigation of the excited-state structural parameters detd. for a large set of small compds. with the dual goals of defining ref. values for further works and assessing the quality of the geometries obtained with relatively cheap computational approaches. In the first stage, we compare the excited-state geometries obtained with ADC(2), CC2, CCSD, CCSDR(3), CC3, and CASPT2 and large at. basis sets. It is found that CASPT2 and CC3 results are generally in very good agreement with one another (typical differences of ca. 3 × 10-3 Å) when all electrons are correlated and when the aug-cc-pVTZ at. basis set is employed with both methods. In a second stage, a statistical anal. reveals that, on the one hand, the excited-state (ES) bond lengths are much more sensitive to the selected level of theory than their ground-state (GS) counterparts and, on the other hand, that CCSDR(3) is probably the most cost-effective method delivering accurate structures. Indeed, CCSD tends to provide too compact multiple bond lengths on an almost systematic basis, whereas both CC2 and ADC(2) tend to exaggerate these bond distances, with more erratic error patterns, esp. for the latter method. The deviations are particularly marked for the polarized CO and CN bonds, as well as for the puckering angle in formaldehyde homologues. In the last part of this contribution, we provide a series of CCSDR(3) GS and ES geometries of medium-sized mols. to be used as refs. in further investigations.
- 52Ongari, D.; Yakutovich, A. V.; Talirz, L.; Smit, B. Building a Consistent and Reproducible Database for Adsorption Evaluation in Covalent–Organic Frameworks. ACS Cent. Sci. 2019, 5, 1663– 1675, DOI: 10.1021/acscentsci.9b00619Google Scholar52https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXhvVeqsL%252FM&md5=8f47a5e3d00532e7064d5259e73ab1f1Building a Consistent and Reproducible Database for Adsorption Evaluation in Covalent-Organic FrameworksOngari, Daniele; Yakutovich, Aliaksandr V.; Talirz, Leopold; Smit, BerendACS Central Science (2019), 5 (10), 1663-1675CODEN: ACSCII; ISSN:2374-7951. (American Chemical Society)We present a workflow that traces the path from the bulk structure of a cryst. material to assessing its performance in carbon capture from coal's postcombustion flue gases. This workflow is applied to a database of 324 covalent-org. frameworks (COFs) reported in the literature, to characterize their CO2 adsorption properties using the following steps: (1) optimization of the crystal structure (at. positions and unit cell) using d. functional theory, (2) fitting at. point charges based on the electron d., (3) characterizing the pore geometry of the structures before and after optimization, (4) computing carbon dioxide and nitrogen isotherms using grand canonical Monte Carlo simulations with an empirical interaction potential, and finally, (5) assessing the CO2 parasitic energy via process modeling. The full workflow has been encoded in the Automated Interactive Infrastructure and Database for Computational Science (AiiDA). Both the workflow and the automatically generated provenance graph of our calcns. are made available on the Materials Cloud, allowing peers to inspect every input parameter and result along the workflow, download structures and files at intermediate stages, and start their research right from where this work has left off. In particular, our set of CURATED (Clean, Uniform, and Refined with Automatic Tracking from Exptl. Database) COFs, having optimized geometry and high-quality DFT-derived point charges, are available for further investigations of gas adsorption properties. We plan to update the database as new COFs are being reported. An automated and reproducible computational workflow is proposed, to systematically optimize the geometry of covalent-org. frameworks and evaluate their performances for carbon capture and storage.
- 53Ongari, D.; Talirz, L.; Smit, B. Too Many Materials and Too Many Applications: An Experimental Problem Waiting for a Computational Solution. ACS Cent. Sci. 2020, 6, 1890– 1900, DOI: 10.1021/acscentsci.0c00988Google Scholar53https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXhvFOgtL3N&md5=56a120a637317936ec763bc398505969Too Many Materials and Too Many Applications: An Experimental Problem Waiting for a Computational SolutionOngari, Daniele; Talirz, Leopold; Smit, BerendACS Central Science (2020), 6 (11), 1890-1900CODEN: ACSCII; ISSN:2374-7951. (American Chemical Society)A review. Finding the best material for a specific application is the ultimate goal of materials discovery. However, there is also the reverse problem: when exptl. groups discover a new material, they would like to know all the possible applications this material would be promising for. Computational modeling can aim to fulfill this expectation, thanks to the sustained growth of computing power and the collective engagement of the scientific community in developing more efficient and accurate workflows for predicting materials' performances. We discuss the impact that reproducibility and automation of the modeling protocols have on the field of gas adsorption in nanoporous crystals. We envision a platform that combines these tools and enables effective matching between promising materials and industrial applications. We identify the opportunity for a computational platform for matching nanoporous materials and gas-related applications, motivating the development of automated and reproducible computational workflows.
- 54Hirata, S.; Head-Gordon, M. Time-dependent density functional theory within the Tamm-Dancoff approximation. Chem. Phys. Lett. 1999, 314, 291– 299, DOI: 10.1016/S0009-2614(99)01149-5Google Scholar54https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK1MXnsl2rsr0%253D&md5=1b62f410de6c2a2193f1011d42f389c5Time-dependent density functional theory within the Tamm-Dancoff approximationHirata, S.; Head-Gordon, M.Chemical Physics Letters (1999), 314 (3,4), 291-299CODEN: CHPLBC; ISSN:0009-2614. (Elsevier Science B.V.)A computationally simple method for mol. excited states, namely, the Tamm-Dancoff approxn. to time-dependent d. functional theory, is proposed and implemented. This method yields excitation energies for several closed- and open-shell mols. that are essentially of the same quality as those obtained from time-dependent d. functional theory itself, when the same exchange-correlation functional is used.
- 55Casida, M. E. In Recent Advances in Density Functional Methods, Part I; Chong, D. P., Ed.; World Scientific: Singapore, 1995; p 155.Google ScholarThere is no corresponding record for this reference.
- 56Guidon, M.; Hutter, J.; VandeVondele, J. Auxiliary Density Matrix Methods for Hartree-Fock Exchange Calculations. J. Chem. Theory Comput. 2010, 6, 2348– 2364, DOI: 10.1021/ct1002225Google Scholar56https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXot1als7Y%253D&md5=7a84c086d49f17bab239fd3f8d004c6eAuxiliary Density Matrix Methods for Hartree-Fock Exchange CalculationsGuidon, Manuel; Hutter, Jurg; Vande Vondele, JoostJournal of Chemical Theory and Computation (2010), 6 (8), 2348-2364CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)The calcn. of Hartree-Fock exchange (HFX) is computationally demanding for large systems described with high-quality basis sets. In this work, we show that excellent performance and good accuracy can nevertheless be obtained if an auxiliary d. matrix is employed for the HFX calcn. Several schemes to derive an auxiliary d. matrix from a high-quality d. matrix are discussed. Key to the accuracy of the auxiliary d. matrix methods (ADMM) is the use of a correction based on std. generalized gradient approxns. for HFX. ADMM integrates seamlessly in existing HFX codes and, in particular, can be employed in linear scaling implementations. Demonstrating the performance of the method, the effect of HFX on the structure of liq. water is investigated in detail using Born-Oppenheimer mol. dynamics simulations (300 ps) of a system of 64 mols. Representative for large systems are calcns. on a solvated protein (Rubredoxin), for which ADMM outperforms the corresponding std. HFX implementation by approx. a factor 20.
- 57Ghosh, D. C.; Islam, N. Whether electronegativity and hardness are manifest two different descriptors of the one and the same fundamental property of atoms? @ TA quest. Int. J. Quantum Chem. 2011, 111, 40– 51, DOI: 10.1002/qua.22415Google Scholar57https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXhtlaqt7%252FE&md5=e601f3a812a1e9f147bfeb376c699815Whether electronegativity and hardness are manifest two different descriptors of the one and the same fundamental property of atoms-A questGhosh, Dulal C.; Islam, NazmulInternational Journal of Quantum Chemistry (2011), 111 (1), 40-51CODEN: IJQCB2; ISSN:0020-7608. (John Wiley & Sons, Inc.)In this report, we have tried to reveal that there is much conceptual commonality between the two fundamental theor. descriptors of chem. and physics-the electronegativity and the hardness. The phys. hardness was introduced and theorized by condensed matter physicists. The chem. hardness was introduced by chemists to generalize and rationalize the HSAB principle. We have tried to establish that the phys. hardness and the chem. hardness with evolution of time have converged to one and the same general principle-the hardness. We have also tried to understand the phys. basis and operational significance of another very important descriptor arising out of theor. constructs of chem.-the electronegativity. We have relied upon the fact that, since these descriptors are not observables, there is no possibility of their quantum mech. evaluation. These descriptors, therefore, should be and must be reified before suggesting ansatz for their evaluation. We have dwelt at length upon the effort of d. functional definition and evaluation of electronegativity and hardness and discovered the inherent inner contradiction of the theory and measurement. We have also noted that a good no. of scientists hold the opinion that the d. functional formula of electronegativity is χ = I and that of hardness is η = I, where I is the ionization potential of the chem. system. This study concludes that the two fundamental descriptors-the hardness and the electronegativity originate from the same source-the electron attracting power of the screened nucleus upon valence electrons and discovers the surprising result that if one measures hardness, the electronegativity is simultaneously measured and vice versa. We have also explored the ansatz of semiempirical evaluation of electronegativity and hardness of at. systems involving the radius of the atoms. We have noted that the ansatz for electronegativity and hardness is the same. To justify our statement that if one measures hardness, the electronegativity is simultaneously measured and vice versa, we have used the evaluated set of hardness for 103 elements of the periodic table as a scale of electronegativity and found that such set of hardness data satisfies the sine qua non of a satisfactory scale of electronegativity. The electronegativity and the hardness are two different appearances of the one and the same fundamental property of atoms. They are different and also nondifferent. They are different in their fields of application. They are nondifferent when we discuss the basic philosophical structures of their origin and development. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2011.
- 58Marx, D.; Hutter, J. In Modern Methods and Algorithms of Quantum Chemistry; Grotendorst, J., Ed.; John von Neumann Institute for Computing, Forschungszentrum Jülich: Jülich, Germany, 2000; pp 301– 449 (first edition, paperback, ISBN 3-00-005618-1) or pp 329–477 (second edition, hardcover, ISBN 3-00-005834-6), see http://www.theochem.rub.de/go/cprev.html.Google ScholarThere is no corresponding record for this reference.
- 59Essmann, U.; Perera, L.; Berkowitz, M. L.; Darden, T.; Lee, H.; Pedersen, L. G. A smooth particle mesh Ewald method. J. Chem. Phys. 1995, 103, 8577– 8593, DOI: 10.1063/1.470117Google Scholar59https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK2MXptlehtrw%253D&md5=092a679dd3bee08da28df41e302383a7A smooth particle mesh Ewald methodEssmann, Ulrich; Perera, Lalith; Berkowitz, Max L.; Darden, Tom; Lee, Hsing; Pedersen, Lee G.Journal of Chemical Physics (1995), 103 (19), 8577-93CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The previously developed particle mesh Ewald method is reformulated in terms of efficient B-spline interpolation of the structure factors. This reformulation allows a natural extension of the method to potentials of the form 1/rp with p ≥ 1. Furthermore, efficient calcn. of the virial tensor follows. Use of B-splines in the place of Lagrange interpolation leads to analytic gradients as well as a significant improvement in the accuracy. The authors demonstrate that arbitrary accuracy can be achieved, independent of system size N, at a cost that scales as N log(N). For biomol. systems with many thousands of atoms and this method permits the use of Ewald summation at a computational cost comparable to that of a simple truncation method of 10 Å or less.
- 60Mayhall, N. J.; Raghavachari, K.; Hratchian, H. P. ONIOM-based QM:QM electronic embedding method using Löwdin atomic charges: Energies and analytic gradients. J. Chem. Phys. 2010, 132, 114107, DOI: 10.1063/1.3315417Google Scholar60https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXjslWhsb4%253D&md5=2ac557eea7ce3b5887cf31841ac6508dONIOM-based QM:QM electronic embedding method using Lowdin atomic charges: Energies and analytic gradientsMayhall, Nicholas J.; Raghavachari, Krishnan; Hratchian, Hrant P.Journal of Chemical Physics (2010), 132 (11), 114107/1-114107/6CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We report a new quantum mech.:quantum mech. (QM:QM) method which provides explicit electronic polarization of the high-level region by using the Lowdin at. charges from the low-level region. This provides an embedding potential which naturally evolves with changes in nuclear geometry. However, this coupling of the high-level and low-level regions introduces complications in the energy gradient evaluation. Following previous work, we derive and implement efficient gradients where a single set of SCF response equations is solved. We provide results for the calcn. of deprotonation energies of a hydroxylated spherosiloxane cluster (Si8O12H7OH) and the dissocn. energy of a water mol. from a [ZnIm3(H2O)]2+ complex. The Lowdin charge embedding model provides results which are not only an improvement over mech. embedding (no electronic embedding) but which are also resistant to large overpolarization effects which occur more often with Mulliken charge embedding. Finally, a scaled-Loewdin charge embedding method is also presented which provides a method for fine tuning the extent of electronic polarization. (c) 2010 American Institute of Physics.
- 61Davidson, E. R. The iterative calculation of a few of the lowest eigenvalues and corresponding eigenvectors of large real-symmetric matrices. J. Comput. Phys. 1975, 17, 87– 94, DOI: 10.1016/0021-9991(75)90065-0Google ScholarThere is no corresponding record for this reference.
- 62Chassaing, T. Excitation Energy Calculations with TD-DFT. PhD thesis, University of Zurich, Zurich, Switzerland, 2005.Google ScholarThere is no corresponding record for this reference.
- 63Ryabinkin, I. G.; Nagesh, J.; Izmaylov, A. F. Fast Numerical Evaluation of Time-Derivative Nonadiabatic Couplings for Mixed Quantum–Classical Methods. J. Phys. Chem. Lett. 2015, 6, 4200– 4203, DOI: 10.1021/acs.jpclett.5b02062Google Scholar63https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2MXhs1Sis7%252FF&md5=2a46a6400734c540a5a72ed2c8b24057Fast Numerical Evaluation of Time-Derivative Nonadiabatic Couplings for Mixed Quantum-Classical MethodsRyabinkin, Ilya G.; Nagesh, Jayashree; Izmaylov, Artur F.Journal of Physical Chemistry Letters (2015), 6 (21), 4200-4203CODEN: JPCLCD; ISSN:1948-7185. (American Chemical Society)We have developed a numerical differentiation scheme that eliminates evaluation of overlap determinants in calcg. the time-deriv. nonadiabatic couplings (TDNACs). Evaluation of these determinants was the bottleneck in previous implementations of mixed quantum-classical methods using numerical differentiation of electronic wave functions in the Slater determinant representation. The central idea of our approach is, first, to reduce the analytic time derivs. of Slater determinants to time derivs. of MOs and then to apply a finite-difference formula. Benchmark calcns. prove the efficiency of the proposed scheme showing impressive several-order-of-magnitude speedups of the TDNAC calcn. step for midsize mols.
- 64Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Petersson, G. A.; Nakatsuji, H.; Li, X.; Caricato, M.; Marenich, A. V.; Bloino, J.; Janesko, B. G.; Gomperts, R.; Mennucci, B.; Hratchian, H. P.; Ortiz, J. V.; Izmaylov, A. F.; Sonnenberg, J. L.; Williams-Young, D.; Ding, F.; Lipparini, F.; Egidi, F.; Goings, J.; Peng, B.; Petrone, A.; Henderson, T.; Ranasinghe, D.; Zakrzewski, V. G.; Gao, J.; Rega, N.; Zheng, G.; Liang, W.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Throssell, K.; Montgomery, J. A., Jr.; Peralta, J. E.; Ogliaro, F.; Bearpark, M. J.; Heyd, J. J.; Brothers, E. N.; Kudin, K. N.; Staroverov, V. N.; Keith, T. A.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A. P.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Millam, J. M.; Klene, M.; Adamo, C.; Cammi, R.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Farkas, O.; Foresman, J. B.; Fox, D. J. Gaussian 16, Revision C.01; Gaussian Inc.: Wallingford, CT, 2016.Google ScholarThere is no corresponding record for this reference.
- 65Weigend, F.; Ahlrichs, R. Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. Phys. Chem. Chem. Phys. 2005, 7, 3297, DOI: 10.1039/b508541aGoogle Scholar65https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2MXpsFWgu7o%253D&md5=a820fb6055c993b50c405ba0fc62b194Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracyWeigend, Florian; Ahlrichs, ReinhartPhysical Chemistry Chemical Physics (2005), 7 (18), 3297-3305CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)Gaussian basis sets of quadruple zeta valence quality for Rb-Rn are presented, as well as bases of split valence and triple zeta valence quality for H-Rn. The latter were obtained by (partly) modifying bases developed previously. A large set of more than 300 mols. representing (nearly) all elements-except lanthanides-in their common oxidn. states was used to assess the quality of the bases all across the periodic table. Quantities investigated were atomization energies, dipole moments and structure parameters for Hartree-Fock, d. functional theory and correlated methods, for which we had chosen Moller-Plesset perturbation theory as an example. Finally recommendations are given which type of basis set is used best for a certain level of theory and a desired quality of results.
- 66Balasubramani, S. G.; Chen, G. P.; Coriani, S.; Diedenhofen, M.; Frank, M. S.; Franzke, Y. J.; Furche, F.; Grotjahn, R.; Harding, M. E.; Hättig, C.; Hellweg, A.; Helmich-Paris, B.; Holzer, C.; Huniar, U.; Kaupp, M.; Marefat Khah, A.; Karbalaei Khani, S.; Müller, h.; Mack, F.; Nguyen, B. D.; Parker, S. M.; Perlt, E.; Rappoport, D.; Reiter, K.; Roy, S.; Rückert, M.; Schmitz, G.; Sierka, M.; Tapavicza, E.; Tew, D. P.; van Wüllen, C.; Voora, V. K.; Weigend, F.; Wodyński, A.; Yu, J. M. TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations. J. Chem. Phys. 2020, 152, 184107, DOI: 10.1063/5.0004635Google Scholar66https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXps1Ogs7s%253D&md5=836d350d0cdc56aeea644742a098b23bTURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulationsBalasubramani, Sree Ganesh; Chen, Guo P.; Coriani, Sonia; Diedenhofen, Michael; Frank, Marius S.; Franzke, Yannick J.; Furche, Filipp; Grotjahn, Robin; Harding, Michael E.; Hattig, Christof; Hellweg, Arnim; Helmich-Paris, Benjamin; Holzer, Christof; Huniar, Uwe; Kaupp, Martin; Marefat Khah, Alireza; Karbalaei Khani, Sarah; Muller, Thomas; Mack, Fabian; Nguyen, Brian D.; Parker, Shane M.; Perlt, Eva; Rappoport, Dmitrij; Reiter, Kevin; Roy, Saswata; Ruckert, Matthias; Schmitz, Gunnar; Sierka, Marek; Tapavicza, Enrico; Tew, David P.; van Wullen, Christoph; Voora, Vamsee K.; Weigend, Florian; Wodynski, Artur; Yu, Jason M.Journal of Chemical Physics (2020), 152 (18), 184107CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)A review. TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chem. simulations of mols., clusters, periodic systems, and solns. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy-cost ratio, such as d. functional theory including local hybrids and the RPA, GW-Bethe-Salpeter methods, second-order Moller-Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resoln.-of-the-identity approxn., imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent addns. to TURBOMOLE's functionality, including excited-state methods, RPA and Green's function methods, relativistic approaches, high-order mol. properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE's current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE's development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted. (c) 2020 American Institute of Physics.
- 67Weigend, F.; Furche, F.; Ahlrichs, R. Gaussian basis sets of quadruple zeta valence quality for atoms H-Kr. J. Chem. Phys. 2003, 119, 12753– 12762, DOI: 10.1063/1.1627293Google Scholar67https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3sXpvFOns7s%253D&md5=9d1664e0004c3a900b7efcb5da10dd6dGaussian basis sets of quadruple zeta valence quality for atoms H-KrWeigend, Florian; Furche, Filipp; Ahlrichs, ReinhartJournal of Chemical Physics (2003), 119 (24), 12753-12762CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We present Gaussian basis sets of quadruple zeta valence quality with a segmented contraction scheme for atoms H to Kr. This extends earlier work on segmented contracted split valence (SV) and triple zeta valence (TZV) basis sets. Contraction coeffs. and orbital exponents are fully optimized in at. Hartree-Fock (HF) calcns. As opposed to other quadruple zeta basis sets, the basis set errors in at. ground-state HF energies are less than 1 mEh and increase smoothly across the Periodic Table, while the no. of primitives is comparably small. Polarization functions are taken partly from previous work, partly optimized in at. MP2 calcns., and for a few cases detd. at the HF level for excited at. states nearly degenerate with the ground state. This leads to basis sets denoted QZVP for HF and d. functional theory (DFT) calcns., and for some atoms to a larger basis recommended for correlated treatments, QZVPP. We assess the performance of the basis sets in mol. HF, DFT, and MP2 calcns. for a sample of diat. and small polyat. mols. by a comparison of energies, bond lengths, and dipole moments with results obtained numerically or using very large basis sets. It is shown that basis sets of quadruple zeta quality are necessary to achieve an accuracy of 1 kcal/mol per bond in HF and DFT atomization energies. For compds. contg. third row as well as alk. and earth alk. metals it is demonstrated that the inclusion of high-lying core orbitals in the active space can be necessary for accurate correlated treatments. The QZVPP basis sets provide sufficient flexibility to polarize the core in those cases. All test calcns. indicate that the new basis sets lead to consistent accuracies in HF, DFT, or correlated treatments even in crit. cases where other basis sets may show deficiencies.
- 68VandeVondele, J.; Hutter, J. Gaussian basis sets for accurate calculations on molecular systems in gas and condensed phases. J. Chem. Phys. 2007, 127, 114105, DOI: 10.1063/1.2770708Google Scholar68https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2sXhtFSrsLvM&md5=d7fdb937efb88cf3fca85792bb49ec27Gaussian basis sets for accurate calculations on molecular systems in gas and condensed phasesVandeVondele, Joost; Hutter, JurgJournal of Chemical Physics (2007), 127 (11), 114105/1-114105/9CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We present a library of Gaussian basis sets that has been specifically optimized to perform accurate mol. calcns. based on d. functional theory. It targets a wide range of chem. environments, including the gas phase, interfaces, and the condensed phase. These generally contracted basis sets, which include diffuse primitives, are obtained minimizing a linear combination of the total energy and the condition no. of the overlap matrix for a set of mols. with respect to the exponents and contraction coeffs. of the full basis. Typically, for a given accuracy in the total energy, significantly fewer basis functions are needed in this scheme than in the usual split valence scheme, leading to a speedup for systems where the computational cost is dominated by diagonalization. More importantly, binding energies of hydrogen bonded complexes are of similar quality as the ones obtained with augmented basis sets, i.e., have a small (down to 0.2 kcal/mol) basis set superposition error, and the monomers have dipoles within 0.1 D of the basis set limit. However, contrary to typical augmented basis sets, there are no near linear dependencies in the basis, so that the overlap matrix is always well conditioned, also, in the condensed phase. The basis can therefore be used in first principles mol. dynamics simulations and is well suited for linear scaling calcns.
- 69Brémond, E.; Savarese, M.; Adamo, C.; Jacquemin, D. Accuracy of TD-DFT Geometries: A Fresh Look. J. Chem. Theory Comput. 2018, 14, 3715– 3727, DOI: 10.1021/acs.jctc.8b00311Google Scholar69https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXhtFSrtL%252FE&md5=bb253fe791ccd2ba3b30042b9b9db4c0Accuracy of TD-DFT Geometries: A Fresh LookBremond, Eric; Savarese, Marika; Adamo, Carlo; Jacquemin, DenisJournal of Chemical Theory and Computation (2018), 14 (7), 3715-3727CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)We benchmark a panel of 48 DFT exchange-correlation functionals in the framework of TD-DFT optimizations of the geometry of valence singlet excited states. To this end, we use a set of 41 small- and medium-sized org. mols. for which ref. geometries were obtained at high level of theory, typically, CC3 or CCSDR(3), with the aug-cc-pVTZ at. basis set. For the ground-state parameters, the tested functionals provide av. deviations that are small (0.010 Å and 0.5° for bond lengths and valence angles) and not very sensitive to the selected (hybrid) functional, but the errors are larger for the most polarized bonds (CO, CN, and so on). Nevertheless, DFT has a tendency to provide too compact distances, a trend slightly enhanced for functionals including a large amt. of exact exchange. The av. errors largely increase when going to the excited-state for most bond types, i.e., TD-DFT delivers less accurate excited-state distances than DFT for ground state. In particular TD-DFT combined with hybrid functionals provides significantly too short CO and CS/CSe bonds with resp. av. errors in the -0.026/-0.052 Å and -0.015/-0.082 Å ranges, depending on the selected hybrid functional. For the carbonyl bonds, the sizes of the TD-DFT deviations obtained when selecting std. hybrid functionals are of the same order of magnitude as the EOM-CCSD ones.
- 70Chan, G. K.-L.; Handy, N. C. C8H8: a density functional theory study of molecular geometries introducing the localised bond density. J. Chem. Soc., Faraday Trans. 1996, 92, 3015– 3021, DOI: 10.1039/ft9969203015Google Scholar70https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK28XlsVGhtbo%253D&md5=9c09739d70a3d81a2a7c13aee9c25355C8H8: A density functional theory study of molecular geometries introducing the localized bond densityChan, Garnet K-L.; Handy, Nicholas C.Journal of the Chemical Society, Faraday Transactions (1996), 92 (17), 3015-3021CODEN: JCFTEV; ISSN:0956-5000. (Royal Society of Chemistry)D. functional theory was used with all the common exchange-correlation functionals to investigate the structures of three isomers of C8H8 found in F.A. Cotton's text, barrelene, cyclooctatetraene, tetramethylenecyclobutane and also ethane and ethene. All calcns. were performed with TZ2P basis sets and large quadrature. The results are compared with expt. and those obtained with Hartree-Fock theory. Delocalization in the three mols. is discussed. A localized bond d. is introduced to explain the transferability of the trends in the predictions of the functionals between different mols. Three-parameter adiabatic connection functionals are examd. and their usefulness in geometry prediction questioned. Finally a phys. picture of the correlation as modeled by d. functional theory is presented and used to explain trends in the overestimation or underestimation of bond lengths.
- 71Neumann, R.; Nobes, R. H.; Handy, N. C. Exchange functionals and potentials. Mol. Phys. 1996, 87, 1– 36, DOI: 10.1080/00268979600100011Google Scholar71https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK28XhvFegs7c%253D&md5=14030dab7e55c7c6faaf9737dbcbcb37Exchange functionals and potentialsNeumann, Ralf; Nobes, Ross H.; Handy, Nicholas C.Molecular Physics (1996), 87 (1), 1-36CODEN: MOPHAM; ISSN:0026-8976. (Taylor & Francis)The commonly used exchange-correlation functionals of d. functional theory and their potentials are examd. numerically following the first such investigation of Perdew. They are also investigated for Ne and Kr. Their behavior for large gradients of the d. and for large distances is not satisfactory. In particular, the correct asymptotic r-1 behavior is difficult to achieve. Following van Leeuwen and Baerends, this is linked to the energy εmax of the highest occupied orbital arising from the Kohn-Sham equations. This deficiency is linked also with the poor prediction of mol. polarizabilities. The Becke-Roussel (BR) exchange functional is examd., which is derived by assuming a hydrogen-like exchange hole at all spatial points, and it has the attraction of being dependent on both the kinetic energy d. and the Laplacian of the d. and has no adjustable parameters. Becke has presented encouraging results using this functional in a hybrid manner. Fully self-consistent Kohn-sham calcns. are performed using it in combination with Perdew's 1986 correlation functional. The results are very encouraging indeed, so much so that this exchange functional is the best generalized gradient approxn. (GGA) yet discovered. In particular, bond lengths of many mol. show a substantial improvement over results from other GGAs. For example, many CH bonds are now within exptl. accuracy, instead of being typically 0.02 Å too long. Our ab initio understanding of non-dynamic correlation and dynamic correlation is then linked with d. functional theory. It is argued that correlation functionals should pick up the local dynamic correlation, whereas exchange functionals should include non-dynamic correlation effects. For these reasons it is considered that exchange functionals are best modeled on a system for which there is effectively no non-dynamic correlation, for which the optimum example is the Ne atom. Thus, again following Becke and Roussel, the spherically averaged Hartree-Fock exchange hole for Ne is examd., compared with the BR model functional hole. An excellent overlap is found, and thus the above good results are explained. As a final contribution, the dissocn. of the H2 mol. is re-examd., looking at it in terms of the exchange hole. For a ref. electron near one proton A, the RHF model has half an exchange electron near it, and half the exchange electron near the other proton B, whereas the BR functional has one electron near the other proton b, whereas the Br functional has one electron near A, which is the correct picture. For this reason the (restricted) BR functional gives a greatly improved dissocn. curve for H2 when compared with the Hartree-Fock curve. In summary, the Becke-Roussel functional is found to be a most attractive exchange functional.
- 72Sancho-García, J. C.; Cornil, J. Assessment of recently developed exchange-correlation functionals for the description of torsion potentials in π-conjugated molecules. J. Chem. Phys. 2004, 121, 3096– 3101, DOI: 10.1063/1.1774976Google Scholar72https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2cXmt1Slsrg%253D&md5=e1684f5d68a4a19ab35897b8e1f18dd0Assessment of recently developed exchange-correlation functionals for the description of torsion potentials in π-conjugated moleculesSancho-Garcia, Juan Carlos; Cornil, JeromeJournal of Chemical Physics (2004), 121 (7), 3096-3101CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Newly developed exchange-correlation functionals in d. functional theory (DFT) have been applied to describe conjugation effects in org. mols. The performance of the various approaches is assessed through the calcn. of torsion energy profiles and their crit. comparison with available exptl. data. Our results indicate that the OPTX-B95 exchange-correlation functional as well as its corresponding hybrid versions perform better than the well-established BLYP or B3LYP schemes when dealing with π-conjugated mols. In contrast, the recently introduced VSXC functional is not as reliable as other DFT methods for the systems examd. here.
- 73Vikramaditya, T.; Lin, S.-T. Limitations of Global Hybrids in Predicting the Geometries and Torsional Energy Barriers of Dimeric Systems and the Role of Hartree Fock and DFT Exchange. J. Comput. Chem. 2019, 40, 2810– 2818, DOI: 10.1002/jcc.26056Google Scholar73https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXhs1CitrvE&md5=9dc520039c1dbe19cee658bdef98f2f0Limitations of Global Hybrids in Predicting the Geometries and Torsional Energy Barriers of Dimeric Systems and the Role of Hartree Fock and DFT ExchangeVikramaditya, Talapunur; Lin, Shiang-TaiJournal of Computational Chemistry (2019), 40 (32), 2810-2818CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)Prediction of accurate geometries is a prerequisite for accurate prediction of mol. properties. Impact of Hartree Fock (HF) exchange (a0) on geometry in the framework of DFT is investigated by monitoring dihedral angles, bond length alternations, and torsional energy barriers of 10 dimeric systems against CCSD (ADZ/ATZ) benchmarks. A strong correlation is obsd. between the fraction of HF exchange, equil. dihedral angles, and the potential energy barriers in global hybrids. Full HF exchange is crit. to accurately predict the nonplanarity. Lower fractions of (a0)/larger DFT exchange (1-a0) results in overestimation of torsional energy barriers at 900 and underestimation at 00. Large contributions of (1-a0) in global hybrid functionals tend to overestimate torsional energy barriers (900) and are biased toward planar geometries. However, inclusion of larger fractions of (a0)/lower (1-a0) also overestimate the torsional energy barriers in syn-conformations due to the localization errors assocd. with HF exchange in global hybrids. Hence, irresp. of the fraction of HF/DFT exchange incorporated, global hybrids fail to accurately predict torsional energy barriers at 00 and 900 simultaneously. Long-range cor. (LC) functionals, which employ full HF exchange at longer regions, outperform global hybrid functionals in predicting geometries and torsional energy barriers of the dimeric mols. The distance dependence of (a0) thus provides a balanced fraction of HF exchange as the dihedral torsion varies. Impact of range sepn. parameter on geometries is marginal in altering the planarity/nonplanarity. However, range sepn. parameter within 0.20-0.40 bohr-1 predicts more reliable torsional energies and geometries. © 2019 Wiley Periodicals, Inc.
- 74Fang, C.; Oruganti, B.; Durbeej, B. How Method-Dependent Are Calculated Differences between Vertical, Adiabatic, and 0─0 Excitation Energies?. J. Phys. Chem. A 2014, 118, 4157– 4171, DOI: 10.1021/jp501974pGoogle Scholar74https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXotlWht7Y%253D&md5=ca7e486acb318eb2955f3b03fd91675cHow Method-Dependent Are Calculated Differences between Vertical, Adiabatic, and 0-0 Excitation Energies?Fang, Changfeng; Oruganti, Baswanth; Durbeej, BoJournal of Physical Chemistry A (2014), 118 (23), 4157-4171CODEN: JPCAFH; ISSN:1089-5639. (American Chemical Society)Through a large no. of benchmark studies, the performance of different quantum chem. methods in calcg. vertical excitation energies is today quite well established. These efforts have in recent years been complemented by a few benchmarks focusing instead on adiabatic excitation energies. However, it is much less well established how calcd. differences between vertical, adiabatic and 0-0 excitation energies vary between methods, which may be due to the cost of evaluating zero-point vibrational energy corrections for excited states. To fill this gap, we have calcd. vertical, adiabatic, and 0-0 excitation energies for a benchmark set of mols. covering both org. and inorg. systems. Considering in total 96 excited states and using both TD-DFT with a variety of exchange-correlation functionals and the ab initio CIS and CC2 methods, it is found that while the vertical excitation energies obtained with the various methods show an av. (over the 96 states) std. deviation of 0.39 eV, the corresponding std. deviations for the differences between vertical, adiabatic, and 0-0 excitation energies are much smaller: 0.10 (difference between adiabatic and vertical) and 0.02 eV (difference between 0-0 and adiabatic). These results provide a quant. measure showing that the calcn. of such quantities in photochem. modeling is well amenable to low-level methods. In addn., we also report on how these energy differences vary between chem. systems and assess the performance of TD-DFT, CIS, and CC2 in reproducing exptl. 0-0 excitation energies.
- 75Peach, M. J. G.; Benfield, P.; Helgaker, T.; Tozer, D. J. Excitation energies in density functional theory: An evaluation and a diagnostic test. J. Chem. Phys. 2008, 128, 044118, DOI: 10.1063/1.2831900Google Scholar75https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXhslOmsLo%253D&md5=6360f5b16a0fac8535963503c82d96d4Excitation energies in density functional theory: An evaluation and a diagnostic testPeach, Michael J. G.; Benfield, Peter; Helgaker, Trygve; Tozer, David J.Journal of Chemical Physics (2008), 128 (4), 044118/1-044118/8CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Electronic excitation energies are detd. using the CAM-B3LYP Coulomb-attenuated functional, together with a std. generalized gradient approxn. (GGA) and hybrid functional. The degree of spatial overlap between the occupied and virtual orbitals involved in an excitation is measured using a quantity Λ, and the extent to which excitation energy errors correlate with Λ is quantified. For a set of 59 excitations of local, Rydberg, and intramol. charge-transfer character in 18 theor. challenging main-group mols., CAM-B3LYP provides by far the best overall performance; no correlation is obsd. between excitation energy errors and Λ, reflecting the good quality, balanced description of all three categories of excitation. By contrast, a clear correlation is obsd. for the GGA and, to a lesser extent, the hybrid functional, allowing a simple diagnostic test to be proposed for judging the reliability of a general excitation from these functionals-when Λ falls below a prescribed threshold, excitations are likely to be in very significant error. The study highlights the ambiguous nature of the term "charge transfer," providing insight into the observation that while many charge-transfer excitations are poorly described by GGA and hybrid functionals, others are accurately reproduced. (c) 2008 American Institute of Physics.
- 76Send, R.; Kühn, M.; Furche, F. Assessing Excited State Methods by Adiabatic Excitation Energies. J. Chem. Theory Comput. 2011, 7, 2376– 2386, DOI: 10.1021/ct200272bGoogle Scholar76https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXptVOiurk%253D&md5=06f4095c2f002c0081e026db938f9fa1Assessing Excited State Methods by Adiabatic Excitation EnergiesSend, Robert; Kuhn, Michael; Furche, FilippJournal of Chemical Theory and Computation (2011), 7 (8), 2376-2386CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)We compile a 109-membered benchmark set of adiabatic excitation energies (AEEs) from high-resoln. gas-phase expts. Our data set includes a variety of org. chromophores with up to 46 atoms, radicals, and inorg. transition metal compds. Many of the 91 mols. in our set are relevant to atm. chem., photovoltaics, photochem., and biol. The set samples valence, Rydberg, and ionic states of various spin multiplicities. As opposed to vertical excitation energies, AEEs are rigorously defined by energy differences of vibronic states, directly observable, and insensitive to errors in equil. structures. We supply optimized ground state and excited state structures, which allows fast and convenient evaluation of AEEs with two single-point energy calcns. per system. We apply our benchmark set to assess the performance of time-dependent d. functional theory using common semilocal functionals and the CI singles method. Hybrid functionals such as B3LYP and PBE0 yield the best results, with mean abs. errors around 0.3 eV. We also investigate basis set convergence and correlations between different methods and between the magnitude of the excited state relaxation energy and the AEE error. A smaller, 15-membered subset of AEEs is introduced and used to assess the correlated wave function methods CC2 and ADC(2). These methods improve upon hybrid TDDFT for systems with single-ref. ground states but perform less well for radicals and small-gap transition metal compds. None of the investigated methods reaches "chem. accuracy" of 0.05 eV in AEEs.
- 77Côté, A. P.; Benin, A. I.; Ockwig, N. W.; O’Keeffe, M.; Matzger, A. J.; Yaghi, O. M. Porous, Crystalline, Covalent Organic Frameworks. Science 2005, 310, 1166– 1170, DOI: 10.1126/science.1120411Google Scholar77https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2MXht1WgtLnF&md5=e2e033fd4af37870fc73cbb18c1eaed1Porous, Crystalline, Covalent Organic FrameworksCote, Adrien P.; Benin, Annabelle I.; Ockwig, Nathan W.; O'Keeffe, Michael; Matzger, Adam J.; Yaghi, Omar M.Science (Washington, DC, United States) (2005), 310 (5751), 1166-1170CODEN: SCIEAS; ISSN:0036-8075. (American Association for the Advancement of Science)Covalent org. frameworks (COFs) have been designed and successfully synthesized by condensation reactions of Ph diboronic acid {C6H4[B(OH)2]2} and hexahydroxytriphenylene [C18H6(OH)6]. Powder x-ray diffraction studies of the highly cryst. products (C3H2BO)6•(C9H12)1 (COF-1) and C9H4BO2 (COF-5) revealed expanded porous graphitic layers that are either staggered (COF-1, P63/mmc) or eclipsed (COF-5, P6/mmm). Their crystal structures are entirely held by strong bonds between B, C, and O atoms to form rigid porous architectures with pore sizes ranging from 7 to 27 angstroms. COF-1 and COF-5 exhibit high thermal stability (to temps. up to 500° to 600°C), permanent porosity, and high surface areas (711 and 1590 square meters per g, resp.).
- 78Das, G.; Biswal, B. P.; Kandambeth, S.; Venkatesh, V.; Kaur, G.; Addicoat, M.; Heine, T.; Verma, S.; Banerjee, R. Chemical sensing in two dimensional porous covalent organic nanosheets. Chem. Sci. 2015, 6, 3931– 3939, DOI: 10.1039/C5SC00512DGoogle Scholar78https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2MXntlWktLg%253D&md5=f2f12867a8edaafb56789f5184445ce6Chemical sensing in two dimensional porous covalent organic nanosheetsDas, Gobinda; Biswal, Bishnu P.; Kandambeth, Sharath; Venkatesh, V.; Kaur, Gagandeep; Addicoat, Matthew; Heine, Thomas; Verma, Sandeep; Banerjee, RahulChemical Science (2015), 6 (7), 3931-3939CODEN: CSHCCN; ISSN:2041-6520. (Royal Society of Chemistry)Two new imide-based cryst., porous, and chem. stable covalent org. frameworks (COFs) (TpBDH and TfpBDH) have been successfully synthesized employing solvothermal crystn. route. Furthermore, thin layered covalent org. nanosheets (CONs) were derived from these bulk COFs by the simple liq. phase exfoliation method. These 2D CONs showcase increased luminescence intensity compared to their bulk counterparts (COFs). Notably, TfpBDH-CONs showcase good selectivity and prominent, direct visual detection towards different nitroarom. analytes over TpBDH-CONs. Quite interestingly, TfpBDH-CONs exhibit a superior "turn-on" detection capability for 2,4,6-trinitrophenol (TNP) in the solid state, but conversely, they also show a "turn-off" detection in the dispersion state. These findings describe a new approach towards developing an efficient, promising fluorescence chemosensor material for both visual and spectroscopic detection of nitroarom. compds. with very low [10-5 (M)] analyte concns.
- 79Das, P.; Chakraborty, G.; Tyagi, S.; Mandal, S. K. Design of Fluorescent and Robust Covalent Organic Framework Host Matrices for Illuminating Mechanistic Insight into Solvatochromic Decoding. ACS Appl. Mater. Interfaces 2020, 12, 52527– 52537, DOI: 10.1021/acsami.0c12785Google Scholar79https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXit12rsrbE&md5=550f96e9c2920b9eaa975ee6e6fe1045Design of Fluorescent and Robust Covalent Organic Framework Host Matrices for Illuminating Mechanistic Insight into Solvatochromic DecodingDas, Prasenjit; Chakraborty, Gouri; Tyagi, Sparsh; Mandal, Sanjay K.ACS Applied Materials & Interfaces (2020), 12 (47), 52527-52537CODEN: AAMICK; ISSN:1944-8244. (American Chemical Society)Two functional covalent org. frameworks (COFs), constructed from 3-connected triazine-based amine or hydrazine with linear dialdehyde, are decorated with mol. docking sites to showcase their solvatochromic decoding behavior toward volatile solvent mols. (VSMs). These luminescent and cryst. COFs, namely, COF-N and COF-NN, are characterized by numerous anal. techniques. After accommodation of different VSMs as guests, the inclusion compds. of COF-N and COF-NN display solvatochromism. More fascinatingly, the singlet energy, band gaps, and lifetime of these VSM-accommodated COF-N and COF-NN are linearly correlated with the properties of VSMs. D. functional theory (DFT) and Monte Carlo simulation studies further support the interaction of VSMs with COF-N and COF-NN. The presence of extra amine functionality in COF-NN leads to the better interaction with VSMs and, therefore, results in different modes of interaction and correlation. Considering their inestimable chem. diversity, this study introduces a new path for finely tuned solvatochromic properties by COFs.
Cited By
Smart citations by scite.ai include citation statements extracted from the full text of the citing article. The number of the statements may be higher than the number of citations provided by ACS Publications if one paper cites another multiple times or lower if scite has not yet processed some of the citing articles.
This article is cited by 23 publications.
- Huiling Ji, Linlin Sun, Jigang Wang, Nan Jiang. Van der Waals Effect-Driven Charge–Transfer Excitation in 2D Polymeric Carbon Nitrides: A Vertical Excitation Model. The Journal of Physical Chemistry C 2025, 129
(2)
, 1282-1292. https://doi.org/10.1021/acs.jpcc.4c06175
- Xin He, Jiazheng Diao, Qi Yao, Hui Wang, Xia Zhong, Wen Si, Qiang Fu, Long Yuan, Xiaodong Zhang, Yi Xie. Activating Intrinsic Self-Trapped Exciton Emission in Bismuth Oxyhalides by Edge Iodine Doping. ACS Photonics 2025, 12
(1)
, 384-391. https://doi.org/10.1021/acsphotonics.4c01770
- Wen-Kai Chen, Xi Zhao, Xiang-Yang Liu, Xiao-Ying Xie, Yanli Zeng, Ganglong Cui. Photoinduced Nonadiabatic Dynamics of a Single-Walled Carbon Nanotube-Porphyrin Complex. The Journal of Physical Chemistry A 2024, 128
(40)
, 8803-8815. https://doi.org/10.1021/acs.jpca.4c04544
- Beliz Sertcan Gökmen, Jürg Hutter, Anna-Sophia Hehn. Excited-State Forces with the Gaussian and Augmented Plane Wave Method for the Tamm–Dancoff Approximation of Time-Dependent Density Functional Theory. Journal of Chemical Theory and Computation 2024, 20
(19)
, 8494-8504. https://doi.org/10.1021/acs.jctc.4c00614
- Xiaojuan Ni, Hong Li, Jean-Luc Brédas. Enhanced Organic–Inorganic Electronic Coupling in Two-Dimensional Hybrid Perovskites through Molecular Engineering of Dipolar Pyrene-Based Cations. ACS Materials Letters 2024, 6
(8)
, 3436-3442. https://doi.org/10.1021/acsmaterialslett.4c00995
- Eva Vandaele, Momir Mališ, Sandra Luber. The Role of Aqueous Solvation on the Intersystem Crossing of Nitrophenols. Journal of Chemical Theory and Computation 2024, 20
(8)
, 3258-3272. https://doi.org/10.1021/acs.jctc.3c01400
- Marco Caricato. A Perspective on the Simulation of Electronic Circular Dichroism and Circularly Polarized Luminescence Spectra in Chiral Solid Materials. The Journal of Physical Chemistry A 2024, 128
(7)
, 1197-1206. https://doi.org/10.1021/acs.jpca.3c08095
- Lili Rassouli, Michel Dupuis. Electronic Structure of Excitons in Hematite Fe2O3. The Journal of Physical Chemistry C 2024, 128
(2)
, 743-758. https://doi.org/10.1021/acs.jpcc.3c06808
- Hung H. Nguyen, Vyacheslav S. Bryantsev, Claudio J. Margulis. Are High-Temperature Molten Salts Reactive with Excess Electrons? Case of ZnCl2. The Journal of Physical Chemistry B 2023, 127
(42)
, 9155-9164. https://doi.org/10.1021/acs.jpcb.3c04210
- Sheng-Rui Wang, Qiu Fang, Xiang-Yang Liu, Wei-Hai Fang, Ganglong Cui. Machine learning accelerated nonadiabatic dynamics simulations of materials with excitonic effects. The Journal of Chemical Physics 2025, 162
(2)
https://doi.org/10.1063/5.0248228
- Jinggang Lan, Majed Chergui, Alfredo Pasquarello. Dynamics of the charge transfer to solvent process in aqueous iodide. Nature Communications 2024, 15
(1)
https://doi.org/10.1038/s41467-024-46772-0
- Stefano Battaglia, Max Rossmannek, Vladimir V. Rybkin, Ivano Tavernelli, Jürg Hutter. A general framework for active space embedding methods with applications in quantum computing. npj Computational Materials 2024, 10
(1)
https://doi.org/10.1038/s41524-024-01477-2
- Andrea Cocco, Sara Paniziutti, Chiara Olla, Riccardo Corpino, Carlo Maria Carbonaro, Pier Carlo Ricci, Nicola Melis, Paola Caria, Giuseppina Sanna, Eli Zysman‐Colman, Francesco Secci. Design, Synthesis, and Photophysical Characterization of Biocompatible Thermally Activated Delayed Fluorescent Carbazole‐Coumarins for Sensing Applications. Chemistry – A European Journal 2024, 30
(55)
https://doi.org/10.1002/chem.202401263
- Sebastian A. Suarez. The master key: structural science in unlocking functional materials advancements. Journal of Applied Crystallography 2024, 57
(3)
, 606-622. https://doi.org/10.1107/S1600576724003674
- Marc de Wergifosse, Stefan Grimme. The eXact integral simplified time-dependent density functional theory (XsTD-DFT). The Journal of Chemical Physics 2024, 160
(20)
https://doi.org/10.1063/5.0206380
- Raphaël Rullan, Pauline Colinet, Quentin Desdion, Stephan N. Steinmann, Tangui Le Bahers. Modeling the polychromism of oxide minerals: The case of alexandrite and cordierite. Journal of Computational Chemistry 2024, 45
(12)
, 834-842. https://doi.org/10.1002/jcc.27288
- Michael Ingham, Alex Aziz, Devis Di Tommaso, Rachel Crespo-Otero. Simulating excited states in metal organic frameworks: from light-absorption to photochemical CO
2
reduction. Materials Advances 2023, 4
(22)
, 5388-5419. https://doi.org/10.1039/D3MA00518F
- Jun‐Hua Wei, Jian‐Bin Luo, Zi‐Lin He, Zhi‐Zhong Zhang, Dai‐Bin Kuang. A New Class of Donor‐
π
‐Acceptor Phosphonium Salt Exhibiting Intriguing Thermally Activated Delayed Fluorescence. Advanced Optical Materials 2023, 11
(17)
https://doi.org/10.1002/adom.202300328
- Vikram Gavini, Stefano Baroni, Volker Blum, David R Bowler, Alexander Buccheri, James R Chelikowsky, Sambit Das, William Dawson, Pietro Delugas, Mehmet Dogan, Claudia Draxl, Giulia Galli, Luigi Genovese, Paolo Giannozzi, Matteo Giantomassi, Xavier Gonze, Marco Govoni, François Gygi, Andris Gulans, John M Herbert, Sebastian Kokott, Thomas D Kühne, Kai-Hsin Liou, Tsuyoshi Miyazaki, Phani Motamarri, Ayako Nakata, John E Pask, Christian Plessl, Laura E Ratcliff, Ryan M Richard, Mariana Rossi, Robert Schade, Matthias Scheffler, Ole Schütt, Phanish Suryanarayana, Marc Torrent, Lionel Truflandier, Theresa L Windus, Qimen Xu, Victor W-Z Yu, D Perez. Roadmap on electronic structure codes in the exascale era. Modelling and Simulation in Materials Science and Engineering 2023, 31
(6)
, 063301. https://doi.org/10.1088/1361-651X/acdf06
- Fabian Belleflamme, Jürg Hutter. Radicals in aqueous solution: assessment of density-corrected SCAN functional. Physical Chemistry Chemical Physics 2023, 25
(31)
, 20817-20836. https://doi.org/10.1039/D3CP02517A
- Jhonathan R. Souza, Carles Curutchet, Yuri A. Aoto, Paula Homem-De-Mello. Benchmarking DFT functionals for photophysics of pyranoflavylium cations. Journal of Molecular Graphics and Modelling 2023, 122 , 108460. https://doi.org/10.1016/j.jmgm.2023.108460
- Manasi R. Mulay, Natalia Martsinovich. Interaction of organic pollutants with TiO
2
: a density functional theory study of carboxylic acids on the anatase (101) surface. Molecular Physics 2023, 121
(7-8)
https://doi.org/10.1080/00268976.2023.2165981
- Haibei Huang, Lijuan Xue, Gang Lu, Shibo Cheng, Yuxiang Bu. Hydrated electrons as nodes in porous clathrate hydrates. The Journal of Chemical Physics 2023, 158
(11)
https://doi.org/10.1063/5.0135335
Article Views are the COUNTER-compliant sum of full text article downloads since November 2008 (both PDF and HTML) across all institutions and individuals. These metrics are regularly updated to reflect usage leading up to the last few days.
Citations are the number of other articles citing this article, calculated by Crossref and updated daily. Find more information about Crossref citation counts.
The Altmetric Attention Score is a quantitative measure of the attention that a research article has received online. Clicking on the donut icon will load a page at altmetric.com with additional details about the score and the social media presence for the given article. Find more information on the Altmetric Attention Score and how the score is calculated.
Recommended Articles
Abstract
Figure 1
Figure 1. Error (in pm) for selected bond lengths for (a) optimized ground- and (b) first singlet excited-state geometries comparing conventional PBE0 and ADMM-PBE0 computations for different auxiliary (ABS) and primary basis sets (PBS) indicated as ABS/PBS.
Figure 2
Figure 2. Error (in deg) for selected angles for optimized first singlet excited-state geometries comparing conventional PBE0 and ADMM-PBE0 computations.
Figure 3
Figure 3. Normal distributions based on the mean errors (ME) and standard deviations (STD) wrt EOM-CCSD reference geometries in (pm)/(deg)/(deg) for selected (a) bond lengths/(b) angles/(c) dihedral angles for optimized first singlet excited-state geometries comparing conventional PBE0, ADMM-PBE0, and sTDA computations using triple-ζ basis sets.
Figure 4
Figure 4. Normal distributions based on the MEs and STDs for vertical excitation energies (in eV) depicting (a) the error introduced by ADMM-PBE0 in comparison to conventional PBE0 computations (upper left plot), (b) the performance of PBE0, ADMM-PBE0, and sTDA kernel computations using double- and triple-ζ basis sets in comparison to a PBE0/TZV2P/CP2K reference (upper right plot) as well as (c) an analogous assessment of the kernels with respect to a PBE0/def2-QZVPP/Turbomole reference (lower plot).
Figure 5
Figure 5. ADMM error distribution in (a) adiabatic excitation and (b) fluorescence energies comparing approximated ADMM-PBE0 and conventional PBE0 results for different auxiliary basis set sizes (in eV).
Figure 6
Figure 6. Comparison of (a) adiabatic excitation energies with respect to EOM-CCSD reference data and (b) fluorescence energies with respect to PBE0/def2-QZVPP/Turbomole reference data for PBE0, ADMM-PBE0, and sTDA kernels (in eV).
Figure 7
Figure 7. Comparison of ADMM-PBE0 and sTDA absorption spectra for a series of fluorescent COFs.
References
This article references 79 other publications.
- 1González, L.; Escudero, D.; Serrano-Andrés, L. Progress and Challenges in the Calculation of Electronic Excited States. ChemPhysChem 2012, 13, 28– 51, DOI: 10.1002/cphc.2011002001https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC38XjvVGntw%253D%253D&md5=710b04d12866b33b409e5377a8d3be72Progress and Challenges in the Calculation of Electronic Excited StatesGonzalez, Leticia; Escudero, Daniel; Serrano-Andres, LuisChemPhysChem (2012), 13 (1), 28-51CODEN: CPCHFT; ISSN:1439-4235. (Wiley-VCH Verlag GmbH & Co. KGaA)A review. A detailed understanding of the properties of electronic excited states and the reaction mechanisms that mols. undergo after light irradn. is a fundamental ingredient for following light-driven natural processes and for designing novel photonic materials. The aim of this review is to present an overview of the ab initio quantum chem. and time-dependent d. functional theory methods that can be used to model spectroscopy and photochem. in mol. systems. The applicability and limitations of the different methods as well as the main frontiers are discussed. To illustrate the progress achieved by excited-state chem. in the recent years as well as the main challenges facing computational chem., three main applications that reflect the authors' experience are addressed: the UV/Vis spectroscopy of org. mols., the assignment of absorption and emission bands of organometallic complexes, and finally, the obtainment of non-adiabatic photoinduced pathways mediated by conical intersections. In the latter case, special emphasis is put on the photochem. of DNA. These applications show that the description of electronically excited states is a rewarding but challenging area of research.
- 2Dreuw, A.; Head-Gordon, M. Single-Reference ab Initio Methods for the Calculation of Excited States of Large Molecules. Chem. Rev. 2005, 105, 4009– 4037, DOI: 10.1021/cr05056272https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2MXhtVKmtLrK&md5=31bf0c10a174b6d5f2aa914bb0118957Single-Reference ab Initio Methods for the Calculation of Excited States of Large MoleculesDreuw, Andreas; Head-Gordon, MartinChemical Reviews (Washington, DC, United States) (2005), 105 (11), 4009-4037CODEN: CHREAY; ISSN:0009-2665. (American Chemical Society)A review of single-ref. ab initio excited state methods, which are applicable to large mols. and do not explicitly include correlation through the ground-state wave function. CIS, TDHF, and TDDFT are rigorously introduced by outlining their derivations and theor. footing, where special emphasis is put on their relations to each other. Different methods for the anal. of complicated electronically excited states are reviewed and comparatively discussed. The applicability of the presented methods is outlined, their limitations are show, and out their strengths and weaknesses are pointed out.
- 3Laurent, A. D.; Jacquemin, D. TD-DFT benchmarks: A review. Int. J. Quantum Chem. 2013, 113, 2019– 2039, DOI: 10.1002/qua.244383https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3sXlsVGgt74%253D&md5=2d81e67d518499c1a10f72ec369ec6aeTD-DFT benchmarks: A reviewLaurent, Adele D.; Jacquemin, DenisInternational Journal of Quantum Chemistry (2013), 113 (17), 2019-2039CODEN: IJQCB2; ISSN:0020-7608. (John Wiley & Sons, Inc.)A review. Time-Dependent D. Functional Theory (TD-DFT) has become the most widely-used theor. approach to simulate the optical properties of both org. and inorg. mols. In this contribution, we review TD-DFT benchmarks that have been performed during the last decade. The aim is often to pinpoint the most accurate or adequate exchange-correlation functional(s). We present both the different strategies used to assess the functionals and the main results obtained in terms of accuracy. In particular, we discuss both vertical and adiabatic benchmarks and comparisons with both exptl. and theor. ref. transition energies. More specific benchmarks (oscillator strengths, excited-state geometries, dipole moments, vibronic shapes, etc.) are summarized as well. © 2013 Wiley Periodicals, Inc.
- 4de Wergifosse, M.; Grimme, S. Perspective on Simplified Quantum Chemistry Methods for Excited States and Response Properties. J. Phys. Chem. A 2021, 125, 3841– 3851, DOI: 10.1021/acs.jpca.1c023624https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3MXpsl2ltLo%253D&md5=279421848df210a3de326d8bd50e783cPerspective on Simplified Quantum Chemistry Methods for Excited States and Response PropertiesDe wergifosse, Marc; Grimme, StefanJournal of Physical Chemistry A (2021), 125 (18), 3841-3851CODEN: JPCAFH; ISSN:1089-5639. (American Chemical Society)A review. We review recent developments in the framework of simplified quantum chem. for excited state and optical response properties (sTD-DFT) and present future challenges for new method developments to improve accuracy and extend the range of application. In recent years, the scope of sTD-DFT was extended to mol. response calcns. of the polarizability, optical rotation, first hyperpolarizability, two-photon absorption (2PA), and excited-state absorption for large systems with hundreds to thousands of atoms. The recently introduced spin-flip simplified time-dependent d. functional theory (SF-sTD-DFT) variant enables an ultrafast treatment for diradicals and related strongly correlated systems. A few drawbacks were also identified, specifically for the computation of 2PA cross sections. We propose solns. to this problem and how to generally improve the accuracy of simplified schemes. New possible simplified schemes are also introduced for strongly correlated systems, e.g., with a second-order perturbative correlation correction. Interpretation tools that can ext. chem. structure-property relationships from excited state or response calcns. are also discussed. In particular, the recently introduced method-agnostic RespA approach based on natural response orbitals (NROs) as the key concept is employed.
- 5Grimme, S.; Schreiner, P. R. Computational Chemistry: The Fate of Current Methods and Future Challenges. Angew. Chem., Int. Ed. 2018, 57, 4170– 4176, DOI: 10.1002/anie.2017099435https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhvFGitrzO&md5=4890cef4031e2238fd8db2f01049cbedComputational Chemistry: The Fate of Current Methods and Future ChallengesGrimme, Stefan; Schreiner, Peter R.Angewandte Chemie, International Edition (2018), 57 (16), 4170-4176CODEN: ACIEF5; ISSN:1433-7851. (Wiley-VCH Verlag GmbH & Co. KGaA)In this essay, we attempt to make predictions about the fate and development of the computational mol. sciences. Of course, it is not the first time that the future challenges for computational org. chem. and biochem. are considered; these were outlined in complementary contexts recently. In this article, the authors take a somewhat different perspective and emphasize the changes expected for chem. that are triggered by the rapid developments and increasingly stronger influences from theory, algorithms, and data-driven technologies. The authors of this essay are about the same age and have a general overview of a period of about 25 years during which they have actively contributed to the field of computational chem. Hence, it appears sensible to make predictions extending 25 years into the future, roughly to the year 2043 (when both authors will long be retired).
- 6Dierksen, M.; Grimme, S. The Vibronic Structure of Electronic Absorption Spectra of Large Molecules: A Time-Dependent Density Functional Study on the Influence of Exact Hartree-Fock Exchange. J. Phys. Chem. A 2004, 108, 10225– 10237, DOI: 10.1021/jp047289h6https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2cXotlKktLk%253D&md5=b85a6bb8357480b91217f853af23c6a3The Vibronic Structure of Electronic Absorption Spectra of Large Molecules: A Time-Dependent Density Functional Study on the Influence of "Exact" Hartree-Fock ExchangeDierksen, Marc; Grimme, StefanJournal of Physical Chemistry A (2004), 108 (46), 10225-10237CODEN: JPCAFH; ISSN:1089-5639. (American Chemical Society)The functional dependence of excited-state geometries and normal modes calcd. with time-dependent d. functional theory (TDDFT) is studied from vibronic structure calcns. of the absorption spectra of large mols. For a set of mols. covering a wide range of different structures including org. dyes, biol. chromophores, and mols. of importance in material science, quantum mech. simulations of the vibronic structure are performed. In total over 40 singlet-singlet transitions of neutral closed-shell compds. and doublet-doublet transitions of neutral radicals, radical cations, and anions are considered. Calcns. with different std. d. functionals show that the predicted vibronic structure critically depends on the fraction of the exact Hartree-Fock exchange (EEX) included in hybrid functionals. The effect can been traced back to a large influence of EEX on the geometrical displacement upon excitation. On the contrary, the dependence of the results on the choice of the local exchange-correlation functional is rather small. From detailed comparisons with exptl. spectra conclusions are drawn concerning the optimum amt. of EEX mixing for a proper description of the excited-state properties. The relation of the quality of the simulated spectra with the errors for 0-0 transition energies is discussed. For the studied singlet-singlet π → π* transitions and the 1st strongly dipole-allowed transitions of PAH radical cations some rules of thumb concerning the optimum portion of EEX are derived. However, in general no universal amt. of EEX seems to exist that gives a uniformly good description for all systems and states. Nevertheless an inclusion of ∼30-40% of EEX in the functional is found empirically to yield in most cases simulated spectra that compare very well with those from expt. and thus seems to be necessary for an accurate description of the excited-state geometry. Pure d. functionals that are computationally more efficient provide less accurate spectra in most cases and their application is recommended solely for comparison purposes to obtain ests. for the reliability of the theor. predictions.
- 7Merlot, P.; Izsák, R.; Borgoo, A.; Kjærgaard, T.; Helgaker, T.; Reine, S. Charge-constrained auxiliary-density-matrix methods for the Hartree-Fock exchange contribution. J. Chem. Phys. 2014, 141, 094104, DOI: 10.1063/1.48942677https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXhsFSjtrzP&md5=0d138a12fe1e543edb43a3e92af05f10Charge-constrained auxiliary-density-matrix methods for the Hartree-Fock exchange contributionMerlot, Patrick; Izsak, Robert; Borgoo, Alex; Kjaergaard, Thomas; Helgaker, Trygve; Reine, SimenJournal of Chemical Physics (2014), 141 (9), 094104/1-094104/11CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Three new variants of the auxiliary-d.-matrix method (ADMM) of Guidon, Hutter, and VandeVondele [J. Chem. Theory Comput. 6, 2348 (2010)] are presented with the common feature that they have a simplified constraint compared with the full orthonormality requirement of the earlier ADMM1 method. All ADMM variants are tested for accuracy and performance in all-electron B3LYP calcns. with several commonly used basis sets. The effect of the choice of the exchange functional for the ADMM exchange-correction term is also investigated. (c) 2014 American Institute of Physics.
- 8Guidon, M.; Hutter, J.; VandeVondele, J. Auxiliary Density Matrix Methods for Hartree-Fock Exchange Calculations. J. Chem. Theory Comput. 2010, 6, 2348– 2364, DOI: 10.1021/ct10022258https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXot1als7Y%253D&md5=7a84c086d49f17bab239fd3f8d004c6eAuxiliary Density Matrix Methods for Hartree-Fock Exchange CalculationsGuidon, Manuel; Hutter, Jurg; Vande Vondele, JoostJournal of Chemical Theory and Computation (2010), 6 (8), 2348-2364CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)The calcn. of Hartree-Fock exchange (HFX) is computationally demanding for large systems described with high-quality basis sets. In this work, we show that excellent performance and good accuracy can nevertheless be obtained if an auxiliary d. matrix is employed for the HFX calcn. Several schemes to derive an auxiliary d. matrix from a high-quality d. matrix are discussed. Key to the accuracy of the auxiliary d. matrix methods (ADMM) is the use of a correction based on std. generalized gradient approxns. for HFX. ADMM integrates seamlessly in existing HFX codes and, in particular, can be employed in linear scaling implementations. Demonstrating the performance of the method, the effect of HFX on the structure of liq. water is investigated in detail using Born-Oppenheimer mol. dynamics simulations (300 ps) of a system of 64 mols. Representative for large systems are calcns. on a solvated protein (Rubredoxin), for which ADMM outperforms the corresponding std. HFX implementation by approx. a factor 20.
- 9Kumar, C.; Fliegl, H.; Jensen, F.; Teale, A. M.; Reine, S.; Kjærgaard, T. Accelerating Kohn-Sham response theory using density fitting and the auxiliary-density-matrix method. Int. J. Quantum Chem. 2018, 118, e25639 DOI: 10.1002/qua.25639There is no corresponding record for this reference.
- 10Rebolini, E.; Izsák, R.; Reine, S. S.; Helgaker, T.; Pedersen, T. B. Comparison of Three Efficient Approximate Exact-Exchange Algorithms: The Chain-of-Spheres Algorithm, Pair-Atomic Resolution-of-the-Identity Method, and Auxiliary Density Matrix Method. J. Chem. Theory Comput. 2016, 12, 3514– 3522, DOI: 10.1021/acs.jctc.6b0007410https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28Xos1Sku7s%253D&md5=dccdeec570871ad4ed29846f8f6d70a4Comparison of Three Efficient Approximate Exact-Exchange Algorithms: The Chain-of-Spheres Algorithm, Pair-Atomic Resolution-of-the-Identity Method, and Auxiliary Density Matrix MethodRebolini, Elisa; Izsak, Robert; Reine, Simen Sommerfelt; Helgaker, Trygve; Pedersen, Thomas BondoJournal of Chemical Theory and Computation (2016), 12 (8), 3514-3522CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)We compare the performance of three approx. methods for speeding up evaluation of the exchange contribution in Hartree-Fock and hybrid Kohn-Sham calcns.: the chain-of-spheres algorithm, and the auxiliary d. matrix method. Both the efficiency relative to that of a conventional linear-scaling algorithm and the accuracy of total, atomization, and orbital energies are compared for a subset contg. 25 of the 200 mols. in the Rx200 set using double-, triple-, and quadruple-ζ basis sets. The accuracy of relative energies is further compared for small alkane conformers (ACONF test set) and Diels-Alder reactions (DARC test set). Overall, we find that the COSX method provides good accuracy for orbital energies as well as total and relative energies, and the method delivers a satisfactory speedup. The PARI-K and in particular ADMM algorithms require further development and optimization to fully exploit their indisputable potential.
- 11Merlot, P.; Kjærgaard, T.; Helgaker, T.; Lindh, R.; Aquilante, F.; Reine, S.; Pedersen, T. B. Attractive electron–electron interactions within robust local fitting approximations. J. Comput. Chem. 2013, 34, 1486– 1496, DOI: 10.1002/jcc.2328411https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3sXltlOntrg%253D&md5=d51f5c5ca0807fb8dd24244873b2a0d2Attractive electron-electron interactions within robust local fitting approximationsMerlot, Patrick; Kjaergaard, Thomas; Helgaker, Trygve; Lindh, Roland; Aquilante, Francesco; Reine, Simen; Pedersen, Thomas BondoJournal of Computational Chemistry (2013), 34 (17), 1486-1496CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)An anal. of Dunlap's robust fitting approach reveals that the resulting two-electron integral matrix is not manifestly pos. semidefinite when local fitting domains or non-Coulomb fitting metrics are used. We present a highly local approx. method for evaluating four-center two-electron integrals based on the resoln.-of-the-identity (RI) approxn. and apply it to the construction of the Coulomb and exchange contributions to the Fock matrix. In this pair-at. resoln.-of-the-identity (PARI) approach, at.-orbital (AO) products are expanded in auxiliary functions centered on the two atoms assocd. with each product. Numerical tests indicate that in 1% or less of all Hartree-Fock and Kohn-Sham calcns., the indefinite integral matrix causes nonconvergence in the self-consistent-field iterations. In these cases, the two-electron contribution to the total energy becomes neg., meaning that the electronic interaction is effectively attractive, and the total energy is dramatically lower than that obtained with exact integrals. In the vast majority of our test cases, however, the indefiniteness does not interfere with convergence. The total energy accuracy is comparable to that of the std. Coulomb-metric RI method. The speed-up compared with conventional algorithms is similar to the RI method for Coulomb contributions; exchange contributions are accelerated by a factor of up to eight with a triple-zeta quality basis set. A pos. semidefinite integral matrix is recovered within PARI by introducing local auxiliary basis functions spanning the full AO product space, as may be achieved by using Cholesky-decompn. techniques. Local completion, however, slows down the algorithm to a level comparable with or below conventional calcns. © 2013 Wiley Periodicals, Inc.
- 12Manzer, S. F.; Epifanovsky, E.; Head-Gordon, M. Efficient Implementation of the Pair Atomic Resolution of the Identity Approximation for Exact Exchange for Hybrid and Range-Separated Density Functionals. J. Chem. Theory Comput. 2015, 11, 518– 527, DOI: 10.1021/ct500858612https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXitFehtrrN&md5=bc28289b8fa1caf00b97bdda95c12c88Efficient Implementation of the Pair Atomic Resolution of the Identity Approximation for Exact Exchange for Hybrid and Range-Separated Density FunctionalsManzer, Samuel F.; Epifanovsky, Evgeny; Head-Gordon, MartinJournal of Chemical Theory and Computation (2015), 11 (2), 518-527CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)An efficient new MO basis algorithm is reported implementing the pair at. resoln. of the identity approxn. (PARI) to evaluate the exact exchange contribution (K) to SCF methods, such as hybrid and range-sepd. hybrid d. functionals. The PARI approxn., in which AO basis function pairs are expanded using auxiliary basis functions centered only on their two resp. atoms, was recently investigated by Merlot et al. Our algorithm is significantly faster than quartic scaling RI-K, with an asymptotic exchange speedup for hybrid functionals of (1 + X/N), where N and X are the AO and auxiliary basis dimensions. The asymptotic speedup is 2 + 2X/N for range sepd. hybrids such as CAM-B3LYP, ωB97X-D, and ωB97X-V which include short- and long-range exact exchange. The obsd. speedup for exchange in ωB97X-V for a C68 graphene fragment in the cc-pVTZ basis is 3.4 relative to RI-K. Like conventional RI-K, our method greatly outperforms conventional integral evaluation in large basis sets; a speedup of 19 is obtained in the cc-pVQZ basis on a C54 graphene fragment. Negligible loss of accuracy relative to exact integral evaluation is demonstrated on databases of bonded and nonbonded interactions. We also demonstrate both anal. and numerically that the PARI-K approxn. is variationally stable.
- 13Neese, F.; Wennmohs, F.; Hansen, A.; Becker, U. Efficient, approximate and parallel Hartree–Fock and hybrid DFT calculations. A ‘chain-of-spheres’ algorithm for the Hartree–Fock exchange. Chem. Phys. 2009, 356, 98– 109, DOI: 10.1016/j.chemphys.2008.10.03613https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1MXhs1eru7g%253D&md5=0a7b02cc3eeb4cad044b9b726b8c64ddEfficient, approximate and parallel Hartree-Fock and hybrid DFT calculations. A 'chain-of-spheres' algorithm for the Hartree-Fock exchangeNeese, Frank; Wennmohs, Frank; Hansen, Andreas; Becker, UteChemical Physics (2009), 356 (1-3), 98-109CODEN: CMPHC2; ISSN:0301-0104. (Elsevier B.V.)In this paper, the possibility is explored to speed up Hartree-Fock and hybrid d. functional calcns. by forming the Coulomb and exchange parts of the Fock matrix by different approxns. For the Coulomb part the previously introduced Split-RI-J variant of the well-known d. fitting' approxn. is used. The exchange part is formed by semi-numerical integration techniques that are closely related to Friesner's pioneering pseudo-spectral approach. Our potentially linear scaling realization of this algorithm is called the 'chain-of-spheres exchange' (COSX). A combination of semi-numerical integration and d. fitting is also proposed. Both Split-RI-J and COSX scale very well with the highest angular momentum in the basis sets. It is shown that for extended basis sets speed-ups of up to two orders of magnitude compared to traditional implementations can be obtained in this way. Total energies are reproduced with an av. error of <0.3 kcal/mol as detd. from extended test calcns. with various basis sets on a set of 26 mols. with 20-200 atoms and up to 2000 basis functions. Reaction energies agree to within 0.2 kcal/mol (Hartree-Fock) or 0.05 kcal/mol (hybrid DFT) with the canonical values. The COSX algorithm parallelizes with a speedup of 8.6 obsd. for 10 processes. Min. energy geometries differ by less than 0.3 pm in the bond distances and 0.5° in the bond angles from their canonical values. These developments enable highly efficient and accurate SCF calcns. including nonlocal Hartree-Fock exchange for large mols. In combination with the RI-MP2 method and large basis sets, second-order many body perturbation energies can be obtained for medium sized mols. with unprecedented efficiency. The algorithms are implemented into the ORCA electronic structure system.
- 14Izsák, R.; Neese, F. An overlap fitted chain of spheres exchange method. J. Chem. Phys. 2011, 135, 144105, DOI: 10.1063/1.364692114https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXht12hsr3P&md5=9c8a746be8f60e71deff75d4029477f9An overlap fitted chain of spheres exchange methodIzsak, Robert; Neese, FrankJournal of Chemical Physics (2011), 135 (14), 144105/1-144105/11CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The "chain of spheres" (COS) algorithm, as part of the RIJCOSX SCF procedure, approximates the exchange term by performing analytic integration with respect to the coordinates of only one of the two electrons, whereas for the remaining coordinates, integration is carried out numerically. In the present work, we attempt to enhance the efficiency of the method by minimizing numerical errors in the COS procedure. The main idea is based on the work of Friesner and consists of finding a fitting matrix, Q, which leads the numerical and anal. evaluated overlap matrixes to coincide. Using Q, the evaluation of exchange integrals can indeed be improved. Improved results and timings are obtained with the present default grid setup for both single point calcns. and geometry optimizations. The fitting procedure results in a redn. of grid sizes necessary for achieving chem. accuracy. We demonstrate this by testing a no. of grids and comparing results to the fully analytic and the earlier COS approxns. This turns out to be favorable for total and reaction energies, for which chem. accuracy can now be reached with a corresponding ∼30% speedup over the original RIJCOSX procedure for single point energies. Results are slightly less favorable for the accuracy of geometry optimizations, but the procedure is still shown to yield geometries with errors well below the method inherent errors of the employed theor. framework. (c) 2011 American Institute of Physics.
- 15Laqua, H.; Thompson, T. H.; Kussmann, J.; Ochsenfeld, C. Highly Efficient, Linear-Scaling Seminumerical Exact-Exchange Method for Graphic Processing Units. J. Chem. Theory Comput. 2020, 16, 1456– 1468, DOI: 10.1021/acs.jctc.9b0086015https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXivVejs7Y%253D&md5=606be25e20a70cc545ff71ebabf28152Highly Efficient, Linear-Scaling Seminumerical Exact-Exchange Method for Graphic Processing UnitsLaqua, Henryk; Thompson, Travis H.; Kussmann, Joerg; Ochsenfeld, ChristianJournal of Chemical Theory and Computation (2020), 16 (3), 1456-1468CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)We present a highly efficient and asymptotically linear-scaling graphic processing unit accelerated seminumerical exact-exchange method (sn-LinK). We go beyond our previous central processing unit-based method (H. Laqua et al., 2018) by employing our recently developed integral bounds (T.H. Thomson and C. Ochsenfeld, 2019) and high-accuracy numerical integration grid (H. Laqua et al., 2018). The accuracy is assessed for several established test sets, providing errors significantly below 1mEh for the smallest grid. Moreover, a comprehensive performance anal. for large mols. between 62 and 1347 atoms is provided, revealing the outstanding performance of our method, in particular, for large basis sets such as the polarized quadruple-zeta level with diffuse functions.
- 16Holzer, C. An improved seminumerical Coulomb and exchange algorithm for properties and excited states in modern density functional theory. J. Chem. Phys. 2020, 153, 184115, DOI: 10.1063/5.002275516https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXitlelur7N&md5=cad44a4ffbef63e25cf320bbcd316688An improved seminumerical Coulomb and exchange algorithm for properties and excited states in modern density functional theoryHolzer, ChristofJournal of Chemical Physics (2020), 153 (18), 184115CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)A seminumerical algorithm capable of performing large-scale (time-dependent) d. functional theory (TD-DFT) calcns. to ext. excitation energies and other ground-state and excited-state properties is outlined. The algorithm uses seminumerical integral techniques for evaluating Coulomb and exchange parts for a set of d. matrixes as occurring in std. TD-DFT or similar methods for the evaluation of vibrational frequencies. A suitable optimized de-aliasing procedure is introduced. The latter does not depend on further auxiliary quantities and retains the symmetry of a given d. matrix. The algorithm is self-contained and applicable to any orbital basis set available without the need for further auxiliary basis sets or optimized de-aliasing grids. Relativistic two-component excited-state TD-DFT calcns. are reported for the first time using the developed seminumerical algorithm for std. and local hybrid d. functional approxns. Errors are compared with the widely used "resoln. of the identity" (RI) approxns. for Coulomb (RI-J) and exchange integrals (RI-K). The fully seminumerical algorithm does not exhibit an enlarged error for std. DFT functionals compared to the RI approxn. For the more involved local hybrid functionals and within strong external fields, accuracy is even considerably improved. (c) 2020 American Institute of Physics.
- 17Ohno, K. Some remarks on the Pariser-Parr-Pople method. Theor. Chim. Acta 1964, 2, 219– 227, DOI: 10.1007/BF0052828117https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaF2cXkt12htLk%253D&md5=b652ff406687369fbbdb1d53b19870e4The Pariser-Parr-Pople methodOhno, KimioTheoretica Chimica Acta (1964), 2 (3), 219-27CODEN: TCHAAM; ISSN:0040-5744.Basic assumptions which characterize the Pariser-Parr (CA 49, 10725f)-Pople (CA 48, 8639g) method of computing mol. electronic wave functions are examd. crit. By restricted variational calcn. of the valence state of C and N atoms and ions, it is demonstrated that the usual methods of evaluation of center Coulomb integrals and at. core energies are rather good. A semitheoretical means of estg. the core resonance integral is proposed and shown to give fair agreement with the empirical values for C-C, O-O, C-N, and C-O bonds.
- 18Klopman, G. A Semiempirical Treatment of molecular Structures. II. Molecular Terms and Application to diatomic Molecules. J. Am. Chem. Soc. 1964, 86, 4550– 4557, DOI: 10.1021/ja01075a00818https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaF2cXkvVygsr4%253D&md5=be4431018f6ae6937415a7bdf8370bc6A semiempirical treatment of molecular structures. II. Molecular terms and application to diatomic moleculesKlopman, G.Journal of the American Chemical Society (1964), 86 (21), 4550-7CODEN: JACSAT; ISSN:0002-7863.cf. CA 60, 12679e. A self-consistent semiempirical method which is designed for the calcn. of heats of formation and charge distribution of nonconjugated mols. is outlined. The method is based on an antisymmetrized product of mol. orbitals, simplified in such a way as to make the values of all involved integrals directly available from at. spectra (loc. cit.) and mol. bond distances. In a preliminary study, this method was used to calc. satisfactory values of bond energies and reasonable values of charge distributions in 80 diat. mols. (σ-bonded).
- 19Nishimoto, K.; Mataga, N. Z. Electronic Structure and Spectra of Some Nitrogen Heterocycles. Z. Phys. Chem. 1957, 12, 335, DOI: 10.1524/zpch.1957.12.5_6.33519https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaG2sXps1ymsg%253D%253D&md5=95bf1fdd3345cb66b0f8af10bf781f0bElectronic structure and spectra of some nitrogen heretocyclesNishimoto, Kitisuke; Mataga, NoboruZeitschrift fuer Physikalische Chemie (Muenchen, Germany) (1957), 12 (), 335-8CODEN: ZPCFAX; ISSN:0044-3336.The energies of lower excited states of pyridine, pyrazine, and sym-triazine and the oscillator strengths of the transitions to those excited states have been calcd. The results are listed and discussed.
- 20Bannwarth, C.; Grimme, S. A simplified time-dependent density functional theory approach for electronic ultraviolet and circular dichroism spectra of very large molecules. Comput. Theor. Chem. 2014, 1040–1041, 45– 53, DOI: 10.1016/j.comptc.2014.02.02320https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXksFylt7o%253D&md5=a51265489cf75ee7284bf9cd31ece615A simplified time-dependent density functional theory approach for electronic ultraviolet and circular dichroism spectra of very large moleculesBannwarth, Christoph; Grimme, StefanComputational & Theoretical Chemistry (2014), 1040-1041 (), 45-53CODEN: CTCOA5; ISSN:2210-271X. (Elsevier B.V.)We present a simplified time-dependent d. functional theory approach (sTD-DFT) that allows fast computation of electronic UV or CD spectra of mols. with 500-1000 atoms. The matrix elements are treated in the same way as in the recently proposed simplified Tamm-Dancoff approach but instead of applying the Tamm-Dancoff approxn., the std. linear-response d. functional theory problem is solved. Compared to sTDA, the method leads to an increase in computation time (typically a factor of 2-5 compared to the corresponding sTDA) which is justified since the resulting transition dipole moments are in general of higher quality. This becomes important if spectral intensities (e.g. single-photon oscillator and rotatory transition strengths) are of interest. Comparison of UV and CD spectra obtained from sTD-DFT and sTDA for some typical systems employing std. hybrid functionals shows that both yield very similar excitation energies but the advantage of using the former approach for transition moments. In order to show the applicability of sTD-DFT to systems which are far beyond the scope of conventional TD-DFT, we present the CD spectrum of a substituted, chiral fullerene over a range of almost 1200 excited states. We propose this method as a more reliable alternative for the prediction esp. of the more challenging CD spectra.
- 21Grimme, S. A simplified Tamm-Dancoff density functional approach for the electronic excitation spectra of very large molecules. J. Chem. Phys. 2013, 138, 244104, DOI: 10.1063/1.481133121https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3sXpvFCit7s%253D&md5=1c591b50b78ab4df676b394c5754f9d3A simplified Tamm-Dancoff density functional approach for the electronic excitation spectra of very large moleculesGrimme, StefanJournal of Chemical Physics (2013), 138 (24), 244104/1-244104/14CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Two approxns. in the Tamm-Dancoff d. functional theory approach (TDA-DFT) to electronically excited states are proposed which allow routine computations for electronic UV- or CD spectra of mols. with 500-1000 atoms. Speed-ups compared to conventional time-dependent DFT (TD-DFT) treatments of about two to three orders of magnitude in the excited state part at only minor loss of accuracy are obtained. The method termed sTDA ("s" for simplified) employs atom-centered Loewdin-monopole based two-electron repulsion integrals with the asymptotically correct 1/R behavior and perturbative single excitation configuration selection. It is formulated generally for any std. global hybrid d. functional with given Fock-exchange mixing parameter ax. The method performs well for two std. benchmark sets of vertical singlet-singlet excitations for values of ax in the range 0.2-0.6. The mean abs. deviations from ref. data are only 0.2-0.3 eV and similar to those from std. TD-DFT. In three cases (two dyes and one polypeptide), good mutual agreement between the electronic spectra (up to 10-11 eV excitation energy) from the sTDA method and those from TD(A)-DFT is obtained. The computed UV- and CD-spectra of a few typical systems (e.g., C60, two transition metal complexes, [7]helicene, polyalanine, a supramol. aggregate with 483 atoms and about 7000 basis functions) compare well with corresponding exptl. data. The method is proposed together with medium-sized double- or triple-zeta type at.-orbital basis sets as a quantum chem. tool to investigate the spectra of huge mol. systems at a reliable DFT level. (c) 2013 American Institute of Physics.
- 22Grimme, S.; Bannwarth, C. Ultra-fast computation of electronic spectra for large systems by tight-binding based simplified Tamm-Dancoff approximation (sTDA-xTB). J. Chem. Phys. 2016, 145, 054103, DOI: 10.1063/1.495960522https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28Xht1Klt7nK&md5=2ff48457d9d88c019084d3fdf6fee0e9Ultra-fast computation of electronic spectra for large systems by tight-binding based simplified Tamm-Dancoff approximation (sTDA-xTB)Grimme, Stefan; Bannwarth, ChristophJournal of Chemical Physics (2016), 145 (5), 054103/1-054103/20CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The computational bottleneck of the extremely fast simplified Tamm-Dancoff approximated (sTDA) time-dependent d. functional theory procedure [S. Grimme, J. Chem. Phys. 138, 244104 (2013)] for the computation of electronic spectra for large systems is the detn. of the ground state Kohn-Sham orbitals and eigenvalues. This limits such treatments to single structures with a few hundred atoms and hence, e.g., sampling along mol. dynamics trajectories for flexible systems or the calcn. of chromophore aggregates is often not possible. The aim of this work is to solve this problem by a specifically designed semiempirical tight binding (TB) procedure similar to the well established self-consistent-charge d. functional TB scheme. The new special purpose method provides orbitals and orbital energies of hybrid d. functional character for a subsequent and basically unmodified sTDA procedure. Compared to many previous semiempirical excited state methods, an advantage of the ansatz is that a general eigenvalue problem in a non-orthogonal, extended AO basis is solved and therefore correct occupied/virtual orbital energy splittings as well as Rydberg levels are obtained. A key idea for the success of the new model is that the detn. of at. charges (describing an effective electron-electron interaction) and the one-particle spectrum is decoupled and treated by two differently parametrized Hamiltonians/basis sets. The three-diagonalization-step composite procedure can routinely compute broad range electronic spectra (0-8 eV) within minutes of computation time for systems composed of 500-1000 atoms with an accuracy typical of std. time-dependent d. functional theory (0.3-0.5 eV av. error). An easily extendable parametrization based on coupled-cluster and d. functional computed ref. data for the elements H-Zn including transition metals is described. The accuracy of the method termed sTDA-xTB is first benchmarked for vertical excitation energies of open- and closed-shell systems in comparison to other semiempirical methods and applied to exemplary problems in electronic spectroscopy. As side products of the development, a robust and efficient valence electron TB method for the accurate detn. of at. charges as well as a more accurate calcn. scheme of dipole rotatory strengths within the Tamm-Dancoff approxn. is proposed. (c) 2016 American Institute of Physics.
- 23de Wergifosse, M.; Bannwarth, C.; Grimme, S. A Simplified Spin-Flip Time-Dependent Density Functional Theory Approach for the Electronic Excitation Spectra of Very Large Diradicals. J. Phys. Chem. A 2019, 123, 5815– 5825, DOI: 10.1021/acs.jpca.9b0317623https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXhtFyjsbzE&md5=e9769d3847d1c74864f2d04259fd5a9dA Simplified Spin-Flip Time-Dependent Density Functional Theory Approach for the Electronic Excitation Spectra of Very Large Diradicalsde Wergifosse, Marc; Bannwarth, Christoph; Grimme, StefanJournal of Physical Chemistry A (2019), 123 (27), 5815-5825CODEN: JPCAFH; ISSN:1089-5639. (American Chemical Society)Experimentalists working with diradicals are often facing the question of what kind of species among singlet or triplet diradicals or closed-shell mols. are obsd. To treat large diradicals with a high d. of electronic states, we propose a simplified version of the spin-flip time-dependent d. functional theory (SF-TD-DFT) method for a fast computation of their state energies and absorption spectra with an accuracy similar to the nonsimplified scheme. An ultrafast tight-binding variant called SF-sTD-DFT-xTB is also developed to treat even larger systems. For a benchmark set of nine diradicals, good agreement between simplified and conventional SF excitation energies for std. functionals is found. This shows that the proposed parameterization is robust for a wide range of Fock exchange mixing values. With the asymptotically correct response integrals used in SF-sTD-DFT and a correction factor of √(2) for the transition moments, the SF-sTD-DFT/B5050LYP/cc-pVDZ method even outperforms the nonsimplified scheme at drastically reduced computational effort when comparing to the exptl. absorption spectra for this set of diradicals. To showcase the actual performance of the method, absorption spectra of two μ-hydroxo-bridged dimers of corrole tape Ga(III) complex derivs. were computed and compared to the expt., providing good qual. agreement. Finally, a comparison with the high-spin triplet spectrum of a perylene bisimide biradical and the one detd. at the SF-sTD-DFT level showed that at room temp., mostly triplet diradicals are present in soln.
- 24De Wergifosse, M.; Grimme, S. Nonlinear-response properties in a simplified time-dependent density functional theory (sTD-DFT) framework: Evaluation of excited-state absorption spectra. J. Chem. Phys. 2019, 150, 094112, DOI: 10.1063/1.508019924https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXktlersLY%253D&md5=8f0cea7b5ce0b88cc73e98a90cce4240Nonlinear-response properties in a simplified time-dependent density functional theory (sTD-DFT) framework: Evaluation of excited-state absorption spectrade Wergifosse, Marc; Grimme, StefanJournal of Chemical Physics (2019), 150 (9), 094112/1-094112/9CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The energy conversion efficiency of org. solar cells seems crucial for a clean future. The design of new light-harvesting devices needs an in-depth understanding of their optical properties, including the excited-state absorption (ESA). In biol., the optical characterization of photochem./phys. processes happening in photosynthetic pigments and proteins can be difficult to interpret due to their structural complexities. Exptl., an ultrafast transient absorption expt. can probe the excited state interaction with light. Quantum chem. could play an important role to model the transient absorption spectrum of excited states. However, systems that need to be investigated can be way too large for existent software implementations. In this contribution, we present the first sTDA/sTD-DFT (simplified time-dependent d. functional theory with and without Tamm Dancoff approxn.) implementation to evaluate the ESA of mols. The ultrafast ESA evaluation presents a negligible extra cost with respect to sTDA/sTD-DFT original schemes for std. ground state absorption. The sTD-DFT method shows ability to assign ESA spectra to the correct excited state. We showed that in the literature, wrong assignments were proposed as for the L34/L44 mixt. and N-methylfulleropyrrolidine. In addn., sTDA/sTD-DFT-xTB tight-binding variants are also available, allowing the evaluation of ESA for systems of a few thousands of atoms, e.g., the spectrum of the photoactive yellow protein composed of 1931 atoms. (c) 2019 American Institute of Physics.
- 25De Wergifosse, M.; Grimme, S. Nonlinear-response properties in a simplified time-dependent density functional theory (sTD-DFT) framework: Evaluation of the first hyperpolarizability. J. Chem. Phys. 2018, 149, 024108, DOI: 10.1063/1.503766525https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXhtlWqurzI&md5=8a2bb42755257c0761582ab77e2efdc5Nonlinear-response properties in a simplified time-dependent density functional theory (sTD-DFT) framework: Evaluation of the first hyperpolarizabilityde Wergifosse, Marc; Grimme, StefanJournal of Chemical Physics (2018), 149 (2), 024108/1-024108/13CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Recent developments in nonlinear imaging microscopy show the need to implement new theor. tools, which are able to characterize nonlinear optical properties in an efficient way. For second-harmonic imaging microscopy (SHIM), quantum chem. could play an important role to design new exogenous dyes with enhanced first hyperpolarizabilities or to characterize the response origin in large endogenous biol. systems. Such methods should be able to screen a large no. of compds. while reproducing their trends and to treat large systems in reasonable computation times. To fulfill these requirements, we present a new simplified time-dependent d. functional theory (sTD-DFT) implementation to evaluate the first hyperpolarizability where the Coulomb and exchange integrals are approximated by short-range damped Coulomb interactions of transition d. monopoles. For an ultra-fast computation of the first hyperpolarizability, a tight-binding version (sTD-DFT-xTB) is also proposed. In our implementation, a sTD-DFT calcn. is more than 600 time faster with respect to a regular TD-DFT treatment, while the xTB version speeds up the entire calcn. further by at least two orders of magnitude. We challenge our implementation on three test cases: typical push-pull π-conjugated compds., fluorescent proteins, and a collagen model, which were selected to model requirements for SHIM applications. (c) 2018 American Institute of Physics.
- 26Bannwarth, C.; Seibert, J.; Grimme, S. Electronic Circular Dichroism of [16]Helicene With Simplified TD-DFT: Beyond the Single Structure Approach. Chirality 2016, 28, 365– 369, DOI: 10.1002/chir.2259426https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XlvV2mtrc%253D&md5=40b2daee5ca75233db6ffaa58bccbb58Electronic Circular Dichroism of [16]Helicene With Simplified TD-DFT: Beyond the Single Structure ApproachBannwarth, Christoph; Seibert, Jakob; Grimme, StefanChirality (2016), 28 (5), 365-369CODEN: CHRLEP; ISSN:0899-0042. (Wiley-Liss, Inc.)The electronic CD (ECD) spectrum of the recently synthesized [16]helicene and a deriv. comprising two triisopropylsilyloxy protection groups was computed by means of the very efficient simplified time-dependent d. functional theory (sTD-DFT) approach. Different from many previous ECD studies of helicenes, nonequil. structure effects were accounted for by computing ECD spectra on "snapshots" obtained from a mol. dynamics (MD) simulation including solvent mols. The trajectories are based on a mol. specific classical potential as obtained from the recently developed quantum chem. derived force field (QMDFF) scheme. The reduced computational cost in the MD simulation due to the use of the QMDFF (compared to ab-initio MD) as well as the sTD-DFT approach make realistic spectral simulations feasible for these compds. that comprise more than 100 atoms. While the ECD spectra of [16]helicene and its deriv. computed vertically on the resp. gas phase, equil. geometries show noticeable differences, these are "washed" out when nonequil. structures are taken into account. The computed spectra with two recommended d. functionals (ωB97X and BHLYP) and extended basis sets compare very well with the exptl. one. In addn. we provide an est. for the missing abs. intensities of the latter. The approach presented here could also be used in future studies to capture nonequil. effects, but also to systematically av. ECD spectra over different conformations in more flexible mols. Chirality 00:000-000, 2016. 2016 Wiley Periodicals, Inc.
- 27Risthaus, T.; Hansen, A.; Grimme, S. Excited states using the simplified Tamm-Dancoff-Approach for range-separated hybrid density functionals: development and application. Phys. Chem. Chem. Phys. 2014, 16, 14408– 14419, DOI: 10.1039/C3CP54517B27https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXhtVKgs7bM&md5=978384902d7b53a6ef2eefa57ef407feExcited states using the simplified Tamm-Dancoff-Approach for range-separated hybrid density functionals: development and applicationRisthaus, Tobias; Hansen, Andreas; Grimme, StefanPhysical Chemistry Chemical Physics (2014), 16 (28), 14408-14419CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)The recently introduced sTDA methodol. [S. Grimme, J. Chem. Phys., 2013, 138, 244104] to compute excitation spectra of huge mol. systems is extended to range-sepd. hybrid (RSH) d. functionals. The three empirical parameters of the method which describe a screened two-electron interaction are obtained for some common RSH functionals (ωB97 family, CAM-B3LYP, LC-BLYP) from a fit to theor. SCS-CC2 ref. vertical excitation energies for a set of small to medium-sized chromophores. The method is cross-validated on a set of inter- and intramol. charge transfer states and a set composed of typical valence transitions. Overall small deviations from ref. data of only about 0.2-0.4 eV are found with best performance for CAM-B3LYP and ωB97X-D3. To demonstrate versatility and robustness of the new methodol., applications (the UV/Vis spectrum of the pyridine polymer and the ECD spectrum of (P)-[11]helicene) and frequently used charge transfer examples are discussed. In one case, 11 000 + excited electronic states of a system contg. 330 atoms were calcd. We show that the asymptotically correct sTDA-RSH combination yields results often superior to those based on global hybrids and that it opens up new possibilities for the computation of excited states in materials science and bio-mol. systems.
- 28Elstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, T.; Suhai, S.; Seifert, G. Self-consistent-charge density-functional tight-binding method for simulations of complex materials properties. Phys. Rev. B 1998, 58, 7260– 7268, DOI: 10.1103/PhysRevB.58.726028https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK1cXmtVamsb0%253D&md5=e3b3143c475c70203badf7db6ef668c3Self-consistent-charge density-functional tight-binding method for simulations of complex materials propertiesElstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, Th.; Suhai, S.; Seifert, G.Physical Review B: Condensed Matter and Materials Physics (1998), 58 (11), 7260-7268CODEN: PRBMDO; ISSN:0163-1829. (American Physical Society)We outline details about an extension of the tight-binding (TB) approach to improve total energies, forces, and transferability. The method is based on a second-order expansion of the Kohn-Sham total energy in d.-functional theory (DFT) with respect to charge-d. fluctuations. The zeroth-order approach is equiv. to a common std. non-self-consistent (TB) scheme, while at second-order a transparent, parameter-free, and readily calculable expression for generalized Hamiltonian matrix elements may be derived. These are modified by a self-consistent redistribution of Mulliken charges (SCC). Besides the usual "band structure" and short-range repulsive terms the final approx. Kohn-Sham energy addnl. includes a Coulomb interaction between charge fluctuations. At large distances this accounts for long-range electrostatic forces between two point charges and approx. includes self-interaction contributions of a given atom if the charges are located at one and the same atom. We apply the new SCC scheme to problems where deficiencies within the non-SCC std. TB approach become obvious. We thus considerably improve transferability.
- 29Seibert, J.; Pisarek, J.; Schmitz, S.; Bannwarth, C.; Grimme, S. Extension of the element parameter set for ultra-fast excitation spectra calculation (sTDA-xTB). Mol. Phys. 2019, 117, 1104– 1116, DOI: 10.1080/00268976.2018.151014129https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXhsFCjs7bK&md5=57617dc9f5fcc3316f1edabcb4581b83Extension of the element parameter set for ultra-fast excitation spectra calculation (sTDA-xTB)Seibert, Jakob; Pisarek, Jana; Schmitz, Sarah; Bannwarth, Christoph; Grimme, StefanMolecular Physics (2019), 117 (9-12), 1104-1116CODEN: MOPHAM; ISSN:0026-8976. (Taylor & Francis Ltd.)The extension of the parameter set for an ultra-fast electronic excitation spectra calcn. is presented. The semiempirical theory based on a tight-binding approach, called extended tight-binding (xTB) in combination with the simplified Tamm-Dancoff approxn. (sTDA) shows remarkable accuracy at very low computational cost for the calcn. of vertical excitation energies of mols. It enables the possibility for computing even large systems up to thousands of atoms or sampling along mol. dynamic (MD) trajectories. The original publication of the sTDA-xTB method included parameters for the most important elements (H-Zn,Br,I). In this work, element parameters for 4d and 5d metals, and the missing ones in 4p, 5p and 6p element blocks are presented and analyzed for their quality. Comparisons to theory and expt. show that sTDA-xTB provides similar good results as for the elements in the original publication with an av. deviation of excitation energies of 0.3-0.5 eV.
- 30Cho, Y.; Bintrim, S. J.; Berkelbach, T. C. A simplified GW/BSE approach for charged and neutral excitation energies of large molecules and nanomaterials ; 2021, arXiv:2109.04421. DOI: 10.48550/arXiv.2109.04421 .There is no corresponding record for this reference.
- 31Rüger, R.; van Lenthe, E.; Heine, T.; Visscher, L. Tight-binding approximations to time-dependent density functional theory ─ A fast approach for the calculation of electronically excited states. J. Chem. Phys. 2016, 144, 184103, DOI: 10.1063/1.494864731https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XnslWru7g%253D&md5=0efdefe5a46f1be35cbf902b0cfc3726Tight-binding approximations to time-dependent density functional theory - A fast approach for the calculation of electronically excited statesRueger, Robert; van Lenthe, Erik; Heine, Thomas; Visscher, LucasJournal of Chemical Physics (2016), 144 (18), 184103/1-184103/12CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We propose a new method of calcg. electronically excited states that combines a d. functional theory based ground state calcn. with a linear response treatment that employs approxns. used in the time-dependent d. functional based tight binding (TD-DFTB) approach. The new method termed time-dependent d. functional theory TD-DFT+TB does not rely on the DFTB parametrization and is therefore applicable to systems involving all combinations of elements. We show that the new method yields UV/Vis absorption spectra that are in excellent agreement with computationally much more expensive TD-DFT calcns. Errors in vertical excitation energies are reduced by a factor of two compared to TD-DFTB. (c) 2016 American Institute of Physics.
- 32Asadi-Aghbolaghi, N.; Pototschnig, J.; Jamshidi, Z.; Visscher, L. Effects of ligands on (de-)enhancement of plasmonic excitations of silver, gold and bimetallic nanoclusters: TD-DFT+TB calculations. Phys. Chem. Chem. Phys. 2021, 23, 17929– 17938, DOI: 10.1039/D1CP03220H32https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3MXhslGms7%252FN&md5=00c70d22cdec2f0da8c7a1a1de065a96Effects of ligands on (de-)enhancement of plasmonic excitations of silver, gold and bimetallic nanoclusters: TD-DFT+TB calculationsAsadi-Aghbolaghi, Narges; Pototschnig, Johann; Jamshidi, Zahra; Visscher, LucasPhysical Chemistry Chemical Physics (2021), 23 (33), 17929-17938CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)Metal nanoclusters can be synthesized in various sizes and shapes and are typically protected with ligands to stabilize them. These ligands can also be used to tune the plasmonic properties of the clusters as the absorption spectrum of a protected cluster can be significantly altered compared to the bare cluster. In this paper, we computationally investigate the influence of thiolate ligands on the plasmonic intensity for silver, gold and alloy clusters. Using time-dependent d. functional theory with tight-binding approxns., TD-DFT+TB, we show that this level of theory can reproduce the broad exptl. spectra of Au144(SR)60 and Ag53Au91(SR)60 (R = CH3) compds. with satisfactory agreement. As TD-DFT+TB does not depend on atom-type parameters we were able to apply this approach on large ligand-protected clusters with various compns. With these calcns. we predict that the effect of ligands on the absorption can be a quenching as well as an enhancement. We furthermore show that it is possible to unambiguously identify the plasmonic peaks by the scaled Coulomb kernel technique and explain the influence of ligands on the intensity (de-)enhancement by analyzing the plasmonic excitations in terms of the dominant orbital contributions.
- 33Van Caillie, C.; Amos, R. D. Geometric derivatives of excitation energies using SCF and DFT. Chem. Phys. Lett. 1999, 308, 249– 255, DOI: 10.1016/S0009-2614(99)00646-633https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK1MXksFWku70%253D&md5=b0a42882739b7c4d495c632cc8205f5aGeometric derivatives of excitation energies using SCF and DFTVan Caillie, Carole; Amos, Roger D.Chemical Physics Letters (1999), 308 (3,4), 249-255CODEN: CHPLBC; ISSN:0009-2614. (Elsevier Science B.V.)There is increasing interest in using the methods of time-dependent d. functional theory to calc. electronic excitation energies. We have implemented an analytic gradient method to find the geometric derivs. of the excitation energies. When added to the gradient for the ground state, this yields the excited-state energy derivs. This enables the efficient generation and searching of excited-state potential energy surfaces to obtain excited-state geometries and other properties. The initial implementation is for SCF methods and for the local d. approxn. Some examples of excited-state geometry optimizations are given.
- 34Van Caillie, C.; Amos, R. D. Geometric derivatives of density functional theory excitation energies using gradient-corrected functionals. Chem. Phys. Lett. 2000, 317, 159– 164, DOI: 10.1016/S0009-2614(99)01346-934https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3cXotlyksw%253D%253D&md5=fb54eacfa34d1288afe2075376db3e3cGeometric derivatives of density functional theory excitation energies using gradient-corrected functionalsVan Caillie, C.; Amos, R. D.Chemical Physics Letters (2000), 317 (1,2), 159-164CODEN: CHPLBC; ISSN:0009-2614. (Elsevier Science B.V.)D. functional theory (DFT) is having increasing success in predicting excitation energies using the methods of time-dependent DFT. As a result, it should be possible to generate potential energy surfaces for excited states by adding the excitation energy, as a function of geometry, to the ground-state energy. It is easier to find stationary points such as min. and transition states if the gradient of the energy is known. The present Letter extends earlier work on the gradients on excited-state surfaces using SCF and LDA (local d. approxn.) methods, to use gradient-cor. and hybrid functionals. Some examples of geometry optimizations are given.
- 35Furche, F.; Ahlrichs, R. Adiabatic time-dependent density functional methods for excited state properties. J. Chem. Phys. 2002, 117, 7433– 7447, DOI: 10.1063/1.150836835https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD38XnvVWrurY%253D&md5=061f0264a5993d772854715400d3d189Adiabatic time-dependent density functional methods for excited state propertiesFurche, Filipp; Ahlrichs, ReinhartJournal of Chemical Physics (2002), 117 (16), 7433-7447CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)This work presents theory, implementation, and validation of excited state properties obtained from time-dependent d. functional theory (TDDFT). Based on a fully variational expression for the excited state energy, a compact derivation of first order properties is given. We report an implementation of analytic excited state gradients and charge moments for local, gradient cor., and hybrid functionals, as well as for the CI singles (CIS) and time-dependent Hartree-Fock (TDHF) methods. By exploiting analogies to ground state energy and gradient calcns., efficient techniques can be transferred to excited state methods. Benchmark results demonstrate that, for low-lying excited states, geometry optimizations are not substantially more expensive than for the ground state, independent of the mol. size. We assess the quality of calcd. adiabatic excitation energies, structures, dipole moments, and vibrational frequencies by comparison with accurate exptl. data for a variety of excited states and mols. Similar trends are obsd. for adiabatic excitation energies as for vertical ones. TDDFT is more robust than CIS and TDHF, in particular, for geometries differing significantly from the ground state min. The TDDFT excited state structures, dipole moments, and vibrational frequencies are of a remarkably high quality, which is comparable to that obtained in ground state d. functional calcns. Thus, yielding considerably more accurate results at similar computational cost, TDDFT rivals CIS as a std. method for calcg. excited state properties in larger mols.
- 36Furche, F.; Ahlrichs, R. Erratum: ”Time-dependent density functional methods for excited state properties” [J. Chem. Phys. 117, 7433 (2002)]. J. Chem. Phys. 2004, 121, 12772– 12773, DOI: 10.1063/1.182490336https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2cXhtFSnsrnE&md5=c2047481400c80542918b0de9189c0baAdiabatic time-dependent density functional methods for excited state properties. [Erratum to document cited in CA138:044999]Furche, Filipp; Ahlrichs, ReinhartJournal of Chemical Physics (2004), 121 (24), 12772-12773CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)A factor 1/2 is missing in the expression for the diagonal elements of the energy weighted difference d. matrix W throughout the paper. The cor. versions of Equations (24), (A4), (A10), and (A13) and Tables V, VI, and VIII are given. The implementation and results presented in the paper were based on the correct expressions. Some of the calcd. spectroscopic consts. of the 1 1.sum.u- state of N2, of the 1 1.sum.u+ state of Mg2, and of the 2 1.sum.+ state of CuH in Tables V, VI, and VII are incorrect and should be replaced by the values in the cor. equations. The conclusions are not affected by these changes.
- 37Hutter, J. Excited state nuclear forces from the Tamm-Dancoff approximation to time-dependent density functional theory within the plane wave basis set framework. J. Chem. Phys. 2003, 118, 3928– 3934, DOI: 10.1063/1.154010937https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3sXht1CltLc%253D&md5=4099596d811bae697319ed324e79ce8cExcited state nuclear forces from the Tamm-Dancoff approximation to time-dependent density functional theory within the plane wave basis set frameworkHutter, JurgJournal of Chemical Physics (2003), 118 (9), 3928-3934CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)An efficient formulation of time-dependent linear response d. functional theory for the use within the plane wave basis set framework is presented. The method avoids the transformation of the Kohn-Sham matrix into the canonical basis and refs. virtual orbitals only through a projection operator. Using a Lagrangian formulation nuclear derivs. of excited state energies within the Tamm-Dancoff approxn. are derived. The algorithms were implemented into a pseudo potential/plane wave code and applied to the calcn. of adiabatic excitation energies, optimized geometries and vibrational frequencies of three low lying states of formaldehyde. An overall good agreement with other time-dependent d. functional calcns., multireference CI calcns. and exptl. data was found.
- 38Petrenko, T.; Kossmann, S.; Neese, F. Efficient time-dependent density functional theory approximations for hybrid density functionals: Analytical gradients and parallelization. J. Chem. Phys. 2011, 134, 054116, DOI: 10.1063/1.353344138https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXhsFSqtrs%253D&md5=1eae10d380aea0e3920f1870fc2a181fEfficient time-dependent density functional theory approximations for hybrid density functionals: Analytical gradients and parallelizationPetrenko, Taras; Kossmann, Simone; Neese, FrankJournal of Chemical Physics (2011), 134 (5), 054116/1-054116/14CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We present the implementation of efficient approxns. to time-dependent d. functional theory (TDDFT) within the Tamm-Dancoff approxn. (TDA) for hybrid d. functionals. For the calcn. of the TDDFT/TDA excitation energies and anal. gradients, we combine the resoln. of identity (RI-J) algorithm for the computation of the Coulomb terms and the recently introduced "chain of spheres exchange" (COSX) algorithm for the calcn. of the exchange terms. It is shown that for extended basis sets, the RIJCOSX approxn. leads to speedups of up to 2 orders of magnitude compared to traditional methods, as demonstrated for hydrocarbon chains. The accuracy of the adiabatic transition energies, excited state structures, and vibrational frequencies is assessed on a set of 27 excited states for 25 mols. with the CI singles and hybrid TDDFT/TDA methods using various basis sets. Compared to the canonical values, the typical error in transition energies is of the order of 0.01 eV. Similar to the ground-state results, excited state equil. geometries differ by less than 0.3 pm in the bond distances and 0.5° in the bond angles from the canonical values. The typical error in the calcd. excited state normal coordinate displacements is of the order of 0.01, and relative error in the calcd. excited state vibrational frequencies is less than 1%. The errors introduced by the RIJCOSX approxn. are, thus, insignificant compared to the errors related to the approx. nature of the TDDFT methods and basis set truncation. For TDDFT/TDA energy and gradient calcns. on Ag-TB2-helicate (156 atoms, 2732 basis functions), it is demonstrated that the COSX algorithm parallelizes almost perfectly (speedup ∼26-29 for 30 processors). The exchange-correlation terms also parallelize well (speedup ∼27-29 for 30 processors). The soln. of the Z-vector equations shows a speedup of ∼24 on 30 processors. The parallelization efficiency for the Coulomb terms can be somewhat smaller (speedup ∼15-25 for 30 processors), but their contribution to the total calcn. time is small. Thus, the parallel program completes a Becke3-Lee-Yang-Parr energy and gradient calcn. on the Ag-TB2-helicate in less than 4 h on 30 processors. We also present the necessary extension of the Lagrangian formalism, which enables the calcn. of the TDDFT excited state properties in the frozen-core approxn. The algorithms described in this work are implemented into the ORCA electronic structure system. (c) 2011 American Institute of Physics.
- 39Grotjahn, R.; Furche, F.; Kaupp, M. Development and Implementation of Excited-State Gradients for Local Hybrid Functionals. J. Chem. Theory Comput. 2019, 15, 5508– 5522, DOI: 10.1021/acs.jctc.9b0065939https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXhs12jsbbE&md5=5943c43f7445f4cf7546917346d65f36Development and Implementation of Excited-State Gradients for Local Hybrid FunctionalsGrotjahn, Robin; Furche, Filipp; Kaupp, MartinJournal of Chemical Theory and Computation (2019), 15 (10), 5508-5522CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)Local hybrid functionals are a relatively recent class of exchange-correlation functionals that use a real-space dependent admixt. of exact exchange. Here, we present the first implementation of time-dependent d. functional theory excited-state gradients for these functionals. Based on the ansatz of a fully variational auxiliary Lagrangian of the excitation energy, the working equations for the case of a local hybrid functional are deduced. For the implementation, we derive the third-order functional derivs. used in the hyper-kernel and kernel-gradients following a seminumerical integration scheme. The first assessment for a test set of small mols. reveals competitive performance for excited-state structural parameters with typical mean abs. errors (MAEs) of 1.2 pm (PBE0: 1.4 pm) for bond lengths as well as for vibrational frequencies with typical MAEs of 81 cm-1 (PBE0: 76 cm-1). Excellent performance was found for adiabatic triplet excitation energies with typical MAEs of 0.08 eV (PBE0: 0.32 eV). In a detailed case anal. of the first singlet and triplet excited states of formaldehyde the conceptional (dis-)advantages of the local hybrid scheme for excited-state gradients are exposed.
- 40Kretz, B.; Egger, D. A. Accurate Molecular Geometries in Complex Excited-State Potential Energy Surfaces from Time-Dependent Density Functional Theory. J. Chem. Theory Comput. 2021, 17, 357– 366, DOI: 10.1021/acs.jctc.0c0085840https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXisVyrurjN&md5=74c6717ff8768412998953095ed42cf5Accurate Molecular Geometries in Complex Excited-State Potential Energy Surfaces from Time-Dependent Density Functional TheoryKretz, Bernhard; Egger, David A.Journal of Chemical Theory and Computation (2021), 17 (1), 357-366CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)The interplay of electronic excitations and structural changes in mols. impacts nonradiative decay and charge transfer in the excited state, thus influencing excited-state lifetimes and photocatalytic reaction rates in optoelectronic and energy devices. To capture such effects requires computational methods providing an accurate description of excited-state potential energy surfaces and geometries. We suggest time-dependent d. functional theory using optimally tuned range-sepd. hybrid (OT-RSH) functionals as an accurate approach to obtain excited-state mol. geometries. We show that OT-RSH provides accurate mol. geometries in excited-state potential energy surfaces that are complex and involve an interplay of local and charge-transfer excitations, for which conventional semilocal and hybrid functionals fail. At the same time, the nonempirical OT-RSH approach maintains the high accuracy of parametrized functionals (e.g., B3LYP) for predicting excited-state geometries of small org. mols. showing valence excited states.
- 41Sokolov, M.; Bold, B. M.; Kranz, J. J.; Höfener, S.; Niehaus, T. A.; Elstner, M. Analytical Time-Dependent Long-Range Corrected Density Functional Tight Binding (TD-LC-DFTB) Gradients in DFTB+: Implementation and Benchmark for Excited-State Geometries and Transition Energies. J. Chem. Theory Comput. 2021, 17, 2266– 2282, DOI: 10.1021/acs.jctc.1c0009541https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3MXlvFWit7Y%253D&md5=433eaf26623e19cbabe5f82254b17193Analytical Time-Dependent Long-Range Corrected Density Functional Tight Binding (TD-LC-DFTB) Gradients in DFTB+: Implementation and Benchmark for Excited-State Geometries and Transition EnergiesSokolov, Monja; Bold, Beatrix M.; Kranz, Julian J.; Hoefener, Sebastian; Niehaus, Thomas A.; Elstner, MarcusJournal of Chemical Theory and Computation (2021), 17 (4), 2266-2282CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)The absorption and emission of light is a ubiquitous process in chem. and biol. processes, making a theor. description inevitable for understanding and predicting such properties. Although ab initio and DFT methods are capable of describing excited states with good accuracy in many cases, the investigation of dynamical processes and the need to sample the phase space in complex systems often requires methods with reduced computational costs but still sufficient accuracy. In the present work, we report the derivation and implementation of anal. nuclear gradients for time-dependent long-range cor. d. functional tight binding (TD-LC-DFTB) in the DFTB+ program. The accuracy of the TD-LC-DFTB potential-energy surfaces is benchmarked for excited-state geometries and adiabatic as well as vertical transition energies. The benchmark set consists of more than 100 org. mols. taken as subsets from available benchmark sets. The reported method yields a mean deviation of 0.31 eV for adiabatic excitation energies with respect to CC2. In order to study more subtle effects, seminumerical second derivs. based on the anal. gradients are employed to simulate vibrationally resolved UV/vis spectra. This extensive test exhibits few problematic cases, which can be traced back to the parametrization of the repulsive potential.
- 42Kühne, T. D.; Iannuzzi, M.; Del Ben, M.; Rybkin, V. V.; Seewald, P.; Stein, F.; Laino, T.; Khaliullin, R. Z.; Schütt, O.; Schiffmann, F.; Golze, D.; Wilhelm, J.; Chulkov, S.; Bani-Hashemian, M. H.; Weber, V.; Borštnik, U.; Taillefumier, M.; Jakobovits, A. S.; Lazzaro, A.; Pabst, H.; Müller, T.; Schade, R.; Guidon, M.; Andermatt, S.; Holmberg, N.; Schenter, G. K.; Hehn, A.; Bussy, A.; Belleflamme, F.; Tabacchi, G.; Glöß, A.; Lass, M.; Bethune, I.; Mundy, C. J.; Plessl, C.; Watkins, M.; VandeVondele, J.; Krack, M.; Hutter, J. CP2K: An electronic structure and molecular dynamics software package - Quickstep: Efficient and accurate electronic structure calculations. J. Chem. Phys. 2020, 152, 194103, DOI: 10.1063/5.000704542https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXhtVSgtrzF&md5=b9e5975bc402f0d53e2d99da998adf5fCP2K: An electronic structure and molecular dynamics software package - Quickstep: Efficient and accurate electronic structure calculationsKuehne, Thomas D.; Iannuzzi, Marcella; Del Ben, Mauro; Rybkin, Vladimir V.; Seewald, Patrick; Stein, Frederick; Laino, Teodoro; Khaliullin, Rustam Z.; Schuett, Ole; Schiffmann, Florian; Golze, Dorothea; Wilhelm, Jan; Chulkov, Sergey; Bani-Hashemian, Mohammad Hossein; Weber, Valery; Borstnik, Urban; Taillefumier, Mathieu; Jakobovits, Alice Shoshana; Lazzaro, Alfio; Pabst, Hans; Mueller, Tiziano; Schade, Robert; Guidon, Manuel; Andermatt, Samuel; Holmberg, Nico; Schenter, Gregory K.; Hehn, Anna; Bussy, Augustin; Belleflamme, Fabian; Tabacchi, Gloria; Gloess, Andreas; Lass, Michael; Bethune, Iain; Mundy, Christopher J.; Plessl, Christian; Watkins, Matt; VandeVondele, Joost; Krack, Matthias; Hutter, JuergJournal of Chemical Physics (2020), 152 (19), 194103CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)A review. CP2K is an open source electronic structure and mol. dynamics software package to perform atomistic simulations of solid-state, liq., mol., and biol. systems. It is esp. aimed at massively parallel and linear-scaling electronic structure methods and state-of-the-art ab initio mol. dynamics simulations. Excellent performance for electronic structure calcns. is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2K to perform efficient and accurate electronic structure simulations. The emphasis is put on d. functional theory and multiple post-Hartree-Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension. (c) 2020 American Institute of Physics.
- 43VandeVondele, J.; Krack, M.; Mohamed, F.; Parrinello, M.; Chassaing, T.; Hutter, J. Quickstep: Fast and accurate density functional calculations using a mixed Gaussian and plane waves approach. Comput. Phys. Commun. 2005, 167, 103– 128, DOI: 10.1016/j.cpc.2004.12.01443https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2MXjt1aitb4%253D&md5=8c5393031c9dbd341e0e73fcdacad486QUICKSTEP: fast and accurate density functional calculations using a mixed Gaussian and plane waves approachVandeVondele, Joost; Krack, Matthias; Mohamed, Fawzi; Parrinello, Michele; Chassaing, Thomas; Hutter, JuergComputer Physics Communications (2005), 167 (2), 103-128CODEN: CPHCBZ; ISSN:0010-4655. (Elsevier B.V.)We present the Gaussian and plane waves (GPW) method and its implementation in which is part of the freely available program package CP2K. The GPW method allows for accurate d. functional calcns. in gas and condensed phases and can be effectively used for mol. dynamics simulations. We show how derivs. of the GPW energy functional, namely ionic forces and the Kohn-Sham matrix, can be computed in a consistent way. The computational cost of computing the total energy and the Kohn-Sham matrix is scaling linearly with the system size, even for condensed phase systems of just a few tens of atoms. The efficiency of the method allows for the use of large Gaussian basis sets for systems up to 3000 atoms, and we illustrate the accuracy of the method for various basis sets in gas and condensed phases. Agreement with basis set free calcns. for single mols. and plane wave based calcns. in the condensed phase is excellent. Wave function optimization with the orbital transformation technique leads to good parallel performance, and outperforms traditional diagonalisation methods. Energy conserving Born-Oppenheimer dynamics can be performed, and a highly efficient scheme is obtained using an extrapolation of the d. matrix. We illustrate these findings with calcns. using commodity PCs as well as supercomputers.
- 44Lippert, G.; Hutter, J.; Parrinello, M. A hybrid Gaussian and plane wave density functional scheme. Mol. Phys. 1997, 92, 477– 488, DOI: 10.1080/0026897970948211944https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK2sXntVSrs7s%253D&md5=ece29370fc0b1b66ddd0b82b937bfd6bA hybrid Gaussian and plane wave density functional schemeLippert, Gerald; Hutter, Juerg; Parrinello, MicheleMolecular Physics (1997), 92 (3), 477-487CODEN: MOPHAM; ISSN:0026-8976. (Taylor & Francis)A d.-functional theory-based algorithm for periodic and nonperiodic ab initio calcns. is presented. This scheme uses pseudopotentials in order to integrate out the core electrons from the problem. The valence pseudo wave functions are expanded in Gaussian-type orbitals and the d. is represented in a plane wave auxiliary basis. The Gaussian basis functions make it possible to use the efficient anal. integration schemes and screening algorithms of quantum chem. Novel recursion relations are developed for the calcn. of the matrix elements of the d.-dependent Kohn-Sham self-consistent potential. At the same time the use of a plane wave basis for the electron d. permits efficient calcn. of the Hartree energy using fast Fourier transforms, thus circumventing one of the major bottlenecks of std. Gaussian based calcns. Furthermore, this algorithm avoids the fitting procedures that go along with intermediate basis sets for the charge d. The performance and accuracy of this new scheme are discussed and selected examples are given.
- 45Iannuzzi, M.; Chassaing, T.; Wallman, T.; Hutter, J. Ground and Excited State Density Functional Calculations with the Gaussian and Augmented-Plane-Wave Method. CHIMIA Int. J. Chem. 2005, 59, 499– 503, DOI: 10.2533/00094290577767616445https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2MXpvF2msL4%253D&md5=2ef7b146fbe6660372f872678e1bdc8eGround and excited state density functional calculations with the Gaussian and Augmented-Plane-Wave methodIannuzzi, Marcella; Chassaing, Thomas; Wallman, Thomas; Hutter, JurgChimia (2005), 59 (7-8), 499-503CODEN: CHIMAD; ISSN:0009-4293. (Swiss Chemical Society)The calcn. of the electronic structure of large systems by methods based on d. functional theory has recently gained a central role in mol. simulations. However, the extensive study of quantities like excited states and related properties is still out of reach due to high computational costs. We present a new implementation of a hybrid method, the Gaussian and APW (GAPW) method, where the electronic d. is partitioned in hard and soft contributions. The former are local terms naturally expanded in a Gaussian basis, whereas the soft contributions are expanded in plane-waves by using a low energy cutoff, without loss in accuracy, even for all-electron calcns. For the calcn. of excitation energies a recently developed, time-dependent d. functional response theory (TD-DFRT) technique is joined with the GAPW procedure. We demonstrate the accuracy of the method by comparison with std. quantum chem. calcns. for a set of small mols. To highlight the performance and efficiency of GAPW we show calcns. on systems with several thousands of basis functions.
- 46Strand, J.; Chulkov, S. K.; Watkins, M. B.; Shluger, A. L. First principles calculations of optical properties for oxygen vacancies in binary metal oxides. J. Chem. Phys. 2019, 150, 044702, DOI: 10.1063/1.507868246https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXhvFOhtro%253D&md5=97b0d171cbcf065900b510bd2222a25eFirst principles calculations of optical properties for oxygen vacancies in binary metal oxidesStrand, Jack; Chulkov, Sergey K.; Watkins, Matthew B.; Shluger, Alexander L.Journal of Chemical Physics (2019), 150 (4), 044702/1-044702/10CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Using an advanced computational methodol. implemented in CP2K, a non-local PBE0-TC-LRC d. functional and the recently implemented linear response formulation of the Time-dependent D. Functional Theory equations, we test the interpretation of the optical absorption and photoluminescence signatures attributed by previous exptl. and theor. studies to O-vacancies in two widely used oxides-cubic MgO and monoclinic (m)-HfO2. The results obtained in large periodic cells including up to 1000 atoms emphasize the importance of accurate predictions of defect-induced lattice distortions. They confirm that optical transitions of O-vacancies in 0, +1, and +2 charge states in MgO all have energies close to 5 eV. We test the models of photoluminescence of O-vacancies proposed in the literature. The photoluminescence of V+2O centers in m-HfO2 is predicted to peak at 3.7 eV and originate from radiative tunneling transition between a V+1O center and a self-trapped hole created by the 5.2 eV excitation. (c) 2019 American Institute of Physics.
- 47Poli, E.; Elliott, J. D.; Chulkov, S. K.; Watkins, M. B.; Teobaldi, G. The Role of Cation-Vacancies for the Electronic and Optical Properties of Aluminosilicate Imogolite Nanotubes: A Non-local, Linear-Response TDDFT Study. Front. Chem. 2019, 7, 210, DOI: 10.3389/fchem.2019.0021047https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXhtlOnu7jK&md5=139ba03de31f5c581350d6fd5a16473dCation-vacancies for the electronic and optical properties of aluminosilicate imogolite nanotubes and a non-local, linear-response TDDFT studyPoli, Emiliano; Elliott, Joshua D.; Chulkov, Sergey K.; Watkins, Matthew B.; Teobaldi, GilbertoFrontiers in Chemistry (Lausanne, Switzerland) (2019), 7 (), 210CODEN: FCLSAA; ISSN:2296-2646. (Frontiers Media S.A.)We report a combined non-local (PBE-TC-LRC) D. Functional Theory (DFT) and linear-response time-dependent DFT (LR-TDDFT) study of the structural, electronic, and optical properties of the cation-vacancy based defects in aluminosilicate (AlSi) imogolite nanotubes (Imo-NTs) that have been recently proposed on the basis of NMR (NMR) expts. Following numerical detn. of the smallest AlSi Imo-NT model capable of accommodating the defect-induced relaxation with negligible finite-size errors, we analyze the defect-induced structural deformations in the NTs and ensuing changes in the NTs' electronic structure. The NMR-derived defects are found to introduce both shallow and deep occupied states in the pristine NTs' band gap (BG). These BG states are found to be highly localized at the defect site. No empty defect-state is modeled for any of the considered systems. LR-TDDFT simulation of the defects reveal increased low-energy optical absorbance for all but one defects, with the appearance of optically active excitations at energies lower than for the defect-free NT. These results enable interpretation of the low-energy tail in the exptl. UV-vis spectra for AlSi NTs as being due to the defects. Finally, the PBE-TC-LRC-approximated exciton binding energy for the defects' optical transitions is found to be substantially lower (up to 0.8 eV) than for the pristine defect-free NT's excitations (1.1 eV).
- 48Sternheimer, R. On Nuclear Quadrupole Moments. Phys. Rev. 1951, 84, 244– 253, DOI: 10.1103/PhysRev.84.24448https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaG38XhtlOltw%253D%253D&md5=dccfdcd6ac94c5ac879db40925286c28Nuclear quadrupole momentsSternheimer, R.Physical Review (1951), 84 (), 244-53CODEN: PHRVAO; ISSN:0031-899X.cf. C.A. 45, 443i. The correction to nuclear quadrupole moments on account of the quadrupole moment induced in the electron shells was obtained by solving the Schroedinger equation for the perturbed core wave functions for Li, Al, and Cl. The correction factor by which the average 〈l/r3〉 over the valence electron function in the equation for the quadrupole coupling should be multiplied to take account of the induced effect is 1.11, 0.83, and 0.68, resp. The previously described Thomas-Fermi calcn. of this effect was carried out for F, Cl, Cu, Br, Y, Ag, I, La, Lu, Pt, Tl, At, and U.
- 49Baroni, S.; de Gironcoli, S.; Dal Corso, A.; Giannozzi, P. Phonons and related crystal properties from density-functional perturbation theory. Rev. Mod. Phys. 2001, 73, 515– 562, DOI: 10.1103/RevModPhys.73.51549https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3MXlvFKrtLc%253D&md5=20ea8e1535ceb775168384a30fa2846dPhonons and related crystal properties from density-functional perturbation theoryBaroni, Stefano; De Gironcoli, Stefano; Dal Corso, Andrea; Giannozzi, PaoloReviews of Modern Physics (2001), 73 (2), 515-562CODEN: RMPHAT; ISSN:0034-6861. (American Physical Society)This article reviews with many refs. the current status of lattice-dynamical calcns. in crystals, using d.-functional perturbation theory, with emphasis on the plane-wave pseudopotential method. Several specialized topics are treated, including the implementation for metals, the calcn. of the response to macroscopic elec. fields and their relevance to long-wavelength vibrations in polar materials, the response to strain deformations, and higher-order responses. The success of this methodol. is demonstrated with a no. of applications existing in the literature.
- 50Sternheimer, R. M. Electronic Polarizabilities of Ions from the Hartree-Fock Wave Functions. Phys. Rev. 1954, 96, 951– 968, DOI: 10.1103/PhysRev.96.95150https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaG2MXhsVOktw%253D%253D&md5=629629b68b30bee1c7fe6de9cca98a06Electron polarizabilities of ions from the Hartree-Fock wave functionsSternheimer, R. M.Physical Review (1954), 96 (), 951-68CODEN: PHRVAO; ISSN:0031-899X.cf. ibid. 95, 655. The electronic polarizability α was calcd. for several ions by obtaining the perturbation of the wave functions by an external field from a numerical solution of the differential equation satisfied by the perturbation. For the He-like ions an analytic solution was obtained by using the wave functions of L.ovrddot.owdin (C.A. 47, 9135c). The calcd. values of α are, in general, 1 to 1.5 times the observed values. For several ions values were calcd. for the quadrupole polarizability, which measures the quadrupole moment induced in the ion by an external charge. The effect of the dipole moment induced in the ion on the elec. field at the nucleus is discussed.
- 51Budzak, S.; Scalmani, G.; Jacquemin, D. Accurate Excited-State Geometries: A CASPT2 and Coupled-Cluster Reference Database for Small Molecules. J. Chem. Theory Comput. 2017, 13, 6237– 6252, DOI: 10.1021/acs.jctc.7b0092151https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhvVWlsb3I&md5=9b7945dbdfd31d2b2055a7efbe1eacfbAccurate Excited-State Geometries: A CASPT2 and Coupled-Cluster Reference Database for Small MoleculesBudzak, Simon; Scalmani, Giovanni; Jacquemin, DenisJournal of Chemical Theory and Computation (2017), 13 (12), 6237-6252CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)We present an investigation of the excited-state structural parameters detd. for a large set of small compds. with the dual goals of defining ref. values for further works and assessing the quality of the geometries obtained with relatively cheap computational approaches. In the first stage, we compare the excited-state geometries obtained with ADC(2), CC2, CCSD, CCSDR(3), CC3, and CASPT2 and large at. basis sets. It is found that CASPT2 and CC3 results are generally in very good agreement with one another (typical differences of ca. 3 × 10-3 Å) when all electrons are correlated and when the aug-cc-pVTZ at. basis set is employed with both methods. In a second stage, a statistical anal. reveals that, on the one hand, the excited-state (ES) bond lengths are much more sensitive to the selected level of theory than their ground-state (GS) counterparts and, on the other hand, that CCSDR(3) is probably the most cost-effective method delivering accurate structures. Indeed, CCSD tends to provide too compact multiple bond lengths on an almost systematic basis, whereas both CC2 and ADC(2) tend to exaggerate these bond distances, with more erratic error patterns, esp. for the latter method. The deviations are particularly marked for the polarized CO and CN bonds, as well as for the puckering angle in formaldehyde homologues. In the last part of this contribution, we provide a series of CCSDR(3) GS and ES geometries of medium-sized mols. to be used as refs. in further investigations.
- 52Ongari, D.; Yakutovich, A. V.; Talirz, L.; Smit, B. Building a Consistent and Reproducible Database for Adsorption Evaluation in Covalent–Organic Frameworks. ACS Cent. Sci. 2019, 5, 1663– 1675, DOI: 10.1021/acscentsci.9b0061952https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXhvVeqsL%252FM&md5=8f47a5e3d00532e7064d5259e73ab1f1Building a Consistent and Reproducible Database for Adsorption Evaluation in Covalent-Organic FrameworksOngari, Daniele; Yakutovich, Aliaksandr V.; Talirz, Leopold; Smit, BerendACS Central Science (2019), 5 (10), 1663-1675CODEN: ACSCII; ISSN:2374-7951. (American Chemical Society)We present a workflow that traces the path from the bulk structure of a cryst. material to assessing its performance in carbon capture from coal's postcombustion flue gases. This workflow is applied to a database of 324 covalent-org. frameworks (COFs) reported in the literature, to characterize their CO2 adsorption properties using the following steps: (1) optimization of the crystal structure (at. positions and unit cell) using d. functional theory, (2) fitting at. point charges based on the electron d., (3) characterizing the pore geometry of the structures before and after optimization, (4) computing carbon dioxide and nitrogen isotherms using grand canonical Monte Carlo simulations with an empirical interaction potential, and finally, (5) assessing the CO2 parasitic energy via process modeling. The full workflow has been encoded in the Automated Interactive Infrastructure and Database for Computational Science (AiiDA). Both the workflow and the automatically generated provenance graph of our calcns. are made available on the Materials Cloud, allowing peers to inspect every input parameter and result along the workflow, download structures and files at intermediate stages, and start their research right from where this work has left off. In particular, our set of CURATED (Clean, Uniform, and Refined with Automatic Tracking from Exptl. Database) COFs, having optimized geometry and high-quality DFT-derived point charges, are available for further investigations of gas adsorption properties. We plan to update the database as new COFs are being reported. An automated and reproducible computational workflow is proposed, to systematically optimize the geometry of covalent-org. frameworks and evaluate their performances for carbon capture and storage.
- 53Ongari, D.; Talirz, L.; Smit, B. Too Many Materials and Too Many Applications: An Experimental Problem Waiting for a Computational Solution. ACS Cent. Sci. 2020, 6, 1890– 1900, DOI: 10.1021/acscentsci.0c0098853https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXhvFOgtL3N&md5=56a120a637317936ec763bc398505969Too Many Materials and Too Many Applications: An Experimental Problem Waiting for a Computational SolutionOngari, Daniele; Talirz, Leopold; Smit, BerendACS Central Science (2020), 6 (11), 1890-1900CODEN: ACSCII; ISSN:2374-7951. (American Chemical Society)A review. Finding the best material for a specific application is the ultimate goal of materials discovery. However, there is also the reverse problem: when exptl. groups discover a new material, they would like to know all the possible applications this material would be promising for. Computational modeling can aim to fulfill this expectation, thanks to the sustained growth of computing power and the collective engagement of the scientific community in developing more efficient and accurate workflows for predicting materials' performances. We discuss the impact that reproducibility and automation of the modeling protocols have on the field of gas adsorption in nanoporous crystals. We envision a platform that combines these tools and enables effective matching between promising materials and industrial applications. We identify the opportunity for a computational platform for matching nanoporous materials and gas-related applications, motivating the development of automated and reproducible computational workflows.
- 54Hirata, S.; Head-Gordon, M. Time-dependent density functional theory within the Tamm-Dancoff approximation. Chem. Phys. Lett. 1999, 314, 291– 299, DOI: 10.1016/S0009-2614(99)01149-554https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK1MXnsl2rsr0%253D&md5=1b62f410de6c2a2193f1011d42f389c5Time-dependent density functional theory within the Tamm-Dancoff approximationHirata, S.; Head-Gordon, M.Chemical Physics Letters (1999), 314 (3,4), 291-299CODEN: CHPLBC; ISSN:0009-2614. (Elsevier Science B.V.)A computationally simple method for mol. excited states, namely, the Tamm-Dancoff approxn. to time-dependent d. functional theory, is proposed and implemented. This method yields excitation energies for several closed- and open-shell mols. that are essentially of the same quality as those obtained from time-dependent d. functional theory itself, when the same exchange-correlation functional is used.
- 55Casida, M. E. In Recent Advances in Density Functional Methods, Part I; Chong, D. P., Ed.; World Scientific: Singapore, 1995; p 155.There is no corresponding record for this reference.
- 56Guidon, M.; Hutter, J.; VandeVondele, J. Auxiliary Density Matrix Methods for Hartree-Fock Exchange Calculations. J. Chem. Theory Comput. 2010, 6, 2348– 2364, DOI: 10.1021/ct100222556https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXot1als7Y%253D&md5=7a84c086d49f17bab239fd3f8d004c6eAuxiliary Density Matrix Methods for Hartree-Fock Exchange CalculationsGuidon, Manuel; Hutter, Jurg; Vande Vondele, JoostJournal of Chemical Theory and Computation (2010), 6 (8), 2348-2364CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)The calcn. of Hartree-Fock exchange (HFX) is computationally demanding for large systems described with high-quality basis sets. In this work, we show that excellent performance and good accuracy can nevertheless be obtained if an auxiliary d. matrix is employed for the HFX calcn. Several schemes to derive an auxiliary d. matrix from a high-quality d. matrix are discussed. Key to the accuracy of the auxiliary d. matrix methods (ADMM) is the use of a correction based on std. generalized gradient approxns. for HFX. ADMM integrates seamlessly in existing HFX codes and, in particular, can be employed in linear scaling implementations. Demonstrating the performance of the method, the effect of HFX on the structure of liq. water is investigated in detail using Born-Oppenheimer mol. dynamics simulations (300 ps) of a system of 64 mols. Representative for large systems are calcns. on a solvated protein (Rubredoxin), for which ADMM outperforms the corresponding std. HFX implementation by approx. a factor 20.
- 57Ghosh, D. C.; Islam, N. Whether electronegativity and hardness are manifest two different descriptors of the one and the same fundamental property of atoms? @ TA quest. Int. J. Quantum Chem. 2011, 111, 40– 51, DOI: 10.1002/qua.2241557https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXhtlaqt7%252FE&md5=e601f3a812a1e9f147bfeb376c699815Whether electronegativity and hardness are manifest two different descriptors of the one and the same fundamental property of atoms-A questGhosh, Dulal C.; Islam, NazmulInternational Journal of Quantum Chemistry (2011), 111 (1), 40-51CODEN: IJQCB2; ISSN:0020-7608. (John Wiley & Sons, Inc.)In this report, we have tried to reveal that there is much conceptual commonality between the two fundamental theor. descriptors of chem. and physics-the electronegativity and the hardness. The phys. hardness was introduced and theorized by condensed matter physicists. The chem. hardness was introduced by chemists to generalize and rationalize the HSAB principle. We have tried to establish that the phys. hardness and the chem. hardness with evolution of time have converged to one and the same general principle-the hardness. We have also tried to understand the phys. basis and operational significance of another very important descriptor arising out of theor. constructs of chem.-the electronegativity. We have relied upon the fact that, since these descriptors are not observables, there is no possibility of their quantum mech. evaluation. These descriptors, therefore, should be and must be reified before suggesting ansatz for their evaluation. We have dwelt at length upon the effort of d. functional definition and evaluation of electronegativity and hardness and discovered the inherent inner contradiction of the theory and measurement. We have also noted that a good no. of scientists hold the opinion that the d. functional formula of electronegativity is χ = I and that of hardness is η = I, where I is the ionization potential of the chem. system. This study concludes that the two fundamental descriptors-the hardness and the electronegativity originate from the same source-the electron attracting power of the screened nucleus upon valence electrons and discovers the surprising result that if one measures hardness, the electronegativity is simultaneously measured and vice versa. We have also explored the ansatz of semiempirical evaluation of electronegativity and hardness of at. systems involving the radius of the atoms. We have noted that the ansatz for electronegativity and hardness is the same. To justify our statement that if one measures hardness, the electronegativity is simultaneously measured and vice versa, we have used the evaluated set of hardness for 103 elements of the periodic table as a scale of electronegativity and found that such set of hardness data satisfies the sine qua non of a satisfactory scale of electronegativity. The electronegativity and the hardness are two different appearances of the one and the same fundamental property of atoms. They are different and also nondifferent. They are different in their fields of application. They are nondifferent when we discuss the basic philosophical structures of their origin and development. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2011.
- 58Marx, D.; Hutter, J. In Modern Methods and Algorithms of Quantum Chemistry; Grotendorst, J., Ed.; John von Neumann Institute for Computing, Forschungszentrum Jülich: Jülich, Germany, 2000; pp 301– 449 (first edition, paperback, ISBN 3-00-005618-1) or pp 329–477 (second edition, hardcover, ISBN 3-00-005834-6), see http://www.theochem.rub.de/go/cprev.html.There is no corresponding record for this reference.
- 59Essmann, U.; Perera, L.; Berkowitz, M. L.; Darden, T.; Lee, H.; Pedersen, L. G. A smooth particle mesh Ewald method. J. Chem. Phys. 1995, 103, 8577– 8593, DOI: 10.1063/1.47011759https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK2MXptlehtrw%253D&md5=092a679dd3bee08da28df41e302383a7A smooth particle mesh Ewald methodEssmann, Ulrich; Perera, Lalith; Berkowitz, Max L.; Darden, Tom; Lee, Hsing; Pedersen, Lee G.Journal of Chemical Physics (1995), 103 (19), 8577-93CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The previously developed particle mesh Ewald method is reformulated in terms of efficient B-spline interpolation of the structure factors. This reformulation allows a natural extension of the method to potentials of the form 1/rp with p ≥ 1. Furthermore, efficient calcn. of the virial tensor follows. Use of B-splines in the place of Lagrange interpolation leads to analytic gradients as well as a significant improvement in the accuracy. The authors demonstrate that arbitrary accuracy can be achieved, independent of system size N, at a cost that scales as N log(N). For biomol. systems with many thousands of atoms and this method permits the use of Ewald summation at a computational cost comparable to that of a simple truncation method of 10 Å or less.
- 60Mayhall, N. J.; Raghavachari, K.; Hratchian, H. P. ONIOM-based QM:QM electronic embedding method using Löwdin atomic charges: Energies and analytic gradients. J. Chem. Phys. 2010, 132, 114107, DOI: 10.1063/1.331541760https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXjslWhsb4%253D&md5=2ac557eea7ce3b5887cf31841ac6508dONIOM-based QM:QM electronic embedding method using Lowdin atomic charges: Energies and analytic gradientsMayhall, Nicholas J.; Raghavachari, Krishnan; Hratchian, Hrant P.Journal of Chemical Physics (2010), 132 (11), 114107/1-114107/6CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We report a new quantum mech.:quantum mech. (QM:QM) method which provides explicit electronic polarization of the high-level region by using the Lowdin at. charges from the low-level region. This provides an embedding potential which naturally evolves with changes in nuclear geometry. However, this coupling of the high-level and low-level regions introduces complications in the energy gradient evaluation. Following previous work, we derive and implement efficient gradients where a single set of SCF response equations is solved. We provide results for the calcn. of deprotonation energies of a hydroxylated spherosiloxane cluster (Si8O12H7OH) and the dissocn. energy of a water mol. from a [ZnIm3(H2O)]2+ complex. The Lowdin charge embedding model provides results which are not only an improvement over mech. embedding (no electronic embedding) but which are also resistant to large overpolarization effects which occur more often with Mulliken charge embedding. Finally, a scaled-Loewdin charge embedding method is also presented which provides a method for fine tuning the extent of electronic polarization. (c) 2010 American Institute of Physics.
- 61Davidson, E. R. The iterative calculation of a few of the lowest eigenvalues and corresponding eigenvectors of large real-symmetric matrices. J. Comput. Phys. 1975, 17, 87– 94, DOI: 10.1016/0021-9991(75)90065-0There is no corresponding record for this reference.
- 62Chassaing, T. Excitation Energy Calculations with TD-DFT. PhD thesis, University of Zurich, Zurich, Switzerland, 2005.There is no corresponding record for this reference.
- 63Ryabinkin, I. G.; Nagesh, J.; Izmaylov, A. F. Fast Numerical Evaluation of Time-Derivative Nonadiabatic Couplings for Mixed Quantum–Classical Methods. J. Phys. Chem. Lett. 2015, 6, 4200– 4203, DOI: 10.1021/acs.jpclett.5b0206263https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2MXhs1Sis7%252FF&md5=2a46a6400734c540a5a72ed2c8b24057Fast Numerical Evaluation of Time-Derivative Nonadiabatic Couplings for Mixed Quantum-Classical MethodsRyabinkin, Ilya G.; Nagesh, Jayashree; Izmaylov, Artur F.Journal of Physical Chemistry Letters (2015), 6 (21), 4200-4203CODEN: JPCLCD; ISSN:1948-7185. (American Chemical Society)We have developed a numerical differentiation scheme that eliminates evaluation of overlap determinants in calcg. the time-deriv. nonadiabatic couplings (TDNACs). Evaluation of these determinants was the bottleneck in previous implementations of mixed quantum-classical methods using numerical differentiation of electronic wave functions in the Slater determinant representation. The central idea of our approach is, first, to reduce the analytic time derivs. of Slater determinants to time derivs. of MOs and then to apply a finite-difference formula. Benchmark calcns. prove the efficiency of the proposed scheme showing impressive several-order-of-magnitude speedups of the TDNAC calcn. step for midsize mols.
- 64Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Petersson, G. A.; Nakatsuji, H.; Li, X.; Caricato, M.; Marenich, A. V.; Bloino, J.; Janesko, B. G.; Gomperts, R.; Mennucci, B.; Hratchian, H. P.; Ortiz, J. V.; Izmaylov, A. F.; Sonnenberg, J. L.; Williams-Young, D.; Ding, F.; Lipparini, F.; Egidi, F.; Goings, J.; Peng, B.; Petrone, A.; Henderson, T.; Ranasinghe, D.; Zakrzewski, V. G.; Gao, J.; Rega, N.; Zheng, G.; Liang, W.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Throssell, K.; Montgomery, J. A., Jr.; Peralta, J. E.; Ogliaro, F.; Bearpark, M. J.; Heyd, J. J.; Brothers, E. N.; Kudin, K. N.; Staroverov, V. N.; Keith, T. A.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A. P.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Millam, J. M.; Klene, M.; Adamo, C.; Cammi, R.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Farkas, O.; Foresman, J. B.; Fox, D. J. Gaussian 16, Revision C.01; Gaussian Inc.: Wallingford, CT, 2016.There is no corresponding record for this reference.
- 65Weigend, F.; Ahlrichs, R. Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. Phys. Chem. Chem. Phys. 2005, 7, 3297, DOI: 10.1039/b508541a65https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2MXpsFWgu7o%253D&md5=a820fb6055c993b50c405ba0fc62b194Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracyWeigend, Florian; Ahlrichs, ReinhartPhysical Chemistry Chemical Physics (2005), 7 (18), 3297-3305CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)Gaussian basis sets of quadruple zeta valence quality for Rb-Rn are presented, as well as bases of split valence and triple zeta valence quality for H-Rn. The latter were obtained by (partly) modifying bases developed previously. A large set of more than 300 mols. representing (nearly) all elements-except lanthanides-in their common oxidn. states was used to assess the quality of the bases all across the periodic table. Quantities investigated were atomization energies, dipole moments and structure parameters for Hartree-Fock, d. functional theory and correlated methods, for which we had chosen Moller-Plesset perturbation theory as an example. Finally recommendations are given which type of basis set is used best for a certain level of theory and a desired quality of results.
- 66Balasubramani, S. G.; Chen, G. P.; Coriani, S.; Diedenhofen, M.; Frank, M. S.; Franzke, Y. J.; Furche, F.; Grotjahn, R.; Harding, M. E.; Hättig, C.; Hellweg, A.; Helmich-Paris, B.; Holzer, C.; Huniar, U.; Kaupp, M.; Marefat Khah, A.; Karbalaei Khani, S.; Müller, h.; Mack, F.; Nguyen, B. D.; Parker, S. M.; Perlt, E.; Rappoport, D.; Reiter, K.; Roy, S.; Rückert, M.; Schmitz, G.; Sierka, M.; Tapavicza, E.; Tew, D. P.; van Wüllen, C.; Voora, V. K.; Weigend, F.; Wodyński, A.; Yu, J. M. TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations. J. Chem. Phys. 2020, 152, 184107, DOI: 10.1063/5.000463566https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXps1Ogs7s%253D&md5=836d350d0cdc56aeea644742a098b23bTURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulationsBalasubramani, Sree Ganesh; Chen, Guo P.; Coriani, Sonia; Diedenhofen, Michael; Frank, Marius S.; Franzke, Yannick J.; Furche, Filipp; Grotjahn, Robin; Harding, Michael E.; Hattig, Christof; Hellweg, Arnim; Helmich-Paris, Benjamin; Holzer, Christof; Huniar, Uwe; Kaupp, Martin; Marefat Khah, Alireza; Karbalaei Khani, Sarah; Muller, Thomas; Mack, Fabian; Nguyen, Brian D.; Parker, Shane M.; Perlt, Eva; Rappoport, Dmitrij; Reiter, Kevin; Roy, Saswata; Ruckert, Matthias; Schmitz, Gunnar; Sierka, Marek; Tapavicza, Enrico; Tew, David P.; van Wullen, Christoph; Voora, Vamsee K.; Weigend, Florian; Wodynski, Artur; Yu, Jason M.Journal of Chemical Physics (2020), 152 (18), 184107CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)A review. TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chem. simulations of mols., clusters, periodic systems, and solns. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy-cost ratio, such as d. functional theory including local hybrids and the RPA, GW-Bethe-Salpeter methods, second-order Moller-Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resoln.-of-the-identity approxn., imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent addns. to TURBOMOLE's functionality, including excited-state methods, RPA and Green's function methods, relativistic approaches, high-order mol. properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE's current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE's development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted. (c) 2020 American Institute of Physics.
- 67Weigend, F.; Furche, F.; Ahlrichs, R. Gaussian basis sets of quadruple zeta valence quality for atoms H-Kr. J. Chem. Phys. 2003, 119, 12753– 12762, DOI: 10.1063/1.162729367https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3sXpvFOns7s%253D&md5=9d1664e0004c3a900b7efcb5da10dd6dGaussian basis sets of quadruple zeta valence quality for atoms H-KrWeigend, Florian; Furche, Filipp; Ahlrichs, ReinhartJournal of Chemical Physics (2003), 119 (24), 12753-12762CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We present Gaussian basis sets of quadruple zeta valence quality with a segmented contraction scheme for atoms H to Kr. This extends earlier work on segmented contracted split valence (SV) and triple zeta valence (TZV) basis sets. Contraction coeffs. and orbital exponents are fully optimized in at. Hartree-Fock (HF) calcns. As opposed to other quadruple zeta basis sets, the basis set errors in at. ground-state HF energies are less than 1 mEh and increase smoothly across the Periodic Table, while the no. of primitives is comparably small. Polarization functions are taken partly from previous work, partly optimized in at. MP2 calcns., and for a few cases detd. at the HF level for excited at. states nearly degenerate with the ground state. This leads to basis sets denoted QZVP for HF and d. functional theory (DFT) calcns., and for some atoms to a larger basis recommended for correlated treatments, QZVPP. We assess the performance of the basis sets in mol. HF, DFT, and MP2 calcns. for a sample of diat. and small polyat. mols. by a comparison of energies, bond lengths, and dipole moments with results obtained numerically or using very large basis sets. It is shown that basis sets of quadruple zeta quality are necessary to achieve an accuracy of 1 kcal/mol per bond in HF and DFT atomization energies. For compds. contg. third row as well as alk. and earth alk. metals it is demonstrated that the inclusion of high-lying core orbitals in the active space can be necessary for accurate correlated treatments. The QZVPP basis sets provide sufficient flexibility to polarize the core in those cases. All test calcns. indicate that the new basis sets lead to consistent accuracies in HF, DFT, or correlated treatments even in crit. cases where other basis sets may show deficiencies.
- 68VandeVondele, J.; Hutter, J. Gaussian basis sets for accurate calculations on molecular systems in gas and condensed phases. J. Chem. Phys. 2007, 127, 114105, DOI: 10.1063/1.277070868https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2sXhtFSrsLvM&md5=d7fdb937efb88cf3fca85792bb49ec27Gaussian basis sets for accurate calculations on molecular systems in gas and condensed phasesVandeVondele, Joost; Hutter, JurgJournal of Chemical Physics (2007), 127 (11), 114105/1-114105/9CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We present a library of Gaussian basis sets that has been specifically optimized to perform accurate mol. calcns. based on d. functional theory. It targets a wide range of chem. environments, including the gas phase, interfaces, and the condensed phase. These generally contracted basis sets, which include diffuse primitives, are obtained minimizing a linear combination of the total energy and the condition no. of the overlap matrix for a set of mols. with respect to the exponents and contraction coeffs. of the full basis. Typically, for a given accuracy in the total energy, significantly fewer basis functions are needed in this scheme than in the usual split valence scheme, leading to a speedup for systems where the computational cost is dominated by diagonalization. More importantly, binding energies of hydrogen bonded complexes are of similar quality as the ones obtained with augmented basis sets, i.e., have a small (down to 0.2 kcal/mol) basis set superposition error, and the monomers have dipoles within 0.1 D of the basis set limit. However, contrary to typical augmented basis sets, there are no near linear dependencies in the basis, so that the overlap matrix is always well conditioned, also, in the condensed phase. The basis can therefore be used in first principles mol. dynamics simulations and is well suited for linear scaling calcns.
- 69Brémond, E.; Savarese, M.; Adamo, C.; Jacquemin, D. Accuracy of TD-DFT Geometries: A Fresh Look. J. Chem. Theory Comput. 2018, 14, 3715– 3727, DOI: 10.1021/acs.jctc.8b0031169https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXhtFSrtL%252FE&md5=bb253fe791ccd2ba3b30042b9b9db4c0Accuracy of TD-DFT Geometries: A Fresh LookBremond, Eric; Savarese, Marika; Adamo, Carlo; Jacquemin, DenisJournal of Chemical Theory and Computation (2018), 14 (7), 3715-3727CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)We benchmark a panel of 48 DFT exchange-correlation functionals in the framework of TD-DFT optimizations of the geometry of valence singlet excited states. To this end, we use a set of 41 small- and medium-sized org. mols. for which ref. geometries were obtained at high level of theory, typically, CC3 or CCSDR(3), with the aug-cc-pVTZ at. basis set. For the ground-state parameters, the tested functionals provide av. deviations that are small (0.010 Å and 0.5° for bond lengths and valence angles) and not very sensitive to the selected (hybrid) functional, but the errors are larger for the most polarized bonds (CO, CN, and so on). Nevertheless, DFT has a tendency to provide too compact distances, a trend slightly enhanced for functionals including a large amt. of exact exchange. The av. errors largely increase when going to the excited-state for most bond types, i.e., TD-DFT delivers less accurate excited-state distances than DFT for ground state. In particular TD-DFT combined with hybrid functionals provides significantly too short CO and CS/CSe bonds with resp. av. errors in the -0.026/-0.052 Å and -0.015/-0.082 Å ranges, depending on the selected hybrid functional. For the carbonyl bonds, the sizes of the TD-DFT deviations obtained when selecting std. hybrid functionals are of the same order of magnitude as the EOM-CCSD ones.
- 70Chan, G. K.-L.; Handy, N. C. C8H8: a density functional theory study of molecular geometries introducing the localised bond density. J. Chem. Soc., Faraday Trans. 1996, 92, 3015– 3021, DOI: 10.1039/ft996920301570https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK28XlsVGhtbo%253D&md5=9c09739d70a3d81a2a7c13aee9c25355C8H8: A density functional theory study of molecular geometries introducing the localized bond densityChan, Garnet K-L.; Handy, Nicholas C.Journal of the Chemical Society, Faraday Transactions (1996), 92 (17), 3015-3021CODEN: JCFTEV; ISSN:0956-5000. (Royal Society of Chemistry)D. functional theory was used with all the common exchange-correlation functionals to investigate the structures of three isomers of C8H8 found in F.A. Cotton's text, barrelene, cyclooctatetraene, tetramethylenecyclobutane and also ethane and ethene. All calcns. were performed with TZ2P basis sets and large quadrature. The results are compared with expt. and those obtained with Hartree-Fock theory. Delocalization in the three mols. is discussed. A localized bond d. is introduced to explain the transferability of the trends in the predictions of the functionals between different mols. Three-parameter adiabatic connection functionals are examd. and their usefulness in geometry prediction questioned. Finally a phys. picture of the correlation as modeled by d. functional theory is presented and used to explain trends in the overestimation or underestimation of bond lengths.
- 71Neumann, R.; Nobes, R. H.; Handy, N. C. Exchange functionals and potentials. Mol. Phys. 1996, 87, 1– 36, DOI: 10.1080/0026897960010001171https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK28XhvFegs7c%253D&md5=14030dab7e55c7c6faaf9737dbcbcb37Exchange functionals and potentialsNeumann, Ralf; Nobes, Ross H.; Handy, Nicholas C.Molecular Physics (1996), 87 (1), 1-36CODEN: MOPHAM; ISSN:0026-8976. (Taylor & Francis)The commonly used exchange-correlation functionals of d. functional theory and their potentials are examd. numerically following the first such investigation of Perdew. They are also investigated for Ne and Kr. Their behavior for large gradients of the d. and for large distances is not satisfactory. In particular, the correct asymptotic r-1 behavior is difficult to achieve. Following van Leeuwen and Baerends, this is linked to the energy εmax of the highest occupied orbital arising from the Kohn-Sham equations. This deficiency is linked also with the poor prediction of mol. polarizabilities. The Becke-Roussel (BR) exchange functional is examd., which is derived by assuming a hydrogen-like exchange hole at all spatial points, and it has the attraction of being dependent on both the kinetic energy d. and the Laplacian of the d. and has no adjustable parameters. Becke has presented encouraging results using this functional in a hybrid manner. Fully self-consistent Kohn-sham calcns. are performed using it in combination with Perdew's 1986 correlation functional. The results are very encouraging indeed, so much so that this exchange functional is the best generalized gradient approxn. (GGA) yet discovered. In particular, bond lengths of many mol. show a substantial improvement over results from other GGAs. For example, many CH bonds are now within exptl. accuracy, instead of being typically 0.02 Å too long. Our ab initio understanding of non-dynamic correlation and dynamic correlation is then linked with d. functional theory. It is argued that correlation functionals should pick up the local dynamic correlation, whereas exchange functionals should include non-dynamic correlation effects. For these reasons it is considered that exchange functionals are best modeled on a system for which there is effectively no non-dynamic correlation, for which the optimum example is the Ne atom. Thus, again following Becke and Roussel, the spherically averaged Hartree-Fock exchange hole for Ne is examd., compared with the BR model functional hole. An excellent overlap is found, and thus the above good results are explained. As a final contribution, the dissocn. of the H2 mol. is re-examd., looking at it in terms of the exchange hole. For a ref. electron near one proton A, the RHF model has half an exchange electron near it, and half the exchange electron near the other proton B, whereas the BR functional has one electron near the other proton b, whereas the Br functional has one electron near A, which is the correct picture. For this reason the (restricted) BR functional gives a greatly improved dissocn. curve for H2 when compared with the Hartree-Fock curve. In summary, the Becke-Roussel functional is found to be a most attractive exchange functional.
- 72Sancho-García, J. C.; Cornil, J. Assessment of recently developed exchange-correlation functionals for the description of torsion potentials in π-conjugated molecules. J. Chem. Phys. 2004, 121, 3096– 3101, DOI: 10.1063/1.177497672https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2cXmt1Slsrg%253D&md5=e1684f5d68a4a19ab35897b8e1f18dd0Assessment of recently developed exchange-correlation functionals for the description of torsion potentials in π-conjugated moleculesSancho-Garcia, Juan Carlos; Cornil, JeromeJournal of Chemical Physics (2004), 121 (7), 3096-3101CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Newly developed exchange-correlation functionals in d. functional theory (DFT) have been applied to describe conjugation effects in org. mols. The performance of the various approaches is assessed through the calcn. of torsion energy profiles and their crit. comparison with available exptl. data. Our results indicate that the OPTX-B95 exchange-correlation functional as well as its corresponding hybrid versions perform better than the well-established BLYP or B3LYP schemes when dealing with π-conjugated mols. In contrast, the recently introduced VSXC functional is not as reliable as other DFT methods for the systems examd. here.
- 73Vikramaditya, T.; Lin, S.-T. Limitations of Global Hybrids in Predicting the Geometries and Torsional Energy Barriers of Dimeric Systems and the Role of Hartree Fock and DFT Exchange. J. Comput. Chem. 2019, 40, 2810– 2818, DOI: 10.1002/jcc.2605673https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXhs1CitrvE&md5=9dc520039c1dbe19cee658bdef98f2f0Limitations of Global Hybrids in Predicting the Geometries and Torsional Energy Barriers of Dimeric Systems and the Role of Hartree Fock and DFT ExchangeVikramaditya, Talapunur; Lin, Shiang-TaiJournal of Computational Chemistry (2019), 40 (32), 2810-2818CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)Prediction of accurate geometries is a prerequisite for accurate prediction of mol. properties. Impact of Hartree Fock (HF) exchange (a0) on geometry in the framework of DFT is investigated by monitoring dihedral angles, bond length alternations, and torsional energy barriers of 10 dimeric systems against CCSD (ADZ/ATZ) benchmarks. A strong correlation is obsd. between the fraction of HF exchange, equil. dihedral angles, and the potential energy barriers in global hybrids. Full HF exchange is crit. to accurately predict the nonplanarity. Lower fractions of (a0)/larger DFT exchange (1-a0) results in overestimation of torsional energy barriers at 900 and underestimation at 00. Large contributions of (1-a0) in global hybrid functionals tend to overestimate torsional energy barriers (900) and are biased toward planar geometries. However, inclusion of larger fractions of (a0)/lower (1-a0) also overestimate the torsional energy barriers in syn-conformations due to the localization errors assocd. with HF exchange in global hybrids. Hence, irresp. of the fraction of HF/DFT exchange incorporated, global hybrids fail to accurately predict torsional energy barriers at 00 and 900 simultaneously. Long-range cor. (LC) functionals, which employ full HF exchange at longer regions, outperform global hybrid functionals in predicting geometries and torsional energy barriers of the dimeric mols. The distance dependence of (a0) thus provides a balanced fraction of HF exchange as the dihedral torsion varies. Impact of range sepn. parameter on geometries is marginal in altering the planarity/nonplanarity. However, range sepn. parameter within 0.20-0.40 bohr-1 predicts more reliable torsional energies and geometries. © 2019 Wiley Periodicals, Inc.
- 74Fang, C.; Oruganti, B.; Durbeej, B. How Method-Dependent Are Calculated Differences between Vertical, Adiabatic, and 0─0 Excitation Energies?. J. Phys. Chem. A 2014, 118, 4157– 4171, DOI: 10.1021/jp501974p74https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXotlWht7Y%253D&md5=ca7e486acb318eb2955f3b03fd91675cHow Method-Dependent Are Calculated Differences between Vertical, Adiabatic, and 0-0 Excitation Energies?Fang, Changfeng; Oruganti, Baswanth; Durbeej, BoJournal of Physical Chemistry A (2014), 118 (23), 4157-4171CODEN: JPCAFH; ISSN:1089-5639. (American Chemical Society)Through a large no. of benchmark studies, the performance of different quantum chem. methods in calcg. vertical excitation energies is today quite well established. These efforts have in recent years been complemented by a few benchmarks focusing instead on adiabatic excitation energies. However, it is much less well established how calcd. differences between vertical, adiabatic and 0-0 excitation energies vary between methods, which may be due to the cost of evaluating zero-point vibrational energy corrections for excited states. To fill this gap, we have calcd. vertical, adiabatic, and 0-0 excitation energies for a benchmark set of mols. covering both org. and inorg. systems. Considering in total 96 excited states and using both TD-DFT with a variety of exchange-correlation functionals and the ab initio CIS and CC2 methods, it is found that while the vertical excitation energies obtained with the various methods show an av. (over the 96 states) std. deviation of 0.39 eV, the corresponding std. deviations for the differences between vertical, adiabatic, and 0-0 excitation energies are much smaller: 0.10 (difference between adiabatic and vertical) and 0.02 eV (difference between 0-0 and adiabatic). These results provide a quant. measure showing that the calcn. of such quantities in photochem. modeling is well amenable to low-level methods. In addn., we also report on how these energy differences vary between chem. systems and assess the performance of TD-DFT, CIS, and CC2 in reproducing exptl. 0-0 excitation energies.
- 75Peach, M. J. G.; Benfield, P.; Helgaker, T.; Tozer, D. J. Excitation energies in density functional theory: An evaluation and a diagnostic test. J. Chem. Phys. 2008, 128, 044118, DOI: 10.1063/1.283190075https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXhslOmsLo%253D&md5=6360f5b16a0fac8535963503c82d96d4Excitation energies in density functional theory: An evaluation and a diagnostic testPeach, Michael J. G.; Benfield, Peter; Helgaker, Trygve; Tozer, David J.Journal of Chemical Physics (2008), 128 (4), 044118/1-044118/8CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Electronic excitation energies are detd. using the CAM-B3LYP Coulomb-attenuated functional, together with a std. generalized gradient approxn. (GGA) and hybrid functional. The degree of spatial overlap between the occupied and virtual orbitals involved in an excitation is measured using a quantity Λ, and the extent to which excitation energy errors correlate with Λ is quantified. For a set of 59 excitations of local, Rydberg, and intramol. charge-transfer character in 18 theor. challenging main-group mols., CAM-B3LYP provides by far the best overall performance; no correlation is obsd. between excitation energy errors and Λ, reflecting the good quality, balanced description of all three categories of excitation. By contrast, a clear correlation is obsd. for the GGA and, to a lesser extent, the hybrid functional, allowing a simple diagnostic test to be proposed for judging the reliability of a general excitation from these functionals-when Λ falls below a prescribed threshold, excitations are likely to be in very significant error. The study highlights the ambiguous nature of the term "charge transfer," providing insight into the observation that while many charge-transfer excitations are poorly described by GGA and hybrid functionals, others are accurately reproduced. (c) 2008 American Institute of Physics.
- 76Send, R.; Kühn, M.; Furche, F. Assessing Excited State Methods by Adiabatic Excitation Energies. J. Chem. Theory Comput. 2011, 7, 2376– 2386, DOI: 10.1021/ct200272b76https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXptVOiurk%253D&md5=06f4095c2f002c0081e026db938f9fa1Assessing Excited State Methods by Adiabatic Excitation EnergiesSend, Robert; Kuhn, Michael; Furche, FilippJournal of Chemical Theory and Computation (2011), 7 (8), 2376-2386CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)We compile a 109-membered benchmark set of adiabatic excitation energies (AEEs) from high-resoln. gas-phase expts. Our data set includes a variety of org. chromophores with up to 46 atoms, radicals, and inorg. transition metal compds. Many of the 91 mols. in our set are relevant to atm. chem., photovoltaics, photochem., and biol. The set samples valence, Rydberg, and ionic states of various spin multiplicities. As opposed to vertical excitation energies, AEEs are rigorously defined by energy differences of vibronic states, directly observable, and insensitive to errors in equil. structures. We supply optimized ground state and excited state structures, which allows fast and convenient evaluation of AEEs with two single-point energy calcns. per system. We apply our benchmark set to assess the performance of time-dependent d. functional theory using common semilocal functionals and the CI singles method. Hybrid functionals such as B3LYP and PBE0 yield the best results, with mean abs. errors around 0.3 eV. We also investigate basis set convergence and correlations between different methods and between the magnitude of the excited state relaxation energy and the AEE error. A smaller, 15-membered subset of AEEs is introduced and used to assess the correlated wave function methods CC2 and ADC(2). These methods improve upon hybrid TDDFT for systems with single-ref. ground states but perform less well for radicals and small-gap transition metal compds. None of the investigated methods reaches "chem. accuracy" of 0.05 eV in AEEs.
- 77Côté, A. P.; Benin, A. I.; Ockwig, N. W.; O’Keeffe, M.; Matzger, A. J.; Yaghi, O. M. Porous, Crystalline, Covalent Organic Frameworks. Science 2005, 310, 1166– 1170, DOI: 10.1126/science.112041177https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2MXht1WgtLnF&md5=e2e033fd4af37870fc73cbb18c1eaed1Porous, Crystalline, Covalent Organic FrameworksCote, Adrien P.; Benin, Annabelle I.; Ockwig, Nathan W.; O'Keeffe, Michael; Matzger, Adam J.; Yaghi, Omar M.Science (Washington, DC, United States) (2005), 310 (5751), 1166-1170CODEN: SCIEAS; ISSN:0036-8075. (American Association for the Advancement of Science)Covalent org. frameworks (COFs) have been designed and successfully synthesized by condensation reactions of Ph diboronic acid {C6H4[B(OH)2]2} and hexahydroxytriphenylene [C18H6(OH)6]. Powder x-ray diffraction studies of the highly cryst. products (C3H2BO)6•(C9H12)1 (COF-1) and C9H4BO2 (COF-5) revealed expanded porous graphitic layers that are either staggered (COF-1, P63/mmc) or eclipsed (COF-5, P6/mmm). Their crystal structures are entirely held by strong bonds between B, C, and O atoms to form rigid porous architectures with pore sizes ranging from 7 to 27 angstroms. COF-1 and COF-5 exhibit high thermal stability (to temps. up to 500° to 600°C), permanent porosity, and high surface areas (711 and 1590 square meters per g, resp.).
- 78Das, G.; Biswal, B. P.; Kandambeth, S.; Venkatesh, V.; Kaur, G.; Addicoat, M.; Heine, T.; Verma, S.; Banerjee, R. Chemical sensing in two dimensional porous covalent organic nanosheets. Chem. Sci. 2015, 6, 3931– 3939, DOI: 10.1039/C5SC00512D78https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2MXntlWktLg%253D&md5=f2f12867a8edaafb56789f5184445ce6Chemical sensing in two dimensional porous covalent organic nanosheetsDas, Gobinda; Biswal, Bishnu P.; Kandambeth, Sharath; Venkatesh, V.; Kaur, Gagandeep; Addicoat, Matthew; Heine, Thomas; Verma, Sandeep; Banerjee, RahulChemical Science (2015), 6 (7), 3931-3939CODEN: CSHCCN; ISSN:2041-6520. (Royal Society of Chemistry)Two new imide-based cryst., porous, and chem. stable covalent org. frameworks (COFs) (TpBDH and TfpBDH) have been successfully synthesized employing solvothermal crystn. route. Furthermore, thin layered covalent org. nanosheets (CONs) were derived from these bulk COFs by the simple liq. phase exfoliation method. These 2D CONs showcase increased luminescence intensity compared to their bulk counterparts (COFs). Notably, TfpBDH-CONs showcase good selectivity and prominent, direct visual detection towards different nitroarom. analytes over TpBDH-CONs. Quite interestingly, TfpBDH-CONs exhibit a superior "turn-on" detection capability for 2,4,6-trinitrophenol (TNP) in the solid state, but conversely, they also show a "turn-off" detection in the dispersion state. These findings describe a new approach towards developing an efficient, promising fluorescence chemosensor material for both visual and spectroscopic detection of nitroarom. compds. with very low [10-5 (M)] analyte concns.
- 79Das, P.; Chakraborty, G.; Tyagi, S.; Mandal, S. K. Design of Fluorescent and Robust Covalent Organic Framework Host Matrices for Illuminating Mechanistic Insight into Solvatochromic Decoding. ACS Appl. Mater. Interfaces 2020, 12, 52527– 52537, DOI: 10.1021/acsami.0c1278579https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXit12rsrbE&md5=550f96e9c2920b9eaa975ee6e6fe1045Design of Fluorescent and Robust Covalent Organic Framework Host Matrices for Illuminating Mechanistic Insight into Solvatochromic DecodingDas, Prasenjit; Chakraborty, Gouri; Tyagi, Sparsh; Mandal, Sanjay K.ACS Applied Materials & Interfaces (2020), 12 (47), 52527-52537CODEN: AAMICK; ISSN:1944-8244. (American Chemical Society)Two functional covalent org. frameworks (COFs), constructed from 3-connected triazine-based amine or hydrazine with linear dialdehyde, are decorated with mol. docking sites to showcase their solvatochromic decoding behavior toward volatile solvent mols. (VSMs). These luminescent and cryst. COFs, namely, COF-N and COF-NN, are characterized by numerous anal. techniques. After accommodation of different VSMs as guests, the inclusion compds. of COF-N and COF-NN display solvatochromism. More fascinatingly, the singlet energy, band gaps, and lifetime of these VSM-accommodated COF-N and COF-NN are linearly correlated with the properties of VSMs. D. functional theory (DFT) and Monte Carlo simulation studies further support the interaction of VSMs with COF-N and COF-NN. The presence of extra amine functionality in COF-NN leads to the better interaction with VSMs and, therefore, results in different modes of interaction and correlation. Considering their inestimable chem. diversity, this study introduces a new path for finely tuned solvatochromic properties by COFs.
Supporting Information
Supporting Information
The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jctc.2c00144.
Treatment of the virtual space, structural information for molecular benchmark set, geometrical data on optimized molecular geometries, statistical analysis for ADMM and sTDA vertical excitation energies, adiabatic excitation energies and statistical analysis, fluorescence energies and statistical analysis, and structural information for covalent organic frameworks (PDF)
Terms & Conditions
Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.