1
|
Kussmann J, Lemke Y, Weinbrenner A, Ochsenfeld C. A Constraint-Based Orbital-Optimized Excited State Method (COOX). J Chem Theory Comput 2024; 20:8461-8473. [PMID: 39345090 DOI: 10.1021/acs.jctc.4c00467] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/01/2024]
Abstract
In this work, we present a novel method to directly calculate targeted electronic excited states within a self-consistent field calculation based on constrained density functional theory (cDFT). The constraint is constructed from the static occupied-occupied and virtual-virtual parts of the excited state difference density from (simplified) linear-response time-dependent density functional theory calculations (LR-TDDFT). Our new method shows a stable convergence behavior, provides an accurate excited state density adhering to the Aufbau principle, and can be solved within a restricted SCF for singlet excitations to avoid spin contamination. This also allows the straightforward application of post-SCF electron-correlation methods like MP2 or direct RPA methods. We present the details of our constraint-based orbital-optimized excited state method (COOX) and compare it to similar schemes. The accuracy of excitation energies will be analyzed for a benchmark of systems, while the quality of the resulting excited state densities is investigated by evaluating excited state nuclear forces and excited state structure optimizations. We also investigate the performance of the proposed COOX method for long-range charge transfer excitations and conical intersections with the ground-state.
Collapse
Affiliation(s)
- Jörg Kussmann
- Chair of Theoretical Chemistry, Department of Chemistry, Ludwig-Maximilians-Universität in Munich (LMU), München D-81377, Germany
| | - Yannick Lemke
- Chair of Theoretical Chemistry, Department of Chemistry, Ludwig-Maximilians-Universität in Munich (LMU), München D-81377, Germany
| | - Anthea Weinbrenner
- Chair of Theoretical Chemistry, Department of Chemistry, Ludwig-Maximilians-Universität in Munich (LMU), München D-81377, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, Ludwig-Maximilians-Universität in Munich (LMU), München D-81377, Germany
- Max-Planck-Institute for Solid State Research, Stuttgart D-70659, Germany
| |
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
|
Hashemi A, Peljo P, Laasonen K. Understanding Electron Transfer Reactions Using Constrained Density Functional Theory: Complications Due to Surface Interactions. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2023; 127:3398-3407. [PMID: 36865990 PMCID: PMC9969872 DOI: 10.1021/acs.jpcc.2c06537] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/13/2022] [Revised: 01/25/2023] [Indexed: 06/18/2023]
Abstract
The kinetic rates of electrochemical reactions depend on electrodes and molecules in question. In a flow battery, where the electrolyte molecules are charged and discharged on the electrodes, the efficiency of the electron transfer is of crucial importance for the performance of the device. The purpose of this work is to present a systematic atomic-level computational protocol for studying electron transfer between electrolyte and electrode. The computations are done by using constrained density functional theory (CDFT) to ensure that the electron is either on the electrode or in the electrolyte. The ab initio molecular dynamics (AIMD) is used to simulate the movement of the atoms. We use the Marcus theory to predict electron transfer rates and the combined CDFT-AIMD approach to compute the parameters for the Marcus theory where it is needed. We model the electrode with a single layer of graphene and methylviologen, 4,4'-dimethyldiquat, desalted basic red 5, 2-hydroxy-1,4-naphthaquinone, and 1,1-di(2-ethanol)-4,4-bipyridinium were selected for the electrolyte molecules. All of these molecules undergo consecutive electrochemical reactions with one electron being transferred at each stage. Because of significant electrode-molecule interactions, it is not possible to evaluate outer-sphere ET. This theoretical study contributes toward the development of a realistic-level prediction of electron transfer kinetics suitable for energy storage applications.
Collapse
Affiliation(s)
- Arsalan Hashemi
- Research
Group of Computational Chemistry, Department of Chemistry and Materials
Science, Aalto University, FI-00076 Aalto, Finland
| | - Pekka Peljo
- Research
Group of Battery Materials and Technologies, Department of Mechanical
and Materials Engineering, Faculty of Technology, University of Turku, 20014 Turun Yliopisto, Finland
| | - Kari Laasonen
- Research
Group of Computational Chemistry, Department of Chemistry and Materials
Science, Aalto University, FI-00076 Aalto, Finland
| |
Collapse
|
4
|
Ahart CS, Rosso KM, Blumberger J. Implementation and Validation of Constrained Density Functional Theory Forces in the CP2K Package. J Chem Theory Comput 2022; 18:4438-4446. [PMID: 35700315 PMCID: PMC9281399 DOI: 10.1021/acs.jctc.2c00284] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
Constrained density
functional theory (CDFT) is a powerful tool
for the prediction of electron transfer parameters in condensed phase
simulations at a reasonable computational cost. In this work we present
an extension to CDFT in the popular mixed Gaussian/plane wave electronic
structure package CP2K, implementing the additional force terms arising
from a constraint based on Hirshfeld charge partitioning. This improves
upon the existing Becke partitioning scheme, which is prone to give
unphysical atomic charges. We verify this implementation for a variety
of systems: electron transfer in (H2O)2+ in a vacuum, electron tunnelling
between oxygen vacancy centers in solid MgO, and electron self-exchange
in aqueous Ru2+–Ru3+. We find good agreement
with previous plane-wave CDFT results for the same systems, but at
a significantly lower computational cost, and we discuss the general
reliability of condensed phase CDFT calculations.
Collapse
Affiliation(s)
- Christian S Ahart
- Department of Physics and Astronomy and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
| | - Kevin M Rosso
- Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Jochen Blumberger
- Department of Physics and Astronomy and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
| |
Collapse
|
5
|
Stella M, Thapa K, Genovese L, Ratcliff LE. Transition-Based Constrained DFT for the Robust and Reliable Treatment of Excitations in Supramolecular Systems. J Chem Theory Comput 2022; 18:3027-3038. [PMID: 35471972 PMCID: PMC9097287 DOI: 10.1021/acs.jctc.1c00548] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2021] [Indexed: 11/28/2022]
Abstract
Despite the variety of available computational approaches, state-of-the-art methods for calculating excitation energies, such as time-dependent density functional theory (TDDFT), are computationally demanding and thus limited to moderate system sizes. Here, we introduce a new variation of constrained DFT (CDFT), wherein the constraint corresponds to a particular transition (T), or a combination of transitions, between occupied and virtual orbitals, rather than a region of the simulation space as in traditional CDFT. We compare T-CDFT with TDDFT and ΔSCF results for the low-lying excited states (S1 and T1) of a set of gas-phase acene molecules and OLED emitters and with reference results from the literature. At the PBE level of theory, T-CDFT outperforms ΔSCF for both classes of molecules, while also proving to be more robust. For the local excitations seen in the acenes, T-CDFT and TDDFT perform equally well. For the charge transfer (CT)-like excitations seen in the OLED molecules, T-CDFT also performs well, in contrast to the severe energy underestimation seen with TDDFT. In other words, T-CDFT is equally applicable to both local excitations and CT states, providing more reliable excitation energies at a much lower computational cost than TDDFT cost. T-CDFT is designed for large systems and has been implemented in the linear-scaling BigDFT code. It is therefore ideally suited for exploring the effects of explicit environments on excitation energies, paving the way for future simulations of excited states in complex realistic morphologies, such as those which occur in OLED materials.
Collapse
Affiliation(s)
- Martina Stella
- Department
of Materials, Imperial College London, London SW7 2AZ, U.K.
- The
Abdus Salam International Centre for Theoretical Physics, Condensed Matter and Statistical Physics, Trieste 34151, Italy
| | - Kritam Thapa
- Department
of Materials, Imperial College London, London SW7 2AZ, U.K.
| | - Luigi Genovese
- Université
Grenoble Alpes, CEA, IRIG-MEM-L_Sim, Grenoble 38000, France
| | - Laura E. Ratcliff
- Department
of Materials, Imperial College London, London SW7 2AZ, U.K.
- Centre
for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K.
| |
Collapse
|
6
|
Dawson W, Degomme A, Stella M, Nakajima T, Ratcliff LE, Genovese L. Density functional theory calculations of large systems: Interplay between fragments, observables, and computational complexity. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1574] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
| | | | - Martina Stella
- Department of Materials Imperial College London London UK
| | | | | | - Luigi Genovese
- Université Grenoble Alpes, INAC‐MEM, L_Sim Grenoble France
| |
Collapse
|
7
|
Zotti LA, Dednam W, Lombardi EB, Palacios JJ. Constrained DFT for Molecular Junctions. NANOMATERIALS 2022; 12:nano12071234. [PMID: 35407352 PMCID: PMC9002544 DOI: 10.3390/nano12071234] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/18/2022] [Revised: 03/31/2022] [Accepted: 04/04/2022] [Indexed: 02/01/2023]
Abstract
We have explored the use of constrained density functional theory (cDFT) for molecular junctions based on benzenediamine. By elongating the junction, we observe that the energy gap between the ionization potential and the electronic affinity increases with the stretching distance. This is consistent with the trend expected from the electrostatic screening. A more detailed analysis shows how this influences the charge distribution of both the individual metal layers and the molecular atoms. Overall, our work shows that constrained DFT is a powerful tool for studying screening effects in molecular junctions.
Collapse
Affiliation(s)
- Linda Angela Zotti
- Departamento de Física Teórica de la Materia Condensada, Universidad Autónoma de Madrid, E-28049 Madrid, Spain
- Condensed Matter Physics Center (IFIMAC), Universidad Autónoma de Madrid, E-28049 Madrid, Spain;
- Correspondence:
| | - Wynand Dednam
- Department of Physics, Science Campus, University of South Africa, Private Bag X6, Florida Park 1710, South Africa; (W.D.); (E.B.L.)
| | - Enrico B. Lombardi
- Department of Physics, Science Campus, University of South Africa, Private Bag X6, Florida Park 1710, South Africa; (W.D.); (E.B.L.)
| | - Juan Jose Palacios
- Condensed Matter Physics Center (IFIMAC), Universidad Autónoma de Madrid, E-28049 Madrid, Spain;
- Departamento de Física de la Materia Condensada and Instituto Nicolás Cabrera (INC), Universidad Autónoma de Madrid, E-28049 Madrid, Spain
| |
Collapse
|
8
|
Shao X, Mi W, Pavanello M. Efficient DFT Solver for Nanoscale Simulations and Beyond. J Phys Chem Lett 2021; 12:4134-4139. [PMID: 33887132 DOI: 10.1021/acs.jpclett.1c00716] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
We present the one-orbital ensemble self-consistent field (OE-SCF), an alternative orbital-free DFT solver that extends the applicability of DFT to beyond nanoscale system sizes, retaining the accuracy required to be predictive. OE-SCF treats the Pauli potential as an external potential updating it iteratively, dramatically outperforming current solvers because only few iterations are needed to reach convergence. OE-SCF enabled us to carry out the largest ab initio simulations for silicon-based materials to date by employing only 1 CPU. We computed the energy of bulk-cut Si nanoparticles as a function of their diameter up to 16 nm, and the polarization and interface charge transfer when a Si slab is sandwiched between two metal slabs where lattice matching mandated a large contact area. Additionally, OE-SCF opens the door to adopting even more accurate functionals in orbital-free DFT simulations while still tackling large system sizes.
Collapse
|
9
|
Romero-Muñiz C, Nakata A, Pou P, Bowler DR, Miyazaki T, Pérez R. High-accuracy large-scale DFT calculations using localized orbitals in complex electronic systems: the case of graphene-metal interfaces. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2018; 30:505901. [PMID: 30468156 DOI: 10.1088/1361-648x/aaec4c] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Over many years, computational simulations based on density functional theory (DFT) have been used extensively to study many different materials at the atomic scale. However, its application is restricted by system size, leaving a number of interesting systems without a high-accuracy quantum description. In this work, we calculate the electronic and structural properties of a graphene-metal system significantly larger than in previous plane-wave calculations with the same accuracy. For this task we use a localised basis set with the Conquest code, both in their primitive, pseudo-atomic orbital form, and using a recent multi-site approach. This multi-site scheme allows us to maintain accuracy while saving computational time and memory requirements, even in our exemplar complex system of graphene grown on Rh(1 1 1) with and without intercalated atomic oxygen. This system offers a rich scenario that will serve as a benchmark, demonstrating that highly accurate simulations in cells with over 3000 atoms are feasible with modest computational resources.
Collapse
Affiliation(s)
- Carlos Romero-Muñiz
- Departamento de Física Teórica de la Materia Condensada, Universidad Autónoma de Madrid, E-28049 Madrid, Spain. First-Principles Simulation Group, Nano-Theory Field, International Center for Materials Nanoarchitectonics (WPI-MANA), National Institute for Materials Science (NIMS), 1-1 Namiki, Tsukuba, Ibaraki 305-0044, Japan
| | | | | | | | | | | |
Collapse
|
10
|
Nakata A, Futamura Y, Sakurai T, Bowler DR, Miyazaki T. Efficient Calculation of Electronic Structure Using O(N) Density Functional Theory. J Chem Theory Comput 2017; 13:4146-4153. [DOI: 10.1021/acs.jctc.7b00385] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ayako Nakata
- First-Principles
Simulation Group, Nano-Theory Field, International Center for Materials
Nanoarchitectonics (WPI-MANA), National Institute for Materials Science
(NIMS), 1-1 Namiki, Tsukuba, Ibaraki 305-0044, Japan
| | - Yasunori Futamura
- Department
of Computer Science, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-8573, Japan
| | - Tetsuya Sakurai
- Department
of Computer Science, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-8573, Japan
- CREST, Japan Science
and Technology Agency, 4-1-8 Hon-cho, Kawaguchi, Saitama 332-0012, Japan
| | - David R Bowler
- Department of Physics & Astronomy, University College London, Gower Street, London WC1E 6BT, U.K
- WPI-MANA, National
Institute for Materials Science (NIMS), 1-1 Namiki, Tsukuba, Ibaraki 305-0044, Japan
- London Centre for Nanotechnology, University College London, 17-19 Gordon Street, London WC1H 0AH, U.K
| | - Tsuyoshi Miyazaki
- First-Principles
Simulation Group, Nano-Theory Field, International Center for Materials
Nanoarchitectonics (WPI-MANA), National Institute for Materials Science
(NIMS), 1-1 Namiki, Tsukuba, Ibaraki 305-0044, Japan
- Department of Physics & Astronomy, University College London, Gower Street, London WC1E 6BT, U.K
| |
Collapse
|
11
|
Goldey MB, Brawand NP, Vörös M, Galli G. Charge Transport in Nanostructured Materials: Implementation and Verification of Constrained Density Functional Theory. J Chem Theory Comput 2017; 13:2581-2590. [DOI: 10.1021/acs.jctc.7b00088] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Matthew B. Goldey
- Institute
for Molecule Engineering, University of Chicago, Chicago, Illinois 60637, United States
| | - Nicholas P. Brawand
- Institute
for Molecule Engineering, University of Chicago, Chicago, Illinois 60637, United States
| | - Márton Vörös
- Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Giulia Galli
- Institute
for Molecule Engineering, University of Chicago, Chicago, Illinois 60637, United States
- Argonne National Laboratory, Lemont, Illinois 60439, United States
| |
Collapse
|
12
|
Pelzer KM, Vázquez-Mayagoitia Á, Ratcliff LE, Tretiak S, Bair RA, Gray SK, Van Voorhis T, Larsen RE, Darling SB. Molecular dynamics and charge transport in organic semiconductors: a classical approach to modeling electron transfer. Chem Sci 2017; 8:2597-2609. [PMID: 28553494 PMCID: PMC5431633 DOI: 10.1039/c6sc04547b] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2016] [Accepted: 01/03/2017] [Indexed: 11/21/2022] Open
Abstract
Organic photovoltaics (OPVs) are a promising carbon-neutral energy conversion technology, with recent improvements pushing power conversion efficiencies over 10%. A major factor limiting OPV performance is inefficiency of charge transport in organic semiconducting materials (OSCs). Due to strong coupling with lattice degrees of freedom, the charges form polarons, localized quasi-particles comprised of charges dressed with phonons. These polarons can be conceptualized as pseudo-atoms with a greater effective mass than a bare charge. We propose that due to this increased mass, polarons can be modeled with Langevin molecular dynamics (LMD), a classical approach with a computational cost much lower than most quantum mechanical methods. Here we present LMD simulations of charge transfer between a pair of fullerene molecules, which commonly serve as electron acceptors in OSCs. We find transfer rates consistent with experimental measurements of charge mobility, suggesting that this method may provide quantitative predictions of efficiency when used to simulate materials on the device scale. Our approach also offers information that is not captured in the overall transfer rate or mobility: in the simulation data, we observe exactly when and why intermolecular transfer events occur. In addition, we demonstrate that these simulations can shed light on the properties of polarons in OSCs. Much remains to be learned about these quasi-particles, and there are no widely accepted methods for calculating properties such as effective mass and friction. Our model offers a promising approach to exploring mass and friction as well as providing insight into the details of polaron transport in OSCs.
Collapse
Affiliation(s)
- Kenley M Pelzer
- Center for Nanoscale Materials , Argonne National Laboratory , 9700 Cass Ave., Lemont , IL 60439 , USA . ; Tel: +1-630-252-7020
- Materials Science Division , Argonne National Laboratory , 9700 Cass Ave, Lemont , IL 60439 , USA
| | - Álvaro Vázquez-Mayagoitia
- Argonne Leadership Computing Facility , Argonne National Laboratory , 9700 Cass Ave. , Lemont , IL 60439 , USA
| | - Laura E Ratcliff
- Argonne Leadership Computing Facility , Argonne National Laboratory , 9700 Cass Ave. , Lemont , IL 60439 , USA
| | - Sergei Tretiak
- Theoretical Division , Center for Nonlinear Studies , Center for Integrated Nanotechnologies , Los Alamos National Laboratory , Los Alamos , NM 87545 , USA
| | - Raymond A Bair
- Mathematics and Computer Science Division , Argonne National Laboratory , 9700 Cass Ave. , Argonne , IL 60439 , USA
- Computation Institute , University of Chicago , 5735 S. Ellis Ave. , Chicago , IL 60637 , USA
- Computer, Environment, and Life Sciences , Argonne National Laboratory , 9700 Cass Ave. , Lemont , IL 60439 , USA
| | - Stephen K Gray
- Center for Nanoscale Materials , Argonne National Laboratory , 9700 Cass Ave., Lemont , IL 60439 , USA . ; Tel: +1-630-252-7020
- Computation Institute , University of Chicago , 5735 S. Ellis Ave. , Chicago , IL 60637 , USA
| | - Troy Van Voorhis
- Department of Chemistry , Massachusetts Institute of Technology , 77 Massachusetts Ave , Cambridge , MA 02139 , USA
| | - Ross E Larsen
- Computational Science Center , National Renewable Energy Laboratory , 15301 Denver W. Parkway, Golden , CO 80401 , USA
| | - Seth B Darling
- Center for Nanoscale Materials , Argonne National Laboratory , 9700 Cass Ave., Lemont , IL 60439 , USA . ; Tel: +1-630-252-7020
- Institute for Molecular Engineering , University of Chicago , 5747 S. Ellis Ave. , Chicago , IL 60637 , USA
| |
Collapse
|
13
|
Holmberg N, Laasonen K. Efficient Constrained Density Functional Theory Implementation for Simulation of Condensed Phase Electron Transfer Reactions. J Chem Theory Comput 2017; 13:587-601. [PMID: 28009515 DOI: 10.1021/acs.jctc.6b01085] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
Constrained density functional theory (CDFT) is a versatile tool for probing the kinetics of electron transfer (ET) reactions. In this work, we present a well-scaling parallel CDFT implementation relying on a mixed basis set of Gaussian functions and plane waves, which has been specifically tailored to investigate condensed phase ET reactions using an explicit, quantum chemical representation of the solvent. The accuracy of our implementation is validated against previous theoretical results for predicting electronic couplings and charge transfer energies. Subsequently, we demonstrate the efficiency of our method by studying the intramolecular ET reaction of an organic mixed-valence compound in water using a CDFT based molecular dynamics simulation.
Collapse
Affiliation(s)
- Nico Holmberg
- COMP Centre of Excellence in Computational Nanoscience, Department of Chemistry, Aalto University , P.O. Box 16100, FI-00076 Aalto, Finland
| | - Kari Laasonen
- COMP Centre of Excellence in Computational Nanoscience, Department of Chemistry, Aalto University , P.O. Box 16100, FI-00076 Aalto, Finland
| |
Collapse
|
14
|
Akimov AV. Nonadiabatic Molecular Dynamics with Tight-Binding Fragment Molecular Orbitals. J Chem Theory Comput 2016; 12:5719-5736. [DOI: 10.1021/acs.jctc.6b00955] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Alexey V. Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260-3000, United States
| |
Collapse
|
15
|
Ratcliff LE, Mohr S, Huhs G, Deutsch T, Masella M, Genovese L. Challenges in large scale quantum mechanical calculations. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2016. [DOI: 10.1002/wcms.1290] [Citation(s) in RCA: 70] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Laura E. Ratcliff
- Argonne Leadership Computing Facility Argonne National Laboratory Lemon IL USA
| | - Stephan Mohr
- Department of Computer Applications in Science and Engineering Barcelona Supercomputing Center (BSC‐CNS) Barcelona Spain
| | - Georg Huhs
- Department of Computer Applications in Science and Engineering Barcelona Supercomputing Center (BSC‐CNS) Barcelona Spain
| | - Thierry Deutsch
- University Grenoble Alpes INAC‐MEM Grenoble France
- CEA, INAC‐MEM Grenoble France
| | - Michel Masella
- Laboratoire de Biologie Structurale et Radiologie, Service de Bioénergétique, Biologie Structurale et Mécanisme Institut de Biologie et de Technologie de Saclay, CEA Saclay Gif‐sur‐Yvette Cedex France
| | - Luigi Genovese
- University Grenoble Alpes INAC‐MEM Grenoble France
- CEA, INAC‐MEM Grenoble France
| |
Collapse
|
16
|
Melander M, Jónsson EÖ, Mortensen JJ, Vegge T, García Lastra JM. Implementation of Constrained DFT for Computing Charge Transfer Rates within the Projector Augmented Wave Method. J Chem Theory Comput 2016; 12:5367-5378. [DOI: 10.1021/acs.jctc.6b00815] [Citation(s) in RCA: 40] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Marko Melander
- Department
of Energy Conversion and Storage, Technical University of Denmark, DK-4000 Roskilde, Denmark
| | - Elvar Ö. Jónsson
- COMP,
Applied Physics Department, Aalto University FI-00076 Aalto, Espoo, Finland
| | - Jens J. Mortensen
- Department
of Physics, Center for Atomic-scale Materials Design, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark
| | - Tejs Vegge
- Department
of Energy Conversion and Storage, Technical University of Denmark, DK-4000 Roskilde, Denmark
| | - Juan Maria García Lastra
- Department
of Energy Conversion and Storage, Technical University of Denmark, DK-4000 Roskilde, Denmark
| |
Collapse
|
17
|
Blumberger J. Recent Advances in the Theory and Molecular Simulation of Biological Electron Transfer Reactions. Chem Rev 2015; 115:11191-238. [DOI: 10.1021/acs.chemrev.5b00298] [Citation(s) in RCA: 238] [Impact Index Per Article: 26.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Jochen Blumberger
- Department of Physics and
Astronomy, University College London, Gower Street, London WC1E 6BT, U.K
| |
Collapse
|
18
|
Ratcliff LE, Genovese L, Mohr S, Deutsch T. Fragment approach to constrained density functional theory calculations using Daubechies wavelets. J Chem Phys 2015; 142:234105. [DOI: 10.1063/1.4922378] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Laura E. Ratcliff
- Argonne Leadership Computing Facility, Argonne National Laboratory, Lemont, Illinois 60439, USA
- Université de Grenoble Alpes, CEA, INAC-SP2M, L_Sim, F-38000 Grenoble, France
| | - Luigi Genovese
- Université de Grenoble Alpes, CEA, INAC-SP2M, L_Sim, F-38000 Grenoble, France
| | - Stephan Mohr
- Université de Grenoble Alpes, CEA, INAC-SP2M, L_Sim, F-38000 Grenoble, France
| | - Thierry Deutsch
- Université de Grenoble Alpes, CEA, INAC-SP2M, L_Sim, F-38000 Grenoble, France
| |
Collapse
|
19
|
Ratcliff LE, Grisanti L, Genovese L, Deutsch T, Neumann T, Danilov D, Wenzel W, Beljonne D, Cornil J. Toward Fast and Accurate Evaluation of Charge On-Site Energies and Transfer Integrals in Supramolecular Architectures Using Linear Constrained Density Functional Theory (CDFT)-Based Methods. J Chem Theory Comput 2015; 11:2077-86. [DOI: 10.1021/acs.jctc.5b00057] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Laura E. Ratcliff
- Université
Grenoble Alpes, INAC-SP2M, L_Sim, F-38000 Grenoble, France
- CEA, INAC-SP2M,
L_Sim, F-38000 Grenoble, France
| | - Luca Grisanti
- Laboratory
for Chemistry of Novel Materials, University of Mons, Place du Parc
20, B-7000 Mons, Belgium
- The Abdus Salam
International Centre for Theoretical Physics (ICTP), Strada Costiera 11, I-34151 Trieste, Italy
| | - Luigi Genovese
- Université
Grenoble Alpes, INAC-SP2M, L_Sim, F-38000 Grenoble, France
- CEA, INAC-SP2M,
L_Sim, F-38000 Grenoble, France
| | - Thierry Deutsch
- Université
Grenoble Alpes, INAC-SP2M, L_Sim, F-38000 Grenoble, France
- CEA, INAC-SP2M,
L_Sim, F-38000 Grenoble, France
| | - Tobias Neumann
- Institute
of Nanotechnology, Karlsruhe Institute of Technology (KIT), 76344 Eggenstein-Leopoldshafen, Germany
| | - Denis Danilov
- Institute
of Nanotechnology, Karlsruhe Institute of Technology (KIT), 76344 Eggenstein-Leopoldshafen, Germany
| | - Wolfgang Wenzel
- Institute
of Nanotechnology, Karlsruhe Institute of Technology (KIT), 76344 Eggenstein-Leopoldshafen, Germany
| | - David Beljonne
- Laboratory
for Chemistry of Novel Materials, University of Mons, Place du Parc
20, B-7000 Mons, Belgium
| | - Jérôme Cornil
- Laboratory
for Chemistry of Novel Materials, University of Mons, Place du Parc
20, B-7000 Mons, Belgium
| |
Collapse
|
20
|
Akimov AV, Prezhdo OV. Large-Scale Computations in Chemistry: A Bird’s Eye View of a Vibrant Field. Chem Rev 2015; 115:5797-890. [DOI: 10.1021/cr500524c] [Citation(s) in RCA: 159] [Impact Index Per Article: 17.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Alexey V. Akimov
- Department
of Chemistry, University of South California, Los Angeles, California 90089, United States
| | - Oleg V. Prezhdo
- Department
of Chemistry, University of South California, Los Angeles, California 90089, United States
| |
Collapse
|
21
|
Nishimoto Y, Fedorov DG, Irle S. Density-Functional Tight-Binding Combined with the Fragment Molecular Orbital Method. J Chem Theory Comput 2014; 10:4801-12. [DOI: 10.1021/ct500489d] [Citation(s) in RCA: 73] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Affiliation(s)
| | - Dmitri G. Fedorov
- Nanosystem
Research Institute (NRI), National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba 305-8565, Japan
| | | |
Collapse
|
22
|
Recent advances in QM/MM free energy calculations using reference potentials. Biochim Biophys Acta Gen Subj 2014; 1850:954-965. [PMID: 25038480 PMCID: PMC4547088 DOI: 10.1016/j.bbagen.2014.07.008] [Citation(s) in RCA: 51] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2014] [Revised: 07/06/2014] [Accepted: 07/07/2014] [Indexed: 01/02/2023]
Abstract
Background Recent years have seen enormous progress in the development of methods for modeling (bio)molecular systems. This has allowed for the simulation of ever larger and more complex systems. However, as such complexity increases, the requirements needed for these models to be accurate and physically meaningful become more and more difficult to fulfill. The use of simplified models to describe complex biological systems has long been shown to be an effective way to overcome some of the limitations associated with this computational cost in a rational way. Scope of review Hybrid QM/MM approaches have rapidly become one of the most popular computational tools for studying chemical reactivity in biomolecular systems. However, the high cost involved in performing high-level QM calculations has limited the applicability of these approaches when calculating free energies of chemical processes. In this review, we present some of the advances in using reference potentials and mean field approximations to accelerate high-level QM/MM calculations. We present illustrative applications of these approaches and discuss challenges and future perspectives for the field. Major conclusions The use of physically-based simplifications has shown to effectively reduce the cost of high-level QM/MM calculations. In particular, lower-level reference potentials enable one to reduce the cost of expensive free energy calculations, thus expanding the scope of problems that can be addressed. General significance As was already demonstrated 40 years ago, the usage of simplified models still allows one to obtain cutting edge results with substantially reduced computational cost. This article is part of a Special Issue entitled Recent developments of molecular dynamics. We present some of the advances to accelerate high-level QM/MM calculations. Quantitative limitations of low-level methods can be overcome by these approaches. Reference potentials make free energy simulations feasible for large systems. Automated fitting reduces the need of expensive sampling of high-level approaches. Application of reference potentials can be extended to a wide range of processes.
Collapse
|
23
|
Giese T, Chen H, Huang M, York DM. Parametrization of an Orbital-Based Linear-Scaling Quantum Force Field for Noncovalent Interactions. J Chem Theory Comput 2014; 10:1086-1098. [PMID: 24803856 PMCID: PMC3985928 DOI: 10.1021/ct401035t] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2013] [Indexed: 01/22/2023]
Abstract
We parametrize a linear-scaling quantum mechanical force field called mDC for the accurate reproduction of nonbonded interactions. We provide a new benchmark database of accurate ab initio interactions between sulfur-containing molecules. A variety of nonbond databases are used to compare the new mDC method with other semiempirical, molecular mechanical, ab initio, and combined semiempirical quantum mechanical/molecular mechanical methods. It is shown that the molecular mechanical force field significantly and consistently reproduces the benchmark results with greater accuracy than the semiempirical models and our mDC model produces errors twice as small as the molecular mechanical force field. The comparisons between the methods are extended to the docking of drug candidates to the Cyclin-Dependent Kinase 2 protein receptor. We correlate the protein-ligand binding energies to their experimental inhibition constants and find that the mDC produces the best correlation. Condensed phase simulation of mDC water is performed and shown to produce O-O radial distribution functions similar to TIP4P-EW.
Collapse
Affiliation(s)
- Timothy
J. Giese
- BioMaPS
Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854-8087, United States
| | - Haoyuan Chen
- BioMaPS
Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854-8087, United States
| | - Ming Huang
- BioMaPS
Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854-8087, United States
- Scientific
Computation, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455−0431, United States
| | - Darrin M. York
- BioMaPS
Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854-8087, United States
| |
Collapse
|
24
|
Terranova U, Bowler DR. Δ Self-Consistent Field Method for Natural Anthocyanidin Dyes. J Chem Theory Comput 2013; 9:3181-8. [PMID: 26583995 DOI: 10.1021/ct400356k] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
We present an application of the Δ self-consistent field (ΔSCF) method, which we have implemented and tested in the DFT code CONQUEST, on the study of excited states of natural anthocyanidin dyes. We show that ΔSCF allows relaxation of the atomic structure for systems in excited states by following gradients on the excited Born-Oppenheimer surface. We compare the vertical excitation energies of some anthocyanidins in gas-phase to results from time-dependent density functional theory (TDDFT) and experiments. To reproduce a typical dye-sensitized solar cell interface, we adsorb cyanidin on TiO2 anatase (101), focusing on the shift of the lowest excitation energy due to the adsorption. We have found that important modifications occur in the excited state geometry of the adsorbed cyanidin.
Collapse
Affiliation(s)
- U Terranova
- Department of Physics and Astronomy, University College London , London WC1E 6BT, U.K.,International Center for Materials Nanoarchitectonics (MANA) , 1-1 Namiki, Tsukuba, Ibaraki 305-0044, Japan
| | - D R Bowler
- Department of Physics and Astronomy, University College London , London WC1E 6BT, U.K.,International Center for Materials Nanoarchitectonics (MANA) , 1-1 Namiki, Tsukuba, Ibaraki 305-0044, Japan
| |
Collapse
|
25
|
Bowler DR, Miyazaki T. O(N) methods in electronic structure calculations. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2012; 75:036503. [PMID: 22790422 DOI: 10.1088/0034-4885/75/3/036503] [Citation(s) in RCA: 210] [Impact Index Per Article: 17.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/21/2023]
Abstract
Linear-scaling methods, or O(N) methods, have computational and memory requirements which scale linearly with the number of atoms in the system, N, in contrast to standard approaches which scale with the cube of the number of atoms. These methods, which rely on the short-ranged nature of electronic structure, will allow accurate, ab initio simulations of systems of unprecedented size. The theory behind the locality of electronic structure is described and related to physical properties of systems to be modelled, along with a survey of recent developments in real-space methods which are important for efficient use of high-performance computers. The linear-scaling methods proposed to date can be divided into seven different areas, and the applicability, efficiency and advantages of the methods proposed in these areas are then discussed. The applications of linear-scaling methods, as well as the implementations available as computer programs, are considered. Finally, the prospects for and the challenges facing linear-scaling methods are discussed.
Collapse
Affiliation(s)
- D R Bowler
- London Centre for Nanotechnology, UCL, 17-19 Gordon St, London WC1H 0AH, UK.
| | | |
Collapse
|
26
|
Řezáč J, Lévy B, Demachy I, de la Lande A. Robust and Efficient Constrained DFT Molecular Dynamics Approach for Biochemical Modeling. J Chem Theory Comput 2012; 8:418-27. [DOI: 10.1021/ct200570u] [Citation(s) in RCA: 42] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- Jan Řezáč
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic and Center for Biomolecules and Complex Molecular Systems, Flemingovo nam. 2, 166 10 Prague 6, Czech Republic
| | - Bernard Lévy
- Laboratoire de Chimie Physique, CNRS UMR 8000, Université Paris-Sud. Bât. 349, Campus d’Orsay, 15 rue Jean Perrin, 91405 Orsay Cedex, France
| | - Isabelle Demachy
- Laboratoire de Chimie Physique, CNRS UMR 8000, Université Paris-Sud. Bât. 349, Campus d’Orsay, 15 rue Jean Perrin, 91405 Orsay Cedex, France
| | - Aurélien de la Lande
- Laboratoire de Chimie Physique, CNRS UMR 8000, Université Paris-Sud. Bât. 349, Campus d’Orsay, 15 rue Jean Perrin, 91405 Orsay Cedex, France
| |
Collapse
|
27
|
Zhang J, Valeev EF. Hybrid one-electron/many-electron methods for ionized states of molecular clusters. Phys Chem Chem Phys 2012; 14:7863-71. [DOI: 10.1039/c2cp40222j] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|