Next Article in Journal
Recent Developments in Plasmonic Nanostructures for Metal Enhanced Fluorescence-Based Biosensing
Next Article in Special Issue
Terahertz Broadband Polarization Conversion for Transmitted Waves Based on Graphene Plasmon Resonances
Previous Article in Journal
Modification of Graphitic Carbon Nitride with Hydrogen Peroxide
 
 
Comment published on 26 December 2022, see Nanomaterials 2023, 13(1), 108.
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Introducing the Novel Mixed Gaussian-Lorentzian Lineshape in the Analysis of the Raman Signal of Biochar

1
Department of Applied Science and Technology, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129 Turin, Italy
2
Faculty of Science, Ontario Tech University, 2000 Simcoe Street North, Oshawa, ON L1G 0C5 T, Canada
3
Consorzio Interuniversitario Nazionale per la Scienza e Tecnologia dei Materiali (INSTM), Via G. Giusti 9, 50121 Florence, Italy
*
Author to whom correspondence should be addressed.
Nanomaterials 2020, 10(9), 1748; https://doi.org/10.3390/nano10091748
Submission received: 28 July 2020 / Revised: 28 August 2020 / Accepted: 31 August 2020 / Published: 3 September 2020
(This article belongs to the Special Issue Carbon-Based Materials: Growth, Characterization, and Applications)

Abstract

:
In this research, an innovative procedure is proposed to elaborate Raman spectra obtained from nanostructured and disordered solids. As a challenging case study, biochar, a bio-derived carbon based material, was selected. The complex structure of biochar (i.e., channeled surface, inorganic content) represents a serious challenge for Raman characterization. As widely reported, the Raman spectra are closely linked to thermal treatments of carbon material. The individual contributions to the Raman spectra are difficult to identify due to the numerous peaks that contribute to the spectra. To tackle this problem, we propose a brand new approach based on the introduction, on sound theoretical grounds, of a mixed Gaussian-–Lorentzian lineshape. As per the experimental part, biochar samples were carbonized in an inert atmosphere at various temperatures and their respective spectra were successfully decomposed using the new lineshape. The evolution of the structure with carbonization temperature was investigated by Raman and XRD analysis. The results of the two techniques fairly well agree. Compared to other approaches commonly reported in the literature this method (i) gives a sounder basis to the lineshape used in disordered materials, and (ii) appears to reduce the number of components, leading to an easier understanding of their origin.

1. Introduction

From the beginning of human history, carbon materials have played a crucial role in human life. In more recent times, new forms of carbon such as fullerene, carbon nanotubes and graphene have become widely available. These new forms have provided new perspectives of application to carbon materials. As the structural quality of a material is often critical to determine its properties, the structural investigation of carbon materials is of paramount importance. Among the characterization techniques used to investigate the local nanostructure of carbon materials Raman microscopy is the most widespread [1]. Raman spectroscopy offers a couple of interesting advantages. On the one hand it is a non destructive technique and as such can also be placed online as a quality control tool in production systems. On the other hand, each ordered carbon material has a distinctive Raman spectrum that constitutes its fingerprint [2] and relevant information can be additionally extracted from the spectra. Two specific spectral ranges are of interest for sp2 coordinated carbon materials, which will form the basis for later discussion: (i) the region in the range 1000 to 2000 cm−1 where the main peaks are D and G [2] seldom accompanied by less intense peaks, and (ii) the range 2400 to 3300 cm−1 where the overtones of the D and G peaks give their contribution [3]. In highly crystalline ordered materials the peaks are sharp and well separated so that the identification of the peaks and their quantification in terms of intensity, centers, and linewidth is straightforward, and their lineshape can be approximated to a Lorentzian curve [4]. However, when the carbon materials are disordered things become more complicated and blurred as the width of the peaks increases, and neighboring peaks start to merge. In the extreme situation of diamond like carbon, featureless bumps appear in both ranges and the identification of the specific components becomes rather challenging [5]. Lorentzian curves can also be used for carbon materials which are mainly sp2 coordinated, but in which the D and G peaks overlap to form a featured single signal in which two peaks are evident. Among the various carbon materials Biochar, a charcoal obtained by pyrolysis of vegetal feedstock in an oxygen-poor environment [6], is an optimum case study for our work, since its structure evolves with carbonization temperature and it is far more complex than carbon films or carbon soot due the massively presence of inorganics. As reported by Franklin [7], graphitization process is affected by the mobility of crystallites during the thermal treatment and it is influenced by C/H and C/O ratio of the material. During the years, the carbonization process has been modelized by using several approached reaching the actual distorted graphene triad proposed by McDonald-Wharry et al. [8]. Due to inorganic content and high oxygen content, biochar could be classified into non-graphitizable carbon according to Franklin definition but by increasing pyrolysis or carbonization temperature, biochar undergoes a turbostratic rearrangement with massive formation of sp2 cluster. Those domains are composed by randomly distributed disordered and nanographitic domains that can reached up to 80 nm in size and eventually evolve to graphite-like materials for temperature above 2000 °C as in the case of carbonization of cellulose nanocrystals [9].
A number of attempts have been made to decompose the 1000–2000 cm−1 range peak into its components [10,11]. To the best of our knowledge, attempts reported in the literature used two components with Gaussian lineshapes. Only Ghosh et al. [12] approached the problem by using a set of Gaussian and Lorentzian curves for Raman spectrum decomposition. However, they did not provide a physical basis for their choice. By fixing the constraint that the lineshape of all peaks should be Gaussian, the number of individual peaks needed to reasonably fit biochar Raman spectra turned out to be quite large [13]. Several authors used Voigt [14] and pseudo-Voigt [15,16] functions to take into account the presence of disordered nanographitic domains. Despite being basically a mix of Lorentzian and Gaussian, in their case the mixing occurs over the whole range of the signal, amounting to assume that two different types of regions (one more ordered, one more disordered) are present in the material, an assumption often not supported by other information. Even if nanographitic domains were identified by microscopical techniques, related Raman signals processing remain largely arbitrary [17]. Therefore, a problem in correlating each peak with a meaningful vibrating mode subsequently arises also in the Voigt and pseudo-Voigt case. Other attempts have been made in this direction; one of the most sophisticated studies is based on first principles computational simulation [18]. However, computation was made for ordered systems and no support from microscopy or other structural investigations was given. In summary, the very complex nature of nanographitic domains and local structure of biochar lead to the problem of properly decomposing Raman spectra being still an open challenge.
In the present paper we will show that by relaxing, on strong physical bases, the constraint of having a Gaussian lineshape, biochar Raman spectra can be decomposed by using two main peaks that can be readily identified as the usual D and G contributions. In order to verify that this approach is consistent we have taken a standardized biochar, produced by pyrolysis at 550 °C, and carbonized it at various temperatures in the range 900–2200 °C. The Raman spectra of all samples were recorded using a laser line at 514 nm. All spectra were satisfactorily decomposed with a single set of peaks. This result was obtained by introducing a lineshape in which three wavenumber ranges are identified for each peak. The central (i.e., closer to the peak highest value) region is fitted with a Lorentzian lineshape, while the two lateral regions are fitted with a Gaussian lineshape. The two lineshapes are matched at the range boundaries providing a unique lineshape for each component used. A thorough discussion in support of this assumption is given the text.

2. Materials and Methods

2.1. Materials

In order to test our new approach, we selected a standard biochar produced by UK biochar research centre (Edinburg, UK) using a rotary kiln pyrolytic unit [19]. The biochar was produced using oil seed rape (OSR 550) at a pyrolytic temperature of 550 °C [20].
In Table 1 we reported the main physiochemical characteristic of OSR 550.
OSR 550 was further investigated by using a field emission scanning electrical microscope (FE-SEM, Zeis SupraTM40, Oberkochen, Germany) equipped with an energy dispersive X-ray detector (EDX, Oxford Inca Energy 450, Oberkochen, Germany) that was used to explore the inorganic composition of OSR 550 as reported in Table 2.
OSR 550 was investigated in three different sites two of them (Site 1 and Site 2) rich in carbon and one rich in inorganic.
The biochar was further thermally treated (see next paragraph for details) in order to modify the structure from a disordered one to a well ordered one as reported by Lehman [22] so that our approach had to withstand a more tough challenge.

2.2. Methods

OSR 550 biochar was subjected to further thermal treatment at temperatures of 900 °C, 1100 °C, 1300 °C, 1500 °C and 2200 °C in order to modify the local structure of the material from disordered to ordered. A vacuum electric furnace from Pro.Ba., Cambiano, Italy was used to carry out the treatments at 900, 1100, 1300 and 1500 °C; the treatment at the highest temperature (2200 °C) was performed using a 448 TAV Vacuum Furnaces oven, Petroceramics, Stezzano, Italy. All treatments were performed under argon atmosphere (99.99% purity, controlled pressure 550 mbar) using a heating rate of 150 °C/h, a dwell at the maximum temperature for 30 min and a cooling to room temperature with the same thermal gradient used for heating. Samples labelling is summarized in Table 3.
Raman spectra were collected using a Renishaw inVia (H43662 model, Gloucestershire, UK) equipped with a green laser line (514 nm) with a 50× objective. Raman spectra were recorded in the range from 250 cm−1 to 3500 cm−1. Decomposition of Raman spectra, the core of the present paper, was focused on the range 1000–2000 cm−1 and performed with a homemade software developed using Matlab® (version R2020a).
X-ray diffraction (XRD) analyses were performed by using Panalytical X’PERT PRO PW3040/60 diffractometer, with Cu Kα radiation at 40 kV and 40 mA, Panalytical BV, Almelo, The Netherlands. The spectra were obtained from biochar powder in the 2θ range from 10 to 90° with a step size of 0.013°. For the decomposition of XRD peaks in the 2θ regions 15–35° and 40–50° two different software were employed in order to verify the reproducibility and consistency of the results: Panalytical HighScore Plus v3.0d and a homemade compiled software developed using Matlab® (version R2020a).
In the XRD analysis, the two peaks under investigation, which were assigned to 002 and 100 reflection of graphite respectively, have been used to determine the degree of biochar crystallinity as a function of temperature. Three main parameters were calculated: d002 which represents the interlayers spacing of aromatic layers; La and Lc on their turn represent the in-plane size and stacking height (size along c axis) of crystallites respectively. The values of d002, La and Lc reported below are the mean values obtained using the two software previously mentioned.
The distance between graphitic layers was estimated using the Bragg equation:
n λ = 2 d002 sin(θ)
where λ is the X-ray wavelength (0.154 nm as the Cu Kα radiation was used), and θ is the diffraction angle. La and Lc were calculated using Debye–Scherrer equation [23]
L i = K i λ B i   c o s θ
where Bi is the full width at half maximum (FWHM) values of the relevant peak, Ka is 1.84 and Kc is 0.89.

3. Results

3.1. New Lineshape Definition: Theoretical Bacground and Mathematical Details

As the key point of this paper consists in the use of a peculiar lineshape for fitting biochar Raman spectra, before presenting and discussing the results it is worth to focus on the origin of such a lineshape. In a continuous wave excitation process the local equilibrium polarizability is perturbed by the external wave while it is restored by the decay processes. This leads to a memory effect that is mapped into the time correlation function (TCF) of the dipole moments [24]. Given a certain resonant vibration its lineshape and linewidth in the frequency domain are determined by the Fourier transform of such TCF [25]. The exact lineshape of the peak related to a given vibration in the Raman spectra is hence related to the mechanisms that dissipate the energy associated to the vibration mode activated by the interaction with the incoming photons. If a single decay mechanism is at work and all the oscillators are in the same environment the typical Debye relaxation occurs, characterized by a single lifetime leading to a Lorentzian shape in the Raman spectra [24]. On the other hand, if the local environment of the same oscillator type at different locations of the solid is different, the sum up of all contributions leads to a different lineshape. In solid state physics it is quite customary [26,27] to assume that this decay relaxation is of a Gaussian type, and its Fourier transform in the Raman spectrum has a Gaussian shape as well. However, this can be an oversimplification, depending on the type and amount of local disorder. Strictly speaking Gaussian distribution can occur only if the various lifetimes are independent while in a partially disordered or in an amorphous solid the situation can be trickier.
A quick look at Figure 1 where a typical Raman spectrum belonging to our set is compared to the spectrum of a highly disordered amorphous carbon (DLC: Diamond-Like Carbon) and of crystalline graphite shows that biochar, as well as many other carbon materials, such as carbon fibers, cannot be associated to either of them. As a consequence, while the use of Gaussian lineshapes for the D and G peaks is appropriate for DLC [5] and that of Lorentzian lineshape is appropriate for graphite our case is worth of a deeper discussion.
Lineshapes other than Lorentzian and Gaussian have been proposed and justified on physical grounds. As the basic phenomena occurs in the time domain it is worth to discuss decay processes in such a domain. Among the various types of relaxation, the stretched exponential one, besides being of interest in liquids [28] is also found quite often in disordered solids [29]. Hence, it appears a good candidate for our material too. The resulting TCF was first discussed by Rothschild, Perrot and Guillaume [30,31,32] and is appropriate to describe strongly interacting liquids [33,34] as well as disordered solids [35]. However, we must note that:
  • we need to fit each Raman spectra in the frequency domain with a number of peaks by using a computational optimization algorithm;
  • the Fourier transform of the TCF arising from the stretched exponential relaxation cannot be calculated in closed form;
These points will lead for every step of the optimization process to numerically compute the Fourier transform of all peaks for every frequency value. As a consequence, the computation time will be beyond any reasonable timescale for a typical PC or workstation.
Taking into account this limitation and also that:
  • biochar materials appear to be only ‘partially’ disordered as shown in Figure 1
  • the effect of disorder arises more strikingly at short and large times [36]
we introduce a specific lineshape in the frequency domain that is able to account for all physical facts although not being yet fully justified on theoretical grounds: a mixed Gaussian-Lorentzian shape (GauLor), symmetrical with respect to its central wavenumber ω0, in which the central part is Lorentzian and the tails are Gaussians. We apply the further constraint that the curve must be of class C1 (i.e., both the functions and the derivatives must be continuous) in the transition points. Being ωth,up (>ω0) the onset wavenumber of the upper Gaussian tail, each resonant peak is described by a two valued function:
I L o r ( ω ) =   I L o r ( ω 0 )   Γ / π 4   ( ω ω 0 ) 2 + ( Γ ) 2
I G a u ( ω ) = I G a u ( ω 0 )   e ( ω     ω 0 ) 2 2   σ 2 σ   2 π
where σ is the width of the Gaussian and Γ is the dissipation term of the Lorentzian. By applying the C1 constraint at ωth,up we have:
σ = 4 ( ω t h , u p     ω 0 ) 2   +   Γ 2 8
To summarize, our function is defined as follows:
I GauLor ( ω ) = { I L o r ( ω ) for   ( ω t h , l o w ω ω t h . u p ) I G a u ( ω )   for   ( ω   ω t h , l o w )   A N D   ( ω   ω t h , u p )
where ωth,low is the threshold wavenumber for transition from Gaussian to Lorentzian shapes at wavenumbers in the low wavenumbers region.
A typical GauLor shape where the two transition points are symmetrical with respect to the peak wavenumber is reported in Figure 2.

3.2. Check Versus Stretched Exponential

As the stretched exponential is the most meaningful time decay in disordered materials, the first test for our GauLor is to check if it is mimicking well enough the Fourier transform arising from the TCF of a stretched exponential decay, that is described by:
I ( t ) = I 0   e ( t τ ) β
where τ is the relaxation time and β is a positive number lower than 1. Hence, we proceed as follows:
  • we fix a set of parameters for a stretched exponential relaxation (β = 0.8)
  • we compute the TCF Gν(t) by using equation [25]
    ln ( G ν ( t ) ) M 2   τ 2 = n = 0 ( 1 ) n ( t / τ ) 2 + n β n ! ( 1 + n   β ) ( 2 + n   β )
    where M2 (a constant value for a given peak) is the vibrational second moment [24]
  • we compute the TCF Fourier transform numerically for each wavenumber
  • we fit this Fourier transform with a single GauLor selecting its parameters with the fitting procedure detailed in the next paragraph
The results of this procedure are shown in Figure 3. We clearly see that:
-
the stretched exponential based Raman peak is symmetrical with respect to its peak wavenumber [28],
-
the GauLor curve mimics very well the stretched exponential based peak,
Taking into account (i) the physical significance of stretched relaxation for disordered materials, (ii) the symmetry of the Raman peak arising from it, and (iii) the excellent mimicking that GauLor provides for the stretched exponential based Raman peak, from now on we confidently use GauLor shapes to fit Raman spectra, assuming for each of them that the two transition points are symmetrical with respect to the peak position.

3.3. Details of the Fitting Procedure and Evaluation of the Contribution of the Lorentzian and Gaussian Parts to the Overall Intensity of a Peak

Given that the recorded Raman spectra are discrete and not continuous we have M wavenumber values in the wavenumber range of interest for the fit. For each wavenumber ωi the measured signal has an intensity S(ωi).
As a first step we define, on meaningful physical basis, the number N of GauLor components we use for the spectra decomposition. For each component we have to give initial tentative values to the following parameters: ω0, ωth,up, σ.
For each wavenumber ωi we compute the fitting curve intensity as follows
S F i t ( ω i ) = j = 1 N I G a u L o r ( ω i , ω 0 , j ,   ω t h , u p , j , σ j )
The optimization algorithm targets the maximum similarity between the measured and fitting curves. Such similarity is obtained applying the Levenberg–Marquardt algorithm [37] targeting the minimization parameter defined as:
Χ = i = 1 M [ S ( ω i ) S F i t ( ω i ) ] 2
The algorithm is run by varying the ω0, ωth,up, σ parameters for each GauLor until the value of X for two subsequent optimizing steps is below a given threshold.
As the role of disorder is to lead to a departure from the Lorentzian shape it is worthy to evaluate the relative contributions of Gaussian and Lorentzian parts to each GauLor. To this purpose we define I L o r t o t   as the total area of the GauLor between ωth,low and ωth,up where the GauLor has a lorentzian shape:
I L o r t o t = ω t h , l o w ω t h , u p I G a u L o r ( ω , ω 0 ,   ω t h , u p , σ ) d ω
By computing as well the total area I G a u L o r t o t under the GauLor peak:
I G a u L o r t o t = + I G a u L o r ( ω , ω 0 ,   ω t h , u p , σ ) d ω
we are able to evaluate the fraction of the GauLor area described by its Lorentzian shaped part:
f L = I L o r t o t I G a u L o r t o t
The fraction of the area for the Gaussian part is:
f G = 1 f L

3.4. Approach to the Fit of the Set of Samples

As the Raman signal of OSR 2200 show well resolved and sharp peaks (see below) we started by fitting it as the uniqueness of the fit is better guaranteed. With the set of GauLor components obtained we fitted the Raman spectra of all other samples. Any additional peak needed to fit such spectra with respect to OSR 2200 was justified on physical basis.

3.5. XRD Analysis

In order to better understand and support the results of our Raman analysis we have performed a thorough XRD analysis with the purpose to investigate the evolution of the nanographitc domains of biochar structure.
Figure 4 compares the X-ray diffraction patterns of biochar coming from oilseed rape after thermal treatment at increasing temperature from 550 °C to 2200 °C. Sample OSR 550 shows the presence of SiO2 in quartz form as the main crystalline phase; starting from OSR 900 cristobalite is present as well. Silica is a common component of biochar especially when it comes from grass, straw and grain husks [38]. Some low intense peaks appear in both the spectra at 2θ~21.8°, 28.0°, 29.3° and 30.2°. The presence of few peaks with very low intensity make difficult to state to which crystalline phase they belong. However, it appears that they cannot be assigned to a single compound. EDS analysis performed on OSR 550 [20] revealed the presence of magnesium, potassium, calcium, sulfur, silicon and aluminum as inorganic components. Moreover, the 2θ position of the peaks under investigation is consistent with the presence of Mg–Al silicate, CaCO3 and K2SO4. Therefore, we can hypothesize that these low intense peaks come from salts present as inorganic components of original biomass. This is in agreement with the literature: Xu et al. [39] reported the presence of a fraction of mineral compounds such as alkali metals or alkaline earth metals in form of carbonates, phosphates or oxides in many kinds of biochar, in addition to the carbon main component. At higher temperature these peaks disappeared probably due to thermal degradation phenomena.
The XRD patterns from OSR 550 to OSR 1500 show a broad hump in the 2θ region from 16° to 32° indicating that biochar contains a portion of highly disordered material in the form of amorphous carbon. When the treatment temperature reaches 1500 °C, a peak at 2θ~26°, characteristics of (002) graphite reflection, emerged from the hump.
The asymmetry of this band can be explained by the presence of a second peak at lower angle. Some authors [40,41,42] fitted the broad hump using two Gaussian curves centered at 20° and 25.5° respectively. The first peak called γ-band can be attributed to saturated structures such as aliphatic components linked to char crystallites edges. The second П-band at 2θ~25.5° can be attributed to the spacing between aromatic ring layers [43]. Moreover, the XRD spectrum of OSR 1500 shows (100) and (110) reflections of graphite at 2θ~43° and 77° respectively; these peaks are also present in OSR 550 and OSR 900 spectra, although with lower intensity and width. The presence of these more defined peaks confirms the existence of crystallites in the treated biochar [43] and their broad shape suggests a very small size of the formed crystallites. These observations suggest that the treatment at 1500 °C leads to an evolution of biochar carbon fraction in a structure which is intermediate between amorphous and the graphitic turbostratic structure [44].
A further increment of pyrolysis temperature up to 2200 °C led to the disappearance of the broad hump; the presence of sharper and more defined peaks clearly indicates a more crystalline structure of biochar. The turbostratic structure seems indeed to further evolve in a more ordered graphitic one whose presence is confirmed by a defined peak at 2θ~26.2°. The less intense peaks at 2θ~43° and 77° correspond to (100) and (110) reflections of graphite as stated above. These previously observed peaks become more shaped and defined when treatment temperature reached 2200 °C.
In addition, the spectra of OSR 1300, OSR 1500 and OSR 2200 show the presence of 3C-SiC phase coming from the reaction between SiO2 and carbon. The presence of β-SiC can be justified considering that the main commercial process for its production involves the carbothermal reduction of silica using different kind of carbonaceous materials such as petroleum coke, anthracite and bituminous coals or carbon black. The process is known as Acheson process [45]; it involves temperatures in excess of 2000 °C and it is strongly dependent both on temperature and gas composition. According to Kurunov et al. [46] the initial temperature for the reaction between silica and carbon ranges from 1200 to 1350 °C. Li et al. [47] reported the presence of a weak peak of β-SiC starting from 1300 °C for samples containing a mixture of silica and graphite treated in argon atmosphere. This justifies the presence of traces of SiC from OSR 1300 sample; it shows a cubic structure which is the most common polytype formed during the carbothermal reduction process.
Figure 5 shows the trend of the three parameters d002, La and Lc, which are representative of the average size of graphitic crystallites, as a function of the increasing thermal treatment temperatures.
Figure 5a reveals similar values of interlayers spacing after treatment at 550 °C and 900 °C which were 0.368 nm and 0.367 nm respectively. After thermal treatment at 1500 °C and 2200 °C, d002 steeply decreases reaching values of 0.344 nm and 0.340 nm, which are closer to the interlayer spacing of graphite (0.335 nm as reported by Jeong et al. [48]. Figure 5b,c reveal that the size of crystallites progressively increases with the increasing of treatment temperature: the in-plane size from around 0.40 nm for OSR 550 to 4.09 nm for OSR 2200. The crystallites size along c axis increases from 0.16 nm for OSR 550 to 1.41 nm for OSR 2200.
Furthermore, XRD analysis represent a powerful tool to evaluate the raising of SiC during the increment of treatment temperature. This evidence will be usefully in the next section for the explanation of Gaussian fraction (fG) trend.

3.6. Raman Spectra Analysis

In Figure 6 we report the Raman spectra measured for the set of thermally treated samples. Spectra have been normalized to the height of the right peak in order to render features more evident.
The spectra suggest that there is an evolution related to the treatment temperature and clearly evidence the correlation among the spectra obtained for the biochars treated at the various temperatures. In Figure 6f the Raman spectrum of OSR2200 is reported together with the fitting obtained using two main I GauLor shaped peaks, one (left–lower wavenumbers named D1 or D1) associated to the D-mode and the other (right–higher wavenumbers G1) associated to the G-mode. Two additional peaks are present, one at low wavenumber and one at high wavenumber namely respectively D2 and G2 (particularly evident in OSR 2200,) on which we will comment in detail later. The same set of GauLor was used to satisfactorily fit the Raman spectra obtained for most other temperatures (Figure 6b–e). As a first comment we note that our approach appears to be consistent in analyzing the various spectra as well as based on sound physical grounds as:
  • it is based on two well assessed peaks for carbon sp2 coordinated materials;
  • it takes into account the peculiarity of the impact of structural disorder and TCF on Raman spectra.
If we focus on the evolution of the signal for increasing temperature considering the thermal evolution of the material, we notice that the D-peak and the G-peak both become narrower.
D peak arises from aromatic rings confined in defective graphitic structures and it is related to the size of crystallites [49] and their edges [50], while G peak is related to in-plane carbon sp2 bonds stretching.
This assignation is supported by the trend vs treatment temperature of the various features of the spectra and of the GauLor components. Figure 7 shows the trends of the more relevant parameters as a function of the carbonization temperature.
D and G peaks were fitted by only two components for temperatures with the exception of OSR 550 where an additional component centered at 1702.2 cm−1 was added. This component has a very low intensity with a relative area percentage of 0.1 % and could be attributed to inorganic residues. Tilli et al. [51] attributed to CaCO3·nH2O active Raman modes ranging from 1738 cm−1 to 1746 cm−1. In our specific case, OSR 550 contains calcium and it is reasonable to assume that the calcium carbonates residues embedded into a carbon matrix could have low Raman frequencies due to the local carbonaceous environment. Furthermore, this component disappeared at higher temperature according to degradation of carbonates at temperature beyond 900 °C [52] and it was not detectable by XRD due the low amount of the species.
Comparing Figure 5 and Figure 7, G1 peak appears to be proportional to reduction of interlayers spacing of aromatic layers and increment of graphitic crystallites showing an origin related to the ordering of the carbonaceous material. G2 area percentage trend was more complex showing a constant increment of both area percentage and band width of up to 1300 °C and a decrement afterwards. We hypothesized G2 was related with bond angle disorder as suggest by Shimodaira et al. [53]. Considering graphitic turbostratic rearrangement, it is reasonable to assume a massive layer stacking driven by π–π orbital alignment when the temperature threshold was reached inducing a short-range order. We located this transition at 1300 °C when we observed a decrement of G2 area percentage due to the reorganization of carbon structures [54]. Both G1 and G2 peak positions were quite stable, being centered respectively at 1563 ± 17 cm−1 and 1601 ± 11 cm−1. D2 peak showed an opposite trend compared to G1 peak and we assumed it to be related to the presence of confined regions of small size as reported by Negri et al. [50]. The interpretation of the other D components was more complex and we hypothesized that they rose from two different but similar phenomena. According to Shimodaira et al. [53], the components named D1 and D1′ are both related to angle disorder but we advance the hypothesis that they were due by graphitic domains misalignment induced by inorganic components. Inorganics represent of up to 19.5 wt.% of OSR 550 [20] and are uniformly dispersed in the carbon matrix. When sp2 clusters started to grow inorganics were excluded from the same area inducing at same time a stress on the aromatic domains. This effect generated the D1 component, replaced by D1′ above 1300 °C. At a temperature higher than 1300 °C in fact, graphitic crystallites size rapidly grew simultaneously to the formation of SiC though carbothermal processes. In this environment, SiC region embedded into aromatic domains can perturb the bond angles of carbon planes originating the band D1′ that show a peak center position at higher frequencies than D1 (Figure 7c). The disorder related to D1 and D1′ was attested by fG values of 100% for D1 close to 90% for D1′. Similarly, D2 show a constant increment in fG attesting a disordering phenomenon. The fG trends of G1 and G2 were more complex showing a decrement with the increment of temperature up to 1100–1300 °C and increment with further temperature increase. These observations witness the presence of inclusion of SiC into graphitic domains formed at high temperature with a consequently Gaussian-like disordered phenomena. As clearly emerged, 1300 °C is a key temperature for the thermal evolution of the OSR to a more ordered state. This was supported also by the evaluation of ID/IG ratio reported in Figure 8.
ID/IG ratio was calculated as a ratio between the total area of the D and of the G components. It constantly decreased up to and beyond 1300 °C while between 1100 °C and 1300 °C increased. This behaviour is partially in agreement with Ferrari et al. [5] that reported an increment of ID/IG ratio following the transition from amorphous carbon to nanocrystalline graphite. In our specific case, we observed an additional phenomenon of ordering due to the degradation of all non-carbonaceous structure together with the clustering of the sp2 phase and evolution of inorganic matters with the formation of SiC. Infiltration of SiC into nanographitic domains was detailed described by Padovano et al. [55] proving its incorporation into the a carbon structure. The combination of the processes above could affect the D evolution creating more defects than in the case of the carbonization of pure carbon film generating the complex trend observed for fG.
We observed high values of fG for both components of D peak (Figure 5b) while the fG components of G peak decreased up to 1300 °C and increased afterwards. This could be explained by the raising of SiC formation, reported also by XRD analysis, that slow down the ordering process of graphitic crystallites. Centers of the components increased alongside all temperature range according with literature [5,56] while components width increased up to 1300 °C and decreased afterwards.
These trends confirmed the reorganization of the OSR with increasing of the carbonization temperature accordingly with XRD data from Figure 6. Both techniques proved that the material is still reorganizing itself to graphitic structure. This process cannot complete in biochar samples considered in this paper due to the presence of inorganic species as detected by XRD and also because of the natural porous structure of the material. These two factors are responsible for the somewhat hindered reorganization of nanographitic domains.
For a strong validation of proposed fitting methodology, we evaluated the fit of Raman spectra by using IGau, ILor and Voigt functions for each sample (Supplementary Materials) and we reported the most significant in Figure 9.
Fit of Raman spectra by using IGau (Figure 9a–c) required a high number of components of up to seven for the OSR 550 spectrum, up to five for OSR 900 and four for OSR 1100, OSR 1300 and OSR 1500 (Supplementary Materials, Figure S1). Spectrum of OSR 2200 (Figure 9c) required only three components but due to the particular shape of the IG peak two of them are concentric. This mathematical fitting was quite accurate but lacked physical ground and failed to explained the mode that originated the components.
Fit of Raman spectra by using ILor were challenging. Fit of OSR 550 (Figure 9d) and OSR 900 (Supplementary Materials, Figure S1) and required respectively five and four components. The fit OSR carbonized at 1100 °C–1500 °C degrees showed some problem related to ILor tails drop (see supporting information and Figure 9e for OSR1300). The line shape of did not allow a correct fit (Figure 9k,n,q) of ID and IG peaks with an overestimation of their area. Nonetheless, ILor fitted correctly the centered region of the ID and IG peaks of materials treated in a temperature range from 1100 °C up to 2200 °C.
Voigt functions results were comparable with IGauLor for OSR 550 (Figure 9g), OSR 1100 and OSR 1300 (Supplementary Materials, Figure S1). Nonetheless, they failed to produce an accurate fit for OSR 1300 (Figure 9h) and they showed the same problem of IGau in the fitting of OSR 2200 (Figure 9i) with the two components of IG peaks that were practically concentric. Fit with pseudo-Voigt lead results to the presence of components without any physical features (i.e., multimodal or negative curves) that could eventually (Supplementary Materials Figure S2).
From data reported in Figure 9, it was arguable that traditional established lineshapes cannot provide a solid base for an accurate fitting on a wide range of carbonization temperature contrary to IGauLor that showed a constant number of components with very accurate fit.
Moreover, it has to be noted that the peak position in the case of Voigt and Gaussian shapes do not move with a clear trend with treatment temperature.

4. Conclusions

In this paper, we successfully propose and validate a ground-breaking approach to Raman signal processing of disordered and nanostructured carbons. The GauLor lineshape introduced for the Raman fit is supported both by numerical accuracy and a valid physics explanation. To test this new methodology, Biochar was investigated as a complex case study due the massive presence of inorganics. These components deeply affected the carbonization process as clearly showed by the XRD analysis that enlighten the formation of SiC inclusion into carbonaceous structure. Due its composition and hierarchical structure, biochar Raman signal processing represent a challenging bench test for our new methodology. The biochar Raman signal was successfully fitted using this innovative GauLor lineshape and it was proved as a reasonable approach to describe complex carbon materials by comparing it with established methodologies. The procedure can also be used for all real carbon samples as it promotes a more accurate decomposition of Raman spectra.

Supplementary Materials

The following are available online at https://www.mdpi.com/2079-4991/10/9/1748/s1. Figure S1: Raman spectra and related fit using respectively IGau, ILor and Voigt functions of (a–c) OSR550, (d–f) OSR900, (g–i) OSR1100, (j–l) OSR 1300, (m–o) OSR1500 and (p–r) OSR 2200., Figure S2: Raman spectra and related fit using respectively pseudo-Voigt functions of (a) OSR550, (b) OSR900, (c) OSR1100, (d) OSR 1300, (e) OSR1500 and (f) OSR 2200

Author Contributions

Conceptualization, A.T.; methodology, A.T., M.R.; software, M.R.; formal analysis, M.R., E.P.; investigation, A.T., M.R., M.B.; resources, A.T.; data curation, A.T., M.R.; writing—original draft preparation, M.B.; writing—review and editing, M.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dresselhaus, M.S.; Jorio, A.; Hofmann, M.; Dresselhaus, G.; Saito, R. Perspectives on carbon nanotubes and graphene Raman spectroscopy. Nano Lett. 2010, 10, 751–758. [Google Scholar] [CrossRef] [PubMed]
  2. Merlen, A.; Buijnsters, J.; Pardanaud, C. A guide to and review of the use of multiwavelength Raman spectroscopy for characterizing defective aromatic carbon solids: From graphene to amorphous carbons. Coatings 2017, 7, 153. [Google Scholar] [CrossRef]
  3. Ferrari, A.C. Raman spectroscopy of graphene and graphite: Disorder, electron–phonon coupling, doping and nonadiabatic effects. Solid State Commun. 2007, 143, 47–57. [Google Scholar] [CrossRef]
  4. Frank, O.; Mohr, M.; Maultzsch, J.; Thomsen, C.; Riaz, I.; Jalil, R.; Novoselov, K.S.; Tsoukleri, G.; Parthenios, J.; Papagelis, K. Raman 2D-band splitting in graphene: Theory and experiment. ACS Nano 2011, 5, 2231–2239. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Ferrari, A.; Robertson, J. Resonant Raman spectroscopy of disordered, amorphous, and diamondlike carbon. Phys. Rev. B 2001, 64, 075414. [Google Scholar] [CrossRef] [Green Version]
  6. Lehmann, J.; Joseph, S. Biochar for Environmental Management: Science, Technology and Implementation; Routledge: London, UK, 2015. [Google Scholar]
  7. Franklin, R.E. Crystallite growth in graphitizing and non-graphitizing carbons. Proc. R. Soc. Lond. Ser. A Math. Phys. Sci. 1951, 209, 196–218. [Google Scholar]
  8. McDonald-Wharry, J.S.; Manley-Harris, M.; Pickering, K.L. Reviewing, Combining, and Updating the Models for the Nanostructure of Non-Graphitizing Carbons Produced from Oxygen-Containing Precursors. Energy Fuels 2016, 30, 7811–7826. [Google Scholar] [CrossRef]
  9. Eom, Y.; Son, S.M.; Kim, Y.E.; Lee, J.-E.; Hwang, S.-H.; Chae, H.G. Structure evolution mechanism of highly ordered graphite during carbonization of cellulose nanocrystals. Carbon 2019, 150, 142–152. [Google Scholar] [CrossRef]
  10. Potgieter-Vermaak, S.; Maledi, N.; Wagner, N.; Van Heerden, J.; Van Grieken, R.; Potgieter, J. Raman spectroscopy for the analysis of coal: A review. J. Raman Spectrosc. 2011, 42, 123–129. [Google Scholar] [CrossRef]
  11. Ferrari, A.C.; Robertson, J. Raman spectroscopy of amorphous, nanostructured, diamond–like carbon, and nanodiamond. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 2004, 362, 2477–2512. [Google Scholar] [CrossRef]
  12. Ghosh, A.; Razzino, C.d.A.; Dasgupta, A.; Fujisawa, K.; Vieira, L.H.S.; Subramanian, S.; Costa, R.S.; Lobo, A.O.; Ferreira, O.P.; Robinson, J.; et al. Structural and electrochemical properties of babassu coconut mesocarp-generated activated carbon and few-layer graphene. Carbon 2019, 145, 175–186. [Google Scholar] [CrossRef]
  13. Feng, D.; Zhao, Y.; Zhang, Y.; Sun, S. Effects of H2O and CO2 on the homogeneous conversion and heterogeneous reforming of biomass tar over biochar. Int. J. Hydrogen Energy 2017, 42, 13070–13084. [Google Scholar] [CrossRef]
  14. Claramunt, S.; Varea, A.; López-Díaz, D.; Velázquez, M.M.; Cornet, A.; Cirera, A. The importance of interbands on the interpretation of the Raman spectrum of graphene oxide. J. Phys. Chem. C 2015, 119, 10123–10129. [Google Scholar] [CrossRef]
  15. Sadezky, A.; Muckenhuber, H.; Grothe, H.; Niessner, R.; Pöschl, U. Raman microspectroscopy of soot and related carbonaceous materials: Spectral analysis and structural information. Carbon 2005, 43, 1731–1742. [Google Scholar] [CrossRef]
  16. Zickler, G.A.; Smarsly, B.; Gierlinger, N.; Peterlik, H.; Paris, O. A reconsideration of the relationship between the crystallite size La of carbons determined by X-ray diffraction and Raman spectroscopy. Carbon 2006, 44, 3239–3246. [Google Scholar] [CrossRef]
  17. Sousa, D.V.d.; Guimarães, L.M.; Félix, J.F.; Ker, J.C.; Schaefer, C.E.R.; Rodet, M.J. Dynamic of the structural alteration of biochar in ancient Anthrosol over a long timescale by Raman spectroscopy. PLoS ONE 2020, 15, e0229447. [Google Scholar] [CrossRef]
  18. Smith, M.W.; Dallmeyer, I.; Johnson, T.J.; Brauer, C.S.; McEwen, J.-S.; Espinal, J.F.; Garcia-Perez, M. Structural analysis of char by Raman spectroscopy: Improving band assignments through computational calculations from first principles. Carbon 2016, 100, 678–692. [Google Scholar] [CrossRef] [Green Version]
  19. Mašek, O.; Buss, W.; Roy-Poirier, A.; Lowe, W.; Peters, C.; Brownsort, P.; Mignard, D.; Pritchard, C.; Sohi, S. Consistency of biochar properties over time and production scales: A characterisation of standard materials. J. Anal. Appl. Pyrolysis 2018, 132, 200–210. [Google Scholar] [CrossRef] [Green Version]
  20. UK Biochar Research Centre. Available online: https://www.biochar.ac.uk/standard_materials.php (accessed on 19 March 2020).
  21. OSR550 Standard Biochar Specification Sheet—Version 1.0. Available online: http://www.charchive.org/record.php?record_id=98 (accessed on 13 August 2020).
  22. Lehmann, J.; Joseph, S. Biochar for Environmental Management; Earthscan: London, UK, 2009. [Google Scholar]
  23. Holzwarth, U.; Gibson, N. The Scherrer equation versus the “Debye-Scherrer equation”. Nature Nanotechnol. 2011, 6, 534. [Google Scholar] [CrossRef]
  24. Kirillov, S. Novel approaches in spectroscopy of interparticle inter-actions. Vibrational line profiles and anomalous non-coincidence effects. In NoVel Approaches to the Structure and Dynamics of Liquids: Experiments, Theories and Simulations; Springer: Berlin/Heidelberg, Germany, 2004; pp. 193–227. [Google Scholar]
  25. Kirillov, S.A. Time-correlation functions from band-shape fits without Fourier transform. Chem. Phys. Lett. 1999, 303, 37–42. [Google Scholar] [CrossRef]
  26. Birge, N.; Jeong, Y.; Nagel, S.R.; Bhattacharya, S.; Susman, S. Distribution of relaxation times in (KBr) 0.5 (KCN) 0.5. Phys. Rev. B 1984, 30, 2306. [Google Scholar] [CrossRef]
  27. Kamba, S.; Porokhonskyy, V.; Pashkin, A.; Bovtun, V.; Petzelt, J.; Nino, J.C.; Trolier-McKinstry, S.; Lanagan, M.T.; Randall, C.A. Anomalous broad dielectric relaxation in Bi1.5Zn1.0Nb1.5O7 pyrochlore. Phys. Rev. B 2002, 66, 054106. [Google Scholar] [CrossRef]
  28. Kirillov, S.A. Markovian frequency modulation in liquids. Analytical description and comparison with the stretched exponential approach. Chem. Phys. Lett. 1993, 202, 459–463. [Google Scholar] [CrossRef]
  29. Kakalios, J.; Street, R.; Jackson, W. Stretched-exponential relaxation arising from dispersive diffusion of hydrogen in amorphous silicon. Phys. Rev. Lett. 1987, 59, 1037. [Google Scholar] [CrossRef] [PubMed]
  30. Rothschild, W.G.; Perrot, M.; Guillaume, F. Vibrational dephasing under fractional (“stretched”) exponential modulation. Chem. Phys. Lett. 1986, 128, 591–594. [Google Scholar] [CrossRef]
  31. Rothschild, W.G.; Perrot, M.; Guillaume, F. On the vibrational T2 processes in partially ordered systems. J. Chem. Phys. 1987, 87, 7293–7299. [Google Scholar] [CrossRef]
  32. Burnstein, A. Motion-broadened and motion-narrowed spectra. Chem. Phys. Lett. 1981, 83, 335–340. [Google Scholar] [CrossRef]
  33. Kirillov, S.A.; Yannopoulos, S.N. Charge-current contribution to low-frequency Raman scattering from glass-forming ionic liquids. Phys. Rev. B 2000, 61, 11391. [Google Scholar] [CrossRef]
  34. Pick, R.; Aouadi, A.; Dreyfus, C.; Torre, R.; Bartolini, P. Light scattering studies of high-frequency relaxation processes in organic glass formers. J. Phys. Condens. Matter 1996, 8, 9593. [Google Scholar] [CrossRef]
  35. Karmwar, P.; Graeser, K.; Gordon, K.C.; Strachan, C.J.; Rades, T. Investigation of properties and recrystallisation behaviour of amorphous indomethacin samples prepared by different methods. Int. J. Pharm. 2011, 417, 94–100. [Google Scholar] [CrossRef]
  36. Kalampounias, A.; Yannopoulos, S.; Steffen, W.; Kirillova, L.; Kirillov, S. Short-time dynamics of glass-forming liquids: Phenyl salicylate (salol) in bulk liquid, dilute solution, and confining geometries. J. Chem. Phys. 2003, 118, 8340–8349. [Google Scholar] [CrossRef]
  37. Moré, J.J. The Levenberg-Marquardt algorithm: Implementation and theory. In Numerical Analysis; Springer: Berlin/Heidelberg, Germany, 1978; pp. 105–116. [Google Scholar]
  38. Amonette, J.E.; Joseph, S. Characteristics of biochar: Microchemical properties. In Biochar for Environmental Management; Routledge: London, UK, 2012; pp. 65–84. [Google Scholar]
  39. Xu, X.; Cao, X.; Zhao, L.; Zhou, H.; Luo, Q. Interaction of organic and inorganic fractions of biochar with Pb (II) ion: Further elucidation of mechanisms for Pb (II) removal by biochar. RSC Adv. 2014, 4, 44930–44937. [Google Scholar] [CrossRef]
  40. Manoj, B. Investigation of nanocrystalline structure in selected carbonaceous materials. Int. J. Miner. Met. Mater. 2014, 21, 940–946. [Google Scholar] [CrossRef]
  41. Takagi, H.; Maruyama, K.; Yoshizawa, N.; Yamada, Y.; Sato, Y. XRD analysis of carbon stacking structure in coal during heat treatment. Fuel 2004, 83, 2427–2433. [Google Scholar] [CrossRef]
  42. Li, Z.; Lu, C.; Xia, Z.; Zhou, Y.; Luo, Z. X-ray diffraction patterns of graphite and turbostratic carbon. Carbon 2007, 45, 1686–1695. [Google Scholar] [CrossRef]
  43. Manoj, B.; Kunjomana, A. Study of stacking structure of amorphous carbon by X-ray diffraction technique. Int. J. Electrochem. Sci. 2012, 7, 3127–3134. [Google Scholar]
  44. Gibson, J.; Holohan, M.; Riley, H. 87. Amorphous carbon. J. Chem. Soc. Resumed 1946, 0, 456–461. [Google Scholar] [CrossRef]
  45. Kumar, P.V.; Gupta, G.S. Study of formation of silicon carbide in the Acheson process. Steel Res. 2002, 73, 31–38. [Google Scholar] [CrossRef]
  46. Kurunov, I. Mechanism and theoretical evaluation of the replacement of coke by schungite in blast-furnace smelting. Izv. Vyssh. Uchebn. Zaved. Chern. Metall 2003, 7, 20–24. [Google Scholar]
  47. Li, X.; Chen, X.; Song, H. Synthesis of β-SiC nanostructures via the carbothermal reduction of resorcinol–formaldehyde/SiO2 hybrid aerogels. J. Mater. Sci. 2009, 44, 4661–4667. [Google Scholar] [CrossRef]
  48. Jeong, H.-K.; Noh, H.-J.; Kim, J.-Y.; Jin, M.; Park, C.; Lee, Y. X-ray absorption spectroscopy of graphite oxide. EPL Europhys. Lett. 2008, 82, 67004. [Google Scholar] [CrossRef] [Green Version]
  49. Tuinstra, F.; Koenig, J.L. Raman spectrum of graphite. J. Chem. Phys. 1970, 53, 1126–1130. [Google Scholar] [CrossRef] [Green Version]
  50. Negri, F.; di Donato, E.; Tommasini, M.; Castiglioni, C.; Zerbi, G.; Müllen, K. Resonance Raman contribution to the D band of carbon materials: Modeling defects with quantum chemistry. J. Chem. Phys. 2004, 120, 11889–11900. [Google Scholar] [CrossRef] [PubMed]
  51. Tlili, M.; Amor, M.B.; Gabrielli, C.; Joiret, S.; Maurin, G.; Rousseau, P. Characterization of CaCO3 hydrates by micro-Raman spectroscopy. J. Raman Spectrosc. 2002, 33, 10–16. [Google Scholar] [CrossRef]
  52. Tonsuaadu, K.; Gross, K.A.; Plūduma, L.; Veiderma, M. A review on the thermal stability of calcium apatites. J. Therm. Anal. Calorim. 2012, 110, 647–659. [Google Scholar] [CrossRef]
  53. Shimodaira, N.; Masui, A. Raman spectroscopic investigations of activated carbon materials. J. Appl. Phys. 2002, 92, 902–909. [Google Scholar] [CrossRef]
  54. Dillon, R.; Woollam, J.A.; Katkanant, V. Use of Raman scattering to investigate disorder and crystallite formation in as-deposited and annealed carbon films. Phys. Rev. B 1984, 29, 3482. [Google Scholar] [CrossRef]
  55. Padovano, E.; Giorcelli, M.; Bianchi, G.; Biamino, S.; Rovere, M.; Tagliaferro, A.; Ortona, A. Graphite-Si-SiC ceramics produced by microwave assisted reactive melt infiltration. J. Eur. Ceram. Soc. 2019, 39, 2232–2243. [Google Scholar] [CrossRef]
  56. Matthews, M.; Pimenta, M.; Dresselhaus, G.; Dresselhaus, M.; Endo, M. Origin of dispersive effects of the Raman D band in carbon materials. Phys. Rev. B 1999, 59, R6585–R6588. [Google Scholar] [CrossRef]
Figure 1. Overview on the lineshape of a highly ordered graphite (solid black), a Diamond-Like Carbon (DLC) (solid red) and oil seed rape (OSR) 550 (solid blue).
Figure 1. Overview on the lineshape of a highly ordered graphite (solid black), a Diamond-Like Carbon (DLC) (solid red) and oil seed rape (OSR) 550 (solid blue).
Nanomaterials 10 01748 g001
Figure 2. Superimposition of IGauLor (solid black), IGau (dotted red) and ILor 550 (dotted blue).
Figure 2. Superimposition of IGauLor (solid black), IGau (dotted red) and ILor 550 (dotted blue).
Nanomaterials 10 01748 g002
Figure 3. Superimposition of IGauLor (dotted blue) and stretched exponential (solid red).
Figure 3. Superimposition of IGauLor (dotted blue) and stretched exponential (solid red).
Nanomaterials 10 01748 g003
Figure 4. XRD patterns of biochar after thermal treatment at the various temperatures.
Figure 4. XRD patterns of biochar after thermal treatment at the various temperatures.
Nanomaterials 10 01748 g004
Figure 5. Variation of crystallites parameters as function of thermal treatment temperatures: (a) interlayers spacing of aromatic layers, (b) in-plane size and (c) stacking height of crystallites.
Figure 5. Variation of crystallites parameters as function of thermal treatment temperatures: (a) interlayers spacing of aromatic layers, (b) in-plane size and (c) stacking height of crystallites.
Nanomaterials 10 01748 g005
Figure 6. Raman spectra and related fit using IGauLor of (a) OSR550, (b) OSR900, (c) OSR1100, (d) OSR 1300, (e) OSR1500 and (f) OSR 2200. Original Raman spectra are reported in black, the fitted Raman spectra are reported in red.
Figure 6. Raman spectra and related fit using IGauLor of (a) OSR550, (b) OSR900, (c) OSR1100, (d) OSR 1300, (e) OSR1500 and (f) OSR 2200. Original Raman spectra are reported in black, the fitted Raman spectra are reported in red.
Nanomaterials 10 01748 g006
Figure 7. Components trends for (a) area percentage, (b) fG, (c) peaks centers and (d) widths.
Figure 7. Components trends for (a) area percentage, (b) fG, (c) peaks centers and (d) widths.
Nanomaterials 10 01748 g007
Figure 8. Trend between ID/IG ratio and carbonization temperature of OSR samples.
Figure 8. Trend between ID/IG ratio and carbonization temperature of OSR samples.
Nanomaterials 10 01748 g008
Figure 9. Raman spectra and related fit using IGau (ac), ILor (df) and Voigt functions (gi) of OSR 550, OSR 1300 and OSR 2200 respectively. Original Raman spectra are reported by using solid black lines, the fitted Raman spectra are reported by using solid red lines and each component is reported in dashed black lines.
Figure 9. Raman spectra and related fit using IGau (ac), ILor (df) and Voigt functions (gi) of OSR 550, OSR 1300 and OSR 2200 respectively. Original Raman spectra are reported by using solid black lines, the fitted Raman spectra are reported by using solid red lines and each component is reported in dashed black lines.
Nanomaterials 10 01748 g009
Table 1. Main chemical and physical properties of biochar samples as reported by safety data sheets [21].
Table 1. Main chemical and physical properties of biochar samples as reported by safety data sheets [21].
4Proximate Analysis
[wt%]
Ultimate Analysis [%]Specific Surface Area [m2/g]
MoistureAshVOCsFixed Carbon aCHO b
OSR 5502.6119.5016.3861.5168.851.828.917.3
a Calculated as difference, b calculated as difference considering the ash content.
Table 2. Elemental composition of OSR 550 as emerged by energy dispersive X-ray detector (EDX).
Table 2. Elemental composition of OSR 550 as emerged by energy dispersive X-ray detector (EDX).
ElementConcentration [wt.%]
Site 1Site 2Site 3
C86.7990.6734.45
O10.828.4942.33
Mg0.100.050.29
Si0.330.150.41
P0.120.050.16
Cl0.100.031.88
K1.320.4220.48
Ca0.420.1334.45
Table 3. Samples labelling.
Table 3. Samples labelling.
Samples Label.Treatment Temperature [°C]
OSR 550550 a
OSR 900900 b
OSR 10001000 b
OSR 11001100 b
OSR 13001300 b
OSR 15001500 b
OSR 22002200 b
a pyrolysis process, b carbonization process.

Share and Cite

MDPI and ACS Style

Tagliaferro, A.; Rovere, M.; Padovano, E.; Bartoli, M.; Giorcelli, M. Introducing the Novel Mixed Gaussian-Lorentzian Lineshape in the Analysis of the Raman Signal of Biochar. Nanomaterials 2020, 10, 1748. https://doi.org/10.3390/nano10091748

AMA Style

Tagliaferro A, Rovere M, Padovano E, Bartoli M, Giorcelli M. Introducing the Novel Mixed Gaussian-Lorentzian Lineshape in the Analysis of the Raman Signal of Biochar. Nanomaterials. 2020; 10(9):1748. https://doi.org/10.3390/nano10091748

Chicago/Turabian Style

Tagliaferro, Alberto, Massimo Rovere, Elisa Padovano, Mattia Bartoli, and Mauro Giorcelli. 2020. "Introducing the Novel Mixed Gaussian-Lorentzian Lineshape in the Analysis of the Raman Signal of Biochar" Nanomaterials 10, no. 9: 1748. https://doi.org/10.3390/nano10091748

APA Style

Tagliaferro, A., Rovere, M., Padovano, E., Bartoli, M., & Giorcelli, M. (2020). Introducing the Novel Mixed Gaussian-Lorentzian Lineshape in the Analysis of the Raman Signal of Biochar. Nanomaterials, 10(9), 1748. https://doi.org/10.3390/nano10091748

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop