1
|
Chen X, Cifuentes-Lopez A, Shao X, Lin L, Prokopchuk D, Pavanello M. Unraveling the Hydration Shell Structure and Dynamics of Group 10 Aqua Ions. J Phys Chem Lett 2024; 15:5517-5528. [PMID: 38749061 DOI: 10.1021/acs.jpclett.4c00464] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/24/2024]
Abstract
We present ab initio simulations based on subsystem DFT of group 10 aqua ions accurately compared against experimental data on hydration structure. Our simulations provide insights into the molecular structures and dynamics of hydration shells, offering recalibrated interpretations of experimental results. We observe a soft, but distinct second hydration shell in Palladium (Pd) due to a balance between thermal fluctuations, metal-water interactions, and hydrogen bonding. Nickel (Ni) and platinum (Pt) exhibit more rigid hydration shells. Notably, our simulations align with experimental findings for Pd, showing axial hydration marked by a broad peak at about 3 Å in the Pd-O radial distribution function, revising the previously sharp "mesoshell" prediction. We introduce the "hydrogen bond dome" concept to describe a resilient network of hydrogen-bonded water molecules around the metal, which plays a critical role in the axial hydration dynamics.
Collapse
Affiliation(s)
- Xin Chen
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
- Department of Physics, Rutgers University, Newark, New Jersey 07102, United States
| | | | - Xuecheng Shao
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| | - Lirong Lin
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| | - Demyan Prokopchuk
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| | - Michele Pavanello
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
- Department of Physics, Rutgers University, Newark, New Jersey 07102, United States
| |
Collapse
|
2
|
Mi W, Luo K, Trickey SB, Pavanello M. Orbital-Free Density Functional Theory: An Attractive Electronic Structure Method for Large-Scale First-Principles Simulations. Chem Rev 2023; 123:12039-12104. [PMID: 37870767 DOI: 10.1021/acs.chemrev.2c00758] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2023]
Abstract
Kohn-Sham Density Functional Theory (KSDFT) is the most widely used electronic structure method in chemistry, physics, and materials science, with thousands of calculations cited annually. This ubiquity is rooted in the favorable accuracy vs cost balance of KSDFT. Nonetheless, the ambitions and expectations of researchers for use of KSDFT in predictive simulations of large, complicated molecular systems are confronted with an intrinsic computational cost-scaling challenge. Particularly evident in the context of first-principles molecular dynamics, the challenge is the high cost-scaling associated with the computation of the Kohn-Sham orbitals. Orbital-free DFT (OFDFT), as the name suggests, circumvents entirely the explicit use of those orbitals. Without them, the structural and algorithmic complexity of KSDFT simplifies dramatically and near-linear scaling with system size irrespective of system state is achievable. Thus, much larger system sizes and longer simulation time scales (compared to conventional KSDFT) become accessible; hence, new chemical phenomena and new materials can be explored. In this review, we introduce the historical contexts of OFDFT, its theoretical basis, and the challenge of realizing its promise via approximate kinetic energy density functionals (KEDFs). We review recent progress on that challenge for an array of KEDFs, such as one-point, two-point, and machine-learnt, as well as some less explored forms. We emphasize use of exact constraints and the inevitability of design choices. Then, we survey the associated numerical techniques and implemented algorithms specific to OFDFT. We conclude with an illustrative sample of applications to showcase the power of OFDFT in materials science, chemistry, and physics.
Collapse
Affiliation(s)
- Wenhui Mi
- Key Laboratory of Material Simulation Methods & Software of Ministry of Education, College of Physics, Jilin University, Changchun 130012, PR China
- State Key Laboratory of Superhard Materials, Jilin University, Changchun 130012, PR China
- International Center of Future Science, Jilin University, Changchun 130012, PR China
| | - Kai Luo
- Department of Applied Physics, Nanjing University of Science and Technology, Nanjing 210094, PR China
| | - S B Trickey
- Quantum Theory Project, Department of Physics and Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States
| | - Michele Pavanello
- Department of Physics and Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| |
Collapse
|
3
|
Jansen M, Reinholdt P, Hedegård ED, König C. Theoretical and Numerical Comparison of Quantum- and Classical Embedding Models for Optical Spectra. J Phys Chem A 2023. [PMID: 37399130 DOI: 10.1021/acs.jpca.3c02540] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/05/2023]
Abstract
Quantum-mechanical (QM) and classical embedding models approximate a supermolecular quantum-chemical calculation. This is particularly useful when the supermolecular calculation has a size that is out of reach for present QM models. Although QM and classical embedding methods share the same goal, they approach this goal from different starting points. In this study, we compare the polarizable embedding (PE) and frozen-density embedding (FDE) models. The former is a classical embedding model, whereas the latter is a density-based QM embedding model. Our comparison focuses on solvent effects on optical spectra of solutes. This is a typical scenario where super-system calculations including the solvent environment become prohibitively large. We formulate a common theoretical framework for PE and FDE models and systematically investigate how PE and FDE approximate solvent effects. Generally, differences are found to be small, except in cases where electron spill-out becomes problematic in the classical frameworks. In these cases, however, atomic pseudopotentials can reduce the electron-spill-out issue.
Collapse
Affiliation(s)
- Marina Jansen
- Institute of Physical Chemistry and Electrochemistry, Leibniz University Hannover, Callinstr. 3A, 30167 Hannover, Germany
| | - Peter Reinholdt
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, 5230 Odense, Denmark
| | - Erik D Hedegård
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, 5230 Odense, Denmark
| | - Carolin König
- Institute of Physical Chemistry and Electrochemistry, Leibniz University Hannover, Callinstr. 3A, 30167 Hannover, Germany
| |
Collapse
|
4
|
Belleflamme F, Hehn AS, Iannuzzi M, Hutter J. A variational formulation of the Harris functional as a correction to approximate Kohn-Sham density functional theory. J Chem Phys 2023; 158:054111. [PMID: 36754794 DOI: 10.1063/5.0122671] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023] Open
Abstract
Accurate descriptions of intermolecular interactions are of great importance in simulations of molecular liquids. We present an electronic structure method that combines the accuracy of the Harris functional approach with the computational efficiency of approximately linear-scaling density functional theory (DFT). The non-variational nature of the Harris functional has been addressed by constructing a Lagrangian energy functional, which restores the variational condition by imposing stationarity with respect to the reference density. The associated linear response equations may be treated with linear-scaling efficiency in an atomic orbital based scheme. Key ingredients to describe the structural and dynamical properties of molecular systems are the forces acting on the atoms and the stress tensor. These first-order derivatives of the Harris Lagrangian have been derived and implemented in consistence with the energy correction. The proposed method allows for simulations with accuracies close to the Kohn-Sham DFT reference. Embedded in the CP2K program package, the method is designed to enable ab initio molecular dynamics simulations of molecular solutions for system sizes of several thousand atoms. Available subsystem DFT methods may be used to provide the reference density required for the energy correction at near linear-scaling efficiency. As an example of production applications, we applied the method to molecular dynamics simulations of the binary mixtures cyclohexane-methanol and toluene-methanol, performed within the isobaric-isothermal ensemble, to investigate the hydrogen bonding network in these non-ideal mixtures.
Collapse
Affiliation(s)
- Fabian Belleflamme
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| | - Anna-Sophia Hehn
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| | - Marcella Iannuzzi
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| | - Jürg Hutter
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| |
Collapse
|
5
|
Treß RS, Liu J, Hättig C, Höfener S. Pushing the limits: Efficient wavefunction methods for excited states in complex systems using frozen-density embedding. J Chem Phys 2022; 157:204101. [DOI: 10.1063/5.0100393] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022] Open
Abstract
Frozen density embedding (FDE) is an embedding method for complex environments that is simple for users to set up. It reduces the computation time by dividing the total system into small subsystems and approximating the interaction by a functional of their densities. Its combination with wavefunction methods is, however, limited to small- or medium-sized molecules because of the steep scaling in computation time of these methods. To mitigate this limitation, we present a combination of the FDE approach with pair natural orbitals (PNOs) in the TURBOMOLE software package. It combines the uncoupled FDE (FDEu) approach for excitation energy calculations with efficient implementations of second-order correlation methods in the ricc2 and pnoccsd programs. The performance of this combination is tested for tetraazaperopyrene (TAPP) molecular crystals. It is shown that the PNO truncation error on environment-induced shifts is significantly smaller than the shifts themselves and, thus, that the local approximations of PNO-based wavefunction methods can without the loss of relevant digits be combined with the FDE method. Computational wall times are presented for two TAPP systems. The scaling of the wall times is compared to conventional supermolecular calculations and demonstrates large computational savings for the combination of FDE- and PNO-based methods. Additionally, the behavior of excitation energies with the system size is investigated. It is found that the excitation energies converge quickly with the size of the embedding environment for the TAPPs investigated in the current study.
Collapse
Affiliation(s)
- Robert S. Treß
- Department of Theoretical Chemistry, Ruhr University Bochum, 44801 Bochum, Germany
| | - Jing Liu
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, Germany
| | - Christof Hättig
- Department of Theoretical Chemistry, Ruhr University Bochum, 44801 Bochum, Germany
| | - Sebastian Höfener
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, Germany
| |
Collapse
|
6
|
Shao X, Lopez AC, Khan Musa MR, Nouri MR, Pavanello M. Adaptive Subsystem Density Functional Theory. J Chem Theory Comput 2022; 18:6646-6655. [PMID: 36179128 DOI: 10.1021/acs.jctc.2c00698] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Subsystem density functional theory (DFT) is emerging as a powerful electronic structure method for large-scale simulations of molecular condensed phases and interfaces. Key to its computational efficiency is the use of approximate nonadditive noninteracting kinetic energy functionals. Unfortunately, currently available nonadditive functionals lead to inaccurate results when the subsystems interact strongly such as when they engage in chemical reactions. This work disrupts the status quo by devising a workflow that extends subsystem DFT's applicability also to strongly interacting subsystems. This is achieved by implementing a fully automated adaptive definition of subsystems which is realized during geometry optimizations or ab initio molecular dynamics simulations. The new method prescribes subsystem merging and splitting events redistributing the resources (both for work and data) in an efficient way making use of modern parallelization strategies and object-oriented programming. We showcase the method with examples probing from moderate-to-strong inter-subsystem interactions, opening the door to using subsystem DFT for modeling chemical reactions in molecular condensed phases with a black box computational tool.
Collapse
Affiliation(s)
- Xuecheng Shao
- Department of Chemistry, Rutgers University, Newark, New Jersey07102, United States
| | | | - Md Rajib Khan Musa
- Department of Chemistry, Rutgers University, Newark, New Jersey07102, United States
| | - Mohammad Reza Nouri
- Department of Chemistry, Rutgers University, Newark, New Jersey07102, United States
| | - Michele Pavanello
- Department of Physics, Rutgers University, Newark, New Jersey07102, United States
| |
Collapse
|
7
|
Shao X, Mi W, Pavanello M. Density Embedding Method for Nanoscale Molecule-Metal Interfaces. J Phys Chem Lett 2022; 13:7147-7154. [PMID: 35901490 DOI: 10.1021/acs.jpclett.2c01424] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
In this work, we extend the applicability of standard Kohn-Sham DFT (KS-DFT) to model realistically sized molecule-metal interfaces where the metal slabs venture into the tens of nanometers in size. Employing state-of-the-art noninteracting kinetic energy functionals, we describe metallic subsystems with orbital-free DFT and combine their electronic structure with molecular subsystems computed at the KS-DFT level resulting in a multiscale subsystem DFT method. The method reproduces within a few millielectronvolts the binding energy difference of water and carbon dioxide molecules adsorbed on the top and hollow sites of an Al(111) surface compared to KS-DFT of the combined supersystem. It is also robust for Born-Oppenheimer molecular dynamics simulations. Very large system sizes are approached with standard computing resources thanks to a parallelization scheme that avoids accumulation of memory at the gather-scatter stage. The results as presented are encouraging and open the door to ab initio simulations of realistically sized, mesoscopic molecule-metal interfaces.
Collapse
Affiliation(s)
- Xuecheng Shao
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| | - Wenhui Mi
- International Center for Computational Method and Software, College of Physics, Jilin University, Changchun 130012, China
| | - Michele Pavanello
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
- Department of Physics, Rutgers University, Newark, New Jersey 07102, United States
| |
Collapse
|
8
|
Della Sala F. Orbital-Free Methods for Plasmonics: Linear Response. J Chem Phys 2022; 157:104101. [DOI: 10.1063/5.0100797] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Plasmonic systems, such as metal nanoparticles, are widely used in different application areas, going from biology to photovoltaics.The modeling of the optical response of such systems is of fundamental importance to analyze their behavior and to design new systems with required properties.When the characteristic sizes/distances reach a few nanometers, non-local and spill-out effects become relevant and conventional classical electrodynamics models are no more appropriate. Methods based on the Time-Dependent Density-Functional Theory (TD-DFT) represent the current reference for the description of quantum effects. However, TD-DFT is based on knowledge of all occupied orbitals whose calculation is computationally prohibitive to model large plasmonic systems of interest for applications.On the other hand, methods based on the Orbital-Free (OF) formulation of TD-DFT, can scale linearly with the system size.In this Review, OF methods ranging from semiclassical models to the quantum hydrodynamic theory, will be derived from the linear response TD-DFT, so that the key approximations and properties of each method can be clearly highlighted. The accuracy of the various approximations will be then validated for the linear optical properties of jellium nanoparticles, the most relevant model system in plasmonics. OF methods can describe the collective excitations in plasmonic systems with great accuracy andwithout system-tuned parameters. The accuracy on these methods depends only on the accuracy on the (universal) kinetic energy functional of the ground-state electronic density. Current approximations and future development directions will be indicated.
Collapse
Affiliation(s)
- Fabio Della Sala
- CNR-IMM, IMM CNR Lecce, Italy
- Istituto Italiano di Tecnologia Center for Biomolecular Nanotechnologies
| |
Collapse
|
9
|
Chakravarty C, Aksu H, Martinez B JA, Ramos P, Pavanello M, Dunietz BD. Role of Dielectric Screening in Calculating Excited States of Solvated Azobenzene: A Benchmark Study Comparing Quantum Embedding and Polarizable Continuum Model for Representing the Solvent. J Phys Chem Lett 2022; 13:4849-4855. [PMID: 35617015 DOI: 10.1021/acs.jpclett.2c00982] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
The low energy excited states of the conformational isomers of solvated azobenzene are calculated with several DFT methods accounting for the solute-solvent interaction implicitly with the polarizable continuum model or explicitly with subsystem DFT. For the latter, embedding potentials are calculated for 21 sampled snapshots of the solvent molecules. First, we find that accounting for the solvent implicitly or explicitly has little effect on the predicted cis-trans S1 excitation energy gap. Second, we find that azobenzene's S1 cis and trans energies are accurate as long as a screened range-separated hybrid exchange-correlation functional is employed. Finally, we also tested a simplified workflow whereby a single, averaged, embedding potential is used. Unfortunately, we find larger deviations against the experiment for the simplified workflow. This highlights a basic flaw in the approach, where the time scale of solvent averaging is much longer than that of the solute's electronic polarization.
Collapse
Affiliation(s)
- Chandrima Chakravarty
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44242, United States
| | - Huseyin Aksu
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44242, United States
- Computational Physics Laboratory, Department of Physics, Pamukkale University, 20010 Denizli, Turkey
| | - Jessica A Martinez B
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| | - Pablo Ramos
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| | - Michele Pavanello
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| | - Barry D Dunietz
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44242, United States
| |
Collapse
|
10
|
Revisiting the trapping of noble gases (He-Kr) by the triatomic H 3+ and Li 3+ species: a density functional reactivity theory study. J Mol Model 2022; 28:122. [PMID: 35437635 DOI: 10.1007/s00894-022-05099-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2021] [Accepted: 03/22/2022] [Indexed: 10/18/2022]
Abstract
Small atomic clusters with exotic stability, bonding, aromaticity, and reactivity properties can be made use of for various purposes. In this work, we revisit the trapping of noble gas atoms (He-Kr) by the triatomic H3+ and Li3+ species by using some analytical tools from density functional theory, conceptual density functional theory, and the information-theoretic approach. Our results showcase that though similar in geometry, H3+ and Li3+ exhibit markedly different behavior in bonding, aromaticity, and reactivity properties after the addition of noble gas atoms. Moreover, the exchange-correlation interaction and steric effect are key energy components in stabilizing the clusters. This study also finds that the origin of the molecular stability of these species is due to the spatial delocalization of the electron density distribution. Our work provides an additional arsenal towards a better understanding of small atomic clusters capturing noble gases.
Collapse
|
11
|
Xu Q, Ma C, Mi W, Wang Y, Ma Y. Nonlocal pseudopotential energy density functional for orbital-free density functional theory. Nat Commun 2022; 13:1385. [PMID: 35296665 PMCID: PMC8927098 DOI: 10.1038/s41467-022-29002-3] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2022] [Accepted: 02/22/2022] [Indexed: 11/09/2022] Open
Abstract
Orbital-free density functional theory (OF-DFT) is an electronic structure method with a low computational cost that scales linearly with the number of simulated atoms, making it suitable for large-scale material simulations. It is generally considered that OF-DFT strictly requires the use of local pseudopotentials, rather than orbital-dependent nonlocal pseudopotentials, for the calculation of electron-ion interaction energies, as no orbitals are available. This is unfortunate situation since the nonlocal pseudopotentials are known to give much better transferability and calculation accuracy than local ones. We report here the development of a theoretical scheme that allows the direct use of nonlocal pseudopotentials in OF-DFT. In this scheme, a nonlocal pseudopotential energy density functional is derived by the projection of nonlocal pseudopotential onto the non-interacting density matrix (instead of "orbitals") that can be approximated explicitly as a functional of electron density. Our development defies the belief that nonlocal pseudopotentials are not applicable to OF-DFT, leading to the creation for an alternate theoretical framework of OF-DFT that works superior to the traditional approach.
Collapse
Affiliation(s)
- Qiang Xu
- International Center for Computational Methods and Software & State Key Lab of Superhard Materials, College of Physics, Jilin University, Changchun, 130012, China
| | - Cheng Ma
- International Center for Computational Methods and Software & State Key Lab of Superhard Materials, College of Physics, Jilin University, Changchun, 130012, China
| | - Wenhui Mi
- International Center for Computational Methods and Software & State Key Lab of Superhard Materials, College of Physics, Jilin University, Changchun, 130012, China
| | - Yanchao Wang
- International Center for Computational Methods and Software & State Key Lab of Superhard Materials, College of Physics, Jilin University, Changchun, 130012, China.
| | - Yanming Ma
- International Center for Computational Methods and Software & State Key Lab of Superhard Materials, College of Physics, Jilin University, Changchun, 130012, China.,International Center of Future Science, Jilin University, Changchun, 130012, China
| |
Collapse
|
12
|
Kinetic Energy Density Functionals Based on a Generalized Screened Coulomb Potential: Linear Response and Future Perspectives. COMPUTATION 2022. [DOI: 10.3390/computation10020030] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
We consider kinetic energy functionals that depend, beside the usual semilocal quantities (density, gradient, Laplacian of the density), on a generalized Yukawa potential, that is the screened Coulomb potential of the density raised to some power. These functionals, named Yukawa generalized gradient approximations (yGGA), are potentially efficient real-space semilocal methods that include significant non-local effects and can describe different important exact properties of the kinetic energy. In this work, we focus in particular on the linear response behavior for the homogeneous electron gas (HEG). We show that such functionals are able to reproduce the exact Lindhard function behavior with a very good accuracy, outperforming all other semilocal kinetic functionals. These theoretical advances allow us to perform a detailed analysis of a special class of yGGAs, namely the linear yGGA functionals. Thus, we show how the present approach can generalize the yGGA functionals improving the HEG linear behavior and leading to an extended formula for the kinetic functional. Moreover, testing on several jellium cluster model systems allows highlighting advantages and limitations of the linear yGGA functionals and future perspectives for the development of yGGA kinetic functionals.
Collapse
|
13
|
Treß RS, Hättig C, Höfener S. Employing Pseudopotentials to Tackle Excited-State Electron Spill-Out in Frozen Density Embedding Calculations. J Chem Theory Comput 2022; 18:1737-1747. [PMID: 35107998 DOI: 10.1021/acs.jctc.1c00732] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
In frozen density embedding (FDE), the properties of a target molecule are computed in the presence of an effective embedding potential, which accounts for the attractive and repulsive contributions of the environment. The formally exact embedding potential, however, is in practice calculated using explicit kinetic-energy functionals for which the resulting potentials are in many cases not repulsive enough to account fully for Pauli repulsion by the electrons of the environment and to compensate thereby the strong electron-nuclear attraction. For the excited states on the target molecule, this leads to charge spill-out when diffuse basis functions are included, which allow that valence electrons are excited to those regions of the environment where the strong nuclear attraction is not sufficiently compensated by repulsive contributions. To reduce this insufficiency, we propose in the present work the inclusion of atomic all-electron pseudopotentials for all environment atoms on top of the conventional embedding potential. In the current work, the pseudopotentials are applied for computing vertical excitation energies of local excited states in complex systems employing the second-order algebraic diagrammatic construction (ADC(2)) scheme. The proposed approach leads to significantly reduced charge spill-out and an improved agreement of FDE and supermolecular calculations in the frozen solvent approximation. In particular, when diffuse functions are employed, the mean absolute deviation (MAD) is reduced from 0.27 to 0.05 eV for the investigated cases.
Collapse
Affiliation(s)
- Robert S Treß
- Department of Theoretical Chemistry, Ruhr University Bochum, Bochum 44801, Germany
| | - Christof Hättig
- Department of Theoretical Chemistry, Ruhr University Bochum, Bochum 44801, Germany
| | - Sebastian Höfener
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Karlsruhe 76131, Germany
| |
Collapse
|
14
|
Pauletti M, Rybkin VV, Iannuzzi M. Subsystem Density Functional Theory Augmented by a Delta Learning Approach to Achieve Kohn-Sham Accuracy. J Chem Theory Comput 2021; 17:6423-6431. [PMID: 34505765 DOI: 10.1021/acs.jctc.1c00592] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Simulations based on electronic structure theory naturally include polarization and have no transferability problems. In particular, Kohn-Sham density functional theory (KS-DFT) has become the method of reference for ab initio molecular dynamics simulations of condensed matter systems. However, the high computational cost often poses strict limits on the affordable system size as well as on the extension of sampling (number of configurations). In this work, we propose an improvement to the subsystem density functional theory approach, known as the Kim-Gordon (KG) scheme, thus enabling the sampling of configurations for condensed molecular systems keeping the KS-DFT level accuracy at a fraction of computer time. Our scheme compensates the known KG shortcomings of the electronic kinetic energy term by adding a simple correction and can match KS-DFT accuracy in energies and forces. The computationally cheap correction is determined by means of a machine learning procedure. The proposed KG scheme is applied within a linear scaling self-consistent field formalism and is assessed by a series of molecular dynamics simulations of liquid water under different conditions. Although system-dependent, the correction is transferable between system sizes and temperatures.
Collapse
Affiliation(s)
- Michela Pauletti
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, Zurich 8057, Switzerland
| | - Vladimir V Rybkin
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, Zurich 8057, Switzerland
| | - Marcella Iannuzzi
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, Zurich 8057, Switzerland
| |
Collapse
|
15
|
Giannone G, Śmiga S, D'Agostino S, Fabiano E, Della Sala F. Plasmon Couplings from Subsystem Time-Dependent Density Functional Theory. J Phys Chem A 2021; 125:7246-7259. [PMID: 34403247 DOI: 10.1021/acs.jpca.1c05384] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Many applications in plasmonics are related to the coupling between metallic nanoparticles (MNPs) or between an emitter and a MNP. The theoretical analysis of such a coupling is thus of fundamental importance to analyze the plasmonic behavior and to design new systems. While classical methods neglect quantum and spill-out effects, time-dependent density functional theory (TD-DFT) considers all of them and with Kohn-Sham orbitals delocalized over the whole system. Thus, within TD-DFT, no definite separation of the subsystems (the single MNP or the emitter) and their couplings is directly available. This important feature is obtained here using the subsystem formulation of TD-DFT, which has been originally developed in the context of weakly interacting organic molecules. In subsystem TD-DFT, interacting MNPs are treated independently, thus allowing us to compute the plasmon couplings directly from the subsystem TD-DFT transition densities. We show that subsystem TD-DFT, as well as a simplified version of it in which kinetic contributions are neglected, can reproduce the reference TD-DFT calculations for gap distances greater than about 6 Å or even smaller in the case of hybrid plasmonic systems (i.e., molecules interacting with MNPs). We also show that the subsystem TD-DFT can be also used as a tool to analyze the impact of charge-transfer effects.
Collapse
Affiliation(s)
- Giulia Giannone
- Center for Biomolecular Nanotechnologies, Istituto Italiano di Tecnologia, Via Barsanti 14, Arnesano (LE) 73010, Italy.,Department of Mathematics and Physics "E. De Giorgi", University of Salento, Via Arnesano, Lecce 73100, Italy
| | - Szymon Śmiga
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Grudzia̧dzka 5, Toruń 87-100, Poland
| | - Stefania D'Agostino
- Center for Biomolecular Nanotechnologies, Istituto Italiano di Tecnologia, Via Barsanti 14, Arnesano (LE) 73010, Italy.,Department of Mathematics and Physics "E. De Giorgi", University of Salento, Via Arnesano, Lecce 73100, Italy.,Institute of Nanotechnology, National Research Council (CNR-NANOTEC), c/o Campus Ecotekne, via Monteroni, Lecce 73100, Italy
| | - Eduardo Fabiano
- Center for Biomolecular Nanotechnologies, Istituto Italiano di Tecnologia, Via Barsanti 14, Arnesano (LE) 73010, Italy.,Institute for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, Lecce 73100, Italy
| | - Fabio Della Sala
- Center for Biomolecular Nanotechnologies, Istituto Italiano di Tecnologia, Via Barsanti 14, Arnesano (LE) 73010, Italy.,Institute for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, Lecce 73100, Italy
| |
Collapse
|
16
|
Guo C, He X, Rong C, Lu T, Liu S, Chattaraj PK. Local Temperature as a Chemical Reactivity Descriptor. J Phys Chem Lett 2021; 12:5623-5630. [PMID: 34110155 DOI: 10.1021/acs.jpclett.1c01207] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Using the electron density and its associated quantities in a molecular system to quantify chemical reactivity in density functional theory is of considerable recent interest. Local temperature based on the kinetic energy density is an intrinsic property of a molecular system, which can be employed for this purpose. In this work, we explore such a possibility. To this end, we examine the local behavior of local temperature with a few choices of the kinetic energy density, apply it to determine regioselectivity of nucleophilic and electrophilic compounds, and then investigate its performance in appreciating reactions along the intrinsic reaction pathway for exothermic, endothermic, and thermoneutral transformations. Our results confirm that local temperature can be used as an effective descriptor of molecular reactivity.
Collapse
Affiliation(s)
- Chunna Guo
- Key Laboratory of Chemical Biology and Traditional Chinese Medicine Research (Ministry of Education of China), Hunan Normal University, Changsha, Hunan 410081, China
| | - Xin He
- Key Laboratory of Chemical Biology and Traditional Chinese Medicine Research (Ministry of Education of China), Hunan Normal University, Changsha, Hunan 410081, China
| | - Chunying Rong
- Key Laboratory of Chemical Biology and Traditional Chinese Medicine Research (Ministry of Education of China), Hunan Normal University, Changsha, Hunan 410081, China
| | - Tian Lu
- Beijing Kein Research Center for Natural Sciences, Beijing 100022, China
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, North Carolina 27599-3420, United States
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599-3290, United States
| | | |
Collapse
|
17
|
Sen R, González-Espinoza CE, Zech A, Dreuw A, Wesolowski TA. Benchmark of the Extension of Frozen-Density Embedding Theory to Nonvariational Correlated Methods: The Embedded-MP2 Case. J Chem Theory Comput 2021; 17:4049-4062. [PMID: 34137597 DOI: 10.1021/acs.jctc.1c00228] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The extension of the frozen-density embedding theory for nonvariational methods [J. Chem. Theory Comput. 2020, 16, 6880] was utilized to evaluate intermolecular interaction energies for complexes in the Zhao-Truhlar basis set. In the applied method (FDET-MP2-FAT-LDA), the same auxiliary system is used to evaluate the correlation energy by means of the second-order Møller-Plesset perturbation theory (MP2), as in our previous work [J. Chem. Phys. 2019, 150, 121101]. Local density approximation is used for ExcTnad[ρA,ρB] in both cases. Additionally, the contribution to the energy due to the neglected correlation potential was evaluated and analyzed. The domain of applicability of the local density approximation for ExcTnad[ρA,ρB] was determined based on deviations from the interaction energies from the conventional MP2 calculations. The local density approximation for ExcTnad[ρA,ρB] performs well for hydrogen- or dipole-bound complexes. The relative errors in the interaction energy lie within 3-30%. While for charge-transfer complexes, this approximation fails consistently, and for other types of complexes, the performance of this approximation is not systematic. The sources of error are discussed in detail.
Collapse
Affiliation(s)
- Reena Sen
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University Heidelberg, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | | | - Alexander Zech
- Department of Physical Chemistry, University of Geneva, 30, Quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University Heidelberg, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Tomasz A Wesolowski
- Department of Physical Chemistry, University of Geneva, 30, Quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| |
Collapse
|
18
|
Shao X, Mi W, Pavanello M. GGA-Level Subsystem DFT Achieves Sub-kcal/mol Accuracy Intermolecular Interactions by Mimicking Nonlocal Functionals. J Chem Theory Comput 2021; 17:3455-3461. [PMID: 33983729 DOI: 10.1021/acs.jctc.1c00283] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
The key feature of nonlocal kinetic energy functionals is their ability to reduce to the Thomas-Fermi functional in the regions of high density and to the von Weizsäcker functional in the region of low-density/high reduced density gradient. This behavior is crucial when these functionals are employed in subsystem DFT simulations to approximate the nonadditive kinetic energy. We propose a GGA nonadditive kinetic energy functional which mimics the good behavior of nonlocal functionals, retaining the computational complexity of typical semilocal functionals. Crucially, this functional depends on the inter-subsystem density overlap. The new functional reproduces Kohn-Sham DFT and benchmark CCSD(T) interaction energies of weakly interacting dimers in the S22-5 and S66 test sets with a mean absolute deviation well below 1 kcal/mol.
Collapse
Affiliation(s)
- Xuecheng Shao
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| | - Wenhui Mi
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| | - Michele Pavanello
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States.,Department of Physics, Rutgers University, Newark, New Jersey 07102, United States
| |
Collapse
|
19
|
Francisco HI, Carmona-Espíndola J, Gázquez JL. Analysis of the kinetic energy functional in the generalized gradient approximation. J Chem Phys 2021; 154:084107. [PMID: 33639771 DOI: 10.1063/5.0040973] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A new density functional for the total kinetic energy in the generalized gradient approximation is developed through an enhancement factor that leads to the correct behavior in the limits when the reduced density gradient tends to 0 and to infinity and by making use of the conjoint conjecture for the interpolation between these two limits, through the incorporation, in the intermediate region of constraints that are associated with the exchange energy functional. The resulting functional leads to a reasonable description of the kinetic energies of atoms and molecules when it is used in combination with Hartree-Fock densities. Additionally, in order to improve the behavior of the kinetic energy density, a new enhancement factor for the Pauli kinetic energy is proposed by incorporating the correct behavior into the limits when the reduced density gradient tends to 0 and to infinity, together with the positivity condition, and imposing through the interpolation function that the sum of its integral over the whole space and the Weiszacker energy must be equal to the value obtained with the enhancement factor developed for the total kinetic energy.
Collapse
Affiliation(s)
- Héctor I Francisco
- Departamento de Química, Universidad Autónoma Metropolitana-Iztapalapa, Av. San Rafael Atlixco 186, México, Ciudad de México 09340, Mexico
| | - Javier Carmona-Espíndola
- Departamento de Química, CONACYT-Universidad Autónoma Metropolitana-Iztapalapa, Av. San Rafael Atlixco 186, Ciudad de México 09340, Mexico
| | - José L Gázquez
- Departamento de Química, Universidad Autónoma Metropolitana-Iztapalapa, Av. San Rafael Atlixco 186, México, Ciudad de México 09340, Mexico
| |
Collapse
|
20
|
Rana B, Herbert JM. Role of hemibonding in the structure and ultraviolet spectroscopy of the aqueous hydroxyl radical. Phys Chem Chem Phys 2020; 22:27829-27844. [DOI: 10.1039/d0cp05216g] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The presence of a two-center, three-electron hemibond in the solvation structure of the aqueous hydroxl radical has long been debated, as its appearance can be sensitive to self-interaction error in density functional theory.
Collapse
Affiliation(s)
- Bhaskar Rana
- Department of Chemistry & Biochemistry
- The Ohio State University
- Columbus
- USA
| | - John M. Herbert
- Department of Chemistry & Biochemistry
- The Ohio State University
- Columbus
- USA
| |
Collapse
|
21
|
Mi W, Ramos P, Maranhao J, Pavanello M. Ab Initio Structure and Dynamics of CO 2 at Supercritical Conditions. J Phys Chem Lett 2019; 10:7554-7559. [PMID: 31738567 DOI: 10.1021/acs.jpclett.9b03054] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Green technologies rely on green solvents and fluids. Among them, supercritical CO2 already finds many important applications. The molecular-level understanding of the dynamics and structure of this supercritical fluid is a prerequisite for rational design of future green technologies. Unfortunately, the commonly employed Kohn-Sham density functional theory (DFT) is too computationally demanding to produce meaningfully converged dynamics within a reasonable time and with a reasonable computational effort. Thanks to subsystem DFT, we analyze finite-size effects by considering simulation cells of varying sizes (up to 256 independent molecules in the cell) and finite-time effects by running 100 ps trajectories. We find that the simulations are in reasonable and semiquantitative agreement with the available neutron diffraction experiments and that, as opposed to the gas phase, the CO2 molecules in the fluid are bent with an average OCO angle of 175.8°. Our simulations also confirm that the dimer T-shape is the most prevalent configuration. Our results further strengthen the experiment-simulation agreement for this fluid when comparing radial distribution functions and diffusion coefficient, confirming subsystem DFT as a viable tool for modeling structure and dynamics of condensed-phase systems.
Collapse
Affiliation(s)
- Wenhui Mi
- Department of Chemistry , Rutgers University , Newark , New Jersey 07102 , United States
| | - Pablo Ramos
- Department of Chemistry , Rutgers University , Newark , New Jersey 07102 , United States
| | - Jack Maranhao
- Department of Chemistry , Rutgers University , Newark , New Jersey 07102 , United States
| | - Michele Pavanello
- Department of Chemistry , Rutgers University , Newark , New Jersey 07102 , United States
- Department of Physics , Rutgers University , Newark , New Jersey 07102 , United States
| |
Collapse
|
22
|
The Role of the Reduced Laplacian Renormalization in the Kinetic Energy Functional Development. COMPUTATION 2019. [DOI: 10.3390/computation7040065] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
The Laplacian of the electronic density diverges at the nuclear cusp, which complicates the development of Laplacian-level meta-GGA (LLMGGA) kinetic energy functionals for all-electron calculations. Here, we investigate some Laplacian renormalization methods, which avoid this divergence. We developed two different LLMGGA functionals, which improve the kinetic energy or the kinetic potential. We test these KE functionals in the context of Frozen-Density-Embedding (FDE), for a large palette of non-covalently interacting molecular systems. These functionals improve over the present state-of-the-art LLMGGA functionals for the FDE calculations.
Collapse
|
23
|
Apostolidou C. OH radical in water from ab initio molecular dynamics simulation employing hybrid functionals. J Chem Phys 2019. [DOI: 10.1063/1.5107479] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Affiliation(s)
- Christina Apostolidou
- Mulliken Center for Theoretical Chemistry, Institute of Physical and Theoretical Chemistry, Rheinische Friedrich-Wilhelms-Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| |
Collapse
|
24
|
Liu S, Zhao D, Rong C, Lu T, Liu S. Using Pauli energy to appraise the quality of approximate semilocal non-interacting kinetic energy density functionals. J Chem Phys 2019; 150:204106. [PMID: 31153167 DOI: 10.1063/1.5095072] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
Abstract
It is well-known that the kinetic energy density (KED) functional is the most difficult to approximate in density functional theory (DFT), yet to take full advantage of DFT with its density-based descriptive capability of molecular properties, an accurate account of KED is a must. To have a better idea of how an approximate KED formula behaves and where we should focus in the future development of better approximate KEDs, in this work we propose to employ the Pauli energy to assess their quality. We tested the performance of a total of 22 approximate semilocal noninteracting KED functionals from the literature for 18 neutral atoms and 20 small molecules. We found that generalized gradient approximation formulas of the KED functional can often reasonably accurately predict the total kinetic energy value for atoms and molecules but failed miserably to forecast the integrated values for Pauli energy related properties. The reason behind this is that presently available approximate KED functionals are unable to accurately account for the kinetic energy distribution in the medium range away from nuclei, where the Pauli energy plays a crucial role. Our results strongly suggest that the key information missing in approximate KED functionals comes from the medium regions, not nuclear cusps nor asymptotic areas, and the Pauli energy is a reliable measure of the quality of approximate KED functionals. Future efforts in developing better KED approximations should be invested in the regions of molecules where chemical bonds are formed in order to accurately account for the Pauli energy.
Collapse
Affiliation(s)
- Siyuan Liu
- Key Laboratory of Chemical Biology and Traditional Chinese Medicine Research (Ministry of Education of China), Synergetic Innovation Center for Quantum Effects and Applications (SICQEA), College of Chemistry and Chemical Engineering, Hunan Normal University, Changsha, Hunan 410081, China
| | - Dongbo Zhao
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Chunying Rong
- Key Laboratory of Chemical Biology and Traditional Chinese Medicine Research (Ministry of Education of China), Synergetic Innovation Center for Quantum Effects and Applications (SICQEA), College of Chemistry and Chemical Engineering, Hunan Normal University, Changsha, Hunan 410081, China
| | - Tian Lu
- Beijing Kein Research Center for Natural Sciences, Beijing 100022, China
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, North Carolina 27599-3420, USA
| |
Collapse
|
25
|
Constantin LA, Fabiano E, Della Sala F. Performance of Semilocal Kinetic Energy Functionals for Orbital-Free Density Functional Theory. J Chem Theory Comput 2019; 15:3044-3055. [DOI: 10.1021/acs.jctc.9b00183] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Lucian A. Constantin
- Center for Biomolecular Nanotechnologies @UNILE, Istituto Italiano di Tecnologia, Via Barsanti, I-73010 Arnesano, Italy
| | - Eduardo Fabiano
- Center for Biomolecular Nanotechnologies @UNILE, Istituto Italiano di Tecnologia, Via Barsanti, I-73010 Arnesano, Italy
- Institute for Microelectronics and Microsystems (CNR-IMM), Campus Unisalento, Via Monteroni, 73100 Lecce, Italy
| | - Fabio Della Sala
- Center for Biomolecular Nanotechnologies @UNILE, Istituto Italiano di Tecnologia, Via Barsanti, I-73010 Arnesano, Italy
- Institute for Microelectronics and Microsystems (CNR-IMM), Campus Unisalento, Via Monteroni, 73100 Lecce, Italy
| |
Collapse
|
26
|
Jiang K, Nafziger J, Wasserman A. Constructing a non-additive non-interacting kinetic energy functional approximation for covalent bonds from exact conditions. J Chem Phys 2018; 149:164112. [PMID: 30384751 DOI: 10.1063/1.5051455] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Kaili Jiang
- Department of Physics and Astronomy, Purdue University, 525 Northwestern Ave., West Lafayette, Indiana 47907, USA
| | - Jonathan Nafziger
- Department of Chemistry, Purdue University, 560 Oval Dr., West Lafayette Indiana 47907, USA
| | - Adam Wasserman
- Department of Physics and Astronomy, Purdue University, 525 Northwestern Ave., West Lafayette, Indiana 47907, USA
- Department of Chemistry, Purdue University, 560 Oval Dr., West Lafayette Indiana 47907, USA
| |
Collapse
|
27
|
Seino J, Kageyama R, Fujinami M, Ikabata Y, Nakai H. Erratum: "Semi-local machine-learned kinetic energy density functional with third-order gradients of electron density" [J. Chem. Phys. 148, 241705 (2018)]. J Chem Phys 2018; 149:079901. [PMID: 30134665 DOI: 10.1063/1.5049687] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Junji Seino
- Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Ryo Kageyama
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Mikito Fujinami
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Yasuhiro Ikabata
- Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Hiromi Nakai
- Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| |
Collapse
|
28
|
Constantin LA, Fabiano E, Della Sala F. Semilocal Pauli-Gaussian Kinetic Functionals for Orbital-Free Density Functional Theory Calculations of Solids. J Phys Chem Lett 2018; 9:4385-4390. [PMID: 30019904 DOI: 10.1021/acs.jpclett.8b01926] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Kinetic energy (KE) approximations are key elements in orbital-free density functional theory. To date, the use of nonlocal functionals, possibly employing system-dependent parameters, has been considered mandatory in order to obtain satisfactory accuracy for different solid-state systems, whereas semilocal approximations are generally regarded as unfit to this aim. Here, we show that, instead, properly constructed semilocal approximations, the Pauli-Gaussian (PG) KE functionals, especially at the Laplacian level of theory, can indeed achieve similar accuracy as nonlocal functionals and can be accurate for both metals and semiconductors, without the need for system-dependent parameters.
Collapse
Affiliation(s)
- Lucian A Constantin
- Center for Biomolecular Nanotechnologies @UNILE , Istituto Italiano di Tecnologia , Via Barsanti , I-73010 Arnesano , Italy
| | - Eduardo Fabiano
- Center for Biomolecular Nanotechnologies @UNILE , Istituto Italiano di Tecnologia , Via Barsanti , I-73010 Arnesano , Italy
- Institute for Microelectronics and Microsystems (CNR-IMM) , Via Monteroni, Campus Unisalento , 73100 Lecce , Italy
| | - Fabio Della Sala
- Center for Biomolecular Nanotechnologies @UNILE , Istituto Italiano di Tecnologia , Via Barsanti , I-73010 Arnesano , Italy
- Institute for Microelectronics and Microsystems (CNR-IMM) , Via Monteroni, Campus Unisalento , 73100 Lecce , Italy
| |
Collapse
|
29
|
Seino J, Kageyama R, Fujinami M, Ikabata Y, Nakai H. Semi-local machine-learned kinetic energy density functional with third-order gradients of electron density. J Chem Phys 2018; 148:241705. [DOI: 10.1063/1.5007230] [Citation(s) in RCA: 47] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Junji Seino
- Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Ryo Kageyama
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Mikito Fujinami
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Yasuhiro Ikabata
- Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Hiromi Nakai
- Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
- CREST, Japan Science and Technology Agency, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan
- Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Katsura, Kyoto 615-8520, Japan
| |
Collapse
|
30
|
Jiang K, Nafziger J, Wasserman A. Non-additive non-interacting kinetic energy of rare gas dimers. J Chem Phys 2018; 148:104113. [PMID: 29544312 DOI: 10.1063/1.5016308] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Approximations of the non-additive non-interacting kinetic energy (NAKE) as an explicit functional of the density are the basis of several electronic structure methods that provide improved computational efficiency over standard Kohn-Sham calculations. However, within most fragment-based formalisms, there is no unique exact NAKE, making it difficult to develop general, robust approximations for it. When adjustments are made to the embedding formalisms to guarantee uniqueness, approximate functionals may be more meaningfully compared to the exact unique NAKE. We use numerically accurate inversions to study the exact NAKE of several rare-gas dimers within partition density functional theory, a method that provides the uniqueness for the exact NAKE. We find that the NAKE decreases nearly exponentially with atomic separation for the rare-gas dimers. We compute the logarithmic derivative of the NAKE with respect to the bond length for our numerically accurate inversions as well as for several approximate NAKE functionals. We show that standard approximate NAKE functionals do not reproduce the correct behavior for this logarithmic derivative and propose two new NAKE functionals that do. The first of these is based on a re-parametrization of a conjoint Perdew-Burke-Ernzerhof (PBE) functional. The second is a simple, physically motivated non-decomposable NAKE functional that matches the asymptotic decay constant without fitting.
Collapse
Affiliation(s)
- Kaili Jiang
- Department of Physics and Astronomy, Purdue University, 525 Northwestern Ave., West Lafayette, Indiana 47907, USA
| | - Jonathan Nafziger
- Department of Chemistry, Purdue University, 560 Oval Dr., West Lafayette, Indiana 47907, USA
| | - Adam Wasserman
- Department of Physics and Astronomy, Purdue University, 525 Northwestern Ave., West Lafayette, Indiana 47907, USA
| |
Collapse
|
31
|
Modified Fourth-Order Kinetic Energy Gradient Expansion with Hartree Potential-Dependent Coefficients. J Chem Theory Comput 2017; 13:4228-4239. [DOI: 10.1021/acs.jctc.7b00705] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
32
|
Śmiga S, Fabiano E, Constantin LA, Della Sala F. Laplacian-dependent models of the kinetic energy density: Applications in subsystem density functional theory with meta-generalized gradient approximation functionals. J Chem Phys 2017; 146:064105. [DOI: 10.1063/1.4975092] [Citation(s) in RCA: 39] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
|
33
|
Huang C, Muñoz-García AB, Pavone M. Effective scheme for partitioning covalent bonds in density-functional embedding theory: From molecules to extended covalent systems. J Chem Phys 2016; 145:244103. [DOI: 10.1063/1.4972012] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
34
|
Cancio AC, Redd JJ. Visualisation and orbital-free parametrisation of the large-Z scaling of the kinetic energy density of atoms. Mol Phys 2016. [DOI: 10.1080/00268976.2016.1246757] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Antonio C. Cancio
- Department of Physics and Astronomy, Ball State University , Muncie, IN, USA
| | - Jeremy J. Redd
- Department of Physics, Utah Valley University , Orem, UT, USA
| |
Collapse
|
35
|
Burke K, Cancio A, Gould T, Pittalis S. Locality of correlation in density functional theory. J Chem Phys 2016; 145:054112. [DOI: 10.1063/1.4959126] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Affiliation(s)
- Kieron Burke
- Department of Chemistry, University of California, Irvine, California 92697, USA
| | - Antonio Cancio
- Department of Physics and Astronomy, Ball State University, Muncie, Indiana 47306, USA
| | - Tim Gould
- Qld Micro- and Nanotechnology Centre, Griffith University, Nathan, Qld 4111, Australia
| | - Stefano Pittalis
- CNR-Istituto di Nanoscienze, Via Campi 213A, I-41125 Modena, Italy
| |
Collapse
|
36
|
Genova A, Ceresoli D, Pavanello M. Avoiding fractional electrons in subsystem DFT based ab-initio molecular dynamics yields accurate models for liquid water and solvated OH radical. J Chem Phys 2016; 144:234105. [DOI: 10.1063/1.4953363] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Alessandro Genova
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, USA
| | - Davide Ceresoli
- CNR-ISTM: Institute of Molecular Sciences and Technologies, Milano, Italy
| | - Michele Pavanello
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, USA
| |
Collapse
|
37
|
Huang C. Patching the Exchange-Correlation Potential in Density Functional Theory. J Chem Theory Comput 2016; 12:2224-33. [DOI: 10.1021/acs.jctc.6b00051] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Chen Huang
- Department
of Scientific
Computing, Florida State University, Tallahassee, Florida 32306-4120, United States
| |
Collapse
|
38
|
|
39
|
Genova A, Pavanello M. Exploiting the locality of periodic subsystem density-functional theory: efficient sampling of the Brillouin zone. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2015; 27:495501. [PMID: 26596499 DOI: 10.1088/0953-8984/27/49/495501] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
In order to approximately satisfy the Bloch theorem, simulations of complex materials involving periodic systems are made n(k) times more complex by the need to sample the first Brillouin zone at n(k) points. By combining ideas from Kohn-Sham density-functional theory (DFT) and orbital-free DFT, for which no sampling is needed due to the absence of waves, subsystem DFT offers an interesting middle ground capable of sizable theoretical speedups against Kohn-Sham DFT. By splitting the supersystem into interacting subsystems, and mapping their quantum problem onto separate auxiliary Kohn-Sham systems, subsystem DFT allows an optimal topical sampling of the Brillouin zone. We elucidate this concept with two proof of principle simulations: a water bilayer on Pt[1 1 1]; and a complex system relevant to catalysis-a thiophene molecule physisorbed on a molybdenum sulfide monolayer deposited on top of an α-alumina support. For the latter system, a speedup of 300% is achieved against the subsystem DTF reference by using an optimized Brillouin zone sampling (600% against KS-DFT).
Collapse
|
40
|
|
41
|
Resonance shifts and spill-out effects in self-consistent hydrodynamic nanoplasmonics. Nat Commun 2015; 6:7132. [PMID: 26013263 DOI: 10.1038/ncomms8132] [Citation(s) in RCA: 216] [Impact Index Per Article: 24.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2014] [Accepted: 04/09/2015] [Indexed: 12/23/2022] Open
Abstract
The standard hydrodynamic Drude model with hard-wall boundary conditions can give accurate quantitative predictions for the optical response of noble-metal nanoparticles. However, it is less accurate for other metallic nanosystems, where surface effects due to electron density spill-out in free space cannot be neglected. Here we address the fundamental question whether the description of surface effects in plasmonics necessarily requires a fully quantum-mechanical ab initio approach. We present a self-consistent hydrodynamic model (SC-HDM), where both the ground state and the excited state properties of an inhomogeneous electron gas can be determined. With this method we are able to explain the size-dependent surface resonance shifts of Na and Ag nanowires and nanospheres. The results we obtain are in good agreement with experiments and more advanced quantum methods. The SC-HDM gives accurate results with modest computational effort, and can be applied to arbitrary nanoplasmonic systems of much larger sizes than accessible with ab initio methods.
Collapse
|
42
|
Wesolowski TA, Shedge S, Zhou X. Frozen-Density Embedding Strategy for Multilevel Simulations of Electronic Structure. Chem Rev 2015; 115:5891-928. [DOI: 10.1021/cr500502v] [Citation(s) in RCA: 219] [Impact Index Per Article: 24.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Tomasz A. Wesolowski
- Department of Physical Chemistry, University of Geneva, CH-1211 Geneva 4, Switzerland
| | - Sapana Shedge
- Department of Physical Chemistry, University of Geneva, CH-1211 Geneva 4, Switzerland
| | - Xiuwen Zhou
- Department of Physical Chemistry, University of Geneva, CH-1211 Geneva 4, Switzerland
| |
Collapse
|
43
|
Śmiga S, Fabiano E, Laricchia S, Constantin LA, Della Sala F. Subsystem density functional theory with meta-generalized gradient approximation exchange-correlation functionals. J Chem Phys 2015; 142:154121. [DOI: 10.1063/1.4917257] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Szymon Śmiga
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Grudziadzka 5, 87-100 Torun, Poland
| | - Eduardo Fabiano
- Istituto Nanoscienze-CNR, Euromediterranean Center for Nanomaterial Modelling and Technology (ECMT), Via per Arnesano 16, I-73100 Lecce, Italy
- Center for Biomolecular Nanotechnologies@UNILE, Istituto Italiano di Tecnologia (IIT), Via Barsanti, 73010 Arnesano, LE, Italy
| | - Savio Laricchia
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Lucian A. Constantin
- Center for Biomolecular Nanotechnologies@UNILE, Istituto Italiano di Tecnologia (IIT), Via Barsanti, 73010 Arnesano, LE, Italy
| | - Fabio Della Sala
- Istituto Nanoscienze-CNR, Euromediterranean Center for Nanomaterial Modelling and Technology (ECMT), Via per Arnesano 16, I-73100 Lecce, Italy
- Center for Biomolecular Nanotechnologies@UNILE, Istituto Italiano di Tecnologia (IIT), Via Barsanti, 73010 Arnesano, LE, Italy
| |
Collapse
|
44
|
Borgoo A, Green JA, Tozer DJ. Molecular Binding in Post-Kohn–Sham Orbital-Free DFT. J Chem Theory Comput 2014; 10:5338-45. [DOI: 10.1021/ct500670h] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Affiliation(s)
- Alex Borgoo
- Department
of Chemistry, Centre for Theoretical and Computational Chemistry, University of Oslo, P.O.
Box 1033, Blindern, Oslo N-0315, Norway
| | - James A. Green
- Department
of Chemistry, Durham University, South Road, Durham DH1 3LE, U.K
| | - David J. Tozer
- Department
of Chemistry, Durham University, South Road, Durham DH1 3LE, U.K
| |
Collapse
|
45
|
Fabiano E, Laricchia S, Della Sala F. Frozen density embedding with non-integer subsystems' particle numbers. J Chem Phys 2014; 140:114101. [PMID: 24655166 DOI: 10.1063/1.4868033] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We extend the frozen density embedding theory to non-integer subsystems' particles numbers. Different features of this formulation are discussed, with special concern for approximate embedding calculations. In particular, we highlight the relation between the non-integer particle-number partition scheme and the resulting embedding errors. Finally, we provide a discussion of the implications of the present theory for the derivative discontinuity issue and the calculation of chemical reactivity descriptors.
Collapse
Affiliation(s)
- Eduardo Fabiano
- National Nanotechnology Laboratory (NNL), Istituto Nanoscienze-CNR, Via per Arnesano 16, I-73100 Lecce, Italy
| | - Savio Laricchia
- National Nanotechnology Laboratory (NNL), Istituto Nanoscienze-CNR, Via per Arnesano 16, I-73100 Lecce, Italy
| | - Fabio Della Sala
- National Nanotechnology Laboratory (NNL), Istituto Nanoscienze-CNR, Via per Arnesano 16, I-73100 Lecce, Italy
| |
Collapse
|
46
|
Jacob CR, Neugebauer J. Subsystem density-functional theory. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2014. [DOI: 10.1002/wcms.1175] [Citation(s) in RCA: 229] [Impact Index Per Article: 22.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Christoph R. Jacob
- Center for Functional Nanostructures and Institute of Physical Chemistry; Karlsruhe Institute of Technology (KIT); Karlsruhe Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut; Westfälische Wilhelms-Universität Münster; Münster Germany
| |
Collapse
|
47
|
Fabiano E, Trevisanutto PE, Terentjevs A, Constantin LA. Generalized Gradient Approximation Correlation Energy Functionals Based on the Uniform Electron Gas with Gap Model. J Chem Theory Comput 2014; 10:2016-26. [DOI: 10.1021/ct500073b] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Eduardo Fabiano
- National
Nanotechnology Laboratory (NNL), Istituto Nanoscienze-CNR, Via
per Arnesano 16, I-73100 Lecce, Italy
| | - Paolo E. Trevisanutto
- Center
for Biomolecular Nanotechnologies @UNILE, Istituto Italiano di Tecnologia, Via Barsanti, I-73010 Arnesano, Italy
| | - Aleksandrs Terentjevs
- National
Nanotechnology Laboratory (NNL), Istituto Nanoscienze-CNR, Via
per Arnesano 16, I-73100 Lecce, Italy
| | - Lucian A. Constantin
- Center
for Biomolecular Nanotechnologies @UNILE, Istituto Italiano di Tecnologia, Via Barsanti, I-73010 Arnesano, Italy
| |
Collapse
|
48
|
Laricchia S, Constantin LA, Fabiano E, Della Sala F. Laplacian-Level Kinetic Energy Approximations Based on the Fourth-Order Gradient Expansion: Global Assessment and Application to the Subsystem Formulation of Density Functional Theory. J Chem Theory Comput 2013; 10:164-79. [DOI: 10.1021/ct400836s] [Citation(s) in RCA: 58] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Savio Laricchia
- Center
for Biomolecular Nanotechnologies @UNILE, Istituto Italiano di Tecnologia (IIT), Via Barsanti, 73010 Arnesano, Lecce, Italy
| | - Lucian A. Constantin
- Center
for Biomolecular Nanotechnologies @UNILE, Istituto Italiano di Tecnologia (IIT), Via Barsanti, 73010 Arnesano, Lecce, Italy
| | - Eduardo Fabiano
- National
Nanotechnology Laboratory (NNL), Istituto Nanoscienze-CNR, Via
per Arnesano 16, 73100 Lecce, Italy
| | - Fabio Della Sala
- Center
for Biomolecular Nanotechnologies @UNILE, Istituto Italiano di Tecnologia (IIT), Via Barsanti, 73010 Arnesano, Lecce, Italy
- National
Nanotechnology Laboratory (NNL), Istituto Nanoscienze-CNR, Via
per Arnesano 16, 73100 Lecce, Italy
| |
Collapse
|
49
|
Bodrenko IV, Della Sala F. A periodic charge-dipole electrostatic model. II. A kinetic-exchange-correlation correction. J Chem Phys 2013; 139:144109. [PMID: 24116605 DOI: 10.1063/1.4824189] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We extend the periodic charge-dipole electrostatic model, see I. V. Bodrenko, M. Sierka, E. Fabiano, and F. Della Sala, J. Chem. Phys. 137, 134702 (2012), to include a kinetic-exchange-correlation (KXC) correction. The KXC correction is approximated by means of an extended-Hückel-type formula, it is exact in the infinite jellium model and it is also computationally efficient as it requires only the computation of overlap integrals. Tests on the linear response of silver slabs to an external electrostatic perturbation show that the KXC correction yields a very accurate description of induced dipole and of the whole induced charge density profile. We also show that the KXC parameters are quite transferable and related to the atomic polarizability.
Collapse
Affiliation(s)
- I V Bodrenko
- National Nanotechnology Laboratory (NNL), Istituto Nanoscienze-CNR, Via per Arnesano 16, 73100 Lecce, Italy
| | | |
Collapse
|
50
|
Höfener S, Gomes ASP, Visscher L. Solvatochromic shifts from coupled-cluster theory embedded in density functional theory. J Chem Phys 2013; 139:104106. [DOI: 10.1063/1.4820488] [Citation(s) in RCA: 42] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
|