Azobenzene Adsorption on the MoS2 (0001) Surface: A Density

6 days ago - The interest in light sensitive organic molecules, such as azobenzene, has increase due to their ability to functionalize two-dimensional...
0 downloads 0 Views 5MB Size
Subscriber access provided by Kaohsiung Medical University

C: Surfaces, Interfaces, Porous Materials, and Catalysis

Azobenzene Adsorption on the MoS2 (0001) Surface: A Density Functional Investigation within van der Waals Corrections Luís Cabral, Fernando Pereira Sabino, Matheus Paes Lima, Gilmar Eugenio Marques, Victor Lopez-Richard, and Juarez L. F. Da Silva J. Phys. Chem. C, Just Accepted Manuscript • Publication Date (Web): 24 Jul 2018 Downloaded from http://pubs.acs.org on July 25, 2018

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 10 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

Azobenzene Adsorption on the MoS2 (0001) Surface: A Density Functional Investigation within van der Waals Corrections L. Cabral,† Fernando P. Sabino,‡ Matheus P. Lima,† Gilmar E. Marques,† Victor Lopez-Richard,† and Juarez L. F. Da Silva∗,¶ †Department of Physics, Federal University of S˜ao Carlos, 13565 − 905, S˜ao Carlos, SP, Brazil ‡S˜ao Carlos Institute of Physics, University of S˜ao Paulo, 13560 − 970, S˜ao Carlos, SP, Brazil ¶S˜ao Carlos Institute of Chemistry, University of S˜ao Paulo, PO Box 780, 13560 − 970, S˜ao Carlos, SP, Brazil E-mail: juarez [email protected]

Phone: +55 16 3373 9930

Abstract The interest in light sensitive organic molecules, such as azobenzene, has increased due to their ability to functionalize two-dimensional layered systems and to engineer their electronic structure. In this work, we explore the azobenzene trans- and cis-isomers adsorption on molybdenum disulfide MoS2 (0001) layer employing the density functional theory (DFT) within van der Waals corrections to the semilocal exchange-correlation functional. We found that the aromatic rings of the azobenzene lay parallel to the surface (2 in the trans-isomer and 1 in the cis-isomer), which contributes to increase the configuration stability by van der Waals interactions. Furthermore, we found a relatively large work function change upon adsorption due to the electron density rearrangement, and hence, it might affect the electronic transport properties within the single MoS2 (0001) layer. We observed an increasing in the relative DFT+vdW total energy among the azobenzene isomers (trans- and cis-) from 0.51 eV in gas-phase (trans-isomer has the lowest energy) to 0.81 eV for azobenzene supported on the MoS2 (0001) surface, which can be explained by the contact of the two rings of the trans-isomer directly to the surface. Thus, the binding of azobenzene on the single MoS2 (0001) layer affects the isomerization process due to the relative energy increase, and this, in turn, influences the transport properties of the single MoS2 (0001) layer due to the changes in the electrostatic potential.

I

Introduction

The emergence of ultra thin semiconductor systems, such as transition-metal dichalcogenides (TMDs), have paved the way for electronic two-dimensional (2D) architectures of new device concepts, and sensors. 1,2 The 2D-TMD materials have been intensively investigated due to the wide range of possible applications, 3,4 which can be extended by the combination of different 2D-TMDs in van der Waals hetero-structures. 5,6 For

example, the control of the number of layers linked by van der Waals interactions 7 and the incidence of electromagnetic fields allows to manipulate their electronic properties. 8 The 2D-TMDs composed by MQ2 layers, where M is a transition metal atom, such as Mo, and Q is a chalcogen atom (S, Se, Te), have attracted great interest. 9 In particular, molybdenum disulfide, MoS2 , 10 has been used in several applications, such as photovoltaics, 11 lubricants, 12 and etc. These systems are flexible, optically active, and their properties can be tuned through composition variation, or stacking, but also by the character of their substrates 13 or the surrounding chemical environment. 14 In that respect, progress has been attained by combining diverse atomic layers and organic interfaces that allows foreseen the functionalisation of 2D devices and their integration at nanoscales. A wide variety of these organicinorganic hybrids has been studied and the reported results have evidenced not just how the film morphologies can be altered, 15 but also their quantum efficiency, 16 and their transport and thermal properties. 17 The integration of organic molecules and two-dimensional layers of TMDs, such as MoS2 , has made these structures attractive for new hybrid solar cells proposals 18 or phototransistor architectures. 19 In the case of photosensitive switches, the use of photochromic azobenzene molecules as a bulding block in these hybrid architectures is very promising. 19 The photoisomerization of azobenzene has already been successfully used for the modulation of biological processes. 20,21 The reversibility of the cis-trans isomerization in the azobenzene groups has also allowed promoting strategies for wavelength-selective control of a wide variety of chemical, 22 mechanical 23 , and electrical properties 24 of their host systems. Azobenzene is an organic molecule formed by the functional group azo, i.e., – N – N – , 25 and two aromatic benzene rings, as represented in Figure 1. Two isomers have been reported so far for azobenzene, namely, the trans-azobenzene with planar structure, and a second isomer, known as cis-, where the two benzene rings are not in the

ACS Paragon Plus Environment 1

The Journal of Physical Chemistry 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

same plan. 21 The transition between the isomers was reported by Hartley, 26 noting that the process cis-trans can occur in the dark, while the inverse isomerization, i.e., trans-cis, is only feasible under the ultraviolet light incidence. 27 As azobenzene directly photoisomerizes in gas-phase, 26,28 this transformation can be quenched when the molecule is placed on metal surfaces, 29,30 which boosted investigations exploring the reversible azobenzene photo-switch. 31 Previous ab initio investigations characterized the potential energy curves of azobenzene isomers and their intermediate configurations. 32 Comstock et. al. studied the influence of the Au(111) surface in the azobenzene photo-switching. 33 The optical properties of the MoS2 mono- and multi-layers exfoliated on Au surfaces functionalized with the azobenzene molecule was also reported. 34 In turn, the magnetic effect of graphene doped by Co in the trans- cis-azobenzene reversibly was investigated by Nurbawono and Zhang, 35 while the influence of graphene in the optical excitation of azobenzene was analyzed by Fu et. al.. 36 In this work, we report the characterization of the adsorption properties of azobenzene conformations on the molybdenum disulfide MoS2 (0001) surface. Our calculations are based on the density functional theory (DFT) within the van der Waals (vdW) correction to improve the description of the weak molecule-surface interactions. We found that the trans-isomer in gas-phase and the conformations located horizontally on the substrate are the lowest energy configurations. Adsorption energies and work functions are also analyzed. In addition, we were able to determine the mechanism of the charge flow between the molecule and the surface within the effective Bader charge concept.

II Theoretical Approach Computational Details

and

A Total Energy Calculations Our total energy calculations are based on DFT 37,38 within the semilocal Perdew–Burke–Ernzerhof 39 (PBE) formulation for the exchange-correlation energy functional. To improve the description of the weak vdW azobenzene-MoS2 (0001) interactions, we employed the vdW D3 correction proposed by S. Grimme, namely, the DFT-D3 framework. 40 For comparison, we performed also DFT-PBE calculations without the vdW correction. In the DFT-D3 framework, as implemented in the Vienna ab initio simulation package DFT+vdW , (VASP), 41,42 the vdW DFT corrected total energy, Etot is obtained by the sum of the vdW energy D3 correction, D3 , and the DFT-PBE total energy, E DFT-PBE , namely, Edisp tot DFT+vdW DFT-PBE D3 Etot = Etot + Edisp ,

(1)

where ij

D3 Edisp =−

ij M M C C  1 XXX fd,6 66 + fd,8 88 . 2 i=1 j=1 L Ri j,L R i j,L

(2)

Page 2 of 10

The sums are over all atoms M and all translations of the unit cell L = (l1 , l2 , l3 ), while the dispersion coefficients, ij ij namely, C6 and C8 , play a crucial role in the magnitude of the vdW interactions. The damping functions, fd,6 and fd,8 , are employed to avoid singularities for small inter-atomic distances, fd,n =

sn ij

1 + 6(Ri j /( sR,n R0 ))−αn

,

(3)

where the sR,n parameters are the scaling factors of the cutoff ij ij ij radii, and R0 = (C8 /C6 )1/2 . The αn parameters are fixed and depends on the selected exchange-correlation functional. Further details on the DFT+vdW framework can be found elsewhere. 40–43 The Kohn–Sham (KS) equations were solved using the projected augmented-wave (PAW) method 44,45 as implemented in VASP, 46,47 version 5.4.1. We employed the PAW projectors provided within VASP, namely, Mo (5s1 4d5 ), S (3s2 3p4 ), C (2s2 2p2 ), N (2s2 2p3 ) and H (1s1 ), where the valence electrons are indicated in parentheses. The equilibrium volume of the bulk MoS2 in the hexagonal P63 /mmc crystal structure was obtained by minimization of the stress tensor and atomic forces using a plane-wave cutoff energy of 689 eV, while for the total energy calculations for all systems (geometric optimizations, density of states, etc), we employed a cutoff energy of 473 eV, which is 12.5 % higher than the largest recommended PAW cutoff energy for the mentioned chemical species. The differences among the cutoff energies is due to the slow (faster) convergence of the stress tensor (total energy) with the number of plane-waves. For the Brillouin zone integration, we employed a kmesh of 8 × 8 × 2 for the bulk MoS2 in the primitive unit cell, while a k-mesh of 8 × 8 × 1 was used for a single MoS2 (0001) monolayer using a 1 × 1 surface unit cell and a vacuum thickness of about 15 Å. For the density of states, we increased the k-point mesh twice for both systems. For the azobenzene/MoS2 (0001) calculations, we employed the repeated slab geometry using a hexagonal 8 × 8 unit cell, one single MoS2 (0001) layer, and 25 Å for the vacuum region due to the large size of the azobenzene molecule, e.g., the trans-azobenzene has a length of 11.26 Å. For convergence tests, we employed also a 4 × 4 and 6 × 6 surface unit cells, and three MoS2 (0001) layers for the slab. For those calculations, we employed only the Γ-point for the Brillouin zone integration due to the larger size of the surface unit cell. For azobenzene in gas-phase, we employed an orthorhombic 3 box with 24 × 25 × 26 Å , which yields a minimum distance between the azobenzene molecules of about 15 Å, and hence, the Γ-point was also used due to the lack of dispersion in the electronic states. We employed a Gaussian smearing of 0.010 eV for all calculations, except for azobenzene in gasphase, where we used 0.0010 eV, which is required to avoid fractional occupation of the electronic states. For the azobenzene/MoS2 (0001) calculations, we considered the adsorption on only one side of the layer, and hence, dipole corrections were employed for all calculations. 48,49 All structures reached the equilibrium geometry once the atomic forces on every atom are smaller

ACS Paragon Plus Environment 2

Page 3 of 10 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

Figure 1 Optimized atomic configurations for the trans- and cis-azobenzene in gas-phase (left side). The highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) are indicated for both configurations (right side). −1

than 0.010 eV Å . For the analysis of the electron density, we employed the Bader charge analysis, where the grid-points were increased by a factor of three for the fast Fourier transform (FFT), which plays a crucial role to obtain well converged effective Bader charges on each atom.

B Azobenzene/MoS2 (0001) Azobenzene is composed by two benzene rings linked by the – Cb – N – N – Cb – chain, where the Cb atoms belong to the benzene rings and the double – N – N – bonds play a crucial role in the azobenzene configurations. 25,50 Azobenzene adopts two configurations, namely, a ground state planar structure, in which the two benzene rings and the – Cb – N – N – Cb – chain are in the same plane (trans-azobenzene). However, under ultraviolet light, 28 mechanical stress, or even electrostatic stimulation, 51 the trans- configuration changes to the cis-azobenzene, where the two benzene rings change their orientation with respect to each other due to the change of the – Cb – N – N – Cb – torsion angle. 50 The reverse isomerization, cis- to trans-, can occur spontaneously under dark conditions. 21 With the aim to understand the role of the azobenene-MoS2 (0001) interaction in the relative stability of both azobenzene isomers, which can play a crucial role in the electronic switch properties, we considered both configurations for our adsorption studies and several adsorption sites, which are shown in the Supporting Information.

III

Results

A Azobenzene in Gas-Phase According to our calculations, we found that the trans-isomer is 0.51 eV lower in energy than the cis-isomer, which is in excellent agreement with experimental results obtained under dark conditions, i.e., 0.52 eV. 28 The optimized gas-phase

structures for the trans- and cis-azobenzene are depicted in Figure 1, where the structural differences can be easily seen, in particular, the dihedral angle for – Cb – N – N – Cb – are 180◦ (trans-) and 10.8◦ (cis-), which are in good agreement with experimental results. 28 From the KS-orbital analysis, we found that the highest occupied molecular orbital (HOMO) is mainly localized on the – Cb – N – N – Cb – fragment in the trans- configuration, while it extends to the benzene rings in the cis-azobenzene, and hence, it indicates a clear difference among both isomers. In contrast, no substancial differences was found in the lowest unoccupied molecular orbital (LUMO) between both molecules. To improve our understanding of the dihedral angle, we Bader , calculated the effective Bader charge, QB e f f = Zval − Q where Zval is the number of valence electrons of the respective chemical species. As expected, the nitrogen atoms have an effective negative charge, which can be explained by the higher electronegativity of N (3.04) compared with C (2.55) atoms. Furthermore, the magnitude of the effective charge on the N atoms is larger for the trans-isomer, i.e., −0.28 e (trans) and −0.23 e (cis-). This is due to the enhanced Coulomb repulsion between the N atoms, and hence, we expect a larger bond length between the N – N atoms in the trans-azobenzene, which is in fact confirmed by our results, e.g., N – N = 1.26 Å (trans-) and 1.25 Å (cis-). Although the difference is small, it correlates with the nature of the interactions in both isomers.

B Bulk and Layered MoS2 Systems The equilibrium lattice constants, a0 and c0 , for the bulk MoS2 in the hexagonal P63 /mmc structure are 3.16 Å and 12.35 Å, respectively, which are consistent with the experimental results, a0 = 3.16 Å and c0 = 12.29 Å, 52 i.e., corresponding to deviations smaller than 0.49 %. For comparison, DFT-PBE yields c0 = 13.45 Å, i.e., a relative error of 9.39 %, and hence, the vdW D3 correction plays a crucial role to reduce the error from 9.39 % to 0.49 %. From the Bader charge analysis, we obtained an effective cationic

ACS Paragon Plus Environment 3

The Journal of Physical Chemistry 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 4 of 10

Figure 2 Band structure and local density of states for MoS2 in the bulk phase (left panel) and for the single MoS2 (0001) layer (right panel) calculated with the DFT-PBE+D3 framework. The Brillouin zone is depicted on the top of band structure with the k-path showed in green lines.

charge on the Mo atoms of 0.74 e, and hence, an anionic charge of −0.37 e on the S atoms, which changes to −0.36 e in the single MoS2 (0001) layer using the same a0 parameter as in the bulk structure, i.e., a tiny change of 0.01 e due to the vdW interactions between the MoS2 layers. The band structure and density of states (DOS) for the MoS2 bulk and single MoS(0001) layer are shown in 2, where the valence band maximum (VBM) is located at the Γ-point for MoS2 bulk, while the conduction band minimum (CBM) is located between the Γ- and K-points, i.e., an indirect band gap of 1.07 eV, which is 17.05 % (1.29 eV) smaller than the experimental result, 53 which is a well known problem in DFTPBE, which cannot be corrected by the addition of the vdW correction as the D3 framework does not affect the electron density directly, but only indirectly due to the changes in the geometric parameters. From the local DOS, the valence band is composed mainly by the Mo and S states and the VBM, in particular, is mostly derived from Mo states. Compared with MoS2 bulk, the single layer DOS is shifted up in energy, which is related to the confinement of the electronic states within the layer.

C Azobenzene Adsorption on the MoS2 (0001) We analyzed the cis- and trans-azobenzene configurations adsorbed on several positions on the MoS2 (0001) surface, namely, perpendicular and planar orientations in different positions on the surface. All calculated configurations are shown in the Supporting Information, while the lowest energy configurations for both isomers on the MoS2 (0001) layer (top and side views) are depicted in Figure 3. Below, we will discuss the structural, energetics, and electronic properties for the lowest energy configurations, while the same parameters are summarized in the Supporting Information for the high energy configurations.

Structural Configurations: In the lowest energy configurations shown in Figure 3, the benzene rings are parallel to the MoS2 (0001) surface, which is expected to enhance the vdW and Coulomb interactions. In particular, the planar trans- structure contains two benzene rings horizontally placed on the surface, while only one benzene ring is parallel to the surface in the cis-isomer. As mentioned above, the binding energy among the N, C, and H atoms is given by the strong covalent and ionic chemical interactions, and hence, expressively larger than the van der Waals interaction with the MoS2 (0001) surface. Thus, we do not expect considerable changes in the conformation of the azobenzene isomers, which are supported by our DFT-D3 calculations. Both molecules (trans- and cis-isomers) and monolayer MoS2 (0001) structures are composed by hexagonal rings, however, there is no alignment of these hexagonal rings upon the adsorption, i.e., both rings are horizontally displaced from one another. We found a distance of 3.26 Å (3.23 Å) between the hexagonal rings of the cis-isomer (trans-) and the MoS2 (0001) surface, which increases to 3.63 Å (3.79 Å) using DFT-PBE. In addition, as the conformation of the molecule on the monolayer changes, the charge redistribution in the geometry also varies, as depicted at the bottom of Figure 3. However, the effective Bader charge of the lowest energy configurations, depicted at the upper part of Figure 3, is very similar to the cis- and trans- isomers in gas-phase, represented in the left part of Figure 1, indicating the weak magnitude of the vdW forces. Adsorption energy: The adsorption energy, Ead , which measures the magnitude of the azobenzene-substrate interaction, was calculated using the following equation, AZO/MoS2 (0001)

Ead = Etot

ACS Paragon Plus Environment 4

MoS2 (0001)

AZO − Etot − Etot

,

(4)

Page 5 of 10 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

Figure 3 Charge density analysis. (a) Top and side views of the lowest energy configurations for the cis- and trans-azobenzene adsorbed on the MoS2 (0001) surface. The effective Bader charges are indicated in the top view, while the vertical distance between the surface and the closest atom is depicted in the side view. (b) The electron density differences induced by the azobenzene adsorption are indicated using a three-dimensional representation. (c) Average electron density difference along the direction perpendicular to the surface. The position of the cis- (red) and the trans- (blue) isomers (closest atom to the surface), as well as the Mo-layer atoms (green), are indicated by the dashed lines. AZO/MoS (0001)

2 where Etot is the total energy of the AZO and azobenzene/MoS2 (0001) configuration, while Etot

MoS (0001)

Etot 2 are the total energies of the azobenzene and single MoS2 (0001) monolayer, respectively. For both azobenzene conformations, we employed the total energy AZO . Thus, the of the respective isomers to calculate Etot adsorption energy focuses only on the interaction energy, neglecting the energy differences due to the geometry variation of the isomers. We found Ead = −0.98 eV (−0.12 eV using PBE) and −0.68 eV (−0.13 eV using PBE) for the trans- and cisisomers, respectively. Thus, trans- binds 0.30 eV stronger than cis-azobenzene on the MoS2 (0001) surface using the DFT-D3, which is expected due to the larger contact of the trans-isomer with the MoS2 (0001) surface, i.e., the vdW interaction and Coulomb forces are stronger in the trans- than in cis-configuration on the layer. The azobenzene rings are in direct contact with the surface and the interaction increases for the two parallel benzene rings of the trans- in comparison with the one horizontal ring of the cis-geometry. Therefore, the vdW D3 correction plays a crucial role in the magnitude of the adsorption energy and in the lowest energy configuration. The perpendicular orientation of the trans-isomer, where

the molecule interacts with a single H atom, increases the DFT-D3 total energy by 0.78 eV, indicating that the parallel orientation of the benzene rings is preferable. On the MoS2 (0001) surface, the trans-isomer is 0.81 eV lower in energy than the cis-azobenzene, whereas this difference is 0.51 eV in the gas-phase. Thus, the interaction with the surface enhances the energy difference between the transand cis-isomers, however, we should keep in mind that those relative total energies are not the energy barriers for the transcis switch. The isomerization switching under light irradiation depends on the azobenzene electronic states, which can be affected by the interaction with the surface. As mentioned above, we found an energy difference between the trans- and cis-isomers supported on the MoS2 (0001) surface of 0.81 eV, which is near to the adsorption energy calculated above, i.e., −0.98 eV, however, it does not indicate light-induced molecular desorption along the trans-cis switch, which is supported by experimental observations, e.g., azobenzene on Au(111) surface. 33 Work Function: The work function, Φ, is the energy required to remove one electron from the system, given

ACS Paragon Plus Environment 5

The Journal of Physical Chemistry by Φ = Ves − EF , where Ves and EF are the electrostatic potential at infinite (middle of the vacuum region) and the Fermi energy, respectively. For the clean MoS2 (0001) surface, we obtained a work function of 5.88 eV, which is reduced upon the adsorption of the azobenzene molecules. In the lowest energy configurations, we obtained work functions of 5.41 eV (∆Φ = 0.47 eV) and 5.24 eV (∆Φ = 0.64 eV) for the trans- and cis-isomers on the surface, respectively. Those results were obtained using the 8 × 8 surface unit cell, while smaller surface unit cells, e.g., 4 × 4 and 6 × 6, yields work functions than differ by less than 0.05 eV, i.e., orders of magnitude smaller than the work function changes. Thus, there is a clear difference in the magnitude of the work function changes due to the different orientation of the trans- and cis-isomers, and hence, the obtained work function reduction can be explained as follows: (i) there is an electron density rearrangement, and as a result, there is an effective electron density flow towards to the MoS2 (0001) surface, which is supported by our analysis indicated in Figure 3 for both isomers; (ii) the difference between the trans- and cisisomers work functions of 0.17 eV is due to the orientation of the cis-isomer, in which one of the benzene rings is nearly perpendicular to the surface. Below, we will discuss our charge transfer analysis, which can helps to clarify these trends. Charge Transfer Analyses: To obtain a better understanding of the substrate work function change upon the adsorption of azobenzene, we performed a charge density analyses using the Bader charge concept, where we Bader . calculated the effective Bader charge, QB e f f = Zval − Q QBader

Zval is the number of valence electrons and is the Bader charge, respectively. Furthermore, we calculated the electron density difference, ∆ρ = ρmol/sub − ρ sub − ρmol , and the averaged xy-density along the z-direction. Theses analyses were performed for the lowest energy configurations and the results are shown in Figure 3. The effective Bader charges on the chemical species in the trans- and cis-isomers are nearly the same as in the azobenzene in gas-phase, which indicates that the effective charges are not affected by the weak van der Waals interactions. Although the average charges are nearly the same for both isomers, there are differences in the local charges, which can be seen from the electron density difference analyses, in particular, the three-dimensional plot and the averaged xy-density along the z-direction, shown in Figure 3, which indicates larger local effects derived from the trans-isomer. It can be explained by the closer approach to the surface, which enhances the interactions. Thus, the larger work function change for the cis-isomer can be explained by the orientation of the second benzene ring, which is nearly perpendicular to the surface and not by the electron density changes. Density of States: To improve our understanding of the electronic levels of the azobenzene molecule on the MoS2 (0001) surface, we calculated the local density of states (LDOS) for the trans- and cis-isomers on the single

2.00 1.50 1.00 0.50 0.00 1.50

Mo

S

1.00 Local Density of States (State/eV)

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 6 of 10

0.50 0.00 1.50 H

1.00

trans cis

0.50 0.00 1.50 C

1.00 0.50 0.00 1.50

N

1.00 0.50 0.00 -5

-4

-3

-2

-1 0 1 Energy (eV)

2

3

4

5

Figure 4 Local density of states of the trans- and cis-azobenzene isomers on the top of a single MoS2 (0001) layer. The electronic states of the semiconductor are smoothly displaced from the Fermi energy, where the azobenzene levels enter inside the gap energy of the MoS2 (0001).

MoS2 (0001) layer, depicted in Figure 4. In the LDOS of the trans- and cis-configurations on the MoS2 (0001) layer, the HOMO states of the single MoS2 (0001) layer is displaced from the Fermi energy by 0.5 eV, which is mostly composed by Mo d- and S p-states, while the LUMO is composed by Mo d- and S p-states. The states of the azobenzene are displaced from each other, where the highest contribution is given by the trans-isomer on the surface. Thus, as the states of the semiconductor are smoothly displaced from the Fermi energy, the electronic levels of the azobenzene enters inside the gap region, changing the electronic properties of the carriers near the Fermi energy. This is an interesting fact to study in the modulation of the electronic transport and its switching under light incidence. As depicted in Figure 4, the calculated band gap of the single MoS2 (0001) layer is 1.42 eV. A ruling factor in the manipulation of the band gap energy is the effect of the confinement profile through different number of stacking layers. Thus, increasing the number of layers, the electronic properties tend to the MoS2 bulk phase. It is important to note that, although the lowest energy configuration shows the higher work function and an adsorption site of preference, the electronic properties of the MoS2 (0001) are smoothly disturbed by the interaction with the azobenzene isomers. The molecule on the MoS2 (0001) does not practically change the electronic properties of the semiconductor.

ACS Paragon Plus Environment 6

Page 7 of 10 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

IV

Summary

In this work, we performed a DFT-PBE+D3 investigation of the adsorption properties of the azobenzene isomers on the MoS2 (0001) surface. Furthermore, for comparison, DFTPBE calculations were also performed. The aromatic rings of azobenzene lay parallel to the surface, e.g., two rings for transand only one for the cis-isomer. The relative total energy among the trans- and cis-isomers in gas-phase, 0.51 eV, increases to 0.81 eV upon the adsorption on the MoS2 (0001) surface. Thus, the adsorption on the 2D substrate might affect the conversion trans-cis process. The electronic properties of the semiconductor layer are locally disturbed by the presence of the azobenzene molecules, which introduce changes in the substrate work function and also in the behavior of the confined carriers of the semiconductor. Therefore, we expect that the prospects for introducing photo-responsive units in electronics and optoelectronics make this study useful.

(5)

(6)

(7)

(8)

(9)

V

Supporting Information

The following additional data are shown in the Supporting Information: all calculated configurations for azobenzene supported on the MoS2 (0001) surface; geometric and energetic parameters; convergence tests using the different surface unit cells.

VI

(10)

(11)

(12)

Acknowledgements

The authors thanks the Brazilian Agencies CAPES (Coordenac¸a˜ o de Aperfeic¸oamento de Pessoal de N´ıvel Superior), CNPq (Conselho Nacional de Desenvolvimento Cient´ıfico e Tecnol´ogico) and FAPESP (Fundac¸a˜ o de Amparo a` Pesquisa do Estado de S˜ao Paulo), grants 2014/19142 − 2, 2014/02112 − 3, for the financial support.

(13)

(14)

References (1) Chhowalla Manish,; Shin Hyeon Suk,; Eda Goki,; Li Lain-Jong,; Loh Kian Ping,; Zhang Hua, The Chemistry of Two–Dimensional Layered Transition Metal Dichalcogenide Nanosheets. Nature Chem. 2013, 5, 263–275. (2) Wang Qing Hua,; Kalantar-Zadeh Kourosh,; Kis Andras,; Coleman Jonathan N.,; Strano Michael S., Electronics and Optoelectronics of Two-dimensional Transition Metal Dichalcogenides. Nat. Nanotechnol. 2012, 7, 699–712. (3) 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. (4) Choi, W.; Choudhary, N.; Han, G. H.; Park, J.; Akinwande, D.; Lee, Y. H. Recent Development of Two–Dimensional Transition Metal Dichalcogenides

(15)

(16)

(17)

and Their Applications. Mater. Today 2017, 20, 116– 130. Guan, Z.; Lian, C.-S.; Hu, S.; Ni, S.; Li, J.; Duan, W. Tunable Structural, Electronic, and Optical Properties of Layered Two-Dimensional C2 N and MoS2 van der Waals Heterostructure as Photovoltaic Material. J. Phys. Chem. C 2017, 121, 3654–3660. Sangwan, V. K.; Beck, M. E.; Henning, A.; Luo, J.; Bergeron, H.; Kang, J.; Balla, I.; Inbar, H.; Lauhon, L. J.; Hersam, M. C. Self–Aligned van der Waals Heterojunction Diodes and Transistors. Nano Lett. 2018, 18, 1421–1427. Jariwala Deep,; Marks Tobin J.,; Hersam Mark C., Mixed-dimensional van der Waals Heterostructures. Nat. Mater. 2016, 16, 170–181. Sun Zhipei,; Martinez Amos,; Wang Feng, Optical Modulators with 2D Layered Materials. Nature Photon. 2016, 10, 227–238. Leng, K.; Chen, Z.; Zhao, X.; Tang, W.; Tian, B.; Nai, C. T.; Zhou, W.; Loh, K. P. Phase Restructuring in Transition Metal Dichalcogenides for Highly Stable Energy Storage. ACS Nano 2016, 10, 9208–9215. Song, I.; Park, C.; Choi, H. C. Synthesis and Properties of Molybdenum Disulphide: From Bulk to Atomic Layers. RSC Adv. 2015, 5, 7495–7514. Fortin, E.; Sears, W. Photovoltaic Effect and Optical Absorption in MoS2 . J. Phys. Chem. Solids 1982, 43, 881–884. Kim, Y.; Huang, J.; Lieber, C. M. Characterization of Nanometer Scale Wear and Oxidation of Transition Metal Dichalcogenide Lubricants by Atomic Force Microscopy. Appl. Phys. Lett. 1991, 59, 3404–3406. Bessonov Alexander A.,; Kirikova Marina N.,; Petukhov Dmitrii I.,; Allen Mark,; Ryhnen Tapani,; Bailey Marc J. A., Layered Memristive and Memcapacitive Switches for Printable Electronics. Nat. Mater. 2014, 14, 199–204. Lee, G.-H.; Yu, Y.-J.; Cui, X.; Petrone, N.; Lee, C.H.; Choi, M. S.; Lee, D.-Y.; Lee, C.; Yoo, W. J.; Watanabe, K.; Taniguchi, T.; Nuckolls, C.; Kim, P.; Hone, J. Flexible and Transparent MoS2 Field–Effect Transistors on Hexagonal Boron Nitride–Graphene Heterostructures. ACS Nano 2013, 7, 7931–7936. Breuer, T.; Maßmeyer, T.; M¨anz, A.; Zoerb, S.; Harbrecht, B.; Witte, G. Structure of van der Waals Bound Hybrids of Organic Semiconductors and Transition Metal Dichalcogenides: the Case of Acene films on MoS2 . Phys. Status Solidi RRL 2016, 10, 905– 910. 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. Wan Chunlei,; Gu Xiaokun,; Dang Feng,; Itoh Tomohiro,; Wang Yifeng,; Sasaki Hitoshi,; Kondo Mami,; Koga Kenji,; Yabuki Kazuhisa,; Snyder G. Jeffrey,; Yang Ronggui,; Koumoto Kunihito, Flexible n– type Thermoelectric Materials by Organic Intercalation

ACS Paragon Plus Environment 7

The Journal of Physical Chemistry 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

(18)

(19)

(20)

(21)

(22)

(23)

(24)

(25) (26) (27)

(28)

(29)

(30)

(31)

(32)

of Layered Transition Metal Dichalcogenide TiS2 . Nat. Mater. 2015, 14, 622–627. 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. Margapoti, E.; Li, J.; Ceylan, O.; Seifert, M.; Nisic, F.; Anh, T. L.; Meggendorfer, F.; Dragonetti, C.; Palma, C.A.; Barth, J. V.; Finley, J. J. A 2D Semiconductor–SelfAssembled Monolayer Photoswitchable Diode. Adv. Mater. 2015, 27, 1426–1431. Kienzler, M. A.; Reiner, A.; Trautman, E.; Yoo, S.; Trauner, D.; Isacoff, E. Y. A Red–Shifted, Fast– Relaxing Azobenzene Photoswitch for Visible Light Control of an Ionotropic Glutamate Receptor. J. Am. Chem. Soc. 2013, 135, 17683–17686. Beharry, A. A.; Woolley, G. A. Azobenzene Photoswitches for Biomolecules. Chem. Soc. Rev. 2011, 40, 4422–4437. Lerch Michael M.,; Hansen Mickel J.,; Velema Willem A.,; Szymanski Wiktor,; Feringa Ben L., Orthogonal Photoswitching in a Multifunctional Molecular System. Nat. Commun. 2016, 7, 12054. Zhou Hongwei,; Xue Changguo,; Weis Philipp,; Suzuki Yasuhito,; Huang Shilin,; Koynov Kaloian,; Auernhammer Gnter K.,; Berger Rdiger,; Butt HansJrgen,; Wu Si, Photoswitching of Glass Transition Temperatures of Azobenzene–Containing Polymers Induces Reversible Solid-to-Liquid Transitions. Nature Chem. 2016, 9, 145–151. Lee, K. E.; Lee, J. U.; Seong, D. G.; Um, M.-K.; Lee, W. Highly Sensitive Ultraviolet Light Sensor Based on Photoactive Organic Gate Dielectrics with an Azobenzene Derivative. J. Phys. Chem. C 2016, 120, 23172–23179. Griffiths, J. II. Photochemistry of Azobenzene and Its Derivatives. Chem. Soc. Rev. 1972, 1, 481–493. Hartley, G. S. The Cis-form of Azobenzene. Nature 1937, 140, 281. Wei-Guang Diau, E. A New Trans–to–Cis Photoisomerization Mechanism of Azobenzene on the S 1 (n, π∗) Surface. J. Phys. Chem. A 2004, 108, 950–956. Merino, E.; Ribagorda, M. Control Over Molecular Motion Using the Cis Trans Photoisomerization of the Azo Group. Beilstein J. Org. Chem. 2012, 8, 1071– 1090. Zhou, X.-L.; Zhu, X.-Y.; White, J. Photochemistry at Adsorbate/Metal Interfaces. Surf. Sci. Rep. 1991, 13, 73–220. Qiu, X. H.; Nazin, G. V.; Ho, W. Vibrationally Resolved Fluorescence Excited with Submolecular Precision. Science 2003, 299, 542–546. Comstock, M. J.; Cho, J.; Kirakosian, A.; Crommie, M. F. Manipulation of Azobenzene Molecules on Au(111) Using Scanning Tunneling Microscopy. Phys. Rev. B 2005, 72, 153414. Monti, S.; Orlandi, G.; Palmieri, P. Features of the

(33)

(34)

(35)

(36)

(37) (38)

(39)

(40)

(41)

(42)

(43)

(44) (45)

(46)

(47)

Page 8 of 10

Photochemically Active State Surfaces of Azobenzene. Chem. Phys. 1982, 71, 87–99. Comstock, M. J.; Levy, N.; Kirakosian, A.; Cho, J.; Lauterwasser, F.; Harvey, J. H.; Strubbe, D. A.; Fr´echet, J. M. J.; Trauner, D.; Louie, S. G.; Crommie, M. F. Reversible Photomechanical Switching of Individual Engineered Molecules at a Metallic Surface. Phys. Rev. Lett. 2007, 99, 038301. Li, J.; Wierzbowski, J.; Ceylan, O.; Klein, J.; Nisic, F.; Anh, T. L.; Meggendorfer, F.; Palma, C.A.; Dragonetti, C.; Barth, J. V.; Finley, J. J.; Margapoti, E. Tuning the Optical Emission of MoS2 Nanosheets Using Proximal Photoswitchable Azobenzene Molecules. Appl. Phys. Lett. 2014, 105, 241116. Nurbawono, A.; Zhang, C. Reversible Magnetism Switching in Graphene–Based Systems via the Decoration of Photochromic Molecules. Appl. Phys. Lett. 2013, 103, 203110. Fu, Q.; Cocchi, C.; Nabok, D.; Gulans, A.; Draxl, C. Graphene–Modulated Photo–Absorption in Adsorbed Azobenzene Monolayers. Phys. Chem. Chem. Phys. 2017, 19, 6196–6205. Hohenberg, P.; Kohn, W. Inhomogeneous Electron Gas. Phys. Rev. 1964, 136, B864–B871. Kohn, W.; Sham, L. J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133–A1138. Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. 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. ´ Buˇcko, T.; Hafner, J.; Leb`egue, S.; Angy´ an, J. G. Improved Description of the Structure of Molecular and Layered Crystals: Ab Initio DFT Calculations with van der Waals Corrections. J. Phys. Chem. A 2010, 114, 11814–11824. Moellmann, J.; Grimme, S. Importance of London Dispersion Effects for the Packing of Molecular Crystals: A Case Study for Intramolecular Stacking in a Bis-Thiophene Derivative. Phys. Chem. Chem. Phys. 2010, 12, 8500–8504. Reckien, W.; Janetzko, F.; Peintinger, M. F.; Bredow, T. Implementation of Empirical Dispersion Corrections to Density Functional Theory for Periodic Systems. J. Comput. Chem. 2012, 33, 2023–2031. Bl¨ochl, P. E. Projector Augmented–Wave Method. Phys. Rev. B 1994, 50, 17953–17979. Kresse, G.; Joubert, D. From Ultrasoft Pseudopotentials to the Projector Augmented-Wave Method. Phys. Rev. B 1999, 59, 1758–1775. Kresse, G.; Hafner, J. Ab Initio Molecular Dynamics for Open–Shell Transition Metals. Phys. Rev. B 1993, 48, 13115–13118. Kresse, G.; Furthm¨uller, J. Efficient Iterative Schemes

ACS Paragon Plus Environment 8

Page 9 of 10 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

(48)

(49)

(50)

(51)

(52)

(53)

The Journal of Physical Chemistry for Ab Initio Total–Energy Calculations Using a PlaneWave Basis Set. Phys. Rev. B 1996, 54, 11169–11186. Neugebauer, J.; Scheffler, M. Adsorbate-Substrate and Adsorbate-Adsorbate Interactions of Na and K Adlayers on Al(111). Phys. Rev. B 1992, 46, 16067–16080. Makov, G.; Payne, M. C. Periodic Boundary Conditions in Ab Initio Calculations. Phys. Rev. B 1995, 51, 4014– 4022. Bandara, H. M. D.; Burdette, S. C. Photoisomerization in Different Classes of Azobenzene. Chem. Soc. Rev. 2012, 41, 1809–1825. Doki´c, J.; Gothe, M.; Wirth, J.; Peters, M. V.; Schwarz, J.; Hecht, S.; Saalfrank, P. Quantum Chemical Investigation of Thermal Cis-to-Trans Isomerization of Azobenzene Derivatives: Substituent Effects, Solvent Effects, and Comparison to Experimental Data. J. Phys. Chem. 2009, 113, 6763–6773. Young, P. A. Lattice Parameter Measurements on Molybdenum Disulphide. J. Phys. D: Appl. Phys. 1968, 1, 936–938. Ellis, J. K.; Lucero, M. J.; Scuseria, G. E. The Indirect to Direct Band Gap Transition in Multilayered MoS2 as Predicted by Screened Hybrid Density Functional Theory. Appl. Phys. Lett. 2011, 99, 261908.

ACS Paragon Plus Environment 9

The Journal of Physical Chemistry 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

Graphical TOC Entry

ACS Paragon Plus Environment 10

Page 10 of 10