Subscriber access provided by UNIV OF WISCONSIN - MADISON
Surfaces, Interfaces, and Catalysis; Physical Properties of Nanomaterials and Materials
Nonadiabatic Dynamics Simulations Reveal Distinct Effects of the Thickness of PTB7 on Interfacial Electron and Hole Transfer Dynamics in PTB7@MoS Heterostructures 2
Xiangyang Liu, Wen-Kai Chen, Weihai Fang, and Ganglong Cui J. Phys. Chem. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.jpclett.9b01066 • Publication Date (Web): 14 May 2019 Downloaded from http://pubs.acs.org on May 14, 2019
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry Letters
Nonadiabatic Dynamics Simulations Reveal Distinct Effects of the Thickness of PTB7 on Interfacial Electron and Hole Transfer Dynamics in PTB7@MoS2 Heterostructures Xiang-Yang Liu, Wen-Kai Chen, Wei-Hai Fang, and Ganglong Cui*
Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
Email: ganglong.cui@bnu.edu.cn
Abstract Mixed-dimensional hybrid heterostructures have attracted a lot of experimental attention because they can provide an ideal charge-separated interface for optoelectronic and photonic applications. In this contribution we have employed first-principles DFT calculations and nonadiabatic dynamics simulations to explore photoinduced interfacial electron and hole transfer processes in two PTB7-nL@MoS2 models (n=1 and 5). The interfacial electron transfer is found to be ultrafast and completes within ca. 10 fs in both PTB7-1L@MoS2 and PTB75L@MoS2 models, which demonstrates that the electron transfer is not sensitive to the thickness of the PTB7 polymer. Differently, the interfacial hole transfer is sensitive to the thickness of the PTB7 polymer. The transfer time is estimated to be ca. 70 ps in PTB71L@MoS2 while it is significantly accelerated to ca. 1 ps in PTB7-5L@MoS2. Finally, we have found that the electron transfer is mainly controlled by adiabatic electron evolution; whereas, 1 ACS Paragon Plus Environment
The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 2 of 24
in the hole transfer, nonadiabatic hoppings play a dominant role. These findings are useful for the design of excellent charge-separated interfaces of mixed-dimensional TMD-based heterojunctions for a variety of optoelectronic applications.
Electronic and optical properties of two-dimensional transition metal dichalcogenides (TMDs) have received much attention due to their potential applications in field-effect transistors, solar cells, and related optoelectronic devices.1-24 In particular, mixed-dimensional hybrid heterostructures combining TMDs with 0- and 1-dimensional materials through weak van der Waals interaction, for example with conjugated polymers, have attracted a lot of experimental attention because they could provide an ideal charge-separated interface for optoelectronic and photonic applications.15,25-38 Recently, Shastry et al. have presented an experimental realization of a large-area type-II photovoltaic heterojunction in which the MoS2 monolayer runs as the primary adsorber and the PTB7 polymer does as the donor.35 This heterostructure exhibits a tunable photoluminescence intensity as a function of the thickness of the PTB7 layer, which finally enables the quenching of the MoS2 photoluminescence.35 But, 2 ACS Paragon Plus Environment
Page 3 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry Letters
the underlying photoinduced carrier dynamics is not clear, which motivates further experimental studies. Zhong et al. have employed time-resolved transient absorption spectroscopies to investigate the ultrafast electron and hole transfer dynamics of the mixeddimensional PTB7@MoS2 heterojunction and concluded that the photoinduced electron transfer from PTB7 to electronically hot states of MoS2 occurs less than 250 fs.27 Followed are an exciton diffusion-limited electron transfer still from PTB7 to MoS2 and a hole transfer from MoS2 to PTB7 of less than 3 ps.27 Although these experimental studies provide preliminary spectroscopic and dynamic information on the carrier dynamics in PTB7@MoS 2, the photoinduced interfacial electron and hole transfer dynamics at the atomistic level is far from complete understanding and still requires theoretical studies.39-47
In this contribution we have constructed two heterostructure models to explore the photoinduced electron and hole transfer dynamics at the interface between PTB7 and MoS 2 by using both electronic structure calculations and nonadiabatic dynamics simulations at the DFT level.48-52 The first model is the simplest one in that it just includes a PTB7 layer, referred to as PTB7-1L@MoS2 hereinafter, which is used to explore the intrinsic interface properties. In order to explore the effects of the thickness of the PTB7 layers on the interfacial electron and hole transfer dynamics, we have constructed the second PTB7-5L@MoS2 heterojunction model, which includes five PTB7 layers, referred to as PTB7-5L@MoS2.
3 ACS Paragon Plus Environment
The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Figure 1.
Page 4 of 24
(left) PBE+D3 optimized PTB7-1L@MoS2 structures and (right) PTB7-5L@MoS2
structures (0 and 300 K). Also shown are inter-distances of the PTB7 layer to the MoS2 surface.
These initially constructed models are first optimized with the PBE+D3 method implemented in the Quickstep/CP2K package (see supporting information for simulation details).53-61 The resultant stable structures are shown in Fig. 1. In the PTB7-1L@MoS2 heterojunction, the PTB7 layer is physically adsorbed on the monolayer MoS 2 surface in a parallel way without forming any chemical bonds at the interface. The adsorption energy is estimated to be 2.12 eV at the HSE06 level.62-66 At this adsorption configuration, the monolayer MoS2 is essentially planar while the PTB7 layer is a little twisted because the π conjugation is broken at the boundary of the two conjugated moieties. The inter-distance between PTB7 and MoS2 is calculated to be 3.531 Å at 0 K. Upon heated to 300 K, the PTB7-1L@MoS2
4 ACS Paragon Plus Environment
Page 5 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry Letters
heterostructure does not change visibly except a subtly larger inter-distance between PTB7 and MoS2, 3.591 Å . In the PTB7-5L@MoS2 heterostructure, the five PTB7 layers are arranged in a nearly parallel way through bilateral non-covalent interaction and they are also physically adsorbed on the monolayer MoS2 surface. At the HSE06 level the corresponding adsorption energy is calculated to be 2.41 eV. Similarly, there is a small twisting due to the breaking of the π conjugation. Compared with those in the PTB7-1L@MoS2 heterojunction, the interdistance between the lowest PTB7 layer and the MoS2 surface is calculated to be 3.483 Å at 0 K (3.623 Å at 300 K).
Figure 2. DOS and PDOS of (A) PTB7-1L@MoS2 and (B) PTB7-5L@MoS2 calculated by the HSE06 functional. Also shown are the components of each state with respect to PTB7 and MoS2 (C and D). See supporting information for those by the PBE+D3 functional. 5 ACS Paragon Plus Environment
The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The density of states (DOS) and projected DOS (PDOS) are very helpful for qualitatively understanding the interfacial electron and hole transfer dynamics and we have thus employed the HSE06 method to calculate the corresponding DOS and PDOS of both PTB7-1L@MoS2 and PTB7-5L@MoS2 heterostructures. The panel A of Fig. 2 shows the DOS and PDOS of PTB7-1L@MoS2 and the panel C depicts the components of each state with respect to PTB7 and MoS2. To discuss more conveniently CBM and VBM are used to represent the conduction band minimum and the valence band maximum of MoS 2 while HOMO and LUMO are used for the highest occupied molecular orbital and the lowest unoccupied molecular orbital of PTB7. Clearly, the two highest occupied states i.e. HOMO and HOMO-1 are all from the PTB7 layer and the lowest unoccupied state i.e. CBM is from the MoS2 monolayer; thus, PTB71L@MoS2 is a type-II heterojunction. At the HSE06 level, CBM and VBM of MoS2 in this heterojunction is calculated to be 0.72 and -1.76 eV, respectively, from which the energy gap is estimated to be 2.48 eV. LUMO, HOMO, HOMO-1, and HOMO-2 energies of PTB7 are calculated to be 1.15, -0.70, -1.32, and -2.30 eV at the same computational level, respectively. HOMO-LUMO energy gap is predicted to be 1.85 eV. Similarly, the panels B and D of Fig. 2 show the counterparts for PTB7-5L@MoS2. Obviously, PTB7-5L@MoS2 is also a type-II heterojunction. CBM and VBM of MoS2 are calculated to be 0.565 and -1.909 eV at the HSE06 level, from which CBM-VBM energy gap is 2.47 eV, close to 2.48 eV in the above PTB71L@MoS2. HOMO and LUMO are estimated to be -0.548 and 0.813 eV; thus, HOMO-LUMO energy gap of the PTB7 layers, 1.36 eV, becomes smaller than 1.85 eV in PTB7-1L@MoS2. 6 ACS Paragon Plus Environment
Page 6 of 24
Page 7 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry Letters
Furthermore, different from PTB7-1L@MoS2, the DOS from the PTB7 layers is significantly increased in PTB7-5L@MoS2. In the former, there exist only two occupied states of PTB7 between CBM and VBM of MoS2; but, in the latter, there are much more states in-between CBM and VBM of MoS2. Of course, this can be easily understood because of more PTB7 layers involved. One can also find that the lowest occupied state of PTB7 between CBM and VBM of MoS2 is much closer to its VBM in PTB7-5L@MoS2 than that in PTB7-1L@MoS2, which implies that the interfacial hole transfer from MoS2 to PTB7 is much easier for PTB7-5L@MoS2 due to a very smaller energy gap.
Figure 3. Spatial distributions of selected states involved in the interfacial electron and hole transfer dynamics in (left) PTB7-1L@MoS2 and (right) PTB7-5L@MoS2 models. 7 ACS Paragon Plus Environment
The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Next we are dedicated to exploring the initially populated electron and hole states in PTB7-1L@MoS2 and PTB7-5L@MoS2, which are useful for studying the following interfacial electron and hole transfer dynamics. These electron and hole states are generated due to the electron transition; thus, we first explore the most possible electron transitions in the FranckCondon region. First, we focus on the electron transition of PTB7-1L@MoS2. At the HSE06 level, it is found that the electron transitions from state-385 (HOMO) to state-400 (LUMO), state-382 to state-386 (CBM), and state-381 to state-387 have the largest oscillator strengths, whose excitation energies [oscillator strengths] are 1.86 eV [0.56], 2.62 eV [0.53], and 2.62 eV [0.53]. The first electron transition of 1.86 eV corresponds to the excitation within the PTB7 layer, close to the experimental value of 1.77 eV [700 nm];27 whereas, the latter two of 2.62 eV are related to the excitation within MoS2 (exp. 2.34 eV, 530 nm).27 As a result, when an electron is promoted from HOMO to LUMO of PTB7 upon photoexcitation, a hole will be generated in HOMO. In this situation, the interfacial hole transfer is not allowed because HOMO is much higher than VBM of MoS2 in energy. On contrary, the interfacial electron transfer is possible in that LUMO is much higher than CBM of MoS2 in energy (see Fig. 2). In parallel, if an electron is excited into either state-386 or state-387 within MoS2, a hole will occupy either state-381 or state-382. In this case, the interfacial electron transfer is not favorable because both state-386 and state-387 are much lower than LUMO of PTB7; instead, the interfacial hole transfer becomes possible owing to the fact that the hole states 381 and 382 are lower than the two occupied states of PTB7, i.e. states 384 and 385 (see Fig. 2). The left panel of Fig. 3 shows the spatial distribution of selected states involved in the interfacial electron and hole transfer 8 ACS Paragon Plus Environment
Page 8 of 24
Page 9 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry Letters
of PTB7-1L@MoS2. It is clear that both HOMO and LUMO are mainly located on PTB7 and the others are on MoS2.
Second, we look at the electronic transition of PTB7-5L@MoS2. As with the situation of PTB7-1L@MoS2, there are two kinds of excitations, which are separately located on PTB7 and MoS2. In the former the electronic transition from state-677 (HOMO) to state-713 (LUMO+4) has the largest oscillator strength with an excitation energy of 1.92 eV, also close to the experimental value of 1.77 eV [700 nm].27 In the latter the electronic transitions from state-664 to state-678 (CBM) and state-665 to state-679 have the largest oscillator strengths with same 2.63 eV excitation energies [exp. 2.34 eV, 530 nm].27 Clearly, the electronic transition of state-677 (HOMO) to state-713 (LUMO+4) within PTB7 excites an electron into LUMO+4 and leaves a hole in HOMO. Since HOMO of PTB7 is much higher than VBM of MoS2, the interfacial hole transfer is not allowed thermodynamically; instead, an interfacial electron transfer is possible in energy because LUMO+4 is much higher than CBM of MoS2. Alternatively, when the electronic transition takes place within MoS2, e.g. from state-664 to state-678 (CBM), the interfacial electron transfer is unfavorable because all unoccupied states of PTB7 are higher than CBM of MoS2 in energy. On contrary, the interfacial hole transfer from MoS2 to PTB7 is allowed thermodynamically. The right panel of Fig. 3 shows the spatial distribution of selected states involved in the interfacial electron and hole transfer dynamics in PTB7-5L@MoS2. Similarly, both HOMO and LUMO+4 are mainly located on PTB7 in spite of that the latter has some tail on MoS2, while the others are essentially on MoS2 except 9 ACS Paragon Plus Environment
The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 10 of 24
state-664 having some tail on PTB7. These tails enhance nonadiabatic couplings benefitting the interfacial hole transfer in PTB7-5L@MoS2 as will be seen in the following nonadiabatic dynamics simulations.
Figure 4. Time-dependent electron amount that is still located on the PTB7 moiety in (A) PTB7-1L@MoS2 and (B) PTB7-5L@MoS2 models; Time-dependent hole amount that is still located on the MoS2 moiety in (C) PTB7-1L@MoS2 and (D) PTB7-5L@MoS2 models.
To estimate the interfacial electron and hole transfer dynamics between PTB7 to MoS2 in PTB7-1L@MoS2 and PTB7-5L@MoS2, we have carried out nonadiabatic dynamics
10 ACS Paragon Plus Environment
Page 11 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry Letters
simulations under the classical path approximation, which has been widely used for simulating carrier dynamics in a series of organic and inorganic materials.39-42
The interfacial electron transfer in PTB7-1L@MoS2 is mainly from the electronic transition from state-385 (HOMO) to state-400 (LUMO) within PTB7 (see above). As a result, state-400 is chosen as initial electron state. Different from state-385 that is entirely on PTB7, state-400 has some tail on MoS2 (see Fig. 3). This will increase the electronic coupling between PTB7 and MoS2, which drives the interfacial electron transfer. The panel A of Fig. 4 shows the time-dependent population of the excited electron that is still localized on PTB7. Due to the tail of state-400, the initial electron population on PTB7 is about 0.5. After only about 10 fs, this amount of electron is transferred to MoS2. This ultrafast electron transfer is mainly ascribed to the fact that the initially populated electron state i.e. state-400 of PTB7 has very small energy gaps with accepting conduction band states and large density of accepting states of MoS2 (see Fig. 2). For PTB7-5L@MoS2, state-713 is chosen as the initial electron state. As with state-400 in the above PTB7-1L@MoS2, state-713 in PTB7-5L@MoS2 also has a tail on MoS2. The interfacial electron transfer is also ultrafast and completed within the first 10 fs as shown in the panel B of Fig. 4. Similarly, small energy gaps between state-713 of PTB7 with accepting conduction band states and large density of accepting states of MoS2 are responsible for this ultrafast interfacial electron transfer. Finally, it is noteworthy that adiabatic electron transfer is dominant and in charge of the interfacial electron transfer in either PTB7-1L@MoS2 or PTB7-5L@MoS2. 11 ACS Paragon Plus Environment
The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 12 of 24
The interfacial hole transfer dynamics from MoS2 to PTB7 is complicated. In PTB71L@MoS2, as discussed above, the two electronic transitions from state-382 to state-386 and from state-381 to state-387 have the largest oscillator strengths and are mainly responsible for generating hole states within MoS2. Therefore, state-382 and state-381 are respectively chosen as initial hole states in our simulations. As seen in Fig. 3, both state-381 and state-382 are completely located on MoS2 without visible tails on PTB7, which implies a small nonadiabatic coupling between MoS2 and PTB7. In addition, due to the existence of a comparable energy gap between states 381-382 of MoS2 and states 384-385 of PTB7, as shown in the panel C of Fig. 2, the interfacial hole transfer from MoS2 to PTB7 should be slow. The panel C of Fig. 4 shows the time-dependent hole population on MoS2, which is averaged over the two dynamics simulations in which state-381 and state-382 are chosen as initial hole states, respectively (see SI for individuals). It can be found that the interfacial hole transfer is indeed very slow. Through a simple exponential fitting on the time-dependent hole population in the panel C of Fig. 4, the time for the interfacial hole transfer is roughly estimated to be ca. 70 ps.
However, the situation is changed in PTB7-5L@MoS2. An ultrafast interfacial hole transfer is observed as shown in the panel D of Fig. 4. Its hole transfer time is estimated to be about 1 ps through fitting the time-dependent hole population with a simple exponential function. Interestingly, in contrast to the electron transfer, nonadiabatic hole hoppings play a major role in the interfacial hole transfer in PTB7-5L@MoS2. How to understand such ultrafast interfacial hole transfer in particular compared with that in PTB7-1L@MoS2? First, the initial 12 ACS Paragon Plus Environment
Page 13 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry Letters
hole states 664 and 665 have remarkable tails on PTB7, which thereby enhance the hole coupling between MoS2 and PTB7. This is unlike the situation of states 381 and 382 in PTB71L@MoS2 where there is no tail spreading to PTB7 (see Fig. 3). Second, the energy gap between the initial hole state of MoS2 and the energetically nearest accepting state of PTB7 is much smaller in PTB7-5L@MoS2 than that in PTB7-1L@MoS2. For example, it is 0.44 eV between VBM of MoS2 and HOMO-1 of PTB7; while, it is reduced to 0.24 eV between VBM of MoS2 and HOMO-9 of PTB7 as shown in Fig. 5.
Figure 5. Schematic mechanism for the interfacial electron and hole transfer dynamics in (left) PTB7-1L@MoS2 and (right) PTB7-5L@MoS2 models. Also shown are the estimated electron and hole transfer times and the corresponding energies of CBM and VBM and some molecular orbitals located on the PTB7 polymer. See text for discussion. 13 ACS Paragon Plus Environment
The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 14 of 24
Our above simulations not only help to rationalize experimental phenomena but also provide new mechanistic insights on the interfacial electron and hole transfer dynamics at the interface of PTB7 and MoS2. First, it is experimentally inferred that 700 nm photon mainly excites PTB7 while 530 nm photon does MoS2.27 This is also confirmed by our calculated electronic transitions in either PTB7-1L@MoS2 or PTB7-5L@MoS2 models. In addition, the interfacial electron transfer is experimentally assumed to start from delocalized electron state so that electron does not need diffuse to the interface before completely transferring to MoS2localized states.27 In our DFT calculations, these delocalized electron states have been observed, for example, state-400 in PTB7-1L@MoS2 and state-713 in PTB7-5L@MoS2 in Fig. 3.
Second, our estimated time of < 10 fs for the interfacial electron transfer in both PTB71L@MoS2 and PTB7-5L@MoS2 models is reasonably consistent with recent experiments in which the photoinduced electron transfer from PTB7 to MoS2 is found to occur in less than IRF-limited 250 fs.27 The discrepancy between experiments and simulations can be understood very well taking into account that our constructed models are far from realistic PTB7@MoS2 heterojunction in which a variety of defects and trapping states of MoS2 exist and the thickness of the PTB7 layer is at the nanometer scale, etc. Nevertheless, our simulations indeed provide important dynamical information and details about this interfacial electron transfer. In addition, it can be found that the thickness of the PTB7 polymer has negligible effects on the interfacial electron transfer dynamics because both PTB7-1L@MoS2 and PTB7-5L@MoS2 models exhibit similar dynamical feature as shown in the panels A and B of Fig. 4. 14 ACS Paragon Plus Environment
Page 15 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry Letters
Third, through comparing the results of PTB7-5L@MoS2 with those of PTB7-1L@MoS2, one can find that the thickness of the PTB7 polymer have remarkable influence on the interfacial hole transfer dynamics. The estimated time for the hole transfer in PTB7-1L@MoS2 is about 70 ps; while, in PTB7-5L@MoS2, it is decreased to ca. 1 ps, which is close to the recent experimental data, i.e. less than 3 ps.27 The different hole transfer time in two models also demonstrates that the thickness of the PTB7 polymer has remarkable influence on the interfacial hole transfer dynamics. Considering that it takes place from MoS 2 to PTB7, we can thus qualitatively understand why the quenching of the MoS2 photoluminescence varies as a function of the thickness of the PTB7 polymer in experiments.35
In summary we have employed first-principles DFT calculations and related nonadiabatic dynamics simulations to explore the photoinduced interfacial electron and hole transfer dynamics in the mixed-dimensional PTB7@MoS2 heterojunction. The interfacial electron transfer is found to be ultrafast and is in complete within about 10 fs in both PTB71L@MoS2 and PTB7-5L@MoS2 models, which is fully consistent with recent experiments.27 This electron transfer is also found to be not sensitive to the thickness of the PTB7 polymer. Different from the electron transfer, the hole transfer is sensitive to the thickness of the PTB7 polymer, the time is estimated to be about 70 ps in PTB7-1L@MoS2 but is decreased to be about 1 ps in PTB7-5L@MoS2, which is close to experimentally determined time of less than 3 ps. 27 Finally, we have found that the electron transfer is mainly controlled by the adiabatic electron evolution; whereas, in the hole transfer, nonadiabatic hoppings play a dominant role. These 15 ACS Paragon Plus Environment
The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 16 of 24
new insights could be useful for the design of excellent charge-separation interfaces of mixeddimensional TMD-based heterojunctions for optoelectronic applications.
Supporting Information. The following files are available free of charge. The theoretical backgrounds, simulation details, additional figures and tables. (PDF)
AUTHOR INFORMATION Corresponding Author Email: ganglong.cui@bnu.edu.cn (G.C.)
Notes The authors declare no competing financial interests. ACKNOWLEDGMENT This work has been supported by the grants NSFC 21522302 (G.C.) and NSFC 21520102005 (G.C. and W.F.).
References (1) Gholamvand, Z.; McAteer, D.; Harvey, A.; Backes, C.; Coleman, J. N. Electrochemical Applications of Two-Dimensional Nanosheets: The Effect of Nanosheet Length and Thickness. Chem. Mater. 2016, 28, 2641-2651. (2) Xue, Y.; Zhang, Y.; Liu, Y.; Liu, H.; Song, J.; Sophia, J.; Liu, J.; Xu, Z.; Xu, Q.; Wang, Z.; Zheng, J.; Liu, Y.; Li, S.; Bao, Q. Scalable Production of a Few-Layer MoS2/WS2 Vertical 16 ACS Paragon Plus Environment
Page 17 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry Letters
Heterojunction Array and Its Application for Photodetectors. ACS Nano 2016, 10, 573580. (3) Tedstone, A. A.; Lewis, D. J.; O’Brien, P. Synthesis, Properties, and Applications of Transition Metal-Doped Layered Transition Metal Dichalcogenides. Chem. Mater. 2016, 28, 1965-1974. (4) Withers, F.; Del Pozo-Zamudio, O.; Mishchenko, A.; Rooney, A. P.; Gholinia, A.; Watanabe, K.; Taniguchi, T.; Haigh, S. J.; Geim, A. K.; Tartakovskii, A. I.; Novoselov, K. S. Light-Emitting Diodes by Band-Structure Engineering in Van Der Waals Heterostructures. Nat. Mater. 2015, 14, 301-306. (5) Peng, B.; Ang, P. K.; Loh, K. P. Two-Dimensional Dichalcogenides for LightHarvesting Applications. Nano Today 2015, 10, 128-137. (6) Jin, C.; Rasmussen, F. A.; Thygesen, K. S. Tuning the Schottky Barrier at the Graphene/MoS2 Interface by Electron Doping: Density Functional Theory and ManyBody Calculations. J. Phys. Chem. C 2015, 119, 19928-19933. (7) Fan, X.; Chang, C. H.; Zheng, W. T.; Kuo, J.-L.; Singh, D. J. The Electronic Properties of Single-Layer and Multilayer MoS2 under High Pressure. J. Phys. Chem. C 2015, 119, 10189-10196. (8) Wang, X.; Xia, F. Stacked 2D Materials Shed Light. Nat. Mater. 2015, 14, 264-265. (9) Hong, X.; Kim, J.; Shi, S.-F.; Zhang, Y.; Jin, C.; Sun, Y.; Tongay, S.; Wu, J.; Zhang, Y.; Wang, F. Ultrafast Charge Transfer in Atomically Thin MoS2/WS2 Heterostructures.
Nat. Nanotechnol. 2014, 9, 682-686.
17 ACS Paragon Plus Environment
The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 18 of 24
(10) Lee, C.-H.; Lee, G.-H.; van der Zande, A. M.; Chen, W.; Li, Y.; Han, M.; Cui, X.; Arefe, G.; Nuckolls, C.; Heinz, T. F.; Guo, J.; Hone, J.; Kim, P. Atomically Thin p-n Junctions with Van Der Waals Heterointerfaces. Nat. Nanotechnol. 2014, 9, 676-681. (11) Fiori, G.; Bonaccorso, F.; Iannaccone, G.; Palacios, T.; Neumaier, D.; Seabaugh, A.; Banerjee, S. K.; Colombo, L. Electronics Based on Two-Dimensional Materials. Nat.
Nanotechnol. 2014, 9, 768-779. (12) Jariwala, D.; Sangwan, V. K.; Lauhon, L. J.; Marks, T. J.; Hersam, M. C. Emerging Device Applications for Semiconducting Two-Dimensional Transition Metal Dichalcogenides. ACS Nano 2014, 8, 1102-1120. (13) Pan, J.; Wang, Z.; Chen, Q.; Hu, J.; Wang, J., Band Structure Engineering of Monolayer MoS2 by Surface Ligand Functionalization for Enhanced Photoelectrochemical Hydrogen Production Activity. Nanoscale 2014, 6, 13565-13571. (14) Britnell, L.; Ribeiro, R. M.; Eckmann, A.; Jalil, R.; Belle, B. D.; Mishchenko, A.; Kim, Y.-J.; Gorbachev, R. V.; Georgiou, T.; Morozov, S. V.; Grigorenko, A. N.; Geim, A. K.; Casir-aghi, C.; Castro Neto, A. H.; Novoselov, K. S. Strong Light-Matter Interactions in Heterostructures of Atomically Thin Films. Science 2013, 340, 1311-1314. (15) Bernardi, M.; Palummo, M.; Grossman, J. C. Extraordinary Sunlight Absorption and One Nanometer Thick Photovoltaics Using Two-Dimensional Monolayer Materials.
Nano Lett. 2013, 13, 3664-3670. (16) Yu, W. J.; Liu, Y.; Zhou, H.; Yin, A.; Li, Z.; Huang, Y.; Duan, X. Highly Efficient GateTunable Photocurrent Generation in Vertical Heterostructures of Layered Materials.
Nat. Nanotechnol. 2013, 8, 952-958.
18 ACS Paragon Plus Environment
Page 19 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry Letters
(17) Wu, C.-C.; Jariwala, D.; Sangwan, V. K.; Marks, T. J.; Hersam, M. C.; Lauhon, L. J. Elucidating the Photoresponse of Ultrathin MoS2 Field-Effect Transistors by Scanning Photocurrent Microscopy. J. Phys. Chem. Lett. 2013, 4, 2508-2513. (18) Nan, H., et al., Strong Photoluminescence Enhancement of MoS2 through Defect Engineering and Oxygen Bonding. ACS Nano 2014, 8, 5738-5745. (19) Chhowalla, M.; Shin, H. S.; Eda, G.; Li, L.-J.; Loh, K. P.; Zhang, H. The Chemistry of Two-Dimensional Layered Transition Metal Dichalcogenide Nanosheets. Nat. Chem. 2013, 5, 263-275. (20) Feng, J.; Qian, X.; Huang, C.-W.; Li, J. Strain-Engineered Artificial Atom as a BroadSpectrum Solar Energy Funnel. Nat. Photonics 2012, 6, 866-872. (21) Wang, Q. H.; Kalantar-Zadeh, K.; Kis, A.; Coleman, J. N.; Strano, M. S. Electronics and Optoelectronics of Two-Dimensional Transition Metal Dichalcogenides. Nat.
Nanotechnol. 2012, 7, 699-712. (22) Eda, G.; Yamaguchi, H.; Voiry, D.; Fujita, T.; Chen, M.; Chhowalla, M. Photoluminescence from Chemically Exfoliated MoS2. Nano Lett. 2011, 11, 5111-5116. (23) Mak, K. F.; Lee, C.; Hone, J.; Shan, J.; Heinz, T. F. Atomically Thin MoS 2: A New Direct-Gap Semiconductor. Phys. Rev. Lett. 2010, 105, 136805. (24) Splendiani, A.; Sun, L.; Zhang, Y.; Li, T.; Kim, J.; Chim, C.-Y.; Galli, G.; Wang, F. Emerging Photoluminescence in Monolayer MoS2. Nano Lett. 2010, 10, 1271-1275. (25) Huang, Y.; Zhuge, F.; Hou, J.; Lv, L.; Luo, P.; Zhou, N.; Gan, L.; Zhai, T. Van Der Waals Coupled Organic Molecules with Monolayer MoS2 for Fast Response Photodetectors with Gate-Tunable Responsivity. ACS Nano 2018, 12, 4062-4073.
19 ACS Paragon Plus Environment
The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 20 of 24
(26) Huang, Y. L.; Zheng, Y. J.; Song, Z.; Chi, D.; Wee, A. T. S.; Quek, S. Y. The Organic2D Transition Metal Dichalcogenide Heterointerface. Chem. Soc. Rev. 2018, 47, 3241– 3264. (27) Zhong, C.; Sangwan, V. K.; Wang, C.; Bergeron, H.; Hersam, M. C.; Weiss, E. A. Mechanisms of Ultrafast Charge Separation in a PTB7/Monolayer MoS2 Van Der Waals Heterojunction. J. Phys. Chem. Lett. 2018, 9, 2484-2491. (28) Homan, S. B.; Sangwan, V. K.; Balla, I.; Bergeron, H.; Weiss, E. A.; Hersam, M. C. Ultrafast Exciton Dissociation and Long-Lived Charge Separation in a Photovoltaic Pentacene-MoS2 Van Der Waals Heterojunction. Nano Lett. 2017, 17, 164-169. (29) Liu, X.; Gu, J.; Ding, K.; Fan, D.; Hu, X.; Tseng, Y.-W.; Lee, Y.-H.; Menon, V.; Forrest, S. R. Photoresponse of an Organic Semiconductor/Two-Dimensional Transition Metal Dichalcogenide Heterojunction. Nano Lett. 2017, 17, 3176-3181. (30) Kafle, T. R.; Kattel, B.; Lane, S. D.; Wang, T.; Zhao, H.; Chan, W.-L. Charge Transfer Exciton and Spin Flipping at Organic-Transition-Metal Dichalcogenide Interfaces.
ACS Nano 2017, 11, 10184-10192. (31) Choi, J.; Zhang, H.; Choi, J. H. Modulating Optoelectronic Properties of Two Dimensional Transition Metal Dichalcogenide Semiconductors by Photoinduced Charge Transfer. ACS Nano 2016, 10, 1671-1680. (32) Jariwala, D.; Howell, S. L.; Chen, K.-S.; Kang, J.; Sangwan, V. K.; Filippone, S. A.; Turrisi, R.; Marks, T. J.; Lauhon, L. J.; Hersam, M. C. Hybrid, Gate-Tunable, Van Der Waals p-n Heterojunctions from Pentacene and MoS2. Nano Lett. 2016, 16, 497-503. (33) Nguyen, E. P.; Carey, B. J.; Harrison, C. J.; Atkin, P.; Berean, K. J.; Della Gaspera, E.; Ou, J. Z.; Kaner, R. B.; Kalantar-zadeh, K.; Daeneke, T. Excitation Dependent 20 ACS Paragon Plus Environment
Page 21 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry Letters
Bidirectional Electron Transfer in Phthalocyanine-Functionalised MoS2 Nanosheets.
Nanoscale 2016, 8, 16276-16283. (34) Zheng, Y. J.; Huang, Y. L.; Chen, Y.; Zhao, W.; Eda, G.; Spataru, C. D.; Zhang, W.; Chang, Y.- H.; Li, L.-J.; Chi, D.; Quek, S. Y.; Wee, A. T. S. Heterointerface Screening Effects
between
Organic
Monolayers
and
Monolayer
Transition
Metal
Dichalcogenides. ACS Nano 2016, 10, 2476–2484. (35) Shastry, T. A.; Balla, I.; Bergeron, H.; Amsterdam, S. H.; Marks, T. J.; Hersam, M. C. Mutual Photoluminescence Quenching and Photovoltaic Effect in Large-Area SingleLayer MoS2-Polymer Heterojunctions. ACS Nano 2016, 10, 10573-10579. (36) Vélez, S.; Ciudad, D.; Island, J.; Buscema, M.; Txoperena, O.; Parui, S.; Steele, G. A.; Casanova, F.; van der Zant, H. S. J.; Castellanos-Gomez, A.; Hueso, L. E. Gate-Tunable Diode and Photovoltaic Effect in an Organic-2D Layered Material p-n Junction.
Nanoscale 2015, 7, 15442-15449. (37) Furchi, M. M.; Pospischil, A.; Libisch, F.; Burgdörfer, J.; Mueller, T. Photovoltaic Effect in an Electrically Tunable Van Der Waals Heterojunction. Nano Lett. 2014, 14, 47854791. (38) Gan, L.-Y.; Zhang, Q.; Cheng, Y.; Schwingenschlögl, U. Photovoltaic Heterojunctions of Fullerenes with MoS2 and WS2 Monolayers. J. Phys. Chem. Lett. 2014, 5, 1445-1449. (39) Liu, X. Y.; Xie, X. Y.; Fang, W.; Cui, G. L. Theoretical Insights into Interfacial Electron Transfer between Zinc Phthalocyanine and Molybdenum Disulfide. J. Phys. Chem. A 2018, 122, 9587-9596. (40) Akimov, A. V.; Neukirch, A. J.; Prezhdo, O. V. Theoretical Insights into Photoinduced Charge Transfer and Catalysis at Oxide Interfaces. Chem. Rev. 2013, 113, 4496-4565. 21 ACS Paragon Plus Environment
The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 22 of 24
(41) Chu, W.; Saidi, W. A.; Zheng, Q.; Xie, Y.; Lan, Z.; Prezhdo, O. V.; Petek, H.; Zhao, J. Ultrafast Dynamics of Photongenerated Holes at a CH3OH/TiO2 Rutile Interface. J. Am.
Chem. Soc. 2016, 138, 13740-13749. (42) Long, R.; Casanova, D.; Fang, W.-H.; Prezhdo, O. V. Donor Acceptor Interaction Deter- mines the Mechanism of Photoinduced Electron Injection from Graphene Quantum Dots into TiO2: π-Stacking Supersedes Covalent Bonding. J. Am. Chem. Soc. 2017, 139, 2619-2629. (43) Niu, X.; Li, Y.; Zhang, Y.; Li, Q.; Zhou, Q.; Zhao, J.; Wang, J., Photo-Oxidative Degradation and Protection Mechanism of Black Phosphorus: Insights from Ultrafast Dynamics. J. Phys. Chem. Lett. 2018, 9, 5034-5039. (44) Nelson, T.; Fernandez-Alberti, S.; Roitberg, A. E.; Tretiak, S., Electronic Delocalization, Vibrational Dynamics, and Energy Transfer in Organic Chromophores. J. Phys. Chem.
Lett. 2017, 8, 3020-3031. (45) Kilina, S. V.; Tamukong, P. K.; Kilin, D. S., Surface Chemistry of Semiconducting Quantum Dots: Theoretical Perspectives. Acc. Chem. Res. 2016, 49, 2127-2135. (46) Ulbricht, R.; Dong, S.; Chang, I. Y.; Mariserla, B. M. K.; Dani, K. M.; Hyeon-Deuk, K.; Loh, Z.-H., Jahn-Teller-Induced Femtosecond Electronic Depolarization Dynamics of the Nitrogen-Vacancy Defect in Diamond. Nat. Commun. 2016, 7, 13510. (47) Erck, A.; Sapp, W.; Kilina, S.; Kilin, D., Photoinduced Charge Transfer at Interfaces of Carbon Nanotube and Lead Selenide Nanowire. J. Phys. Chem. C 2016, 120, 2319723206. (48) Hohenberg, P.; Kohn, W. Inhomogeneous Electron Gas. Phys. Rev. B 1964, 136, B864B871. 22 ACS Paragon Plus Environment
Page 23 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry Letters
(49) Kohn, W.; Sham, L. J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133-A1138. (50) Parr, R. G.; Yang, W. T. Density-Functional Theory of Atoms and Molecules; Oxford University Press: USA; 1994. (51) Fiolhais, K.; Nogueira, F.; Marques, M. A Primer in Density Functional Theory; Springer: Berlin; 2003. (52) Marques, M. A. L.; Ullrich, C. A.; Nogueira, F.; Rubio, A.; Burke, K.; Gross, E. K. U.
Time-dependent Density Functional Theory; Springer: Berlin; 2006. (53) Hutter, J.; Iannuzzi, M.; Schiffmann, F.; VandeVondele, J. CP2K: Atomistic Simulations of Condensed Matter Systems. WIREs Comput. Mol. Sci. 2014, 4, 15-25. (54) VandeVondele, J.; Krack, M.; Mohamed, F.; Parrinello, M.; Chassaing, T.; Hutter, J. QUICKSTEP: Fast and Accurate Density Functional Calculations Using a Mixed Gaussian and Plane Waves Approach. Comput. Phys. Commun. 2005, 167, 103-128. (55) Perdew, J. P.; Ruzsinszky, A.; Csonka, G. I.; Vydrov, O. A.; Scuseria, G. E.; Constantin, L. A.; Zhou, X.; Burke, K. Restoring the Density-Gradient Expansion for Exchange in Solids and Surfaces. Phys. Rev. Lett. 2008, 100, 136406. (56) VandeVondele, J.; Hutter, J. Gaussian Basis Sets for Accurate Calculations on Molecular Systems in Gas and Condensed Phases. J. Chem. Phys. 2007, 127, 114105. (57) Krack, M. Pseudopotentials for H to Kr Optimized for Gradient-Corrected ExchangeCorrelation Functionals. Theor. Chem. Acc. 2005, 114, 145-152. (58) Hartwigsen, C.; Goedecker, S.; Hutter, J. Relativistic Separable Dual-Space Gaussian Pseudopotentials from H to Rn. Phys. Rev. B 1998, 58, 3641-3662. 23 ACS Paragon Plus Environment
The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 24 of 24
(59) Goedecker, S.; Teter, M.; Hutter, J. Separable Dual-Space Gaussian Pseudopotentials.
Phys. Rev. B 1996, 54, 1703-1710. (60) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865-3868. (61) Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A Consistent and Accurate Ab Initio Parametrization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. J. Chem. Phys. 2010, 132, 154104. (62) Guidon, M.; Hutter, J.; VandeVondele, J. Auxiliary Density Matrix Methods for Hartree-Fock Exchange Calculations. J. Chem. Theory Comput. 2010, 6, 2348-2364. (63) Guidon, M.; Hutter, J.; VandeVondele, J. Robust Periodic Hartree-Fock Exchange for Large-Scale Simulations Using Gaussian Basis Sets. J. Chem. Theory Comput. 2009, 5, 3010-3021. (64) Guidon, M.; Schiffmann, F.; Hutter, J.; VandeVondele, J. Ab Initio Molecular Dynamics Using Hybrid Density Functionals. J. Chem. Phys. 2008, 128, 214104. (65) Heyd, J.; Scuseria, G. E.; Ernzerhof, M. Hybrid Functionals Based on a Screened Coulomb Potential. J. Chem. Phys. 2006, 124, 219906. (66) Heyd, J.; Scuseria, G. E.; Ernzerhof, M. Hybrid Functionals Based on a Screened Coulomb Potential. J. Chem. Phys. 2003, 118, 8207-8215.
24 ACS Paragon Plus Environment