Next Article in Journal
Direct 3D Printing of Clear Orthodontic Aligners: Current State and Future Possibilities
Previous Article in Journal
Investigations on the Effect of Layers’ Thickness and Orientations in the Machining of Additively Manufactured Stainless Steel 316L
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Establishing the Carbonation Profile with Raman Spectroscopy: Effects of Fly Ash and Ground Granulated Blast Furnace Slag

1
College of Materials Science and Engineering, Chongqing University, 174 Shazheng Street, Shapingba, Chongqing 400044, China
2
CRANN and AMBER Research Centers, Trinity College Dublin, D02 PN40 Dublin 2, Ireland
3
School of Civil Engineering, University of Leeds, Woodhouse Lane, Leeds LS2 9JT, UK
4
Department of Civil, Environmental and Geomatic Engineering, University College London, Gower Street, London WC1E 6BT, UK
*
Author to whom correspondence should be addressed.
Submission received: 21 February 2021 / Revised: 31 March 2021 / Accepted: 2 April 2021 / Published: 5 April 2021
(This article belongs to the Section Construction and Building Materials)

Abstract

:
Establishing the carbonation profile is of great significance to the prediction of the service life of reinforced concrete structures. In our previous work, Raman spectroscopy was shown to be an efficient tool for characterizing calcium carbonate (CaCO3) polymorphs and their profile in plain Portland cement (PC) matrices. However, as supplementary cementitious materials (SCMs), particularly fly ash (FA) and ground granulated blast furnace slag (GGBS), are widely used in concrete, establishing the carbonation profile without considering the possible effects of these SCMs could be of little significance to the real world. This paper, thus, investigated the effects of FA and GGBS on the working capacity and reliability of Raman spectroscopy for establishing the carbonation profile in PC blends containing SCMs. The thermogravimetry (TG) analysis was also conducted to verify the results from Raman spectroscopy. The results show that Raman spectroscopy demonstrated a good capacity for differentiating the variation of CaCO3 contents in FA or GGBS blends. However, the incorporation of FA and GGBS into the PC system caused some adverse effects on the quantification of CaCO3 by Raman spectroscopy, which could be attributed to the darker color and weak scatter nature of FA and the high content of glassy phases in GGBS.

1. Introduction

Concrete is inevitably exposed to carbonation attacks during its service life owing to the presence of carbon dioxide (CO2) such as that in the atmosphere. Briefly, carbonation occurs when the CO2 diffuses into concrete through its pore network, dissolves in the pore solution as CO32− ions and then reacts with various phases in the cement matrix [1]. Whilst virtually almost all of the Ca-bearing phases are susceptible to carbonation, calcium hydroxide (Ca(OH)2, CH, portlandite) has been commonly considered to be the phase that reacts most readily with CO2, leading to the formation of calcium carbonate (CaCO3) [1,2]. The consumption of CH in particular causes the drop of the pH in concrete, which subsequently triggers the depassivation of the rebar and ultimately leads to its corrosion in concrete [1,3]. In addition to the reaction with CH, CO2 can also attack calcium silicate hydrate (C-S-H), leading to the decalcification of the C-S-H gel as well as the formation of CaCO3 [4,5]. Evidently, the carbonation of concrete contributes significantly to the formation of CaCO3, which could exist in different polymorphs; namely, calcite, aragonite and vaterite [6]. Compared with the naturally occurring carbonation with a low CO2 concentration (around 0.03%) at a very slow rate, the accelerated carbonation (with a higher CO2 concentration), which is commonly adopted in laboratory evaluations, generally leads to the co-existence of all of these three polymorphs [7]. Furthermore, Ca-modified amorphous silica gel could also be formed from the carbonation of the C-S-H gel [7]. Although different products could be formed from the carbonation reaction, the most important phase that needs particular attention is the CaCO3 polymorph phase as this can cause significant changes to the properties of concrete. Specifically, whilst carbonation can enhance the mechanical properties of concrete through densifying the microstructure, the reduction in pH on the other hand can result in the corrosion of the reinforcing steel and hence the degradation of the performance of concrete. Therefore, accurate information about the carbonation in concrete will, no doubt, benefit the assessment of the durability performance of existing reinforced concrete (RC) structures [8,9].
In general, the amount of CaCO3 formed in carbonated concrete decreases with the increase of the distance from the exposed surface of concrete, which is mainly due to the reduced carbonation over the depth into concrete [5,10]. By plotting the CaCO3 content against the depth into concrete, a visual illustration of the distribution of CaCO3 inside concrete could be produced, forming the so-called ‘carbonation profile’ [5,6]. Establishing an accurate and reliable carbonation profile has been considered important not only because this can allow us to follow the carbonation front, but more importantly it can facilitate the prediction of the initiation of rebar corrosion and the service life of a reinforced concrete structure [11]. Additionally, as the formation of CaCO3 can also change a few other properties of concrete such as refining the pore structure as well as increasing the compressive and tensile strengths [12,13] a good knowledge of the distribution of CaCO3 may also benefit the understanding of the cause of these changes. Traditionally, the phenolphthalein spray test has been widely used to determine the carbonation depth by spraying the phenolphthalein indicator onto the freshly split surface of a concrete prism and taking a visual judgement of the color boundary [14]. Although this test can clearly indicate the fully carbonated region, it has been heavily criticized for its incapability of differentiating the uncarbonated and partially carbonated regions [15]. Moreover, the phenolphthalein test also cannot determine the content of CaCO3. Therefore, efforts have been made in the past to find better alternatives for establishing the carbonation profile. For example, the thermogravimetry (TG)/derivative thermogravimetry (DTG) analysis has been employed to establish CaCO3 content at different depths of specimens by tracing the DTG diagram related to CaCO3 decomposition [5,16]. However, this technique is less satisfactory in differentiating each CaCO3 polymorph, particularly the CaCO3 phase with a low thermal stability, as it decomposes in the near/same temperature range as CH. Although the X-ray diffraction (XRD) analysis could be used to establish the distribution of CaCO3, its sample preparation and data processing are complicated [17]. Other techniques such as the Gamma-ray attenuation method (GRAM)/Gammadensimetry and Fourier transformation infrared spectroscopy (FTIR) all suffer different limitations for quantifying different CaCO3 phases [6,14,15,16].
Raman spectroscopy, one of the most important vibrational spectroscopic analysis technologies, is an advanced analytical tool combining both the fingerprint identification capacity and quantification potential. As the wavelength (frequency) shift between the incident laser light and the Raman scattered light is specific to the chemical bond and the symmetry and vibrations of molecules, the Raman scattered signal (i.e., Raman shift) can thus provide ‘fingerprint’ information to recognize substances [18,19]. Due to its various advantages, researchers have attempted to employ Raman spectroscopy in the area of cement and concrete, e.g., for analyzing bottom ash fired bricks [20], as well as for the carbonation products formed in cementitious materials [21]. Although all of the calcium carbonate polymorphs, i.e., calcite, vaterite and aragonite, have been successfully characterized in previous research, no attempt has been made to quantify the calcium carbonate content in order to establish the carbonation profile. Besides, no research has been carried out on the Raman spectroscopic investigation of the carbonation in blended cement materials containing supplementary cementitious materials (SCMs) (e.g., fly ash (FA), ground granulated blast furnace slag (GGBS)), let alone establishing the carbonation profile in the blended cementitious matrices. In recent years, the authors have carried out extensive research to apply both bench-mounted Raman spectroscopy and the optical fiber Raman technique in cement and concrete, covering the different aspects related to cementitious materials such as the anhydrous phases and hydration products as well as the deteriorated phases/products [22,23,24,25,26]. It is worth highlighting that previous works have successfully developed bespoke optics-tailored fiber optic Raman systems for characterizing the carbonation products of concrete. In particular, the systems not only can qualitatively differentiate the polymorphs of calcium carbonates but can also quantitatively establish their profiles [23,26]. Nonetheless, previous studies were fully performed on plain (pure) Portland cement (PC) pastes without incorporating any SCMs. However, as SCMs, particularly FA and GGBS, are widely used in concrete practice to meet the increased requirements for durability and sustainability, establishing the carbonation profile without considering the possible effects of these SCMs on Raman spectra could be of little practical significance. This is, in particular, due to the fact that FA is an intrinsically weak Raman scatterer, which, coupled with its complex (e.g., iron oxide) and heterogeneous nature, has made it a difficult material to be analyzed by Raman spectroscopy [27]. GGBS on the other hand contains almost non-crystalline phases and has also been considered to be tricky to be analyzed by Raman spectroscopy probably due to its high percentage of disordered phases [27]. Indeed, intense fluorescence has been observed both in our own investigations and the studies reported in the literature when Raman spectroscopy has been employed to characterize FA or GGBS [27]. Depending on the intensity of the fluorescence, it could make a strong disturbance to or even totally hamper the weak Raman signals being obtained. Consequently, concerns could be raised over their potential effects on the quantification of CaCO3 when the Raman spectroscopy technique is used to establish the carbonation profile of PC blends containing FA or GGBS.
The current work thus aims to investigate the working capacity of Raman spectroscopy for establishing the carbonation profile of FA or GGBS blended PC systems and more importantly to explore the potential effects that FA and GGBS may have on the working capacity and reliability of Raman spectroscopy. The CaCO3 polymorphs and the carbonation profile (i.e., the CaCO3 content against the depth) in the FA and GGBS blended cement pastes after being exposed to accelerated carbonation were characterized by a bench-mounted Raman spectroscope. The signal-to-noise ratios (SNR) of the Raman spectra collected from different cement systems were compared and the TG analysis was adopted as a supplementary tool to verify the Raman results and quantitatively assess the effects of FA and GGBS on the carbonation profile established based on the content of CaCO3. It is anticipated that the results obtained from this study will lay a foundation for analyzing the similar effects that FA and GGBS may have on the optical fiber Raman system in the future.

2. Materials and Methods

2.1. Materials

In this study, Portland cement (PC) CEM I 42.5N supplied by Quinn Cement (Derrylin, County Fermanagh, Northern Ireland, UK) (BS EN 197-1:2011) was used to formulate blended cementitious systems. The low calcium fly ash (FA), conforming to the Category S according to BS EN 450-1:2012, was supplied by Connexpo (N.I.) Ltd. from Kilroot Power Station in Antrim Northern Ireland, UK. The specific surface area and specific density of the FA was 340 m2/kg and 2520 kg/m3, respectively. The ground granulated blast furnace slag (GGBS) was supplied by Civil and Marine Ltd., Maidenhead, UK. The GGBS was an off-white color powder with a fineness of 527 m2/kg and a loose bulk density of 1000–1100 kg/m3. The photo indicating the color of the FA, PC and GGBS is shown in Figure 1 and their chemical compositions are presented in Table 1. A Sikaguard 680-s acrylic paint supplied by Sika Switzerland was employed to coat the paste samples before carbonation (as detailed in Section 2.3).

2.2. Manufacture of FA/PC and GGBS/PC Blends

The FA/PC and GGBS/PC pastes were manufactured with a water/binder ratio of 0.5, which enabled both mixes to achieve a similar workability in terms of the mini-slump value in the range of 85 ± 5 mm. The percentage of FA and GGBS in the cement blends was 35% and 50%, respectively. These are the typical cement replacement levels for FA and GGBS, respectively, in practice. Prior to mixing, the PC and FA or GGBS powders were placed into a plastic container and mixed by rotating the container by hand for three minutes in an attempt to achieve a thorough mixing of the powders. The pastes were then manufactured by mixing the water and FA/PC or GGBS/PC powder using a Hobart planetary mixer following the procedures specified in BS EN 196-3:2005+A1. To avoid any variation between batches, the specimens needed for each mix (i.e., FA/PC and GGBS/PC) were cast in one batch. The obtained pastes were grey color slurries with appropriate flowability (i.e., a mini-slump value of 85 ± 5 mm as aforementioned).
Immediately after mixing, the cement pastes were cast in two layers into cylindrical PVC molds (size of Ø50 mm × 80 mm) with each layer being vibrated for several seconds to remove the entrapped air. Following this, the specimen surfaces were finished with a steel scraper. The molds (with specimens) were then sealed with lids and cured in a constant-temperature room at 20 (± 1) °C for 24 h.

2.3. Curing, Conditioning and Carbonation

After 24 h initial curing in the molds, the specimens were demolded and wrapped individually with a water-saturated hessian. The hessian was regularly checked and wetted to ensure the samples were fully saturated. The specimens were then sealed in plastic bags and stored again in the constant-temperature room at 20 (± 1) °C for 55 days.
After a total of 56 days’ curing, the specimens were dried in an oven at 40 (± 1) °C for 14 days. The specimens were then wrapped individually with a polythene sheet and placed again in the 40 (± 1) °C oven for another 14 days to redistribute the moisture and achieve an internal relative humidity (RH) of about 60 (± 5)% for the subsequent accelerated carbonation. Following this, the specimens were cooled down to 20 (± 1) °C for one day and all faces, except the face cast against the molds, were coated with three layers of Sika acrylic paint. After being dried under 20 (± 1) °C for five days, the painted specimens were stored in a carbonation chamber (LEEC, UK) at a constant temperature of 20 (± 1) °C with a CO2 concentration of 5 (± 0.5)% and RH of 60 (± 5)% for six weeks.
After the carbonation, the powder samples were collected at five discrete depths from the exposed surface inwards (namely, 0–2 mm, 5–7 mm, 10–12 mm, 15–17 mm and 20–22 mm) of the specimens using a digital drill (with an accuracy of ± 0.1 mm) with a bit of 8 mm diameter and 4 mm deep slant edge. The collected powders were then passed through a 63 µm sieve. To avoid any contamination to the samples, a vacuum cleaner was employed to suck away the powders from the intervals (i.e., 2–5 mm, 7–10 mm, 12–15 mm and 17–20 mm). The obtained powder samples (i.e., the particles below 63 µm) were then put into airtight plastic bags and stored in a vacuum desiccator until being analyzed by Raman and TG tests. Figure 2 illustrates the sample preparation procedure involved during the experiments.

2.4. Bench-Mounted Raman Spectroscopy

A bench-mounted Renishaw Raman microscope with a charged coupled device (CCD) (Renishaw, Gloucestershire, UK) was used in this study, which worked under a controlled temperature of 20 °C. The laser was a 514.5 nm single-line (Argon ions) laser and its beam was focused into the sample through a LEICA (Wetzlar, Germany) N Plan objective (10×). The measured power at the sampling level was around 3.2 mW. Before each experiment, the Raman shift was calibrated using the sharp peak of the TiO2 (anatase) powder. The analyses were performed with an exposure time of 10 s and accumulations of 10 to achieve an appropriate signal-to-noise ratio (SNR).

2.5. Background Subtraction

During the data process, the background of the Raman spectra was subtracted by adopting the baseline correction using Origin 2018 (OriginLab, Northampton, MA, USA). This procedure can remove the strong fluorescence background, allowing the Raman peaks to be clearly observed.

2.6. Calculation of the Signal-to-Noise Ratio (SNR)

The SNR values of the Raman spectra were calculated using the υ1 CO3 peak of the samples obtained at the depth of 0–2 mm according to the method recommended in American Society for Testing Materials (ASTM ) E579-04 as follows:
SNR = Signal level/Noise (RMS) level
where the signal level is the peak height/intensity after subtracting the background and the noise level is obtained by the root mean square (RMS) method, which is the standard deviation of the intensity values of a selected Raman shift region on the spectrum after subtracting the background.

2.7. Thermogravimetry (TG)

To verify the results obtained from the Raman analysis, the well-established TG analysis was employed to assess the CaCO3 quantity. A NETZSCH STA 449C instrument (NETZSCH, Selb, Germany), which worked under an inert environment (N2) with a flow rate of 1 mL/min, was used. The samples were placed in an alumina (Al2O3) crucible and heated from room temperature to 1000 °C at a heating rate of 20 °C/min. The derivative thermogravimetric curves (DTG) were also recorded simultaneously. The CaCO3 quantity was calculated by following its decomposition to CaO and CO2 at a temperature of around 470 °C~820 °C.

3. Results and Discussion

To identify the possible effects that FA or GGBS might impose on the qualitative and, particularly, the quantitative working capacity of Raman spectroscopy for establishing carbonation profiles in blended cementitious systems, the experiment work was organized into two parts and the results are reported accordingly in this section as follows:
  • Establish the carbonation profiles viz. the variation of CaCO3 content against the depth of FA or GGBS containing PC pastes after being subjected to an accelerated carbonation with bench-mounted Raman spectroscopy.
  • Verify the carbonation profiles as developed by Raman spectroscopy at stage 1 by the well-established TG analysis so that the potential effects that FA or GGBS may have on the reliability of the results obtained from Raman spectroscopy could be recognized.

3.1. Establishing the Carbonation Profiles of FA/PC and GGBS/PC Pastes with Bench-Mounted Raman Spectroscopy

3.1.1. Raman Spectroscopy Analysis

As aforementioned, FA is an intrinsically weak Raman scatterer whilst GGBS is also difficult to analyze by Raman spectroscopy due to its highly disordered amorphous phases. As exemplified in Figure 3 by the Raman spectra of the carbonated FA/PC, GGBS/PC and pure PC pastes at the depth of 0–2 mm, which all clearly showed an intense sharp peak at about 1085 cm−1 due to the υ1 (i.e., symmetrical stretching) vibration of CO3 of the CaCO3 phases [28,29,30,31], the intensity of the υ1 CO3 band of the FA/PC blend was much weaker (only 51% of that of the pure PC paste, as shown in Table 2) than those of both GGBS/PC and pure PC. In addition to the weak Raman scatter nature of FA, the weak intensity from the FA/PC blend could also be attributed to the dark color of FA (as shown and compared with PC and GGBS in Figure 1). This is because, as well-established, the phases with the dark color can absorb more illuminating light energy [32]. As a result, the scattered light intensity could be reduced due to the reduced incident light energy. This could partly explain the much reduced Raman intensity of the FA/PC blend in Figure 3. Although the intensity of the υ1 CO3 band of GGBS/PC was not reduced (104% of PC) compared with that of PC, a stronger fluorescence background could be recognized. As the peak height is normally used to quantify the phases in a Raman analysis (as detailed in Section 3.1.2), the interferences caused by FA to the intensity and by GGBS to the fluorescence background of the Raman spectra in Figure 3 could cause concerns about the quantification of CaCO3 formed in the different mixes before establishing the carbonation profiles. To quantitatively identify the possible effects that FA and GGBS could have imposed on the quality of the Raman spectra, the SNR values of the Raman spectra in Figure 3 were calculated using the υ1 CO3 band according to the method described in Section 2.6 and the results are reported in Table 2. As can be seen from Table 2, under the Raman spectroscopy analysis, the SNR of the spectra obtained from the GGBS/PC paste was lower than that of PC, i.e., 75 vs. 81, while an even lower SNR (47) was obtained from the FA/PC paste, which could be primarily due to the much reduced signal level (51% of PC) of the FA/PC paste. On the other hand, the relatively lower SNR from GGBS/PC compared with the pure PC mix could be mainly caused by the higher noise level of GGBS (113% of PC), which is presumably due to the highly disordered amorphous phases in GGBS. Nonetheless, a strong background could be identified from all the three Raman spectra and in particular the GGBS/PC blend and the pure PC samples. As both PC and GGBS were subjected to a grinding treatment during their manufacture, the higher background from PC and GGBS could also be caused by the grinding agent as well as the grinding itself because grinding could also introduce different defects, promoting the formation of a strong fluorescence background [33]. Therefore, before CaCO3 could be quantified, the interference from the strong fluorescence background has to be removed first. This could also facilitate the identification of the possible effects that FA and GGBS could have introduced to the carbonation profile to be established in the next section. The background of the spectra collected at the different depths was therefore subtracted by following the procedures described in Section 2.5 and the details are presented below.
Figure 4 and Figure 5 present the Raman spectra after subtracting the background from five different depths of the carbonated FA/PC and GGBS/PC paste cylinders, respectively. As shown in Figure 4, for the FA/PC blend, there were sharp peaks located at about 1085 cm−1, 1075 cm−1 and 1091 cm−1 (i.e., υ1 symmetrical stretching bands of CO3) being identified at 0–2 mm, 5–7 mm, 10–12 mm and 15–17 mm, indicating the formation of different CaCO3 polymorphs (viz. calcite, aragonite, vaterite) at these depths [28,29,30,31]. As highlighted before, the co-existence of these three CaCO3 polymorphs could be possible in the current study due to the severe carbonation condition applied, i.e., 5% CO2, during the accelerated carbonation test. In addition, the in-plane bending mode of the CO34 CO3) band of calcite/aragonite was also identified at the two outermost depths, viz. 0–2 mm and 5–7 mm, as could be observed at about 707 cm−1 [28,29,30,31]. The weak band that emerged at about 1349 cm−1 could not be assigned appropriately, which could be attributed to a few impurities in the raw cement. In Figure 5, the υ1 CO3 peaks of the CaCO3 polymorphs at the first four depths were also clearly identified by Raman spectroscopy in the GGBS/PC blend as shown at 1085 cm−1 (calcite/aragonite) and 1075 cm−1/1090 cm−1 (vaterite). Additionally, the υ4 CO3 band of calcite/aragonite was identified at the outmost two depths (i.e., 0–2 mm and 5–7 mm), as indicated at about 707 cm−1 [28,29,30,31]. Most importantly, as clearly shown in both Figure 4 and Figure 5, the intensity of the υ1 CO3 band decreased with the increase of the depth, indicating the reduced amount of CaCO3 formed over the depth into the hardened PC blends of both systems, which could be mainly attributed to the reduced penetration of the CO2 into the FA and GGBS bearing PC pastes, respectively, under the accelerated carbonation.

3.1.2. Establishing the Carbonation Profiles

In the previous section, a decreasing trend of the υ1 CO3 band over the depth in both the FA/PC and GGBS/PC pastes was clearly recognized. To identify the feasibility of establishing the carbonation profiles (in terms of the variation of the CaCO3 content against depth) using Raman spectroscopy, the peak heights (acting as the indicator of CaCO3 quantity) of the υ1 CO3 band at the five depths of the FA/PC and GGBS/PC pastes were retrieved and are presented in Table 3. In addition, the peak heights of the pure PC samples are also presented for comparison purposes [26]. In Table 3, it can be noticed that the carbonation fronts of the FA and GGBS containing pastes were deeper (until 15–17 mm) than that of the PC system (until 10–12 mm), as no CaCO3 was identified at the depth of 15–17 mm onwards in the PC pastes. This phenomenon is consistent with the literature showing that the addition of the mineral additions such as FA and GGBS could reduce the resistance to carbonation attack [34,35]. Additionally, for all of the three systems, CaCO3 was not identified at the depth of 20–22 mm, indicating these most inner regions might have not been carbonated.
To develop a more visual illustration of the carbonation profiles, the peak height was then plotted against the depth as shown in Figure 6. Evidently, for both the FA/PC and GGBS/PC blends, there was an inverse relationship between the height (intensity) of the υ1 CO3 peak and the depth, indicating the decrease of the CaCO3 content with the increase of the sampling depth. Obviously, these two profiles displayed very similar patterns to that of the PC paste, i.e., the peak height of all of the three cementitious materials, viz. FA/PC, GGBS/PC and pure PC pastes, showed a declining trend with the increasing depth. This phenomenon correlated well with the typical carbonation profile, i.e., the CaCO3 content decreased with the depth. However, the peak intensities of the FA/PC blends at the depths of 0–2 mm and 5–7 mm were much lower than that of the GGBS/PC and PC pastes. Although a lower CaCO3 content from the FA/PC mix was somehow expected due to its reduced calcium bearing phases, the reduced incident laser power for the Raman scattering of the FA/PC mix owing to the darker color of FA could also have contributed to this phenomenon, which will be further discussed in Section 3.2.2 below.
Therefore, based on the above results and discussion, it can be seen that introducing FA/GGBS into the pure PC mix has indeed caused some effects on the Raman spectra, which could be mainly attributed to the dark color and weak scatter nature of FA and the highly disordered phases in GGBS. However, although there were differences in the background, noise level and signal intensity between the Raman spectra obtained from the different cementitious systems, the qualitative working capacity of Raman spectroscopy for recognizing different CaCO3 polymorphs was not affected. Nonetheless, due to the differences in the peak heights between the different cementitious systems, the quantitative working capacity of Raman spectroscopy might have been affected and this was to be identified through the TG verification and is described in the next section below.

3.2. Verifying the Carbonation Profiles by TG Analysis

In the above section, using the peak height of the υ1 CO3 band as the indicator of CaCO3 content, the carbonation profiles of the FA/PC and GGBS/PC blends were established by Raman spectroscopy. However, due to the effects of FA and GGBS, the reliability of the quantified information as obtained from the Raman spectroscopy technique was still uncertain. In this section, derivative TG analysis (i.e., DTG) was employed as a well-established quantitative technique to verify the quantitative information and, hence, the carbonation profiles established by Raman spectroscopy.

3.2.1. Thermogravimetry (TG) Analysis

To clearly identify each individual phase, the DTG curves of the FA and GGBS blended cement pastes are presented in Figure 7 and Figure 8, respectively. As shown in Figure 7, the decomposition of CaCO3 in the FA/PC paste occurred at the temperature range about 470 °C~820 °C, with a peak centered at about 780 °C. Based on the literature, this could suggest the co-existence of the three CaCO3 polymorphs; namely, calcite, aragonite and vaterite [5]. An endothermic peak emerged at about 475 °C, which could be due to the dehydration of the portlandite in the cementitious materials. However, this peak disappeared from the first three depths of the FA blended paste, indicating that there was almost no portlandite existing in the first three depths of the FA blends, which could have been consumed by the reactions with CO2 (carbonation) or siliceous/alumino-siliceous in FA (pozzolanic reaction). On the other hand, in Figure 8, the endothermic peak for CaCO3 in the GGBS/PC samples emerged at about 760 °C with a related decomposition temperature ranging from around 470 °C to 810 °C, which corroborated the sharp mass loss due to the decomposition of the CaCO3 phases. Meanwhile, the decomposition temperature for portlandite was also observed at about 420 °C~470 °C with its peak centered at about 450 °C. These peaks became more evident in the deeper depths, i.e., 10–12 mm, 15–17 mm and 20–22 mm, indicating that a considerable amount of portlandite still existed in these areas due to reduced carbonation.
Moreover, in both Figure 7 and Figure 8 it is apparent that the integrated area of CaCO3 during the temperature range, which was linked to the quantity of CaCO3, decreased with the sampling depth, implying the possible decrease of the CaCO3 content. To exactly quantify the variation of the CaCO3 content, the area of the DTG curve at the decomposition temperature range of CaCO3 was determined accordingly and then used to calculate the content of CaCO3. Table 4 reports the obtained CaCO3 content and those of the PC paste samples are again listed for comparison purposes [26]. As can be seen from Table 4, considerable amounts (>8%) of CaCO3 were identified at the first four depths of the FA blended paste cylinders while only small quantity of CaCO3 was identified at the depth of 20–22 mm. For the GGBS/PC blends, the CaCO3 contents within the outermost three depths (0–2 mm, 5–7 mm and 10–12 mm) were relatively high (>9%), but decreased from a depth of 15–17 mm. These results, overall, correlate well with the Raman spectroscopy results (Table 3) in which CaCO3 was characterized at the first four depths of the samples. However, as small amounts of CaCO3 were identified by TG at the innermost depth (20–22 mm) of the carbonated paste cylinders, but not detected by Raman spectroscopy, it seems that further research is still needed to investigate how to improve the sensitivity of Raman spectroscopy.

3.2.2. Verifying the Carbonation Profiles

Based on the CaCO3 content obtained from the TG analysis (shown in Table 4), the carbonation profiles of the FA/PC and GGBS/PC pastes established by Raman spectroscopy using the peak height of υ1 CO3 as an indicator of the CaCO3 content were plotted together with the counterpart obtained from the TG analysis as seen in Figure 9 and Figure 10 respectively. Apparently, for both the FA and GGBS blended cement pastes, the carbonation profiles obtained from these two techniques (i.e., Raman spectroscopy and TG), both of which showed a reducing trend, were highly comparable. As TG is a well-established technique for quantifying CaCO3 in cementitious materials, the above results demonstrated a reasonable working capacity of the Raman spectroscopy for differentiating the variation (i.e., the relative quantity) of the CaCO3 content in carbonated FA/PC or GGBS/PC blends over the depth. However, from Figure 9 and Figure 10, it was not clear whether the CaCO3 content established from Raman spectroscopy could be used as an indicator of the real quantity (i.e., the absolute quantity) of CaCO3. This will be discussed in the sections below.
To quantitatively and, hence, fully verify the working capacity of Raman spectroscopy for developing carbonation profiles, the correlations between the Raman and TG analyses were established and compared in Figure 11. As can be seen from Figure 11, for the FA/PC and GGBS/PC blends, there were reasonably good correlations between the contents of CaCO3 obtained from TG and the peak of the Raman spectroscopy analysis as demonstrated by their correlation coefficient (R2) values, i.e., 0.88 for FA/PC and 0.91 for the GGBS/PC system. Nonetheless, there was no doubt that the addition of FA and GGBS into the PC system clearly reduced these correlations as exemplified by the relatively lower R2 values of the FA/PC (0.88) and the GGBS/PC (0.91) blends compared with that of the pure PC (0.97). This indicated that FA and GGBS could cause some adverse effects on the quantification of CaCO3 by Raman spectroscopy with FA demonstrating more negative effects than GGBS. As previously discussed, the relatively darker color of FA could absorb more energy from the incident laser light leading to reduced Raman scattering and, hence, a reduced peak intensity. This, along with the weak Raman scatter nature of FA, could potentially reduce the intensity of the Raman spectrum even if the same amount of CaCO3 might have existed (as clearly demonstrated by the lower peak intensity reported in Table 2). Additionally, it should be highlighted that, as clearly demonstrated in Figure 1, the darkness of the raw materials increased in the order GGBS < PC < FA. Accordingly, one may expect that the consumption of light energy by the raw materials should also be in the same order of GGBS < PC < FA, which clearly coincided well with the peak intensities in Figure 11 even though the exact nature that caused the difference in intensity still needs further research. Additionally, the reduced correlation from the GGBS/PC blend could also be attributed to the possible interference from the increased noise level. As a result, as demonstrated in Figure 11, if the peak intensity was used to compare the quantity of CaCO3 across different cementitious systems, there was a potential risk that CaCO3 in the FA/PC blend could be underestimated whilst that in the GGBS/PC could be overestimated. The above results would, therefore, indicate that simply using the peak intensity as a quantity indicator might not be able to provide reliable information on the quantity of CaCO3 due to the effects of FA and GGBS. To overcome these effects, a quantified Raman analysis with an internal standard [36,37] similar to quantitative X-ray diffraction (QXRD) might need to be considered in future studies. Alternatively, if the Raman spectroscopy system is to be developed into a future monitoring system for the continued monitoring of the durability of concrete structures, individual calibration curves could be established for each type of concrete.

4. Conclusions

Carbonation initiates the corrosion of reinforcing steel and alters the mineralogical composition as well as the microstructure of the cementitious materials. Thus, characterizing the CaCO3 polymorphs and in particular establishing the carbonation profiles are essential for predicting the service life and the health condition of concrete structures. The current work employed Raman spectroscopy to establish the carbonation profiles (in terms of the variation of CaCO3 content against the depth) of FA and GGBS blended cementitious materials. Additionally, the potential effects that FA and GGBS might impose on the quantitative working capacity of Raman spectroscopy was systematically investigated. From the results obtained in this study, it showed that Raman spectroscopy was adequate to differentiate the variation (i.e., the relative quantity) of the CaCO3 content and hence there is a potential that the carbonation profiles of the FA/PC or GGBS/PC blended cementitious materials can be established with the Raman spectroscopy technique. However, the addition of FA and GGBS clearly reduced the correlations between the contents of CaCO3 obtained from TG and the peak intensity of the Raman spectroscopy analysis as exemplified by the lower R2 values obtained. Therefore, using the peak intensity as the indicator of the CaCO3 content, the CaCO3 content in the FA/PC blend could be underestimated whilst that in the GGBS/PC could be overestimated. These adverse effects could be due to the darker color and weak scatter nature of FA and the high content of glassy phases in GGBS, respectively. To solve this problem, a quantified Raman analysis with an internal standard or individual calibration curve has to be employed before this technique can be used as a reliable alternative approach to monitoring the carbonation attack in concrete structures.

Author Contributions

Methodology, Y.Y., J.J.W. and Y.B.; validation, Y.Y. and Y.B.; formal analysis, Y.Y.; investigation, Y.Y.; resources, J.J.W., P.A.M.B. and Y.B.; data curation, Y.Y.; writing—original draft preparation, Y.Y.; writing—review and editing, Y.Y. and Y.B.; supervision, J.J.W., P.A.M.B. and Y.B.; project administration, Y.Y., J.J.W. and Y.B.; funding acquisition, Y.Y., J.J.W. and Y.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant number 51802033), the State Key Laboratory of High Performance Civil Engineering Materials (Grant number 2019CEM009), the Fundamental Research Funds for the Central Universities (Grant number 2019CDXYCL0031), the European Community Research Infrastructures under the FP7 Capacities Programme (the QualityNano Transnational Access Project, Grant number INFRA-2010-262163) and the UK Royal Society (the International Exchange Scheme with Ireland, Grant number IE131481).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Acknowledgments

The authors would like to acknowledge Quinn and Sika for providing the materials for this research.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Neville, A.M. Properties of Concrete, 5th ed.; Prentice Hall: San Francisco, CA, USA, 1995. [Google Scholar]
  2. Šavija, B.; Luković, M. Carbonation of cement paste: Understanding, challenges, and opportunities. Construct. Build. Mater. 2016, 117, 285–301. [Google Scholar] [CrossRef] [Green Version]
  3. Broomfield, J.P. Carbonation and its effects in reinforced concrete. Mater. Perform. 2000, 39, 64–67. [Google Scholar]
  4. Sevelsted, T.F.; Skibsted, J. Carbonation of C–S–H and C–A–S–H samples studied by 13C, 27Al and 29Si MAS NMR spectroscopy. Cem. Concr. Res. 2015, 71, 56–65. [Google Scholar] [CrossRef]
  5. Villain, G.; Thiery, M.; Platret, G. Measurement methods of carbonation profiles in concrete: Thermogravimetry, chemical analysis and gammadensimetry. Cem. Concr. Res. 2007, 37, 1182–1192. [Google Scholar] [CrossRef]
  6. Morandeau, A.; Thiéry, M.; Dangla, P. Investigation of the carbonation mechanism of CH and C-S-H in terms of kinetics, microstructure changes and moisture properties. Cem. Concr. Res. 2014, 56, 153–170. [Google Scholar] [CrossRef] [Green Version]
  7. Liu, W.; Li, Y.Q.; Tang, L.P.; Dong, Z.J. XRD and 29Si MAS NMR study on carbonated cement paste under accelerated carbonation using different concentration of CO2. Mater. Today Commun. 2019, 19, 464–470. [Google Scholar] [CrossRef]
  8. Ruggieri, S.; Porco, F.; Uva, G.; Vamvatsikos, D. Two frugal options to assess class fragility and seismic safety for low-rise reinforced concrete school buildings in Southern Italy. Bull. Earthq. Eng. 2021, 19, 1415–1439. [Google Scholar] [CrossRef]
  9. Porco, F.; Ruggieri, S.; Uva, G. Seismic assessment of irregular existing building: Appraisal of the influence of compressive strength variation by means of nonlinear conventional and multimodal static analysis. Ing. Sismica 2018, 35. [Google Scholar]
  10. Houst, Y.F.; Wittmann, F.H. Depth profiles of carbonates formed during natural carbonation. Cem. Concr. Res. 2002, 32, 1923–1930. [Google Scholar] [CrossRef]
  11. Ann, K.Y.; Pack, S.W.; Hwang, J.P.; Song, H.W.; Kim, S.H. Service life prediction of a concrete bridge structure subjected to carbonation. Constr. Build. Mater. 2010, 24, 1494–1501. [Google Scholar] [CrossRef]
  12. Johannesson, B.; Utgenannt, P. Microstructural changes caused by carbonation of cement mortar. Cem. Concr. Res. 2001, 31, 925–931. [Google Scholar] [CrossRef]
  13. Taylor, H.F.W. Cement Chemistry, 2nd ed.; Thomas Telford: London, UK, 1997. [Google Scholar]
  14. Lee, H.J.; Kim, D.G.; Lee, J.H.; Cho, M.S. A study for carbonation degree on concrete using a phenolphthalein indicator and fourier-transform infrared spectroscopy. World Acad. Sci. Eng. Technol. 2012, 62, 184–190. [Google Scholar]
  15. Lo, Y.; Lee, H.M. Curing effects on carbonation of concrete using a phenolphthalein indicator and Fourier-transform infrared spectroscopy. Build. Environ. 2002, 37, 507–514. [Google Scholar] [CrossRef]
  16. Thiery, M.; Villain, G.; Dangla, P.; Platret, G. Investigation of the carbonation front shape on cementitious materials: Effects of the chemical kinetics. Cem. Concr. Res. 2007, 37, 1047–1058. [Google Scholar] [CrossRef]
  17. Chang, C.-F.; Chen, J.-W. The experimental investigation of concrete carbonation depth. Cem. Concr. Res. 2006, 36, 1760–1767. [Google Scholar] [CrossRef]
  18. Long, D.A. Raman Spectroscopy; McGraw-Hill International Book Company: New York, NY, USA, 1977. [Google Scholar]
  19. McCreery, R.L. Raman Spectroscopy for Chemical Analysis; John Wiley & Sons Inc.: Hoboken, NJ, USA, 2000. [Google Scholar]
  20. Predeanu, G.; Popescu, L.G.; Abagiu, T.A.; Panaitescu, C.; Valentim, B.; Guedes, A. Characterization of bottom ash of Pliocene lignite as ceramic composites raw material by petrographic, SEM/EDS and Raman microspectroscopical methods. Int. J. Coal Geol. 2016, 168, 131–145. [Google Scholar] [CrossRef]
  21. Corvisier, J.; Brunet, F.; Fabbri, A.; Bernard, S.; Findling, N.; Rimmelé, G.; Barlet-Gouédard, V.; Beyssac, O.; Goffé, B. Raman mapping and numerical simulation of calcium carbonates distribution in experimentally carbonated Portland-cement cores. Eur. J. Miner. 2010, 22, 63–74. [Google Scholar] [CrossRef] [Green Version]
  22. Yue, Y.; Wang, J.J.; Basheer, P.A.M.; Bai, Y. Raman spectroscopic investigation of Friedel’s salt. Cem. Concr. Compos. 2018, 86, 306–314. [Google Scholar] [CrossRef]
  23. Yue, Y.; Wang, J.J.; Basheer, P.A.M.; Boland, J.J.; Bai, Y. Characterisation of carbonated Portland cement paste with optical fibre excitation Raman spectroscopy. Constr. Build. Mater. 2017, 135, 369–376. [Google Scholar] [CrossRef]
  24. Yue, Y.; Wang, J.J.; Bai, Y. Tracing the status of silica fume in cementitious materials with Raman microscope. Constr. Build. Mater. 2018, 159, 610–616. [Google Scholar] [CrossRef]
  25. Yue, Y.; Bai, Y.; Basheer, P.A.M.; Boland, J.J.; Wang, J.J. Monitoring the cementitious materials subjected to sulfate attack with optical fiber excitation Raman spectroscopy. Opt. Eng. 2013, 52, 104107. [Google Scholar] [CrossRef]
  26. Yue, Y.; Wang, J.J.; Basheer, P.A.M.; Boland, J.J.; Bai, Y. A Raman spectroscopy based optical fibre system for detecting carbonation profile of cementitious materials. Sens. Actuator. B Chem. 2018, 257, 635–649. [Google Scholar] [CrossRef] [Green Version]
  27. Potgieter-Vermaak, S.S.; Potgieter, J.H.; Belleil, M.; DeWeerdt, F.; Van Grieken, R. The application of Raman spectrometry to the investigation of cement: Part II: A micro-Raman study of OPC, slag and fly ash. Cem. Concr. Res. 2006, 36, 663–670. [Google Scholar] [CrossRef]
  28. Martinez-Ramirez, S.; Sanchez-Cortes, S.; Garcia-Ramos, J.V.; Domingo, C.; Fortes, C.; Blanco-Varela, M.T. Micro-Raman spectroscopy applied to depth profiles of carbonates formed in lime mortar. Cem. Concr. Res. 2003, 33, 2063–2068. [Google Scholar] [CrossRef]
  29. Black, L.; Breen, C.; Yarwood, J.; Garbev, K.; Gasharova, B. Structural Features of C–S–H(I) and Its Carbonation in Air—A Raman Spectroscopic Study. Part II: Carbonated Phases. J. Am. Ceram. Soc. 2007, 90, 908–917. [Google Scholar] [CrossRef]
  30. Tlili, M.M.; Ben Amor, M.; 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]
  31. Ibanez, J.; Artus, L.; Cusco, R.; Lopez, A.; Menendez, E.; Andrade, M.C. Hydration and carbonation of monoclinic C2S and C3S studied by Raman spectroscopy. J. Raman Spectrosc. 2007, 38, 61–67. [Google Scholar] [CrossRef]
  32. Caggiani, M.C.; Colomban, P. Raman identification of strongly absorbing phases: The ceramic black pigments. J. Raman Spectrosc. 2011, 42, 839–843. [Google Scholar] [CrossRef]
  33. Richardson, I.G.; Skibsted, J.; Black, L.; Kirkpatrick, R.J. Characterisation of cement hydrate phases by TEM, NMR and Raman spectroscopy. Adv. Cem. Res. 2010, 22, 233–248. [Google Scholar] [CrossRef]
  34. Hussain, S.; Bhunia, D.; Singh, S.B. Comparative study of accelerated carbonation of plain cement and fly-ash concrete. J. Build. Eng. 2017, 10, 26–31. [Google Scholar] [CrossRef]
  35. Han-Seung, L.; Wang, X.-Y. Evaluation of compressive strength development and carbonation depth of high volume slag-blended concrete. Constr. Build. Mater. 2016, 124, 45–54. [Google Scholar] [CrossRef]
  36. Sun, Q.; Qin, C. Raman OH stretching band of water as an internal standard to determine carbonate concentrations. Chem. Geol. 2011, 283, 274–278. [Google Scholar] [CrossRef]
  37. Moreno, T.; Morán López, M.A.; Huerta Illera, I.; Piqueras, C.M.; Sanz Arranz, A.; García Serna, J.; Cocero, M.J. Quantitative Raman determination of hydrogen peroxide using the solvent as internal standard: Online application in the direct synthesis of hydrogen peroxide. Chem. Eng. J. 2011, 166, 1061–1065. [Google Scholar] [CrossRef]
Figure 1. Photo of fly ash (FA), Portland cement (PC) and ground granulated blast furnace slag (GGBS).
Figure 1. Photo of fly ash (FA), Portland cement (PC) and ground granulated blast furnace slag (GGBS).
Materials 14 01798 g001
Figure 2. Flowchart illustrating the sample preparation procedure.
Figure 2. Flowchart illustrating the sample preparation procedure.
Materials 14 01798 g002
Figure 3. Raman spectra of the powders collected at a depth of 0–2 mm of the carbonated FA/PC, GGBS/PC and PC paste samples.
Figure 3. Raman spectra of the powders collected at a depth of 0–2 mm of the carbonated FA/PC, GGBS/PC and PC paste samples.
Materials 14 01798 g003
Figure 4. Raman spectra (after subtracting the background) of the powders collected at five depths of carbonated FA/PC paste samples.
Figure 4. Raman spectra (after subtracting the background) of the powders collected at five depths of carbonated FA/PC paste samples.
Materials 14 01798 g004
Figure 5. Raman spectra (after the subtracting background) of the powders collected at five depths of carbonated GGBS/PC paste samples.
Figure 5. Raman spectra (after the subtracting background) of the powders collected at five depths of carbonated GGBS/PC paste samples.
Materials 14 01798 g005
Figure 6. Relationship between the depth and peak height obtained from carbonated FA/PC, GGBS/PC and PC paste samples under Raman spectroscopy.
Figure 6. Relationship between the depth and peak height obtained from carbonated FA/PC, GGBS/PC and PC paste samples under Raman spectroscopy.
Materials 14 01798 g006
Figure 7. Derivative thermogravimetric (DTG) curves of the carbonated FA/PC paste samples at five depths.
Figure 7. Derivative thermogravimetric (DTG) curves of the carbonated FA/PC paste samples at five depths.
Materials 14 01798 g007
Figure 8. DTG curves of the carbonated GGBS/PC paste samples at five depths.
Figure 8. DTG curves of the carbonated GGBS/PC paste samples at five depths.
Materials 14 01798 g008
Figure 9. Carbonation profiles of the carbonated FA/PC paste samples established by Raman spectroscopy and TG analyses.
Figure 9. Carbonation profiles of the carbonated FA/PC paste samples established by Raman spectroscopy and TG analyses.
Materials 14 01798 g009
Figure 10. Carbonation profiles of the carbonated GGBS/PC paste samples established by Raman spectroscopy and TG analyses.
Figure 10. Carbonation profiles of the carbonated GGBS/PC paste samples established by Raman spectroscopy and TG analyses.
Materials 14 01798 g010
Figure 11. Correlation between the calcium carbonate in the carbonated FA/PC, GGBS/PC and PC paste samples identified by Raman spectroscopy and TG analyses.
Figure 11. Correlation between the calcium carbonate in the carbonated FA/PC, GGBS/PC and PC paste samples identified by Raman spectroscopy and TG analyses.
Materials 14 01798 g011
Table 1. Chemical composition of the PC, FA and GGBS.
Table 1. Chemical composition of the PC, FA and GGBS.
Oxides/%SiO2Al2O3Fe2O3CaOMgOK2ONa2OSO3
PC23.006.152.9561.301.800.680.222.50
FA56.3023.504.704.401.801.801.000.90
GGBS34.3015.000.4039.408.000.380.450.80
Table 2. Summary of the signal, noise and signal-to-noise ratio (SNR) of the Raman spectra of the powders collected at a depth of 0–2 mm of the carbonated FA/PC, GGBS/PC and PC paste samples.
Table 2. Summary of the signal, noise and signal-to-noise ratio (SNR) of the Raman spectra of the powders collected at a depth of 0–2 mm of the carbonated FA/PC, GGBS/PC and PC paste samples.
Parameters of SpectraFA/PC PasteGGBS/PC PastePC Paste
Signal (peak intensity)5151.5
(51%)
10,491.3
(104%)
10,055.1
(100%)
Noise110.0
(89%)
139.8
(113%)
124.2
(100%)
SNR477581
Note: the signal and noise were normalized based on that of the PC paste and shown in brackets.
Table 3. Summary of the calcium carbonate content (in terms of peak height) at five depths of carbonated FA/PC, GGBS/PC and PC paste samples under Raman spectroscopy.
Table 3. Summary of the calcium carbonate content (in terms of peak height) at five depths of carbonated FA/PC, GGBS/PC and PC paste samples under Raman spectroscopy.
Depth/mmFA/PC PasteGGBS/PC PastePC Paste
0–25151.510491.310,055.1
5–73461.59734.59741.7
10–122993.94776.63014.8
15–172143.72650.00
20–22000
Table 4. Summary of the calcium carbonate content (%) at five depths of carbonated FA/PC, GGBS/PC and PC paste samples under a thermogravimetry (TG) analysis.
Table 4. Summary of the calcium carbonate content (%) at five depths of carbonated FA/PC, GGBS/PC and PC paste samples under a thermogravimetry (TG) analysis.
Depth/mmFA/PC PasteGGBS/PC PastePC Paste
0–218.411.925.2
5–714.411.821.4
10–1214.19.38.7
15–178.25.75.5
20–226.25.55.1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yue, Y.; Wang, J.J.; Basheer, P.A.M.; Bai, Y. Establishing the Carbonation Profile with Raman Spectroscopy: Effects of Fly Ash and Ground Granulated Blast Furnace Slag. Materials 2021, 14, 1798. https://0-doi-org.brum.beds.ac.uk/10.3390/ma14071798

AMA Style

Yue Y, Wang JJ, Basheer PAM, Bai Y. Establishing the Carbonation Profile with Raman Spectroscopy: Effects of Fly Ash and Ground Granulated Blast Furnace Slag. Materials. 2021; 14(7):1798. https://0-doi-org.brum.beds.ac.uk/10.3390/ma14071798

Chicago/Turabian Style

Yue, Yanfei, Jing Jing Wang, P. A. Muhammed Basheer, and Yun Bai. 2021. "Establishing the Carbonation Profile with Raman Spectroscopy: Effects of Fly Ash and Ground Granulated Blast Furnace Slag" Materials 14, no. 7: 1798. https://0-doi-org.brum.beds.ac.uk/10.3390/ma14071798

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