Fabian
Weber
ab,
Jian
Ren
ac,
Tristan
Petit
a and
Annika
Bande
*a
aInstitute of Methods for Material Development, Helmholtz-Zentrum Berlin für Materialien und Energie GmbH, Albert-Einstein-Str. 15, 12489 Berlin, Germany. E-mail: annika.bande@helmholtz-berlin.de; Tel: +49 30 8062-42026
bInstitute of Chemistry and Biochemistry, Freie Universität Berlin, Takustr. 3, 14195 Berlin, Germany
cDepartment of Physics, Freie Universität Berlin, Arnimallee 14, 14195 Berlin, Germany
First published on 6th March 2019
In this work we provide a proof of principle for a theoretical methodology to identify functionalisation patterns in oxidised carbon 2D nanomaterials. The methodology is based on calculating a large number of X-ray absorption spectra of individually excited carbon atoms in different chemical environments using density functional theory. Since each resulting spectrum gives a fingerprint of the local electronic structure surrounding the excited atom, we may relate each spectrum to the functionalisation pattern of that excited atom up to a desired neighbourhood radius. These functionalisation pattern-specific spectra are collected in a database, that allows fast composition of X-ray absorption spectra for arbitrary structures in density functional theory quality. Finally, we present an exemplary application of the database approach to estimate the relative amount of functional groups in two different experimental samples of carbon nanomaterials.
The mission of finding possible local functionalisation patterns is ongoing research in various fields of physical, analytical and theoretical chemistry in which many different structural models have been proposed over the course of time.15–19 Some of the milestones that show the complementary nature of this structural analysis can be summed up as follows: with its first reported syntheses, elemental analyses were conducted that confirmed the non-stoichiometric character of graphite oxide.8 Following the advent of X-ray powder diffraction, the comparison with powder diffraction images of graphite lead to the conclusion that the carbon network of six-membered rings was mostly preserved after oxidation.20–22 Due to the scarce diffraction characteristics, however, no systematic ordering of the oxygen atoms could be detected. These initial measurements already gave rise to an overall structural understanding as ‘stacked two-dimensional macromolecular sheets’ that were assumed to be physically separable into monolayers.22 First infrared spectroscopy measurements supported the existence of carboxyl, aliphatic hydroxyl and epoxy groups.23 With the development of electron microprobe diffraction techniques the monolayer stacking was proven.24 Using this technique, it was also validated that no long-range order for the oxygen functionalisation exists and that the oxygen atoms were most likely bound as hydroxyl or epoxy groups on the sheet surface.24 The basis of the most widely accepted structural model so far was established by Lerf and Klinowski after 13C solid state nuclear magnetic resonance (SSNMR) measurements that are capable of locally probing the atomistic surrounding of nuclei of a specific element.11,19,25 Their findings validated the existence of carbon–carbon double bonds and lead to the conclusion that GO contained graphitic besides oxidized regions proportional to the degree of oxidation. Further improvements to their structural model have been made since, by SSNMR of 13C enriched samples12 and high-resolution transmission electron microscopy images.26
Another field of spectroscopy that provides contributions in terms of element-specific and local structural analysis is X-ray absorption (XA) spectroscopy, where highly energetic X-ray photons are used to excite core electrons into unoccupied orbitals or the ionisation continuum. Since the core electrons are affected by the practically unscreened Coulomb attraction of the spatially close nucleus, each chemical element has specific excitation energies (i.e. absorption edges) at which to probe core-valence transitions. Depending on the excitation energy, different effects may be observed that give rise to three distinguishable energetic regions, namely the pre-edge, X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine structure (EXAFS).
From a theoretical perspective, two main effects contribute to X-ray absorption. The first contribution is connected to the ionisation of the excited core electron and only appears beyond the absorption edge energy. Here, the photoexcited electron is described by continuum wavefunctions that are scattered at several neighbouring atoms. For low ejection energies, (i.e. in the XANES region) several scattering events may occur due to the low mean free path of the ejected electron. Hence, signals in the XANES region contain information on multiple scatterers at once and can be used to identify the nature of functional groups near the excited atom, as well as their relative positioning.27,28 In the EXAFS region, the mean free path of the electron is long and mostly single-scattering events are obtained. Therefore, depending on the periodicity (so called mid-range and long-range order), convoluted oscillatory structures are obtained where the frequency of convolution components can be related to the distance of periodically arranged scattering sites.29 Both effects are covered in the field of multi-scattering theory and have been implemented in programs such as FEFF,30 FDMNES,31 or WIEN2k32 to simulate XANES and EXAFS spectra.
The second contribution is caused by excitation of the core electron into bound and resonance states. While resonance states are difficult to describe, bound final states of excitonic (i.e. extended long-range)33 or molecular (i.e. local short-range) nature can be readily simulated using ab initio quantum chemistry methods such as time-dependent density functional theory (TDDFT) with core-excitation specialized exchange–correlation functionals.34,35
To extract structural information from experimental spectra, there already exist computational procedures based on multi-scattering theory that fit the properties and sites of scatterers to match experimental data in the EXAFS and XANES region. Here, implementations such as IFEFFIT36 use efficient truncations in the scattering expansion length (reciprocal space) or introduce a cutoff radius (real space) in which to optimize the properties of the scattering sites.
However, to the best of our knowledge no similar fitting procedure has been introduced to bound-to-bound state molecular calculations yet, since the non-periodic molecular description of larger systems yields high numbers of bound states, that become prohibitively expensive. To fill this gap, we introduce in this work a methodology that may open up the possibility for such an algorithm. Firstly, we calculate a large number of core-excitations of carbon atoms in different model molecules and then group the resulting spectra with respect to the bonding environment of the excited atom up to a specific radius. These resulting functionalisation specific groups are then used to construct a database, from which we may try to fit the bonding environment of explicit calculations and experimental samples of oxidised 2D carbon nanomaterials.
This paper is structured as follows. First we explain how the database of functionalisation specific spectra is generated in terms of calculation techniques, model structures and how we distinguished the different functional groups. Then, we give a general proof of principle by critically reviewing the choice of a cutoff radius and comparing the resulting spectra of fully-atomistic calculations with compositions from the database. Finally, we show the results of a complementary analysis of two experimental samples of oxidised 2D carbon nanomaterials.
Besides their pristine forms, the three scaffolds were then functionalised with oxygen-containing groups and hydrogen in several ways to obtain different environments for the individual carbon atoms. Whenever not functionalised, the dangling bonds of the edge-positioned carbon atoms were saturated with one hydrogen atom. By saturation of the graphitic scaffold, several sp3 hybridized carbon atoms are obtained that then represent the properties of carbon dots. This way, a total of sixteen model molecules was prepared on the C scaffold, while twelve belong to the CC group. All in all, this gives rise to 1049 carbon atoms from which we construct the database of individual-atom XA spectra. Additionally, one pristine structure was prepared on the CCC scaffold using a smaller basis set for the XA calculations for at least comparing the effects of different model sizes.
Each model molecule was structurally optimised using the density functional theory (DFT) implementation of the ORCA program suite.37 For this step we use the non-empirical meta-GGA TPSS functional,38 as well as Ahlrichs' def2-SVP basis set.39 Dispersion correction is achieved through Grimme's third order atom-pairwise dispersion correction with Becke-Jones (D3BJ) damping.40,41 To speed up the calculation of the exchange term, the RIJCOSX approximation is used42 with the appropriate def2-SVP/J auxiliary basis sets.43 For each structure a numerical frequency analysis was conducted to confirm optimisation of the structure to its minimum energy. All structures and their total energies may be found as attached .xyz files in the ESI.†
(1) |
(2) |
(3) |
In our work, the orbitals for the LR-TDDFT calculation are obtained by performing a single-point calculation with the BH0.57LYP functional35 on the pre-optimised model structures. This hybrid functional with a high portion of Hartree exchange has been shown to give a very good onset energy for carbon K-edge XA spectra.35 Afterwards, the 1s core-space of the resulting ground-state orbitals is localised using the Pipek–Mezey44 criterion. Note that the procedure was applied regardless of whether the core orbitals were actually degenerate or not, to ensure that the resulting XA spectra can be clearly assigned to different individual atoms.
The individual-atom XA spectra are then obtained from standard LR-TDDFT calculations exciting one specific core-electron at a time into the whole virtual space. For each LR-TDDFT calculation, a total of N = 150 roots is calculated. To account for possible relativistic effects, the zeroth-order regular approximation (ZORA) is used.47 Having obtained both the excited states |n〉 and discrete transition energies ω0n from the LR-TDDFT calculation, we construct XA spectra by applying a Gauss-broadening to the sum of the orientationally averaged discrete oscillator strengths (eqn (1)–(3)) with a full width at half maximum of 1.0 eV.
The underestimation of core excitation energies within TDDFT calculations due to approximating the exchange–correlation functional was treated by adding a global shift of 0.25 eV to all calculated XA spectra. The shift was obtained from performing a number of benchmark calculations on smaller molecules like ethene, ethanol and formic acid, and is in the same order of magnitude as shifts from other DFT functionals that were optimized for X-ray calculations and tested on the same molecules.35
To find a balance between accuracy and manageability, the indexing method presented in this work generates group descriptors that consist of the following two components. The first component is based on one of the positions basal, edge and bridge. This determines whether the indexed carbon atom is in the inner region (circles in Fig. 1(a)), at a bridge position (diamonds in Fig. 1(a)) or at the edge (triangles in Fig. 1(a)). To distinguish these positions, the indexing scheme has to take two layers of neighbours into consideration – i.e. the directly bound atoms and their respective bound atoms. We shall refer to this as the next-nearest neighbour (NNN) radius in the following.
The second component of the group descriptor concerns the functionalisation of the indexed carbon atom. To determine this, we chose to use a reduced NNN radius. This means, we nominally look at the nearest and next-nearest neighbours of the indexed atom to be able to distinguish hydroxyl from epoxy groups and such, but then store only the functionalisation of the actual indexed carbon atom, rather than the full surrounding. In other words, we end up with a more simplified group descriptor like “basal C–OH” instead of “basal C–OH bound to three aromatic basal carbon atoms”. The combination of both descriptors then gives rise to a finite number of possible functionalisation patterns, i.e. NNN distinguishable groups.
Using this scheme, we thus may distinguish aromatic atoms (C), hydrogenated atoms (C–H) and hydroxyl atoms (phenolic C–OH and non-phenolic CH–OH). Furthermore, atoms within an epoxy-group (C-Epo and CH-Epo) as well as carbonylic atoms (CHO), ketones (C-Keto) and carboxylic atoms (COOH). Additionally, we make a distinction between the atoms connecting to the carbonylic or carboxylic carbon atoms (C–CHO or C–COOH, respectively) from other aromatic carbon atoms.
A list of how often the different groups are found within the C and CC model structures can be taken from Table 1. Since all functionalised sites were embedded in structurally optimised model structures, atom types that are needed to form the carbon scaffold around functional groups are naturally more often needed to complete a meaningful molecular model structure. After assigning one of the NNN distinguishable groups to each individual-atom XA spectrum, we then construct the catalogue of so-called mean group XA spectra. These are obtained by summing up all individual spectra belonging to one group and divide the result by the number of occurrences. This way, effects caused by electronic transitions beyond the NNN radius are averaged out to some degree (see Section A of the ESI† on conjugation effects), since the individual spectra are in a different embedding environment each time. Nonetheless, it shall be noted that the quality naturally is higher, when more entries are included in the mean spectrum. The mean group XA spectra then reflect the average XA spectrum when exciting a carbon atom in the respective NNN distinguishable environment.
Group | C | CC | Total |
---|---|---|---|
Note, that numbers carrying an asterisk indicate that an additional hydrogen atom is bound to the carbon as well. Also, carboxyl and carbonyl groups were only attached at edge positions, since addition onto the basal plane was not considered sensible for oxidation of a pristine graphene sheet. | |||
CH | 2/181*/3 | 14/192*/— | 16/373*/3 |
C | 76/—/84 | 191/—/125 | 267/—/209 |
C–OH | 4/1/2 | 23/9*/3 | 27/10*/5 |
C-Epo | 14/3*/7 | 60/2*/16 | 74/5*/23 |
C–CHO | —/2/— | —/5/— | —/7/— |
CHO | —/2/— | —/5/— | —/7/— |
C–COOH | —/2/— | —/8/— | —/10/— |
COOH | —/2/— | —/8/— | —/10/— |
C-Keto | —/3/— | —/—/— | —/3/— |
Sum | 388 | 661 | 1049 |
The composition of total XA spectra is then performed by addition of mean group spectra weighted by the number of carbon atoms that belong to the respective groups in a structure. This way, we may predict the theoretical total XA spectrum for molecular structures within the underlying TDDFT quality in very short times. The other way around, this methodology may open up the possibility to identify spectral features on oxidised carbon materials without prior knowledge of their exact structure.
Take note that the central parameter of this methodology is the choice of indexing radius, because by increasing the indexing radius towards including all atoms of an example molecule, the composition of that molecule from the database will eventually converge to the exact theoretical XA spectrum. The reason for this can be understood when considering an example linear molecule ABCD, where each letter shall represent an atom of a distinguishable environment. The exact theoretical spectrum F[ABCD] is obtained when the core-excitation of every atom is calculated individually and all contributions f are summed up.
F[ABCD] = f[ABCD] + f[ABCD] + f[ABCD] + f[ABCD] | (4) |
[ABCD] = [ABC] + f[ABCD] + f[ABCD] + [BCD] | (5) |
To facilitate the nomenclature of transitions, we shall denote transitions as σ*(X–Y) or π*(XY), where the atom labels in parentheses refer to the atoms involved and the greek letter refers to the type of transition under assumption of local symmetry. Note though, that we do not claim a pure π- or σ-character at any time, but rather give interpretations based on what the virtual Kohn–Sham orbitals (KSO) of a given transition look like. Also the bond order in the labelling is only reflecting the local Lewis structure. Hence, transitions of π*(X–Y) type may occur. Additionally, if there were more than one transition of the same type possible, we use numbers in front of the transition, i.e. 2π*(XY) for a second π* transition.
Firstly, it can be seen that all three positions show a π*(CC) transition feature around 285.1 eV which is in accordance with experimental results.49–53 While this is the only noteworthy peak for the pre-edge region of the basal position (uppermost panel of Fig. 1(b)), the other two positions start to develop a signal around 288.4 eV. This peak, which is caused by transitions into σ*(C–H) orbitals becomes more pronounced for the edge position and shifts to 289.3 eV, which is in line with the experimentally observed value of 289.4 eV in small hydrocarbons.49
Since the three positions indeed show a different behaviour, the reduced NNN radius correctly distinguishes these different types of atoms. Also, as there are no major differences among each group for the different scaffold sizes, the C and CC scaffolds are assumed to be large enough to describe the three different positions.
Fig. 2 (a) Structure of a singly epoxy-functionalised CC model molecule and (b) theoretical near-edge carbon K-edge XA spectra of the highlighted atoms. |
The model molecule chosen for this comparison is shown in Fig. 2(a), with the respective carbon K-edge XA spectra in Fig. 2(b). The individual-atom spectra can be connected to different atoms in the structure, corresponding to their color and symbol. Firstly, we shall discuss the circle-marked atoms that would be assigned to the basal C-Epo group by the reduced NNN indexing. Their XA spectra show a pronounced signal just below 287.0 eV (similar to epoxy groups in bulk thin-film graphene oxide at 286.5 eV),50 which can be attributed to π*(C–O) transitions. Both the blue and red spectra are identical, because the atoms are identical with respect to their chemical surrounding as well as molecular symmetry.
The next pair of atoms (highlighted with diamonds) would be assigned as basal C atoms and these atoms are direct neighbours of the basal C-Epo atoms. Their spectra show a signal at 284.7 eV which can be attributed to π*(CC) transitions, which is in qualitative agreement with the experimental value of 285.4 eV in graphite.54 The shift in excitation energy with respect to the experimental reference is likely caused by the proximity to the oxygen atom that is only two bond-lengths away.
Evidence that these direct neighbours are shifted due to the oxygen atom becomes even more apparent, when comparing to a pair of carbon atoms that are three and four bonds away from the oxygen atom, respectively (marked with triangles). In the individual-atom XA spectra, both atoms show a π*-like transition around 285.4 eV, which is in perfect agreement with the experimental value for graphite. However, the two signals are shifted by approximately 0.1 eV with respect to each other, which is due to the fact that one of the atoms (red triangle in Fig. 2(a)) is farther away from the oxygen atom than the other (blue triangle in Fig. 2(a)). Since both of the atoms are otherwise in an identical surrounding with respect to their own next-nearest neighbors, the two carbon atoms seem to have reached a distance from the epoxy-oxygen that is large enough to not affect the individual-atom XA spectrum substantially. Conclusively, in this model molecule the locally excited carbon K-edge XA spectra are sensitive to their surrounding for up to 2–3 Å, i.e. about two bond-lengths.
It has to be noted though, that the triangle-marked atoms would also be assigned as basal C atoms, although their respective π*(CC) transitions lie about 0.7 eV apart from the one of the diamond-marked basal C atoms. While the reduced NNN indexing radius would not be able to distinguish between their environments, already the full NNN indexing radius would indeed be large enough to distinguish these two pairs from each other. With the currently available database size we, however, needed to work with the reduced NNN radius, since otherwise not all possible functionalisation environments would have been covered. Nonetheless, we can use the reduced NNN radius as a first approximation for probing the general usefulness of the overall methodology, since increasing the radius will only improve the quality.
In the first example a comparison to a lowly-functionalised model molecule (see Fig. 3(b)) is made. The solid black line of Fig. 3(a) gives the respective total XA spectrum, when adding up every individual-atom XA response in its exact environment (cf.eqn (4)). The total XA spectrum is dominated by mainly three signals. The first can be identified as a strong π*(CC) transition at 285.0 eV due to a largely conjugated π-electron system that is only slightly disrupted by the central functionalisations. After that, a signal corresponding to σ*(C–H) transitions can be found at 289.1 eV. Finally, for even higher energies a σ*(C–C) signal is obtained at 293.3 eV.
The database composition (ω) for the model molecule in panel (b) at transition energies ω is calculated by addition of the mean group XA spectra i(ω) times the number of occurrences of a specific type of atom. The composition is represented by the cumulative addition of colored planes, where each color represents one of the groups from the different mean group XA spectra. The formula below shows how the resulting composition of Fig. 3(a) has been obtained (from bottom to top), while the colored dots in molecule Fig. 3(b) mark the respective species.
(ω) = basalCH(ω) + 18edgeCH(ω) + 12bridgeC(ω) + 20basalC(ω) + 2basalC-Epo(ω) + basalC–OH(ω) | (6) |
When comparing the black line to the composition, it becomes apparent that the intensities do not match in several regions, while the position of features is well reproduced. The reason for this behaviour is that roughly above 294.0 eV, the intensities begin to smear out for specific signals in some of the mean group XA spectra due to more and more energetic shifts among the multitude of individual-atom spectra. This is probably the case for both signals at 294.8 eV and 298.5 eV, which are not well reproduced by the database composition. While the former one is contained with low intensity in the edge CH group, the latter signal seems to be shifted strongly to higher energies within the database. This is evidence that there are transitions which need a higher indexing radius than the reduced NNN radius for an exact description. However, since these transitions lie beyond the ionisation edge-jump, our approximation using the reduced NNN radius works well enough for the pre-edge region.
Next, we perform the same analysis for a highly-oxidised model molecule (see Fig. 3(d)). The theoretical total XA spectrum (see Fig. 3(c)) is dominated by four transitions. The first transition at 285.2 eV can be identified as a π*(CC) signal that shows a weak shoulder to lower energies caused by π*(CO) transitions of carboxyl and carbonyl groups. The second signal at 287.6 eV is one of σ*(C–O) character and indicates epoxy groups. At 290.5 eV one finds the third peak, which is caused by σ*(C–H) transitions. The last of the dominant signals in the explicit theoretical spectrum can be found at 292.3 eV and is given by σ*(C–C) transitions.
When the spectrum is composed from the database, one can confirm three of these four signals very well in both energetic positioning as well as intensity, with exception of the σ*(C–H) transition. One can furthermore assign more subtle shoulders to specific groups, such as the one at 288.8 eV that is caused mainly from the basal epoxy groups' σ*(C–O) transition as well as the σ*(C–H) transition of the edge CH group. Above around 293.0 eV the intensities again begin to smear out like in the case of the model molecule displayed in Fig. 3(b). Two more signals at 295.9 eV and 299.4 eV can be reproduced qualitatively in the cumulative sum of all group XA spectra. The assignment to one specific type of transition is, however, complicated by the fact that none of the groups is responsible for the signal on its own, but rather all of the groups contribute to the signals to different degrees.
As can be seen, however, both composition XA spectra also yield information beyond qualitative peak assignments, since they can provide the information on how much each functional group relatively contributes to a signal. This is particularly relevant when assigning the near-edge XA spectra of structurally unknown compounds.
Fig. 4 Comparison of the near-edge features of an XA spectrum of graphene oxide (black) with the cumulative sum of four different group XA components. Each component consists of several group XA spectra (see ESI† for detailed composition). |
Fig. 5 Comparison of a XA spectrum of an oxidised CD sample with the cumulative sum of four different group XA components. Each component consists of several group XA spectra (see ESI† for detailed composition). The spectrum was taken from Ren et al.48 |
The compositions provided here are obtained as hand-selected sums of different group XA spectra, which qualitatively reproduce the signals in the pre-edge and low-energy XANES region of the experimental spectra. The selection process was guided along the lines of additional spectroscopic findings,48 like the presence or absence of characteristic vibrations in the infrared spectrum of the respective sample. Hence, the composition presents the result of complementary methods that may allow for estimation of the relative occurrence of functionalisation patterns in the sample. Note that intensities in the region after the ionisation edge-jump around ≈290 eV in the experimental XA spectra have only qualitative meaning. Also, for easier reading we have further grouped the different mean group XA spectra into four main portions, the detailed decomposition of which can be found in Section D of the ESI.† Take note though, that the provided numbers should be seen as instructive examples for these hand-selected compositions. In future work, such compositions shall be obtained by carefully designed data-scientific fitting tools.
The next two signals can be found at 288.9 eV and 290.2 eV. The first is approximately reproduced by the orange component's peak at 289.2 eV in a mixture of π*(C–C), σ*(O–H), σ*(C–H) contributions of mainly the basal CH, basal C–OH and edge CH-Epo groups. The second signal is captured by the uppermost red component's peak at 290.1 eV (18% of carbon atoms) that mainly consists of σ*(C–O) and σ*(O–H) transitions of the edge CH–OH group. It shall be mentioned, though, that especially these two signals allow for a lot of different interpretations and are heavily discussed in literature. In other reports, the signal around 288.9 eV was assigned to CO transitions related to COOH groups,59,60 interlayer states,61 as well as contaminant effects on the GO surface.6 Although our theoretical study cannot rule out the possibility of interlayer effects or contaminants, we can confirm that the signal at 288.9 eV is, in principle, describable as π*(CO) related transition from both carbonylic as well as carboxylic group XA spectra. Nonetheless we believe that in this specific sample, those groups may have limited contribution to the peaks, because adding any of them to the composition XA spectrum would at the same time require a higher absorption intensity around 284 eV to account for π*(CC) transitions. It shall be mentioned that the σ*(C–C) transitions at 293.2 eV and 296.8 eV can be qualitatively reproduced, although with incorrect intensities since they are lying in the region beyond the edge-jump. Overall, based on the coefficients in the sum of individual functional groups in this estimated composition, we also find that the prominent π*-transition signal in the GO sample can be related to a relatively high percentage of 54 at% sp2 hybridized carbon atoms, of which around 75% account for unfunctionalised basal C atoms.
The reason for the energetic shifts when trying to describe the feature at 286.7 eV and 288.9 eV could be the underlying shift for the BH0.57LYP functional, which has been shown to slightly deviate for oxygen rich environments.35 In future work, it shall be investigated whether introducing group-specific shifts during construction of the mean group XA spectra improves the fitting results.
Next in the experimental spectrum, there is a steep rise in intensity around 287.3 eV that is best reproduced by the peak at 287.4 eV in the second lowest component (green, 40% of carbon atoms). This slightly smeared out signal results from mainly σ*(C–O) transitions of basal C-Epo groups. After that, one observes a very intense peak at 288.9 eV. In our composition the signal is described by the second highest and lowest components (orange and green) which peak at 289.2 eV and 289.1 eV, respectively. The orange components' (31% of carbon atoms) peak is dominated by σ*(C–H) transitions of the bridge CH and basal CH groups, while the green components' signal stems from a mixture of π*(C–C) and σ*(O–H) of the basal C–OH group.
Finally, the region between the intense signal at 287.3 eV and the ionisation edge jump is covered by the red component XA spectrum that peaks around 290.3 eV and makes up the remaining 19% of carbon atoms. In this region it contains σ*(C–O), σ*(O–H) and σ*(C–H) transitions from the edge CH2 and edge CH–OH groups.
Although the intensity of the signal at 288.9 eV is not reproduced by the group XA spectra very well, there may be different reasons for that. One reason could be that due to the averaging process in the formulation of group XA spectra, the intensity of specific groups is smeared out at this exact energy. Similar effects were observed when comparing the theoretical total XA spectra with the database compositions (see Section C of the ESI†), which hints at an insufficient indexing radius. Another reason for the very high intensity of this signal could be an excitation into an excitonic state caused by the quantum confinement of the oxidised CD.62
From the XA groups used in the composition XA spectrum of the oxidised CD, we can estimate that about 9% of the carbon atoms are sp2 hybridized, where more than half of these belong to unfunctionalized basal C atoms. This matches with the comparably low absorption intensity in the π* region of the experimental sample.
From this database, one can compose theoretical XA spectra of, in principle, any substance that is describable by the underlying model molecules in density functional theory quality. We confirmed this by successful reproduction of theoretical XA spectra from explicit calculations by simple addition of the respective database entries. Finally, we give an instructive example of how the methodology may be used as a complementary tool to analyse the functionalisation patterns from experimental XANES spectra of micrometer-sized graphene oxide (GO) flakes and oxidised carbon dots (CD).
We anticipate that the quality of the methodology is mainly limited by the complexity of the descriptors available. Hence, a future measure of improvement shall be the incorporation of more neighbour atoms' environments to the descriptors and expanding the number of model molecules. Furthermore, the underlying calculations will be generally improved by increasing the basis set size or the general level of theory. Finally, a data-scientific fitting algorithm shall be set up such that one is able to fit the composition of a measured experimental sample with respect to the carbon, oxygen and nitrogen pre-edge spectra simultaneously. Applying these multiple requirements will substantially improve the results of the composition methodology.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c8cp06620e |
This journal is © the Owner Societies 2019 |