2.1. DsQ Adsorbed on the GrNP Surface
Structure: As a first step, the fully aromatic hydrocarbon structure containing three zig-zag and seven armchair edges as model system for GrNP was fully optimized considering the ωB97X-D3BJ/def2-TZVPP/CPCM level of theory (For details see
Section 3). The geometry optimization leads to an irregular honeycomb lattice for GrNP, where the C–C bond length varies between 1.38–1.47 Å. Hereafter, the geometry of the GrNP was kept fixed, on which the DsQ molecule was superimposed, and the relative position of the molecule with respect to the GrNP was optimized using the same theoretical method. The geometry configuration of the DsQ adsorbed on the GrNP surface can be seen in
Figure 2. In the equilibrium geometry configuration of the binary complex, the DsQ lies almost parallel with the graphene sheet, distances between the graphene plane and the heavy atoms of the catechol fragment are between 3.369 – 3.463 Å, and the N atom of the amine group is at 3.63 Å far from the plane (see
Figure 2b). The adsorption energy, defined as the difference between the energy of the binary complex and the sum of the energies of the individual fragments, is -30.41 kcal/mol, which is more than 50% higher than the adsorption energies obtained for DA (-20.24 kcal/mol) and DoQ (-20.11 kcal/mol), respectively (See Ref. [
33]). This also means that DsQ adheres better to the graphene surface than either DA or DoQ molecules. As for the ground state charge distribution, there is a charge transfer from DsQ to the surface of 0.070
e (so-called molecular-surface charge transfer or MSCT), which is very close to the value obtained for DA (0.068
e, see Ref. [
33]). As for the difference in conformational energy between DA and DsQ adsorbed on graphene, it was found that the GrNP – DA complex is more stable than GrNP – DsQ with an energy value of 10.55 kcal/mol. However, in a very recent study, it was shown that the PCM model incorrectly prefers the neutral form to the zwitterion [
39], which may, for example, give a discrepancy between DA and DsQ results in the present case. In this respect, e.g. the CANDLE solvent model could give more accurate results [
39,
40].
UV absorption: In order to better understand the UV absorption spectrum of the binary system and the fingerprints of the different components, the first 30 electronically excited states were calculated both for the binary system and the components separately (See the black, red and blue curves in
Figure 3). Taking into account that the results for DsQ have already been reported (see results for dopamine zwitterion (or DA
zw) presented in Ref. [
34]), only a brief summary of its photochemical behavior is given. In the present case, a slightly different result was obtained, since instead of the DLPNO-STEOM-CCSD method developed at coupled-cluster theory level, the calculations were performed with TDDFT. Accordingly, DsQ presents spectral characteristics with peak maxima at 259 nm (S
1) and 224 nm (S
2). See the blue spectral curve in
Figure 3. The one-electron molecular as well as the natural difference orbitals involved in the electronic excitations of DsQ molecule can be found in the
Supplementary Material (SM) file (see Tables S1 and S2)
The geometry optimization of the first excited electron state (S1) leads to a shortening of the C–O bonds (from 1.371 Å to 1.367 Å and from 1.298 Å to 1.270 Å) and to a stretching of the carbon ring, with all six bonds being longer than 1.4 Å (the bond values are 1.427, 1.456, 1.400, 1.414, 1.400 and 1.424, respectively in Å). The fluorescence emission wavelength is 308 nm.
Graphene is a 2D solid that has interesting physics due to its unusual electronic band structure [
41]. Of course, this is accompanied by unusual optical properties characterized by a broad absorption spectral range, given by intra-band transitions (π –π *) at low photon energies (0.0-5.0 eV in the far-infrared spectral range) and inter-band transitions ((π – σ*) at higher energies (more than 5.0 eV from the mid-infrared to the ultraviolet) [
42,
43]. As described in more detail in
Section 3, the simplified 0-dimensional GrNP model is able to reproduce most of the spectral features observed in the 2D model. The UV absorption spectrum for GrNP is shown with red line in
Figure 3. In this case, two characteristic peaks are observed, one at 566 nm and the other at 270 nm. The first one appears as an effect of a GrNP finite dimensional model, while the second one is the absorption maximum typical of a 2D graphene sheet [
33,
43].
The UV absorption spectrum of the GrNP – DsQ complex is shown in black in the same
Figure 3, while the electronic excited state energies oscillator strengths, amount of charge transfers between the subsystems and the nature of the electronic states (either localized on one of the subsystems or CT state) for the first thirty electronic excited states are collected in
Table 1. For the spectrum of the binary system, it is observed that both peaks seen at 566 nm and 270 nm for the pristine GrNP system are almost identically preserved. A very small intensity increase is observed for the S
0 → S
2 electronic transition (λ(S
2) = 471 nm). As for the nature of the excited states, ten are localized on the GrNP, eleven are CT-type, while nine are mixed, i.e. containing both localized and CT-type excitation (See data in
Table 1 and the natural difference orbitals (NDO) collected in
Table S3 in the SM file). Compared to the GrNP – DA and GrNP – DoQ cases where seventeen and eighteen electron excitations are localized on GrNP, respectively, DsQ is characterized by fewer localized and more CT-type electron transitions, with no pure DsQ-like excitations. Interestingly, it was observed that excitations localized on GrNP have lower (exception S
0→ S
2), while pure CT-type excitations are characterized by higher energy (upper than S
0→ S
12) transitions. The strongest CT nature of the excited electron states for DsQ is also confirmed by the amount of charge transferred between the two subsystems. While in the case of GrNP–DA there were only three cases where the transferred charge was greater than 0.3
e and only one case greater than 0.7
e, for the present system there are eleven states having CT greater than 0.3
e, three cases greater than 0.7
e, and for the S
2 state even one unitary charge is transferred from the DsQ to the graphene nanoparticle (called as molecule-to-surface charge transfer or MSCT effects). To understand the reason for this significant charge transfer, it is necessary to see what the nature of the highest occupied and lowest unoccupied electron orbitals are. Accordingly, based on the fragment orbital contribution analysis, the molecular orbital scheme of the joined GrNP – DsQ system was built. Its graphics can be seen in
Figure 4. By comparison with the GrNP – DA case, it can be observed that since the HOMO energy of the individual DsQ molecule is higher than that of the DA and thus approaches the HOMO energy of GrNP, not only the HOMO - 1 orbital of the GrNP – DsQ complex will contain a DsQ contribution (more than 80% DsQ), but also the HOMO (18% DsQ). And, of course, not only the HOMO and HOMO - 1 orbitals contain smaller or larger DsQ contributions, but also the other lower lying occupied orbitals. As for the unoccupied orbitals, the first eleven are almost entirely localized on GrNP, but also the subsequent three orbitals contain only a small proportion of DsQ contributions (
Figure 4).
Accordingly, in the molecular orbital excitation scheme the lowest electronic excitation partially localized on DsQ (S
0 → S
25) involves the LUMO + 15 unoccupied orbitals (See
Table S3 in SM file). Overall, it can be concluded that the nature of the electronic transitions observed in GrNP – DsQ is very similar to that seen in GrNP – DA, with the addition that in the first case, charge-transfer electronic excitations appear to be more pronounced.
Excited state relaxation: From the analyses carried out so far, it can be concluded that excited states can be localised on one or the other component of the molecule-surface complex, or charges can easily migrate from the molecule to the surface, or vice versa [
33]. It can therefore be assumed that relaxation processes for different electronic excited states can occur either on the components separately or even inside the binary system. It has already been shown that non-radiative relaxation can occur either on the graphene surface itself [
44] or in the binary complex [
45,
46,
47,
48,
49,
50], while excited states can also decay via fluorescence phenomena in the case of slightly modified GrNP [
51,
52]. Furthermore, it was also shown that for the S
1 excited state, there was no significant geometric change for DA adsorbed on GrNP, whereas for DoQ, significant changes appeared for both C=O and aromatic C–C bonds due to excitation [
33]. To see exactly what effect the excitation has on the DsQ molecule adsorbed on GrNP, the geometry of the molecule was optimized in the first excited state. Similar approaches as for DA and DoQ have been considered in the present case. Namely,
i) non-radiative phenomena were neglected due to the very complex theoretical framework and huge computation resources;
ii) positions of the carbon atoms forming the surface have been frozen. When comparing the ground and first excited state geometries, it can be observed that there is a slightly more pronounced change in the molecular bonds than was seen for DA, but not nearly as significant as was obtained for DoQ. Accordingly, the C–O
- and C–OH bonds stretch from 1.283 Å and 1.359 Å to 1.302 Å and 1.372 Å, respectively, while the aromatic C–C bond lengths oscillate from 1.411-1.383-1.395-1.398-1.374-1.429 Å to 1.405-1.397-1.388-1.404-1.377-1.426 Å. As regarding the relative position of the DsQ, the plane defined by the aromatic ring fragment comes closer to the GrNP surface and appears to be tilted, with the six carbon atoms between 3.08–3.27 Å away from the plane of the surface, the O atoms at 3.00 Å and 3.04 Å, respectively, as well as the N atom of the amine group at 3.36 Å. At the same time, the wavelength of the S
0 → S
1 excitation changes from 570 nm to 627 nm. In conclusion, it can be stated that no significant geometric change occurs in the DsQ molecule during the relaxation of the S
1 excited state, unlike in DoQ (See Ref. [
33]).
2.2. Charge Transfer and Distance from Surface
The effects of changing the separation distance on the charge transfer interactions between luminescent QD and proximal dopamine (in QD–dopamine assemblies) has been already demonstrated by steady-state and time-resolved fluorescence measurements considering different lengths of poly(ethylene glycol) (PEG) as bridge [
53]. In light of this, the question arises as to what extent the perpendicular shift between the plane of the DA, DoQ and DsQ molecules and the GrNP surface affects the charge transfer phenomenon in these assemblies. Accordingly, the molecules have been moved along the z-coordinate axis by different values (Δ
z = -0.3, -0.2, -0.1, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.8 and 1.0 Å) relative to the equilibrium distance (
z0) from the surface. Separately for each of these geometries, the Löwdin electron population was calculated for the first thirty electronic excited states, based on the electron density of the given state. Amounts of charge transfer calculated for different plane distances relative to the equilibrium geometry for the DA, DoQ and DsQ adsorbed on the GrNP surface are presented in
Figure 5. Molecular orbital energy schemes (in eV) built based on the fragment orbital contribution analysis of the individual, GrNP, DA, DoQ and DsQ components and of the mixed GrNP – DA, GrNP – DoQ and GrNP – DsQ binary complexes computed for Δ
z = -0.3, 0.0 and +1.0 relative stacking distances values are shown in
Figure 6.
DA: It was already shown that, the nature of the charge transfer between the DA and the GrNP substrate is MSCT (or metal-to-surface charge transfer) [
33]. For the present case, at the equilibrium geometry configuration (
z0 + Δz, where Δ
z = 0.0) it can be observed three distinct sets of electronic excited states with significant CT character (See
Figure 5a). The first set is formed by S
6, S
7 and S
9, the second group is formed by S
15, S
16 and S
17, while the third set is defined by S
25 – S
29 states. When the DA plane approaches the surface (Δ
z = -0.3, -0.2 and -0.1 Å) it is observed that the number of CT-like electronically excited states increases, but the charge transfer amount of the existing CT states decreases and that of the newly created CT states increases. For Δ
z = -0.3, only S
1, S
3 and S
5 showed a charge transfer lower than 0.1
e, but for S
6, for example, the 0.4
e CT obtained at Δ
z = 0.0 Å for -0.3 Å decreases to 0.18
e. Furthermore, a shift of the CT peaks also can be observed as the distance between the planes decreases. For example, the peak of S
7 at Δ
z = 0.0 Å becomes S
6 at -0.3 Å, S
16 becomes S
15, and S
26 mainly bifurcates into S
25 and S
27. As the distance between the planes increases, the opposite effect is observed. That is, the number of CT states decreases, but the amount of charge transfer increases. The first set of excited states (S
6, S
7 and S
9) collapses as S
14 for Δ
z = +1.0 Å, the second group as S
26, while the third group moves out of the S
1–S
30 interval proposed for the present analysis. However, for S
14 and S
26, the amount of charge transfer is 0.87
e and 0.95
e, respectively. Fragment orbital contribution analysis reveals that the amount of this charge transfer as a function of Δ
z depends on the degree of overlap of the individual molecular orbitals. Indeed, for Δ
z = -0.3, the HOMO and HOMO – 1 orbitals contain larger DA-type contribution, which means stronger CT effects between the DA and the GrNP substrate, while in the case of Δ
z = +1.0 there are hardly any mixed molecular orbitals and therefore the number of CT-type excited electron states is significantly reduced. In the latter case, the S
12 electron excitation is characterized by HOMO –2 → LUMO, while S
26 is characterized by HOMO – 8 → LUMO and HOMO – 9 –> LUMO one-electron transitions, respectively (See
Figure 6a).
DoQ: For the CT states of the complex formed by the DoQ molecule in adsorption on GrNP surface the charge is mostly transferred from the surface to the molecule (called as SMCT or surface-to-molecule charge transfer) [
33], in contrast to the previous DA case. At the equilibrium geometry position, for Δ
z = 0.0 Å, also three distinct sets of electronic excited states with significant CT character can be observed (See
Figure 5b). The first set is made up of the S
3 state alone, the second is formed by the S
16, S
18 and S
21 states, and the third is the S
30 state. In the decreasing distance direction between the planes (Δ
z = -0.3, -0.2 and -0.1 Å), no substantial change in the charge transfer phenomenon is observed. Rather, the transferred charge is slightly dispersed in other states (see S
9 – S
15), but the position of states with high CT values in the excitation spectrum does not change. A similar trend can be observed for the increasing distance direction between the planes. The S
3 excited state still behaves as a CT state even at Δ
z = +1.0 Å, the set of S
16 – S
21 are slightly rearranged, here states S
20 and S
23 will show the strongest charge transfer character, while state S
30 moves out from the interval proposed for the present analysis. Basically, it can be said that for DoQ – GrNP, the CT phenomenon does not show any significant change in the distance interval of -0.3 Å and +0.1 Å. This can also be explained by the fact that, based on the fragment orbital contribution analysis, LUMO + 1 is essentially a DoQ-type molecular orbital and by varying the distance (Δ
z), the orbital nature almost remain the same (See
Figure 6b). For ex. the CT-type S
3 state is defined by the HOMO → LUMO + 1 one-electron transitions.
DsQ: It has already been shown in the previous sub
Section 2.1 that the photochemical behavior of DsQ is similar to that of DA since the nature of the CT effects is the same MSCT-type transition as that observed for DA. The difference is that the number of pure CT-type transitions or their mixing with localized GrNP-type excitation within the first thirty excited states is significantly higher for DsQ than for DA, and on average the magnitude of charge transfer is also higher (Compare
Figure 5a,c). At the equilibrium geometry position of the DsQ, i.e. with Δ
z = 0.0 Å, significant CT transitions are seen for the S
2 state and for all other excited states starting from S
13. For a molecular configuration approaching to the surface, e.g. Δ
z = -0.3 Å, the intensity of charge transfer increases and the CT nature is already observed at S
11. When the molecule starts to move away from the surface, the opposite happens, i.e. intensities decrease and, for example, at Δ
z = +0.4 Å the S
13 – S
15 states hardly contain any CT nature. At Δ
z = +1.0 Å only S
2, S
16 and S
30, and partially the groups of S
20 – S
27, show CT character. Similar to DA, the CT transitions for DsQ are strongly dependent on the distance between the molecule and the surface, but the number of CT states is still larger than that observed for DA. Compared to DA, the change is that for Δ
z = +1.0 Å, the DsQ-type orbital is no longer HOMO – 2, but HOMO – 1 with lower orbital energy (See
Figure 6c).