1
|
Majid A, Raza NZ, Haider S, Alam K, Naeem S. Electronic Transport Properties of Molecular Clusters Sb 4O 6, P 4Se 3, and P 4O 6. J Phys Chem A 2024; 128:4814-4822. [PMID: 38857364 DOI: 10.1021/acs.jpca.4c02757] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2024]
Abstract
Inorganic molecular crystal (IMC) is a trending class of materials in which structural units comprise molecular cages or clusters bonded via van der Waal forces. The structure-property relationship in IMCs is less known due to the unusual assembly of molecular clusters in these materials. In this paper, the density functional theory-calculated electronic transport properties of the molecular clusters of antimony oxide (Sb4O6), phosphorus triselenide (P4Se3), and phosphorus trioxide (P4O6) are described in detail. The calculated values of highest occupied molecular orbital-lowest unoccupied molecular orbital gaps appeared as 5.487, 2.296, and 4.425 eV for Sb4O6, P4Se3, and P4O6, respectively. The work was carried out to explore the charge transport mechanism in IMCs in order to disclose their potential in practical applications. The calculations involved charge-transfer integral based on Marcus theory to compute the electronic coupling (V), reorganization energies (λ), and hopping rate (k) in the structures. The hopping rate for Sb4O6, P4Se3, and P4O6 is found as 8.49 × 10-12, 1.28 × 10-14, and 2.51 × 10-20 s-1, respectively. The transport properties of Sb4O6 are found better, which predicts the application of the relevant IMC for device grade applications. The findings of this study are important for future application of the IMCs in electronic and optoelectronic applications.
Collapse
Affiliation(s)
- Abdul Majid
- Department of Physics, University of Gujrat, Gujrat 50700, Pakistan
| | - Nimra Zaib Raza
- Department of Physics, University of Gujrat, Gujrat 50700, Pakistan
| | - Sajjad Haider
- Chemical Engineering Department, College of Engineering, King Saud University, P.O.Box 800, Riyadh 11421, Saudi Arabia
| | - Kamran Alam
- Department of Chemical Engineering Materials Environment Sapienza, University of Rome, Roma RM 00185, Italy
| | - Samia Naeem
- Department of Physics, Government College Women University Sialkot, Sialkot 51310, Pakistan
| |
Collapse
|
2
|
Domingos SR, Tikhonov DS, Steber AL, Eschenbach P, Gruet S, Hrodmarsson HR, Martin K, Garcia GA, Nahon L, Neugebauer J, Avarvari N, Schnell M. Evolution of the ionisation energy with the stepwise growth of chiral clusters of [4]helicene. Nat Commun 2024; 15:4928. [PMID: 38858352 PMCID: PMC11164862 DOI: 10.1038/s41467-024-48778-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2023] [Accepted: 05/13/2024] [Indexed: 06/12/2024] Open
Abstract
Polycyclic aromatic hydrocarbons (PAHs) are widely established as ubiquitous in the interstellar medium (ISM), but considering their prevalence in harsh vacuum environments, the role of ionisation in the formation of PAH clusters is poorly understood, particularly if a chirality-dependent aggregation route is considered. Here we report on photoelectron spectroscopy experiments on [4]helicene clusters performed with a vacuum ultraviolet synchrotron beamline. Aggregates (up to the heptamer) of [4]helicene, the smallest PAH with helical chirality, were produced and investigated with a combined experimental and theoretical approach using several state-of-the-art quantum-chemical methodologies. The ionisation onsets are extracted for each cluster size from the mass-selected photoelectron spectra and compared with calculations of vertical ionisation energies. We explore the complex aggregation topologies emerging from the multitude of isomers formed through clustering of P and M, the two enantiomers of [4]helicene. The very satisfactory benchmarking between experimental ionisation onsets vs. predicted ionisation energies allows the identification of theoretically predicted potential aggregation motifs and corresponding energetic ordering of chiral clusters. Our structural models suggest that a homochiral aggregation route is energetically favoured over heterochiral arrangements with increasing cluster size, hinting at potential symmetry breaking in PAH cluster formation at the scale of small grains.
Collapse
Affiliation(s)
- Sérgio R Domingos
- Deutsches Elektronen-Synchrotron DESY, Notkestr. 85, 22607, Hamburg, Germany.
- CFisUC, Department of Physics, University of Coimbra, 3004-516, Coimbra, Portugal.
| | - Denis S Tikhonov
- Deutsches Elektronen-Synchrotron DESY, Notkestr. 85, 22607, Hamburg, Germany.
| | - Amanda L Steber
- Deutsches Elektronen-Synchrotron DESY, Notkestr. 85, 22607, Hamburg, Germany.
- Department of Physical Chemistry, Faculty of Science, University of Valladolid, 47011, Valladolid, Spain.
| | - Patrick Eschenbach
- Organisch-Chemisches Institut, University of Münster, 48149, Münster, Germany
- Center for Multiscale Theory and Computation (CMTC), University of Münster, 48149, Münster, Germany
| | - Sebastien Gruet
- Deutsches Elektronen-Synchrotron DESY, Notkestr. 85, 22607, Hamburg, Germany
| | - Helgi R Hrodmarsson
- Synchrotron SOLEIL, L'Orme des Merisiers, 91192, Gif sur Yvette, Cedex, France
- LISA UMR 7583 Université Paris-Est Créteil and Université de Paris, Institut Pierre et Simon Laplace, 61 Avenue du Général de Gaulle, 94010, Créteil, France
| | - Kévin Martin
- Univ Angers, CNRS, MOLTECH-Anjou, SFR MATRIX, 49000, Angers, France
| | - Gustavo A Garcia
- Synchrotron SOLEIL, L'Orme des Merisiers, 91192, Gif sur Yvette, Cedex, France
| | - Laurent Nahon
- Synchrotron SOLEIL, L'Orme des Merisiers, 91192, Gif sur Yvette, Cedex, France
| | - Johannes Neugebauer
- Organisch-Chemisches Institut, University of Münster, 48149, Münster, Germany
- Center for Multiscale Theory and Computation (CMTC), University of Münster, 48149, Münster, Germany
| | - Narcis Avarvari
- Univ Angers, CNRS, MOLTECH-Anjou, SFR MATRIX, 49000, Angers, France
| | - Melanie Schnell
- Deutsches Elektronen-Synchrotron DESY, Notkestr. 85, 22607, Hamburg, Germany.
- Institut für Physikalische Chemie, Christian-Albrechts-Universität zu Kiel, Max-Eyth-Str. 1, 24118, Kiel, Germany.
| |
Collapse
|
3
|
Narita M, Kavungathodi MFM, Dheendayal M, Wagner P, Mori S, Mozer AJ. High Electronic Coupling between Cu Complexes and Oxidized Dyes Confirmed by Measurements of Driving Force Dependent Regeneration Kinetics in Minimal Electrolyte System. J Am Chem Soc 2024; 146:12310-12314. [PMID: 38668078 DOI: 10.1021/jacs.4c02237] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/09/2024]
Abstract
We confirm fast regeneration kinetics between copper complexes and oxidized organic dyes and the major contribution of electronic coupling (HDA). The highest efficiency of dye-sensitized TiO2 solar cells has been shown by employing Cu complex redox couples. Various groups have reported a fast regeneration rate of oxidized dyes by Cu complexes giving a low driving force attributed to low reorganization energy (λ), but the effect of HDA has not been evaluated. The values of HDA and λ can be derived from driving force dependent transient absorption (TA) measurements. However, analyzing TA decay using Cu complexes is not trivial because accelerated recombination by the presence of Cu2+ complexes and biphasic TA decay often complicates the analysis. Here we employ 16 Cu1+ and Co2+ complexes and two dyes. To simplify the system, i.e., making a minimal electrolyte system, Cu2+ and Co3+ complexes and a common additive of 4-tert-butylpyridine are not used. From the driving force dependent TA decays of oxidized dyes by both Cu1+ and Co2+ complexes, λ for the combination of the Cu complexes and dyes is found to be about 0.15 eV lower than that of Co complexes. Approximately 3 to 5 times higher HDA values of Cu complexes than those of Co complexes are obtained, which is the dominant factor for faster rates. The values vary with the structure of the molecules, showing the possibility of increasing the HDA values further. The higher HDA values of a Cu complex than that of a Co complex are also reproduced by quantum chemical calculations.
Collapse
Affiliation(s)
- Mitsuru Narita
- Division of Chemistry and Materials, Faculty of Textile Science, Shinshu University, Ueda, Nagano 386-8567, Japan
- Intelligent Polymer Research Institute, Australian Institute for Innovative Materials, University of Wollongong, Wollongong, New South Wales 2522, Australia
| | - Munavvar Fairoos Mele Kavungathodi
- Intelligent Polymer Research Institute, Australian Institute for Innovative Materials, University of Wollongong, Wollongong, New South Wales 2522, Australia
| | - Mantra Dheendayal
- Intelligent Polymer Research Institute, Australian Institute for Innovative Materials, University of Wollongong, Wollongong, New South Wales 2522, Australia
| | - Pawel Wagner
- Intelligent Polymer Research Institute, Australian Institute for Innovative Materials, University of Wollongong, Wollongong, New South Wales 2522, Australia
| | - Shogo Mori
- Division of Chemistry and Materials, Faculty of Textile Science, Shinshu University, Ueda, Nagano 386-8567, Japan
| | - Attila J Mozer
- Intelligent Polymer Research Institute, Australian Institute for Innovative Materials, University of Wollongong, Wollongong, New South Wales 2522, Australia
| |
Collapse
|
4
|
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
|
5
|
Chakraborty R, Boguslawski K, Tecmer P. Static embedding with pair coupled cluster doubles based methods. Phys Chem Chem Phys 2023; 25:25377-25388. [PMID: 37705409 DOI: 10.1039/d3cp02502k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/15/2023]
Abstract
Quantum embedding methods have recently been significantly developed to model large molecular structures. This work proposes a novel wave function theory in a density functional theory (WTF-in-DFT) embedding scheme based on pair-coupled cluster doubles (pCCD)-type methods. While pCCD can reliably describe strongly-correlated systems with mean-field-like computational cost, the large extent of the dynamic correlation can be accounted for by (linearized) coupled-cluster corrections on top of the pCCD wave function. Here we focus on the linearized coupled-cluster singles and doubles (LCCSD) ansatz for electronic ground states and its extension to excited states within the equation of motion (EOM) formalism. We test our EOM-pCCD-LCCSD-in-DFT approach for the vertical excitation energies of the hydrogen-bonded water-ammonia complex, micro-solvated thymine, and uranyl tetrahalides (UO2X42-, X = F, Cl, Br). Furthermore, we assess the quality of the embedding potential using an orbital entanglement and correlation analysis. The approximate embedding models successfully capture changes in the excitation energies going from bare fragments to supramolecular structures and represent a promising computational method for excited states in large molecular systems.
Collapse
Affiliation(s)
- Rahul Chakraborty
- Institute of Physics, Faculty of Physics, Astronomy, and Informatics, Nicolaus Copernicus University in Toruń, Grudziadzka 5, 87-100 Toruń, Poland.
| | - Katharina Boguslawski
- Institute of Physics, Faculty of Physics, Astronomy, and Informatics, Nicolaus Copernicus University in Toruń, Grudziadzka 5, 87-100 Toruń, Poland.
| | - Paweł Tecmer
- Institute of Physics, Faculty of Physics, Astronomy, and Informatics, Nicolaus Copernicus University in Toruń, Grudziadzka 5, 87-100 Toruń, Poland.
| |
Collapse
|
6
|
Sen S, Visscher L. Towards the description of charge transfer states in solubilised LHCII using subsystem DFT. PHOTOSYNTHESIS RESEARCH 2023; 156:39-57. [PMID: 35988131 PMCID: PMC10070235 DOI: 10.1007/s11120-022-00950-7] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/11/2022] [Accepted: 08/04/2022] [Indexed: 05/05/2023]
Abstract
Light harvesting complex II (LHCII) in plants and green algae have been shown to adapt their absorption properties, depending on the concentration of sunlight, switching between a light harvesting and a non-harvesting or quenched state. In a recent work, combining classical molecular dynamics (MD) simulations with quantum chemical calculations (Liguori et al. in Sci Rep 5:15661, 2015) on LHCII, it was shown that the Chl611-Chl612 cluster of the terminal emitter domain can play an important role in modifying the spectral properties of the complex. In that work the importance of charge transfer (CT) effects was highlighted, in re-shaping the absorption intensity of the chlorophyll dimer. Here in this work, we investigate the combined effect of the local excited (LE) and CT states in shaping the energy landscape of the chlorophyll dimer. Using subsystem Density Functional Theory over the classical [Formula: see text]s MD trajectory we look explicitly into the excitation energies of the LE and the CT states of the dimer and their corresponding couplings. Upon doing so, we observe a drop in the excitation energies of the CT states, accompanied by an increase in the couplings between the LE/LE and the LE/CT states facilitated by a shorter interchromophoric distance upon equilibration. Both these changes in conjunction, effectively produces a red-shift of the low-lying mixed exciton/CT states of the supramolecular chromophore pair.
Collapse
Affiliation(s)
- Souloke Sen
- Division of Theoretical Chemistry, Faculty of Sciences, Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | - Lucas Visscher
- Division of Theoretical Chemistry, Faculty of Sciences, Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| |
Collapse
|
7
|
Zhang J, Tang Z, Zhang X, Zhu H, Zhao R, Lu Y, Gao J. Target State Optimized Density Functional Theory for Electronic Excited and Diabatic States. J Chem Theory Comput 2023; 19:1777-1789. [PMID: 36917687 DOI: 10.1021/acs.jctc.2c01317] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/16/2023]
Abstract
A flexible self-consistent field method, called target state optimization (TSO), is presented for exploring electronic excited configurations and localized diabatic states. The key idea is to partition molecular orbitals into different subspaces according to the excitation or localization pattern for a target state. Because of the orbital-subspace constraint, orbitals belonging to different subspaces do not mix. Furthermore, the determinant wave function for such excited or diabatic configurations can be variationally optimized as a ground state procedure, unlike conventional ΔSCF methods, without the possibility of collapsing back to the ground state or other lower-energy configurations. The TSO method can be applied both in Hartree-Fock theory and in Kohn-Sham density functional theory (DFT). The density projection procedure and the working equations for implementing the TSO method are described along with several illustrative applications. For valence excited states of organic compounds, it was found that the computed excitation energies from TSO-DFT and time-dependent density functional theory (TD-DFT) are of similar quality with average errors of 0.5 and 0.4 eV, respectively. For core excitation, doubly excited states and charge-transfer states, the performance of TSO-DFT is clearly superior to that from conventional TD-DFT calculations. It is shown that variationally optimized charge-localized diabatic states can be defined using TSO-DFT in energy decomposition analysis to gain both qualitative and quantitative insights on intermolecular interactions. Alternatively, the variational diabatic states may be used in molecular dynamics simulation of charge transfer processes. The TSO method can also be used to define basis states in multistate density functional theory for excited states through nonorthogonal state interaction calculations. The software implementing TSO-DFT can be accessed from the authors.
Collapse
Affiliation(s)
- Jun Zhang
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China
| | - Zhen Tang
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China
| | - Xiaoyong Zhang
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China
| | - Hong Zhu
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China.,School of Chemical Biology and Biotechnology, Peking University Shenzhen Graduate School, Shenzhen 518055, P. R. China
| | - Ruoqi Zhao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China.,Institute of Theoretical Chemistry, Jilin University, Changchun, 130023 Jilin, P. R. China
| | - Yangyi Lu
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China
| | - Jiali Gao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China.,School of Chemical Biology and Biotechnology, Peking University Shenzhen Graduate School, Shenzhen 518055, P. R. China.,Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
8
|
Eschenbach P, Artiukhin DG, Neugebauer J. Reliable Isotropic Electron-Paramagnetic-Resonance Hyperfine Coupling Constants from the Frozen-Density Embedding Quasi-Diabatization Approach. J Phys Chem A 2022; 126:8358-8368. [DOI: 10.1021/acs.jpca.2c04959] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Patrick Eschenbach
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| | - Denis G. Artiukhin
- Institut für Chemie und Biochemie, Freie Universität Berlin, Arnimallee 22, 14195 Berlin, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| |
Collapse
|
9
|
Eschenbach P, Neugebauer J. Subsystem density-functional theory: A reliable tool for spin-density based properties. J Chem Phys 2022; 157:130902. [PMID: 36209003 DOI: 10.1063/5.0103091] [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
Subsystem density-functional theory compiles a set of features that allow for efficiently calculating properties of very large open-shell radical systems such as organic radical crystals, proteins, or deoxyribonucleic acid stacks. It is computationally less costly than correlated ab initio wave function approaches and can pragmatically avoid the overdelocalization problem of Kohn-Sham density-functional theory without employing hard constraints on the electron-density. Additionally, subsystem density-functional theory calculations commonly start from isolated fragment electron densities, pragmatically preserving a priori specified subsystem spin-patterns throughout the calculation. Methods based on subsystem density-functional theory have seen a rapid development over the past years and have become important tools for describing open-shell properties. In this Perspective, we address open questions and possible developments toward challenging future applications in connection with subsystem density-functional theory for spin-dependent properties.
Collapse
Affiliation(s)
- Patrick Eschenbach
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| |
Collapse
|
10
|
Alam B, Jiang H, Zimmerman PM, Herbert JM. State-specific solvation for restricted active space spin-flip (RAS-SF) wave functions based on the polarizable continuum formalism. J Chem Phys 2022; 156:194110. [PMID: 35597663 DOI: 10.1063/5.0091636] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The restricted active space spin-flip (RAS-SF) formalism is a particular form of single-reference configuration interaction that can describe some forms of strong correlation at a relatively low cost and which has recently been formulated for the description of charge-transfer excited states. Here, we introduce both equilibrium and nonequilibrium versions of a state-specific solvation correction for vertical transition energies computed using RAS-SF wave functions, based on the framework of a polarizable continuum model (PCM). Ground-state polarization is described using the solvent's static dielectric constant and in the nonequilibrium solvation approach that polarization is modified upon vertical excitation using the solvent's optical dielectric constant. Benchmark calculations are reported for well-studied models of photo-induced charge transfer, including naphthalene dimer, C2H4⋯C2F4, pentacene dimer, and perylene diimide (PDI) dimer, several of which are important in organic photovoltaic applications. For the PDI dimer, we demonstrate that the charge-transfer character of the excited states is enhanced in the presence of a low-dielectric medium (static dielectric constant ɛ0 = 3) as compared to a gas-phase calculation (ɛ0 = 1). This stabilizes mechanistic traps for singlet fission and helps to explain experimental singlet fission rates. We also examine the effects of nonequilibrium solvation on charge-separated states in an intramolecular singlet fission chromophore, where we demonstrate that the energetic ordering of the states changes as a function of solvent polarity. The RAS-SF + PCM methodology that is reported here provides a framework to study charge-separated states in solution and in photovoltaic materials.
Collapse
Affiliation(s)
- Bushra Alam
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Hanjie Jiang
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Paul M Zimmerman
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
11
|
Picconi D. Quantum dynamics of the photoinduced charge separation in a symmetric donor–acceptor–donor triad: The role of vibronic couplings, symmetry and temperature. J Chem Phys 2022; 156:184105. [DOI: 10.1063/5.0089887] [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
The photoinduced charge separation in a symmetric donor–acceptor–donor (D–A–D) triad is studied quantum mechanically using a realistic diabatic vibronic coupling model. The model includes a locally excited DA*D state and two charge-transfer states D+A−D and DA−D+ and is constructed according to a procedure generally applicable to semirigid D–A–D structures and based on energies, forces, and force constants obtained by quantum chemical calculations. In this case, the electronic structure is described by time-dependent density functional theory, and the corrected linear response is used in conjunction with the polarizable continuum model to account for state-specific solvent effects. The multimode dynamics following the photoexcitation to the locally excited state are simulated by the hybrid Gaussian-multiconfigurational time-dependent Hartree method, and temperature effects are included using thermo field theory. The dynamics are connected to the transient absorption spectrum obtained in recent experiments, which is simulated and fully assigned from first principles. It is found that the charge separation is mediated by symmetry-breaking vibrations of relatively low frequency, which implies that temperature should be accounted for to obtain reliable estimates of the charge transfer rate.
Collapse
Affiliation(s)
- David Picconi
- Institute of Chemistry, University of Potsdam, Karl-Liebknecht-Str. 24-25, D-14476 Potsdam-Golm, Germany and Institute of Physics and Astronomy, University of Potsdam, Karl-Liebknecht-Str. 24-25, D-14476 Potsdam-Golm, Germany
| |
Collapse
|
12
|
Eschenbach P, Artiukhin DG, Neugebauer J. Multi-state formulation of the frozen-density embedding quasi-diabatization approach. J Chem Phys 2021; 155:174104. [PMID: 34749503 DOI: 10.1063/5.0070486] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a multi-state implementation of the recently developed frozen-density embedding diabatization (FDE-diab) methodology [D. G. Artiukhin and J. Neugebauer, J. Chem. Phys. 148, 214104 (2018)] in the Serenity program. The new framework extends the original approach such that any number of charge-localized quasi-diabatic states can be coupled, giving an access to calculations of ground and excited state spin-density distributions as well as to excitation energies. We show that it is possible to obtain results similar to those from correlated wave function approaches such as the complete active space self-consistent field method at much lower computational effort. Additionally, we present a series of approximate computational schemes, which further decrease the overall computational cost and systematically converge to the full FDE-diab solution. The proposed methodology enables computational studies on spin-density distributions and related properties for large molecular systems of biochemical interest.
Collapse
Affiliation(s)
- Patrick Eschenbach
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| | - Denis G Artiukhin
- Department of Chemistry, Aarhus Universitet, DK-8000 Aarhus, Denmark
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| |
Collapse
|
13
|
Artiukhin DG, Eschenbach P, Matysik J, Neugebauer J. Theoretical Assessment of Hinge-Type Models for Electron Donors in Reaction Centers of Photosystems I and II as well as of Purple Bacteria. J Phys Chem B 2021; 125:3066-3079. [PMID: 33749260 DOI: 10.1021/acs.jpcb.0c10656] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Hinge-type molecular models for electron donors in reaction centers of Photosystems I and II and purple bacteria were investigated using a two-state computational approach based on frozen-density embedding (FDE). This methodology, dubbed FDE-diab, is known to avoid consequences of the self-interaction error as far as intermolecular phenomena are concerned, which allows a prediction of qualitatively correct spin densities for large biomolecular systems. The calculated spin density distributions are in a good agreement with available experimental results and demonstrated a very high sensitivity to changes in the relative orientation of cofactors and amino acid protonation states. This allows a validation of the previously proposed hinge-type models providing hints on possible protonation states of axial histidine molecules.
Collapse
Affiliation(s)
- Denis G Artiukhin
- Department of Chemistry, Aarhus Universitet, DK-8000 Aarhus, Denmark
| | - Patrick Eschenbach
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| | - Jörg Matysik
- Institut für Analytische Chemie, Universität Leipzig, Linnéstr. 3, 04103 Leipzig, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| |
Collapse
|
14
|
Maldonado AM, Hagiwara S, Choi TH, Eckert F, Schwarz K, Sundararaman R, Otani M, Keith JA. Quantifying Uncertainties in Solvation Procedures for Modeling Aqueous Phase Reaction Mechanisms. J Phys Chem A 2021; 125:154-164. [PMID: 33393781 DOI: 10.1021/acs.jpca.0c08961] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Computational quantum chemistry provides fundamental chemical and physical insights into solvated reaction mechanisms across many areas of chemistry, especially in homogeneous and heterogeneous renewable energy catalysis. Such reactions may depend on explicit interactions with ions and solvent molecules that are nontrivial to characterize. Rigorously modeling explicit solvent effects with molecular dynamics usually brings steep computational costs while the performance of continuum solvent models such as polarizable continuum model (PCM), charge-asymmetric nonlocally determined local-electric (CANDLE), conductor-like screening model for real solvents (COSMO-RS), and effective screening medium method with the reference interaction site model (ESM-RISM) are less well understood for reaction mechanisms. Here, we revisit a fundamental aqueous hydride transfer reaction-carbon dioxide (CO2) reduction by sodium borohydride (NaBH4)-as a test case to evaluate how different solvent models perform in aqueous phase charge migrations that would be relevant to renewable energy catalysis mechanisms. For this system, quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations almost exactly reproduced energy profiles from QM simulations, and the Na+ counterion in the QM/MM simulations plays an insignificant role over ensemble averaged trajectories that describe the reaction pathway. However, solvent models used on static calculations gave much more variability in data depending on whether the system was modeled using explicit solvent shells and/or the counterion. We pinpoint this variability due to unphysical descriptions of charge-separated states in the gas phase (i.e., self-interaction errors), and we show that using more accurate hybrid functionals and/or explicit solvent shells lessens these errors. This work closes with recommended procedures for treating solvation in future computational efforts in studying renewable energy catalysis mechanisms.
Collapse
Affiliation(s)
- Alex M Maldonado
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Satoshi Hagiwara
- National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba 305-8568, Japan
| | - Tae Hoon Choi
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Frank Eckert
- Dassault Systèmes Deutschland GmbH, Imbacher Weg 46, 51379 Leverkusen, Germany
| | - Kathleen Schwarz
- Material Measurement Laboratory, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States
| | - Ravishankar Sundararaman
- Department of Materials Science and Engineering, Rensselaer Polytechnic Institute, Troy, New York 12180, United States
| | - Minoru Otani
- National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba 305-8568, Japan
| | - John A Keith
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| |
Collapse
|
15
|
Tölle J, Cupellini L, Mennucci B, Neugebauer J. Electronic couplings for photo-induced processes from subsystem time-dependent density-functional theory: The role of the diabatization. J Chem Phys 2020; 153:184113. [PMID: 33187428 DOI: 10.1063/5.0022677] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
Subsystem time-dependent density-functional theory (sTDDFT) making use of approximate non-additive kinetic energy (NAKE) functionals is known to be capable of describing excitation energy transfer processes in a variety of applications. Here, we show that sTDDFT, especially when combined with projection-based embedding (PbE), can be employed for the entire range of photo-induced electronic couplings essential for modeling photophysical properties of complex chemical and biological systems and therefore represents a complete toolbox for this class of problems. This means that it is capable of capturing the interaction/coupling associated with local- and charge-transfer (CT) excitons. However, this requires the choice of a reasonable diabatic basis. We therefore propose different diabatization strategies of the virtual orbital space in PbE-sTDDFT and show how CT excitations can be included in sTDDFT using NAKE functionals via a phenomenological approach. Finally, these electronic couplings are compared to couplings from a multistate fragment excitation difference (FED)-fragment charge difference (FCD) diabatization procedure. We show that both procedures, multistate FED-FCD and sTDDFT (with the right diabatization procedure chosen), lead to an overall good agreement for the electronic couplings, despite differences in their general diabatization strategy. We conclude that the entire range of photo-induced electronic couplings can be obtained using sTDDFT (with the right diabatization procedure chosen) in a black-box manner.
Collapse
Affiliation(s)
- Johannes Tölle
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster Corrensstraße 40, 48149 Münster, Germany
| | - Lorenzo Cupellini
- Università di Pisa, Dipartimento di Chimica e Chimica Industriale, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Benedetta Mennucci
- Università di Pisa, Dipartimento di Chimica e Chimica Industriale, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster Corrensstraße 40, 48149 Münster, Germany
| |
Collapse
|
16
|
Jiang H, Zimmerman PM. Charge transfer via spin flip configuration interaction: Benchmarks and application to singlet fission. J Chem Phys 2020; 153:064109. [DOI: 10.1063/5.0018267] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Affiliation(s)
- Hanjie Jiang
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Paul M. Zimmerman
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| |
Collapse
|
17
|
Artiukhin DG, Eschenbach P, Neugebauer J. Computational Investigation of the Spin-Density Asymmetry in Photosynthetic Reaction Center Models from First Principles. J Phys Chem B 2020; 124:4873-4888. [DOI: 10.1021/acs.jpcb.0c02827] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
| | - Patrick Eschenbach
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| |
Collapse
|
18
|
Abstract
By invoking a divide-and-conquer strategy, subsystem DFT dramatically reduces the computational cost of large-scale, ab initio electronic structure simulations of molecules and materials. The central ingredient setting subsystem DFT apart from Kohn-Sham DFT is the nonadditive kinetic energy functional (NAKE). Currently employed NAKEs are at most semilocal (i.e., they only depend on the electron density and its gradient), and as a result of this approximation, so far large-scale simulations only included systems composed of weakly interacting subsystems. In this work, we advance the state-of-the-art by introducing fully nonlocal NAKEs in subsystem DFT simulations for the first time. A benchmark analysis based on the S22-5 test set shows that nonlocal NAKEs considerably improve the computed interaction energies and electron densities compared to commonly employed GGA NAKEs, especially when increasing intersubsystem electron density overlap is considered. Most importantly, we resolve the long-standing problem of too attractive interaction energy curves typically resulting from the use of GGA NAKEs.
Collapse
Affiliation(s)
- Wenhui Mi
- Department of Chemistry , Rutgers University , Newark , New Jersey 07102 , United States
- Department of Physics , 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
|
Periyanagounder D, Wei TC, Li TY, Lin CH, Gonçalves TP, Fu HC, Tsai DS, Ke JJ, Kuo HW, Huang KW, Lu N, Fang X, He JH. Fast-Response, Highly Air-Stable, and Water-Resistant Organic Photodetectors Based on a Single-Crystal Pt Complex. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2020; 32:e1904634. [PMID: 31736151 DOI: 10.1002/adma.201904634] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/19/2019] [Revised: 10/19/2019] [Indexed: 06/10/2023]
Abstract
Organic semiconductors demonstrate several advantages over conventional inorganic materials for novel electronic and optoelectronic applications, including molecularly tunable properties, flexibility, low-cost, and facile device integration. However, before organic semiconductors can be used for the next-generation devices, such as ultrafast photodetectors (PDs), it is necessary to develop new materials that feature both high mobility and ambient stability. Toward this goal, a highly stable PD based on the organic single crystal [PtBr2 (5,5'-bis(CF3 CH2 OCH2 )-2,2'-bpy)] (or "Pt complex (1o)") is demonstrated as the active semiconductor channel-a material that features a lamellar molecular structure and high-quality, intraligand charge transfer. Benefitting from its unique crystal structure, the Pt-complex (1o) device exhibits a field-effect mobility of ≈0.45 cm2 V-1 s-1 without loss of significant performance under ambient conditions even after 40 days without encapsulation, as well as immersion in distilled water for a period of 24 h. Furthermore, the device features a maximum photoresponsivity of 1 × 103 A W-1 , a detectivity of 1.1 × 1012 cm Hz1/2 W-1 , and a record fast response/recovery time of 80/90 µs, which has never been previously achieved in other organic PDs. These findings strongly support and promote the use of the single-crystal Pt complex (1o) in next-generation organic optoelectronic devices.
Collapse
Affiliation(s)
- Dharmaraj Periyanagounder
- Computer, Electrical, and Mathematical Sciences and Engineering Division, King Abdullah University of Science & Technology, Thuwal, 23955-6900, Saudi Arabia
| | - Tzu-Chiao Wei
- Computer, Electrical, and Mathematical Sciences and Engineering Division, King Abdullah University of Science & Technology, Thuwal, 23955-6900, Saudi Arabia
| | - Ting-You Li
- Computer, Electrical, and Mathematical Sciences and Engineering Division, King Abdullah University of Science & Technology, Thuwal, 23955-6900, Saudi Arabia
- Department of Molecular Science and Engineering, National Taipei University of Technology, Taipei, 106, Taiwan, ROC
| | - Chun-Ho Lin
- Computer, Electrical, and Mathematical Sciences and Engineering Division, King Abdullah University of Science & Technology, Thuwal, 23955-6900, Saudi Arabia
| | - Théo Piechota Gonçalves
- KAUST Catalyst Centre, Physical Science and Engineering Division, King Abdullah University of Science & Technology, Thuwal, 23955-6900, Saudi Arabia
| | - Hui-Chun Fu
- Computer, Electrical, and Mathematical Sciences and Engineering Division, King Abdullah University of Science & Technology, Thuwal, 23955-6900, Saudi Arabia
| | - Dung-Sheng Tsai
- Computer, Electrical, and Mathematical Sciences and Engineering Division, King Abdullah University of Science & Technology, Thuwal, 23955-6900, Saudi Arabia
| | - Jr-Jian Ke
- Computer, Electrical, and Mathematical Sciences and Engineering Division, King Abdullah University of Science & Technology, Thuwal, 23955-6900, Saudi Arabia
| | - Hung-Wei Kuo
- Institute of Organic and Polymeric Materials, National Taipei University of Technology, Taipei, 106, Taiwan, ROC
| | - Kuo-Wei Huang
- KAUST Catalyst Centre, Physical Science and Engineering Division, King Abdullah University of Science & Technology, Thuwal, 23955-6900, Saudi Arabia
| | - Norman Lu
- Department of Molecular Science and Engineering, National Taipei University of Technology, Taipei, 106, Taiwan, ROC
- Institute of Organic and Polymeric Materials, National Taipei University of Technology, Taipei, 106, Taiwan, ROC
| | - Xiaosheng Fang
- Department of Materials Science, Fudan University, Shanghai, 200433, China
| | - Jr-Hau He
- Computer, Electrical, and Mathematical Sciences and Engineering Division, King Abdullah University of Science & Technology, Thuwal, 23955-6900, Saudi Arabia
- Department of Materials Science and Engineering, City University of Hong Kong, Kowloon, Hong Kong SAR, 999077, China
| |
Collapse
|
20
|
Tölle J, Böckers M, Niemeyer N, Neugebauer J. Inter-subsystem charge-transfer excitations in exact subsystem time-dependent density-functional theory. J Chem Phys 2019; 151:174109. [DOI: 10.1063/1.5121908] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Johannes Tölle
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| | - Michael Böckers
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| | - Niklas Niemeyer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| |
Collapse
|
21
|
Schulz A, Jacob CR. Description of intermolecular charge transfer with subsystem density-functional theory. J Chem Phys 2019; 151:131103. [PMID: 31594348 DOI: 10.1063/1.5125218] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Efficient quantum-chemical methods that are able to describe intermolecular charge transfer are crucial for modeling organic semiconductors. However, the correct description of intermolecular charge transfer with density-functional theory (DFT) is hampered by the fractional charge error of approximate exchange-correlation (xc) functionals. Here, we investigate the charge transfer induced by an external electric field in a tetrathiafulvalene-tetracyanoquinodimethane (TTF-TCNQ) complex as a test case. For this seemingly simple model system, a supermolecular DFT treatment fails with most conventional xc functionals. Here, we present an extension of subsystem DFT to subsystems with a fractional number of electrons. We show that within such a framework, it becomes possible to overcome the fractional charge error by enforcing the correct dependence of each subsystem's total energy on the subsystem's fractional charge. Such a subsystem DFT approach allows for a correct description of the intermolecular charge transfer in the TTF-TCNQ model complex. The approach presented here can be generalized to larger molecular aggregates and will thus allow for modeling organic semiconductor materials accurately and efficiently.
Collapse
Affiliation(s)
- Anika Schulz
- Technische Universität Braunschweig, Institute of Physical and Theoretical Chemistry, Gaußstr. 17, 38106 Braunschweig, Germany
| | - Christoph R Jacob
- Technische Universität Braunschweig, Institute of Physical and Theoretical Chemistry, Gaußstr. 17, 38106 Braunschweig, Germany
| |
Collapse
|
22
|
Artiukhin DG, Neugebauer J. Frozen-density embedding as a quasi-diabatization tool: Charge-localized states for spin-density calculations. J Chem Phys 2018; 148:214104. [DOI: 10.1063/1.5023290] [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)
- Denis G. Artiukhin
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| |
Collapse
|
23
|
Sato K, Pradhan E, Asahi R, Akimov AV. Charge transfer dynamics at the boron subphthalocyanine chloride/C60 interface: non-adiabatic dynamics study with Libra-X. Phys Chem Chem Phys 2018; 20:25275-25294. [DOI: 10.1039/c8cp03841d] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The Libra-X software for non-adiabatic molecular dynamics is reported. It is used to comprehensively study the charge transfer dynamics at the boron subphtalocyanine chloride (SubPc)/fullerene (C60) interface.
Collapse
Affiliation(s)
- Kosuke Sato
- Toyota Central Research and Development Laboratories, Inc
- Nagakute
- Japan
| | - Ekadashi Pradhan
- Department of Chemistry
- University at Buffalo
- The State University of New York
- New York 14260-3000
- USA
| | - Ryoji Asahi
- Toyota Central Research and Development Laboratories, Inc
- Nagakute
- Japan
| | - Alexey V. Akimov
- Department of Chemistry
- University at Buffalo
- The State University of New York
- New York 14260-3000
- USA
| |
Collapse
|
24
|
Maitra NT. Charge transfer in time-dependent density functional theory. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2017; 29:423001. [PMID: 28766507 DOI: 10.1088/1361-648x/aa836e] [Citation(s) in RCA: 58] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Charge transfer plays a crucial role in many processes of interest in physics, chemistry, and bio-chemistry. In many applications the size of the systems involved calls for time-dependent density functional theory (TDDFT) to be used in their computational modeling, due to its unprecedented balance between accuracy and efficiency. However, although exact in principle, in practise approximations must be made for the exchange-correlation functional in this theory, and the standard functional approximations perform poorly for excitations which have a long-range charge-transfer component. Intense progress has been made in developing more sophisticated functionals for this problem, which we review. We point out an essential difference between the properties of the exchange-correlation kernel needed for an accurate description of charge-transfer between open-shell fragments and between closed-shell fragments. We then turn to charge-transfer dynamics, which, in contrast to the excitation problem, is a highly non-equilibrium, non-perturbative, process involving a transfer of one full electron in space. This turns out to be a much more challenging problem for TDDFT functionals. We describe dynamical step and peak features in the exact functional evolving over time, that are missing in the functionals currently used. The latter underestimate the amount of charge transferred and manifest a spurious shift in the charge transfer resonance position. We discuss some explicit examples.
Collapse
Affiliation(s)
- Neepa T Maitra
- Department of Physics and Astronomy, Hunter College and the Physics Program at the Graduate Center of the City University of New York, 695 Park Avenue, NY 10065, United States of America
| |
Collapse
|
25
|
Oberhofer H, Reuter K, Blumberger J. Charge Transport in Molecular Materials: An Assessment of Computational Methods. Chem Rev 2017. [PMID: 28644623 DOI: 10.1021/acs.chemrev.7b00086] [Citation(s) in RCA: 191] [Impact Index Per Article: 27.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
The booming field of molecular electronics has fostered a surge of computational research on electronic properties of organic molecular solids. In particular, with respect to a microscopic understanding of transport and loss mechanisms, theoretical studies assume an ever-increasing role. Owing to the tremendous diversity of organic molecular materials, a great number of computational methods have been put forward to suit every possible charge transport regime, material, and need for accuracy. With this review article we aim at providing a compendium of the available methods, their theoretical foundations, and their ranges of validity. We illustrate these through applications found in the literature. The focus is on methods available for organic molecular crystals, but mention is made wherever techniques are suitable for use in other related materials such as disordered or polymeric systems.
Collapse
Affiliation(s)
- Harald Oberhofer
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München , Lichtenbergstrasse 4, D-85747 Garching, Germany
| | - Karsten Reuter
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München , Lichtenbergstrasse 4, D-85747 Garching, Germany
| | - Jochen Blumberger
- Department of Physics and Astronomy, University College London , Gower Street, London WC1E 6BT, United Kingdom.,Institute for Advanced Study, Technische Universität München , Lichtenbergstrasse 2 a, D-85748 Garching, Germany
| |
Collapse
|
26
|
Senn F, Park YC. Constricted variational density functional theory for spatially clearly separated charge-transfer excitations. J Chem Phys 2016; 145:244108. [PMID: 28049328 DOI: 10.1063/1.4972231] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Constricted Variational Density Functional Theory (CV-DFT) is known to be one of the successful methods in predicting charge-transfer excitation energies. In this paper, we apply the CV-DFT method to the well-known model systems ethylene-tetrafluoroethylene (C2H4 × C2F4) and the zincbacteriochlorin-bacteriochlorin complex (ZnBC-BC). The analysis of the CV-DFT energies enables us to understand the -1/R charge-transfer behaviour in CV-DFT for large separation distances R. With this we discuss the importance of orbital relaxations using the relaxed version of CV(∞)-DFT, the R-CV(∞)-DFT method. Possible effects of the optimization of the transition matrix for the relaxed self-consistent field version of CV(∞)-DFT, RSCF-CV(∞)-DFT in the case of large fragment separations are shown and we introduce two possible gradient restrictions to avoid the unwanted admixing of other transitions.
Collapse
Affiliation(s)
- Florian Senn
- Department of Chemistry, University of Calgary, 2500 University Drive NorthWest, Calgary, Alberta T2N 1N4, Canada
| | - Young Choon Park
- Department of Chemistry, Korea Advanced Institute of Science and Technology, Daejeon 305-701, South Korea
| |
Collapse
|
27
|
Schlüns D, Franchini M, Götz AW, Neugebauer J, Jacob CR, Visscher L. Analytical gradients for subsystem density functional theory within the slater-function-based amsterdam density functional program. J Comput Chem 2016; 38:238-249. [DOI: 10.1002/jcc.24670] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2016] [Revised: 10/14/2016] [Accepted: 10/19/2016] [Indexed: 11/09/2022]
Affiliation(s)
- Danny Schlüns
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster; Corrensstraße 40 Münster 48149 Germany
| | - Mirko Franchini
- Amsterdam Center for Multiscale Modeling; Vrije Universiteit Amsterdam; De Boelelaan 1083 HV Amsterdam 1081 The Netherlands
- Scientific Computing & Modelling NV; Vrije Universiteit; Theoretical Chemistry, De Boelelaan 1083 HV Amsterdam 1081 The Netherlands
| | - Andreas W. Götz
- San Diego Supercomputer Center; University of California San Diego; 9500 Gilman Drive La Jolla, California 92093-0505 USA
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster; Corrensstraße 40 Münster 48149 Germany
| | - Christoph R. Jacob
- Institute of Physical and Theoretical Chemistry; TU Braunschweig, Hans-Sommer-Straße 10 Braunschweig 38106 Germany
| | - Lucas Visscher
- Amsterdam Center for Multiscale Modeling; Vrije Universiteit Amsterdam; De Boelelaan 1083 HV Amsterdam 1081 The Netherlands
| |
Collapse
|
28
|
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
|
29
|
Abstract
The design of optimal light-harvesting (supra)molecular systems and materials is one of the most challenging frontiers of science. Theoretical methods and computational models play a fundamental role in this difficult task, as they allow the establishment of structural blueprints inspired by natural photosynthetic organisms that can be applied to the design of novel artificial light-harvesting devices. Among theoretical strategies, the application of quantum chemical tools represents an important reality that has already reached an evident degree of maturity, although it still has to show its real potentials. This Review presents an overview of the state of the art of this strategy, showing the actual fields of applicability but also indicating its current limitations, which need to be solved in future developments.
Collapse
Affiliation(s)
- Carles Curutchet
- Departament de Fisicoquímica, Facultat de Farmàcia, Universitat de Barcelona , Av. Joan XXIII s/n, 08028 Barcelona, Spain
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, University of Pisa , via G. Moruzzi 13, 56124 Pisa, Italy
| |
Collapse
|
30
|
Abstract
Constrained Subsystem Density Fucntional Theory (CSDFT) allows to compute diabatic states for charge transfer reactions using the machinery of the constrained DFT method, and at the same time is able to embed such diabatic states in a molecular environment via a subsystem DFT scheme.
Collapse
Affiliation(s)
- Pablo Ramos
- Department of Chemistry
- Rutgers University
- Newark
- USA
| | | |
Collapse
|
31
|
Shrestha K, Jakubikova E. Ground-State Electronic Structure of RC-LH1 and LH2 Pigment Assemblies of Purple Bacteria via the EBF-MO Method. J Phys Chem A 2015. [DOI: 10.1021/acs.jpca.5b05644] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Kushal Shrestha
- Department of Chemistry, North Carolina State University, Raleigh, North Carolina 27695, United States
| | - Elena Jakubikova
- Department of Chemistry, North Carolina State University, Raleigh, North Carolina 27695, United States
| |
Collapse
|
32
|
|
33
|
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
|
34
|
Ramos P, Papadakis M, Pavanello M. Performance of Frozen Density Embedding for Modeling Hole Transfer Reactions. J Phys Chem B 2015; 119:7541-57. [DOI: 10.1021/jp511275e] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Affiliation(s)
- Pablo Ramos
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| | - Markos Papadakis
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| | - Michele Pavanello
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| |
Collapse
|
35
|
Genova A, Ceresoli D, Pavanello M. Periodic subsystem density-functional theory. J Chem Phys 2014; 141:174101. [DOI: 10.1063/1.4897559] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|