Next Article in Journal
Ionizing Radiation Sensing with Functionalized and Copper-Coated SWCNT/PMMA Thin Film Nanocomposites
Previous Article in Journal
Synthesis and Characterization of Environmentally Friendly Chitosan–Arabic Gum Nanoparticles for Encapsulation of Oregano Essential Oil in Pickering Emulsion
Previous Article in Special Issue
Hybrid Porous Nanomaterials for Energy and Environment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Electrochemical Impedance as an Assessment Tool for the Investigation of the Physical and Mechanical Properties of Graphene-Based Cementitious Nanocomposites

by
Eirini Tziviloglou
1,
Zoi S. Metaxa
2,
George Maistros
3,
Stavros K. Kourkoulis
4,
Dionysios S. Karousos
1,5,
Evangelos P. Favvas
5 and
Nikolaos D. Alexopoulos
1,*
1
Research Unit of Advanced Materials, Department of Financial and Management Engineering, University of the Aegean, 82132 Chios, Greece
2
Department of Chemistry, International Hellenic University, 65404 Kavala, Greece
3
ADVISE, 17 Gymnasiarchou Madia Str., 82132 Chios, Greece
4
Department of Mechanics, National Technical University of Athens, 15780 Athens, Greece
5
Institute of Nanoscience and Nanotechnology, N.C.S.R. “Demokritos”, Patr. Gregoriou E & 27 Neapoleos Str., 15341 Agia Paraskevi, Greece
*
Author to whom correspondence should be addressed.
Nanomaterials 2023, 13(19), 2652; https://doi.org/10.3390/nano13192652
Submission received: 1 September 2023 / Revised: 19 September 2023 / Accepted: 24 September 2023 / Published: 27 September 2023
(This article belongs to the Special Issue Hybrid Porous Nanomaterials for Energy and Environment)

Abstract

:
This investigation explores the potential of electrochemical impedance spectroscopy (EIS) in evaluating graphene-based cementitious nanocomposites, focusing on their physical and structural properties, i.e., electrical resistivity, porosity, and fracture toughness. EIS was employed to study cement mixtures with varying graphene nanoplatelet (xGnP) concentrations (0.05–0.40% per dry cement weight), whereas flexural tests assessed fracture toughness and porosimetry analyses investigated the structural characteristics. The research demonstrated that the electrical resistivity initially decreased with increasing xGnP content, leveling off at higher concentrations. The inclusion of xGnPs correlated with an increase in the total porosity of the cement mixtures, which was indicated by both EIS and porosimetry measurements. Finally, a linear correlation emerged between fracture toughness and electrical resistivity, contributing also to underscore the use of EIS as a potent non-destructive tool for evaluating the physical and mechanical properties of conductive nano-reinforced cementitious nanocomposites.

1. Introduction

Cementitious composites have found extensive utility in numerous applications, primarily owing to their ready availability, cost effectiveness, straightforward production processes, and impressive compressive strength characteristics. Nevertheless, their inherent brittleness has posed a significant challenge, as cementitious materials tend to exhibit low tensile strength and fracture toughness. This brittle nature often leads to undesirable consequences, including cracking that can subsequently lead to structural and material degradation. In response to this challenge, substantial scientific efforts have been devoted to mitigating the brittle behavior of cementitious materials. Researchers have explored various concepts, capitalizing on the physical properties of both natural and synthetic macro- and microfibers [1]. These innovative approaches not only seek to enhance mechanical properties but also aim to fulfill the demand for multifunctional smart cementitious materials, offering improved crack-bridging, crack-healing, and sensing capabilities [1,2].
In conjunction with the recent strides in nanomaterial science, there has been a pronounced focus on investigating the impact of incorporating carbon-based nanomaterials (CBNs) into the cementitious matrix [3,4,5]. Owing to their distinctive characteristics, CBNs have gathered attention for their ability to enhance the microstructure of cementitious composites by facilitating the creation of a denser matrix. This enhancement arises from CBNs serving as nucleation sites for the deposition and growth of cement hydration products [4]. Consequently, micro-porosity and structural flaws are typically reduced, resulting in an overall improvement in mechanical performance. Furthermore, CBNs possess noteworthy electrical properties that can be harnessed to craft cementitious nanocomposites with apparent conductivity and piezo-resistive characteristics [2,4]. These tailored electrical attributes hold the potential to furnish valuable insights into the stress–strain conditions of the material or the early detection of potential damage.
Among the various CBNs, a widely used nanofiller in cementitious composites is graphene, a two-dimensional sheet of carbon atoms arranged in a hexagonal lattice, with Young’s modulus and mechanical strength as high as ~1 TPa and ~130 GPa, respectively [6]. Graphene and its derivatives can be produced in large quantities from graphite using top-down (i.e., detachment of the graphene layers from graphite) mechanochemical processes [7,8,9,10]. Graphene nanoplatelets (GnPs) are exfoliated from pristine graphite and consist of several layers of graphene sheets with diameters of several microns and thicknesses up to 100 nm. GnPs are known to contribute to improving the mechanical, thermal and electrical properties of various materials, including cementitious composites [11]. Moreover, GnPs have a reduced production cost compared with other CBNs, such as carbon nanotubes (CNTs) [12]. Hence, their excellent mechanical and electrical properties together with their availability and low cost renders them as an optimal option for incorporation in cementitious composites for improving the mechanical response and providing self-sensing properties.
The homogeneous dispersion of GnPs within the cementitious matrix emerges as a critical concern. The tendency of GnPs to agglomerate, driven by van der Waals interactions, poses a potential detriment to cementitious composites. Such agglomeration can introduce flaws and pores, serving as points for crack initiation [13] and subsequently exerting adverse effects on both the mechanical strength and electrical conductivity. Consequently, various experimental techniques detailed in the literature have been developed to assess the dispersion of GnPs, either prior to their incorporation [14,15,16,17] or subsequent to it [15,16,17,18,19,20,21,22], into the cementitious matrix. One notably effective method, frequently reported in the literature [13,14,15,23,24], involves subjecting graphene to ultrasonication energy in the presence of superplasticizers within a water medium. This procedure ensures the proper dispersion of graphene, mitigating agglomeration and enhancing the performance of cementitious composites.
To explore and understand the ways that the cementitious matrix is affected by the incorporation of GnPs, extensive laboratory studies including mechanical and electrical tests have also been performed on hardened cementitious nanocomposites. These studies are often conducted through compressive and/or flexural strength tests [18,19,20,21,25,26] and direct current (DC) resistivity measurements [27,28,29]. Mechanical testing reveals enhanced compressive and flexural strength up to a specific GnP concentration, which can significantly vary among the different studies. The improved performance was mostly attributed to the denser microstructure of the matrix resulting from the addition of GnPs and was also supported by porosimetry measurements [15,16,22,30]. On the contrary, above a certain nanofiller content, the mechanical response was decreased, presumably due to the excessive agglomeration of the nanoparticles. Correspondingly, the electrical resistivity of the cementitious nanocomposites was also affected in line with the GnPs content. In fact, electrical resistivity, in DC measurements, progressively decreased as the GnP content increased, until a critical volume fraction, above which it remained stable or even increased. This denotes that the continuous conductive network was interrupted by the formation of agglomerates and/or flaws in the cementitious matrix.
Electrical impedance spectroscopy (EIS) has found wide-ranging applications in the realm of cementitious composites, serving as a valuable tool for monitoring various facets of the material. These applications range from studying the cement hydration process [31,32,33] to tracking the evolution of microstructure [34,35] and observing degradation phenomena [36,37,38,39]. Nevertheless, in the context of evolving cementitious nanocomposites, EIS measurements have taken on an additional role. They are employed to investigate the dispersion state of carbon-based nanomaterials (CBNs), offering crucial insights into the percolation threshold [40,41]. Furthermore, EIS proves its worth by facilitating the monitoring of the hydration process and the detection of alterations in the microstructure [42]. Remarkably, it has emerged as a non-destructive method capable of evaluating the mechanical characteristics of cement composites [42,43].
The present research examines the impact of exfoliated GnPs (xGnPs) within the hardened cementitious matrix by EIS measurements alongside flexural tests, mercury intrusion porosimetry, and liquid nitrogen adsorption porosimetry analyses. Beyond the comprehensive examination of the physical, electrical, and mechanical features following xGnP incorporation, this study aims to elucidate a noteworthy correlation between EIS results and experimentally acquired parameters, specifically fracture toughness and porosity. An innovative perspective emerges as electrical measurements demonstrate their potential as a non-destructive tool for assessing these pivotal features of cementitious composites.

2. Materials and Methods

2.1. Preparation of Cementitious Nanocomposites

Cementitious nanocomposites were fabricated using type I ordinary Portland cement (CEM I 42.5 R) and graphene nanoplatelets of grade M (xGnPs) that had an average particle size of 5 μm and a thickness ranging from approximately 6 to 8 nm. These graphene nanoplatelets were sourced from XG Sciences Inc. (East Lansing, MI, USA) in Michigan. A superplasticizer (SP) based on polycarboxylate (ViscoCrete® Ferro 1000, provided by Sika Hellas ABEE - Protomagias 15, 145 68, Kryoneri, Attica, Greece) was used for the current experimental program. The SP selection was based on the results obtained by previous research [23], as ViscoCrete® Ferro 1000 proved to be the most suitable for uniform dispersion of xGnPs in the cementitious matrix among another three types of polycarboxylate-based superplasticizers provided by the same company.
To achieve a homogenous dispersion, the xGnPs and the SP were introduced into water, constituting 6% of the total mixing water volume. Ultrasonic energy was applied using a probe ultrasonicator to facilitate the dispersion of xGnPs within the aqueous solution. After achieving the desired dispersion of xGnPs, additional water was added to attain a specified water-to-cement ratio (w/c) of 0.3. The resulting mixture was manually stirred and combined with the cement, following the procedures outlined in ASTM C305 [44]. Five (5) distinct formulations were prepared, each containing varying xGnP contents: 0.0%, 0.05%, 0.10%, 0.20%, and 0.40%, based on the dry weight of the cement and named M-0.00, M-0.05, M-0.10, M-0.20, and M-0.40 respectively; the number after “M-” identifies the xGnP percentage concentration (mixture). The detailed mixture compositions are presented in Table 1.
As depicted in Figure 1, two types of prismatic specimens, each measuring 20 × 20 × 80 mm, were manufactured for the purpose of investigating (a) the mechanical and (b) the electrical characteristics of the prepared mixtures, respectively. Notched prismatic specimens were utilized for conducting mechanical tests, employing a three-point bending setup. Meanwhile, prismatic specimens were employed for the electrical resistance tests and were equipped with two embedded steel mesh electrodes. These steel mesh electrodes were positioned at an approximate distance of 30 mm from each other. All specimens were released from their moulds 24 h post-casting and were subsequently submerged in a saturated Ca(OH)2 solution (0.82 g/L) to prevent the leaching of calcium ions from the specimens. This immersion process was maintained at room temperature (20 ± 2 °C) until the specimens reached the age of 28 days.

2.2. Experimental Setup

2.2.1. Electrochemical Impedance Spectroscopy Measurements

EIS and equivalent circuits were employed to examine the conductive mechanisms within cementitious matrices reinforced with xGnPs. These analyses were conducted on prismatic specimens containing two embedded steel mesh electrodes, as illustrated in Figure 1b. To ensure nearly complete saturation, all specimens were immersed in tap water for a duration of 7 days under ambient atmospheric conditions prior to conducting the measurements. EIS was carried out using the Dielectric Thermal Analysis System (DETA-SCOPE®, ADVISE), where an alternating voltage signal with an amplitude of 1 V was applied. The frequency of the signal was swept in the range of 1 Hz to 0.391 MHz. The magnitude of impedance and the phase angle were recorded using the data acquisition software of the system. Utilizing the data obtained from these scans, Nyquist plots were generated. These plots depicted the negative imaginary component against the real part of impedance, spanning the corresponding frequency range.
The obtained spectra were analysed using equivalent circuit (EC) models with EIS Spectrum Analyser software [45]. First proposed in the late 1980s, this is a common mechanistic approach to study the electrical characteristics of cementitious materials [46].

2.2.2. Mechanical Tests

Flexural (three-point bending) experiments were carried out on notched prismatic specimens, as depicted in Figure 2. At least three (3) specimens were investigated to obtain a reliable average value of the mechanical behaviour. Following mechanical testing, the fracture toughness (KIc) of each specimen was determined. The bending tests were conducted in the first place for material characterization, to investigate the effect of the inclusion of the xGnPs in the cement matrix. Secondly, the calculated values could also serve as indicators of the quality of the xGnP dispersion in the cement paste and might also reveal a correlation between the mechanical and electrical responses of the cementitious nanocomposites.
The tests were carried out using an MTS (MTS Systems Co., Eden Prairie, MN, USA) testing machine with a 10 kN capacity. During the testing process, the mechanical load was applied at the three points of the specimen, as indicated in Figure 2, with a constant crosshead displacement rate of 0.0015 mm/s. Simultaneously, measurements of applied load, vertical displacement, and crack mouth opening (CMOD) were recorded.
The KIc is given by Equation (1) and calculated by using Equation (2) according to ASTM E399 [47], as:
KIc = PmaxSf(a/W)/(B∙W1.5)
f(a/W) = 3(a/W)0.5∙(1.99 − (a/W)∙(1 − a/W))∙(2.15 − 3.93(a/W) + 2.7(a/W)2)∙(2 + 4∙a/W)∙(1 − a/W)−1.5
where Pmax is the ultimate load, S is the span of the specimen, B is the width of the specimen, W is the height of the specimen, and α is the depth of the notch.

2.2.3. Porosity Measurements

After the completion of flexural test series, the broken parts of each mixture were picked up randomly to measure the porosity and pore size distribution using nitrogen adsorption and mercury porosimetry.

Nitrogen Adsorption Study

Investigations were conducted using N2 adsorption isotherms at 77 K to ascertain the structural attributes of the samples, including pore size distribution, total pore volume, and specific surface area. The experiments were conducted using a Quantachrome Autosorb IQ instrument. Before performing the isotherm measurements, the samples underwent overnight outgassing at 120 °C.

Mercury Porosimetry

Mercury porosimetry measurements were performed on specimens from vacuum up to pressure of 60,000 psig. The experiments were carried out using a PoreMaster 33 instrument from Quantachrome Instruments, Boynton Beach, FL, USA, as a standard procedure for analyzing the pore structure of solids encompassing pore sizes ranging from large mesopores to macropores. The technique is precisely described by Gregg and Sing [48]. The equivalent pore radius r was computed according to the capillary pressure Washburn equation [49]:
R = −3∙(2γcosθ)/PC
where PC is the capillary pressure in Pa, γ is the interfacial tension in Nm−1 (γHg/air = 0.471 Nm−1), and θ is the contact angle (θγHg/air = 140°, as suggested by Cook and Hover [50]). The useful information obtained from the mercury porosimetry intrusion curves mainly concerned the total pore volume of the sample and the pore size distribution.

3. Results and Discussion

3.1. Electrochemical Impedance Spectroscopy

3.1.1. Nyquist Plots and Electrical Resistivity

The EIS experimental findings for cementitious composites with different xGnP concentrations can be illustrated through Nyquist plots, as displayed in Figure 3, where the real and the negative imaginary part of the impedance are presented. The impedance, which expresses the resistance of the AC to flow within the cementitious matrix, is determined as follows:
Z = Z′ + Z″j,
where Z is the complex impedance, Z′ is the real part (ohmic resistance), Z″ is the imaginary part (reactance) of impedance, and j = 1 . Z′ and Z″ have been calculated using Equations (5) and (6), respectively:
Z′ = |Z| ∙ cos(arc{Z}∙π/180)
Z″ = −|Z| ∙ sin(arc{Z}∙π/180)
In each plot, the test frequency was scanned by decreasing from 0.391 MHz to 1 Hz, with corresponding data points spreading from left to right with respect to the horizontal axis. As can be seen in Figure 3, every Nyquist plot comprises of two arcs, with the right-most (lower-frequency) arc corresponding to the geometrical capacitance formed by the two steel electrodes and the second (higher-frequency) arc being associated with the formed electrochemical double layers at the electrolyte/electrode interfaces and the electrolyte ohmic resistance in the pores of the impregnated cementitious composites. Moreover, the point where the two arcs meet (cusp point) and where the reactance, representing the energy storage (capacitive) part [41] of the material, is almost zero, can be described as the total ohmic resistance of the cell, which is the sum of the electrolyte’s ohmic resistance Rs and the charge transfer resistance Rct [46]. Based on the Nyquist plots depicted in Figure 3, it is evident that the reference mixture displays the highest total ohmic resistance, whereas the mixtures enhanced with xGnPs exhibit a noticeable decrease in resistance. It is also noticeable that the cusp points for the mixtures M-0.05 and M-0.10 tend to show comparable values of approximately 2.5 kOhm, whereas the same applies for the other two mixtures (M-0.20 and M-0.40), showing higher total resistance values of nearly 4.0 kOhm.
The electrical resistivity (ρAC) for each specimen was calculated by using the following equation:
ρAC = R∙A/l,
where R is the ohmic resistance at the cusp point, A is the cross section of the specimen, and l is the distance between the internal electrodes.
As can be derived from Figure 4, the three-sample average total electrical resistivity values, obtained from the corresponding cusp points of the Nyquist plots of the investigated mixtures, first decrease with xGnP content and then increase up to an intermediate constant value.

3.1.2. Equivalent Circuit Analysis and Results Discussion

Several EC models have been introduced in the existing literature to effectively match Nyquist plots and interpret the impedance response observed in cementitious composites reinforced with carbon-based nanomaterials (CBNs) [22,29,40,42,43]. Figure 5a presents a Randles circuit, which is often used to describe electrochemical systems in different research fields.
In this EC model, the series resistance, denoted as RS, signifies the ohmic resistance within the electrolyte solution, Rct represents the charge transfer resistance, and Q corresponds to the constant phase element’s admittance Y = 1/|Z| at ω = 1 rad/s in F∙s−1. Q is intricately connected to the double layer capacitance at the interface between the electrode and the liquid electrolyte, addressing their interactions and characteristics. W stands for the Warburg resistance caused by ion diffusion in the electrolyte. Nevertheless, this simple model was not appropriate to describe the electrochemical response of the nano-reinforced specimens of the current research; therefore, the modified Randles circuit of Figure 5b was used instead, which includes three additional elements in a second loop. The second EC model considers the overall specimen’s geometrical capacitance in terms of an xGnP-related capacitor CN in parallel connection with a cementitious matrix-related non-ideal capacitor comprising of Rint and Cint elements in series connection.
As mentioned previously, RS is used to indicate the resistance of the electrolyte (water) that is present in the pores of the cementitious matrix, as shown in Figure 6. Rs mainly depends on three parameters:
(i)
The volume of electrolyte-filled pores in the bulk material far from the electrodes. The higher the pore volume, the easier the passage of ions from electrode to electrode and the smaller the series (electrolyte) resistance [29].
(ii)
The concentration of dissolved ions in the electrolyte, which might be enhanced by oxygen groups at the filler surface; this can be ignored here since the xGnPs used had no such groups.
(iii)
Partial percolation of conductive filler xGnPs, forming electron pathways that extend from the electrodes’ surfaces and shorten the distance between them, causing Rs to drop. This becomes the case when agglomerates form, as is illustrated in a forthcoming figure.
Figure 7 shows the average values of electrolyte resistance for the different mixtures. Interestingly, the pore solution resistance reveals a tendency to decrease as the xGnP concentration increases, which could be partially associated with the increasing pore volume at lower xGnP concentrations and partially with agglomerate formation causing electrode extension and shorter pathways for ions (as will be shown later in this section) at higher xGnP concentrations.
On the other hand, the charge transfer resistance Rct, plotted versus xGnP concentration in Figure 8, is mainly a function of the amount of charge that transfers at the electrode/electrolyte interface. As follows, it is a function of the active surface area of the electrode, i.e., the electrode’s surface area that is in contact with electrolyte. Obviously, the higher the active surface area of the electrode, the smaller the charge transfer resistance. As can be derived from the following sections, with the increase in xGnP content, initially, both microporosity and mesoporosity increase. This increase manifests itself also as an increase in the active electrode surface area and therefore perfectly justifies the initial reduction in Rct up to a 0.10% xGnP concentration. Nevertheless, the increase in Rct as xGnP concentration further increases from 0.10% to 0.40% can only be explained as the result of formation of hydrophobic xGnP agglomerates, which obstruct electrolyte access to parts of the electrodes’ surface, as shown in Figure 9b. It is notable that although an increase in the highly microporous xGnP concentration leads to micropore volume enhancement (Table 2), their hydrophobic nature (oxygen < 1%) practically blocks the passage of water through their slit-like microporous structure and, at the same time, their adhesion with the hydrophilic concrete matrix is not expected to be tight, leaving gaps at the matrix/filler interface. These gaps are estimated to be in the range of mesopores, which fill with electrolyte and, in agreement with mercury porosimetry results (Table 2), become less in terms of mesopore volume upon agglomerate formation, as shown in Figure 9a.

3.2. Porosity Studies

In the current research, porosimetry tests were conducted to examine the effect of the addition of specific concentrations of xGnPs in the total porosity and pore size distribution of the cementitious matrix. According to the literature, CBNs act as nucleation points where C-S-H particles can grow. As a result, the cementitious microstructure is often denser, exhibiting reduced total porosity. In case that there is an excessive clustering of the nanoparticles, the total porosity of the nanocomposites can increase, and the pore size distribution can alter due to pores created by inefficient nanoparticle distribution.
Overall, mercury adsorption porosimetry provided the pore volume and surface area of mesopores, whereas nitrogen adsorption isotherms at 77 K provided the respective parameters for micropores. The porosimetry parameter results versus xGnP content are summed up in Table 2. In general, the amounts of adsorption were low (Figure 10), with LN2-porosimetry-calculated BET surface areas (Brunauer–Emmett–Teller) ranging between 7.6 and 14 m2/g and the measured total adsorbed N2 amounts fluctuating between 33 and 56 cc/g (STP).
In the analysis of pore size distribution, the BJH model [51] was employed, focusing on the desorption branch of the LN2 isotherms. The determined pore sizes (diameters) fell within the range 2 to 100 nm for all five examined samples. These samples exhibited similar isotherms, characterized as type IV with an H3 hysteresis loop, as depicted in Figure 10. The main difference among the five isotherms was the different total absorbed amounts and, in the case of sample “M-0.20”, which has the highest BET area and largest pore volume, also a slightly but evidently wider hysteresis loop. The pores that are added with the increasing xGnP content have pore diameters mainly in the range 4–100 nm. Additionally, the existence of a reproducible, permanent hysteresis loop in all cases is generally associated with capillary condensation and, in this present case, can be attributed to network effects [52]. The type of the presented hysteresis loop H3 corresponds to non-rigid aggregates of plate-like particles (e.g., certain clays) but also to macropores that are not completely filled with pore condensate. In the studied samples, both of the abovementioned cases can co-exist.
As can be derived by comparing the pore volumes given in Table 2, in all samples the main part of the total pore volume is micropores. Nevertheless, upon a xGnP content increase, the micropore volume was also found to increase, whereas the mesopore volume initially increased up to a 0.10 wt.% xGnP content and then dropped to an intermediate value; this was attributed to agglomerate formation at the higher concentrations, as described in Section 3.1.2 (Figure 8).

3.3. Mechanical Properties

Figure 11 shows the average fracture toughness of 28-day-old cement specimens. Overall, the fracture toughness shows elevated average values for the two mixtures reinforced with xGnPs compared with the plain mixture. The increase demonstrates a peak value for the M-0.05 (29% increase); then, for the higher xGnP contents (from 0.10 wt.% to 0.40 wt.%), the average KIc values start to drop, without reaching the reference mixture corresponding value.
The mechanical behaviour observed in the present study, showing an initial enhancement for mixtures with the lowest xGnP concentrations followed by a gradual decline as the nanoparticle content increases, aligns with findings reported in prior studies involving comparable GnP concentrations [16,19,30]. The enhanced toughness observed in the mixtures M-0.05 and M-0.10 can be attributed to the effective inhibition of micro-crack propagation, which is a result of the well-dispersed nanoparticles within the cementitious matrix. Conversely, less efficient dispersion, which may occur at higher nanoparticle concentrations, can compromise the flexural response of the cementitious nanocomposites. In fact, inadequate dispersion of xGnPs may lead to the formation of agglomerations, thereby contributing to an inferior mechanical performance.

3.4. Correlationship between Electrical and Mechanical Properties

Furthermore, from the comparison of the graphical curves of ρAC and KIc as a function of xGnP concentration (Figure 4 and Figure 11), it was observed that the two curves have an inverse response, which was revealed by the linear plot of KIc versus 1/ρAC (Figure 12) regarding the corresponding xGnP concentrations. The equation that describes the functional relationship between 1/ρAC and KIc is given in Equation (8).
KIc = 0.016/ρAC + 0.7237
Nevertheless, the aforementioned equation represents distinct material and experimental factors, including variations in the nanofiller type and concentration, the type of cement and superplasticizer employed, the water-to-cement (w/c) ratio, the specimen dimensions, and various other material and experimental attributes. This indicates that both the 1/ρAC-coefficient (0.016) and the constant value (0.7237) within the regression equation are anticipated to undergo modification when different nanomaterials and experimental conditions are chosen.
Recognizing that the correlation established in the present research article is contingent upon the specific material and experimental parameters employed, it is evident that further research is warranted to assess the potential transferability of this methodology to various types of cements and carbon-based nanofillers (CBNs). This avenue of inquiry holds promise for broadening the scope of fracture toughness assessment in cementitious nanocomposites. Nonetheless, it is evident that EIS can offer a straightforward approach that can be readily utilized for the assessment of fracture toughness in conductive cementitious nanocomposites.

4. Conclusions

In the present research article, electrochemical impedance spectroscopy (EIS) measurements were employed to examine multiple properties of graphene-based cementitious nanocomposites, including resistivity, porosity, and fracture toughness. During our experimental research, several noteworthy conclusions were drawn:
-
The electrical resistivity of the investigated cement mixtures exhibited a distinctive pattern. Initially, as the content of xGnPs increased (ranging from 0.05 wt.% to 0.10 wt.% by cement), the resistivity decreased. Nevertheless, for mixtures incorporating up to 0.40 wt.% xGnPs, a subsequent increase occurred, stabilizing at an intermediate constant value when compared with the reference mixture.
-
There was a visible trend in the reduction of pore solution resistance (Rs) as the xGnP concentration increased. This phenomenon can be attributed, in part, to the expansion of pore volume at lower xGnP concentrations and the formation of agglomerates. These agglomerates extended the electrode interface and consequently led to shorter ion pathways. Porosity measurements further substantiated the increase in total pore volume following the incorporation of xGnPs into the cement mixture.
-
Both EIS and porosimetry analyses provided compelling evidence of agglomerate formation, particularly when xGnP concentrations exceeded 0.10% by dry cement weight.
-
A functional linear relationship was observed between fracture toughness, as assessed through bending tests, and electrical resistivity, determined via EIS measurements.
Consequently, this study establishes the applicability of EIS as a valuable non-destructive tool for the comprehensive assessment of the physical and mechanical properties of cementitious nanocomposites.

Author Contributions

Conceptualization, N.D.A.; methodology, E.T., Z.S.M. and D.S.K..; investigation, E.T., Z.S.M., D.S.K. and E.P.F.; resources, S.K.K., E.P.F. and G.M.; writing—original draft preparation, E.T. and D.S.K.; writing—review and editing, E.P.F., Z.S.M. and N.D.A.; supervision, N.D.A.; project administration, N.D.A.; funding acquisition, N.D.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research is funded by the Research e-Infrastructure “Interregional Digital Transformation for Culture and Tourism in Aegean Archipelagos” (Code Number MIS 5047046), which is implemented within the framework of the “Regional Excellence” Action of the Operational Program “Competitiveness, Entrepreneurship and Innovation”. The action is co-funded by the European Regional Development Fund (ERDF) and the Greek State (Partnership Agreement 2014-2020).

Data Availability Statement

Data are available upon request to the authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Parra-Montesinos, G.J.; Reinhardt, H.W.; Naaman, A.E. High Performance Fiber Reinforced Cement Composites 6, 1st ed.; Springer: Dordrecht, The Netherlands, 2012. [Google Scholar]
  2. Tian, Z.; Li, Y.; Zheng, J.; Wang, S. A state-of-the-art on self-sensing concrete: Materials, fabrication and properties. Compos. Part B Eng. 2019, 177, 107437. [Google Scholar] [CrossRef]
  3. Oxana, V.K.; Leticia, M.T.M.; Boris, I.K. Recent Trends of Reinforcement of Cement with Carbon Nanotubes and Fibers. In Advances in Carbon Nanostructures; Adrian, M.T.S., Sonia, A.C.C., Eds.; IntechOpen: Rijeka, Croatia, 2016; Chapter 7. [Google Scholar]
  4. Du, M.; Jing, H.; Gao, Y.; Su, H.; Fang, H. Carbon nanomaterials enhanced cement-based composites: Advances and challenges. Nanotechnol. Rev. 2020, 9, 115–135. [Google Scholar] [CrossRef]
  5. Saleem, H.; Zaidi, S.J.; Alnuaimi, N.A. Recent Advancements in the Nanomaterial Application in Concrete and Its Ecological Impact. Materials 2021, 14, 6387. [Google Scholar] [CrossRef]
  6. Lee, C.; Xiaoding, W.; Kysarand, J.W.; Hone, J. Measurement of the Elastic Properties and Intrinsic Strength of Monolayer Graphene. Science 2008, 321, 385–388. [Google Scholar] [CrossRef]
  7. Buzaglo, M.; Bar, I.P.; Varenik, M.; Shunak, L.; Pevzner, S.; Regev, O. Graphite-to-Graphene: Total Conversion. Adv. Mater. 2017, 29, 1603528. [Google Scholar] [CrossRef] [PubMed]
  8. Novoselov, K.S.; Fal′ko, V.I.; Colombo, L.; Gellert, P.R.; Schwab, M.G.; Kim, K. A roadmap for graphene. Nature 2012, 490, 192–200. [Google Scholar] [CrossRef] [PubMed]
  9. Marconcini, P.; Macucci, M. Transport Simulation of Graphene Devices with a Generic Potential in the Presence of an Orthogonal Magnetic Field. Nanomaterials 2022, 12, 1087. [Google Scholar] [CrossRef] [PubMed]
  10. Huang, X.; Yin, Z.; Wu, S.; Qi, X.; He, Q.; Zhang, Q.; Yan, Q.; Boey, F.; Zhang, H. Graphene-Based Materials: Synthesis, Characterization, Properties, and Applications. Small 2011, 7, 1876–1902. [Google Scholar] [CrossRef]
  11. Pan, Z.; He, L.; Qiu, L.; Korayem, A.H.; Li, G.; Zhu, J.W.; Collins, F.; Li, D.; Duan, W.H.; Wang, M.C. Mechanical properties and microstructure of a graphene oxide–cement composite. Cem. Concr. Compos. 2015, 58, 140–147. [Google Scholar] [CrossRef]
  12. Lu, J.; Do, I.; Drzal, L.T.; Worden, R.M.; Lee, I. Nanometal-Decorated Exfoliated Graphite Nanoplatelet Based Glucose Biosensors with High Sensitivity and Fast Response. ACS Nano 2008, 2, 1825–1832. [Google Scholar] [CrossRef]
  13. Papanikolaou, I.; Ribeiro de Souza, L.; Litina, C.; Al-Tabbaa, A. Investigation of the dispersion of multi-layer graphene nanoplatelets in cement composites using different superplasticiser treatments. Constr. Build. Mater. 2021, 293, 123543. [Google Scholar] [CrossRef]
  14. Wotring, E.; Mondal, P.; Marsh, C. Characterizing the Dispersion of Graphene Nanoplatelets in Water with Water Reducing Admixture. In Proceedings of the Nanotechnology in Construction; Spring: Cham, Switzerland, 2015; pp. 141–148. [Google Scholar]
  15. Du, H.; Pang, S.D. Enhancement of barrier properties of cement mortar with graphene nanoplatelet. Cem. Concr. Res. 2015, 76, 10–19. [Google Scholar] [CrossRef]
  16. Tao, J.; Wang, X.; Wang, Z.; Zeng, Q. Graphene nanoplatelets as an effective additive to tune the microstructures and piezoresistive properties of cement-based composites. Constr. Build. Mater. 2019, 209, 665–678. [Google Scholar] [CrossRef]
  17. Du, H.; Pang, S.D. Dispersion and stability of graphene nanoplatelet in water and its influence on cement composites. Constr. Build. Mater. 2018, 167, 403–413. [Google Scholar] [CrossRef]
  18. Dela Vega, M.S.D.C.; Vasquez, M.R. Plasma-functionalized exfoliated multilayered graphene as cement reinforcement. Compos. Part B Eng. 2019, 160, 573–585. [Google Scholar] [CrossRef]
  19. Chougan, M.; Marotta, E.; Lamastra, F.R.; Vivio, F.; Montesperelli, G.; Ianniruberto, U.; Bianco, A. A systematic study on EN-998-2 premixed mortars modified with graphene-based materials. Constr. Build. Mater. 2019, 227, 116701. [Google Scholar] [CrossRef]
  20. Dalal, S.P.; Dalal, P. Experimental Investigation on Strength and Durability of Graphene Nanoengineered Concrete. Constr. Build. Mater. 2021, 276, 122236. [Google Scholar] [CrossRef]
  21. Jiang, Z.; Sevim, O.; Ozbulut, O.E. Mechanical properties of graphene nanoplatelets-reinforced concrete prepared with different dispersion techniques. Constr. Build. Mater. 2021, 303, 124472. [Google Scholar] [CrossRef]
  22. Lamastra, F.R.; Chougan, M.; Marotta, E.; Ciattini, S.; Ghaffar, S.H.; Caporali, S.; Vivio, F.; Montesperelli, G.; Ianniruberto, U.; Al-Kheetan, M.J.; et al. Toward a better understanding of multifunctional cement-based materials: The impact of graphite nanoplatelets (GNPs). Ceram. Int. 2021, 47, 20019–20031. [Google Scholar] [CrossRef]
  23. Metaxa, Z.S. Polycarboxylate Based Superplasticizers as Dispersant Agents for Exfoliated Graphene Nanoplatelets Reinforcing Cement Based Materials. J. Eng. Sci. Technol. Rev. 2015, 8, 1–5. [Google Scholar] [CrossRef]
  24. Wang, B.; Jiang, R.; Wu, Z. Investigation of the Mechanical Properties and Microstructure of Graphene Nanoplatelet-Cement Composite. Nanomater 2016, 6, 200. [Google Scholar] [CrossRef] [PubMed]
  25. Metaxa, Z.S. Exfoliated graphene nanoplatelet cement-based nanocomposites as piezoresistive sensors: Influence of nanoreinforcement lateral size on monitoring capability. Cienc. Tecnol. Dos Mater. 2016, 28, 73–79. [Google Scholar] [CrossRef]
  26. Dalla, P.T.; Dassios, K.G.; Tragazikis, I.K.; Exarchos, D.A.; Matikas, T.E. Carbon nanotubes and nanofibers as strain and damage sensors for smart cement. Mater. Today Commun. 2016, 8, 196–204. [Google Scholar] [CrossRef]
  27. Du, H.; Quek, S.T.; Pang, S.D. Smart Multifunctional Cement Mortar Containing Graphite Nanoplatelet; SPIE: Bellingham, WA, USA, 2013; Volume 8692. [Google Scholar]
  28. Al-Dahawi, A.; Sarwary, M.H.; Öztürk, O.; Yıldırım, G.; Akın, A.; Şahmaran, M.; Lachemi, M. Electrical percolation threshold of cementitious composites possessing self-sensing functionality incorporating different carbon-based materials. Smart Mater. Struct. 2016, 25, 105005. [Google Scholar] [CrossRef]
  29. Sun, S.; Ding, S.; Han, B.; Dong, S.; Yu, X.; Zhou, D.; Ou, J. Multi-layer graphene-engineered cementitious composites with multifunctionality/intelligence. Compos. Part B Eng. 2017, 129, 221–232. [Google Scholar] [CrossRef]
  30. Baomin, W.; Shuang, D. Effect and mechanism of graphene nanoplatelets on hydration reaction, mechanical properties and microstructure of cement composites. Constr. Build. Mater. 2019, 228, 116720. [Google Scholar] [CrossRef]
  31. Christensen, B.J.; Coverdale, R.T.; Olson, R.A.; Ford, S.J.; Garboczi, E.J.; Jennings, H.M.; Mason, T.O. Impedance Spectroscopy of Hydrating Cement-Based Materials: Measurement, Interpretation, and Application. J. Am. Ceram. Soc. 1994, 77, 2789–2804. [Google Scholar] [CrossRef]
  32. Xian, X.P.; Wang, Y.S.; Xing, F.; Dong, B.Q. Measuring and Modeling Analysis of Electrochemical Impedance Spectroscopy for Hydration Procedure of Cement Materials. Adv. Mater. Res. 2012, 588–589, 1033–1036. [Google Scholar] [CrossRef]
  33. Dotelli, G.; Mari, C.M. The evolution of cement paste hydration process by impedance spectroscopy. Mater. Sci. Eng. A 2001, 303, 54–59. [Google Scholar] [CrossRef]
  34. Tang, S.W.; Li, Z.J.; Chen, E.; Shao, H.Y. Impedance measurement to characterize the pore structure in Portland cement paste. Constr. Build. Mater. 2014, 51, 106–112. [Google Scholar] [CrossRef]
  35. Long, W.-J.; Gu, Y.-C.; Xing, F.; Khayat, K.H. Microstructure development and mechanism of hardened cement paste incorporating graphene oxide during carbonation. Cem. Concr. Compos. 2018, 94, 72–84. [Google Scholar] [CrossRef]
  36. Díaz, B.; Freire, L.; Merino, P.; Nóvoa, X.R.; Pérez, M.C. Impedance spectroscopy study of saturated mortar samples. Electrochim. Acta 2008, 53, 7549–7555. [Google Scholar] [CrossRef]
  37. Dong, B.; Qiu, Q.; Xiang, J.; Huang, C.; Xing, F.; Han, N. Study on the Carbonation Behavior of Cement Mortar by Electrochemical Impedance Spectroscopy. Materials 2014, 7, 218–231. [Google Scholar] [CrossRef]
  38. Kim, J.; Suryanto, B.; McCarter, W.J. Conduction, relaxation and complex impedance studies on Portland cement mortars during freezing and thawing. Cold Reg. Sci. Technol. 2019, 166, 102819. [Google Scholar] [CrossRef]
  39. Liu, Z.; Beaudoin, J. The permeability of cement systems to chloride ingress and related test methods. Cem. Concr. Aggreg. 2000, 22, 16–23. [Google Scholar]
  40. Wansom, S.; Kidner, N.J.; Woo, L.Y.; Mason, T.O. AC-impedance response of multi-walled carbon nanotube/cement composites. Cem. Concr. Compos. 2006, 28, 509–519. [Google Scholar] [CrossRef]
  41. Danoglidis, P.A.; Konsta-Gdoutos, M.S.; Gdoutos, E.E.; Shah, S.P. Strength, energy absorption capability and self-sensing properties of multifunctional carbon nanotube reinforced mortars. Constr. Build. Mater. 2016, 120, 265–274. [Google Scholar] [CrossRef]
  42. da Silva, G.F.; Martini, S.; Moraes, J.C.B.; Teles, L.K. AC impedance spectroscopy (AC-IS) analysis to characterize the effect of nanomaterials in cement-based mortars. Constr. Build. Mater. 2021, 269, 121260. [Google Scholar] [CrossRef]
  43. Li, W.-W.; Ji, W.-M.; Fang, G.-H.; Liu, Y.-Q.; Xing, F.; Liu, Y.-K.; Dong, B.-Q. Electrochemical impedance interpretation for the fracture toughness of carbon nanotube/cement composites. Constr. Build. Mater. 2016, 114, 499–505. [Google Scholar] [CrossRef]
  44. ASTM C305; Standard Practice for Mechanical Mixing of Hydraulic Cement Pastes and Mortars of Plastic Consistency. ASTM International: West Conshohocken, PA, USA, 2012.
  45. Bondarenko, A.S.; Ragoisha, G.A. Inverse Problem in Potentiodynamic Electrochemical Impedance Spectroscopy. Electrochim. Acta 2005, 50, 89–102. [Google Scholar] [CrossRef]
  46. McCarter, W.J.; Garvin, S.; Bouzid, N. Impedance measurements on cement paste. J. Mater. Sci. Lett. 1988, 7, 1056–1057. [Google Scholar] [CrossRef]
  47. ASTM E399; Standard Test Method for Linear-Elastic Plane-Strain Fracture Toughness of Metallic Materials. ASTM International: West Conshohocken, PA, USA, 2012.
  48. Haul, R.S.J.; Gregg, K.S.W. Berichte der Bunsengesellschaft für physikalische Chemie. In Sing: Adsorption, Surface Area and Porosity; Academic Press: London, UK, 1982; Volume 86, p. 957. [Google Scholar] [CrossRef]
  49. Washburn, E.W. Note on a Method of Determining the Distribution of Pore Sizes in a Porous Material. Proc. Natl. Acad. Sci. USA 1921, 7, 115–116. [Google Scholar] [CrossRef] [PubMed]
  50. Cook, R.A.; Hover, K.C. Experiments on the contact angle between mercury and hardened cement paste. Cem. Concr. Res. 1991, 21, 1165–1175. [Google Scholar] [CrossRef]
  51. Barrett, E.P.; Joyner, L.G.; Halenda, P.P. The Determination of Pore Volume and Area Distributions in Porous Substances. I. Computations from Nitrogen Isotherms. J. Am. Chem. Soc. 1951, 73, 373–380. [Google Scholar] [CrossRef]
  52. Thommes, M.; Kaneko, K.; Neimark, A.V.; Olivier, J.P.; Rodriguez-Reinoso, F.; Rouquerol, J.; Sing, K.S.W. Physisorption of gases, with special reference to the evaluation of surface area and pore size distribution (IUPAC Technical Report). Pure Appl. Chem. 2015, 87, 1051–1069. [Google Scholar] [CrossRef]
Figure 1. Prismatic specimens used for (a) mechanical tests (three-point-bending) and (b) electrochemical impedance spectroscopy measurements.
Figure 1. Prismatic specimens used for (a) mechanical tests (three-point-bending) and (b) electrochemical impedance spectroscopy measurements.
Nanomaterials 13 02652 g001
Figure 2. Flexural (three-point bending) test set-up and specimen key dimensions (W: height, B: width, S: span, and α: depth of the notch of the specimen) used for calculation of fracture toughness (KIc).
Figure 2. Flexural (three-point bending) test set-up and specimen key dimensions (W: height, B: width, S: span, and α: depth of the notch of the specimen) used for calculation of fracture toughness (KIc).
Nanomaterials 13 02652 g002
Figure 3. Typical Nyquist plots of the investigated nanoreinforced cement specimens.
Figure 3. Typical Nyquist plots of the investigated nanoreinforced cement specimens.
Nanomaterials 13 02652 g003
Figure 4. Τotal ohmic resistivity from Nyquist cusp points vs. xGnP concentration.
Figure 4. Τotal ohmic resistivity from Nyquist cusp points vs. xGnP concentration.
Nanomaterials 13 02652 g004
Figure 5. (a) EC model (Randles circuit) describing cement paste without xGnPs and (b) EC model describing cement pastes with xGnPs.
Figure 5. (a) EC model (Randles circuit) describing cement paste without xGnPs and (b) EC model describing cement pastes with xGnPs.
Nanomaterials 13 02652 g005
Figure 6. Set-up representation for EIS measurements on porous electrolyte-impregnated solid.
Figure 6. Set-up representation for EIS measurements on porous electrolyte-impregnated solid.
Nanomaterials 13 02652 g006
Figure 7. Average resistance of electrolyte solution according to EIS for mixtures with various xGnP concentrations, as obtained by fitting Nyquist plots with the EC model of Figure 5b.
Figure 7. Average resistance of electrolyte solution according to EIS for mixtures with various xGnP concentrations, as obtained by fitting Nyquist plots with the EC model of Figure 5b.
Nanomaterials 13 02652 g007
Figure 8. Average resistance (Rct) in the cementitious matrix with varying xGnP concentration, as obtained by the EC model.
Figure 8. Average resistance (Rct) in the cementitious matrix with varying xGnP concentration, as obtained by the EC model.
Nanomaterials 13 02652 g008
Figure 9. Schematics representing (a) mesoporosity reduction upon agglomeration and (b) agglomerate-induced electrode active area reduction, as well as ion path distance reduction. The blue color represents electrolyte (tap water), light grey represents concrete, dark grey represents xGnP filler, and black represents steel (electrodes).
Figure 9. Schematics representing (a) mesoporosity reduction upon agglomeration and (b) agglomerate-induced electrode active area reduction, as well as ion path distance reduction. The blue color represents electrolyte (tap water), light grey represents concrete, dark grey represents xGnP filler, and black represents steel (electrodes).
Nanomaterials 13 02652 g009
Figure 10. N2 adsorption isotherms at 77 K and BJH pore size distribution curves (inset) for different xGnP concentrations: (a) 0–0.1 wt.% and (b) 0–0.40 wt.%. Filled symbols stand for adsorption branch and empty symbols for desorption branch.
Figure 10. N2 adsorption isotherms at 77 K and BJH pore size distribution curves (inset) for different xGnP concentrations: (a) 0–0.1 wt.% and (b) 0–0.40 wt.%. Filled symbols stand for adsorption branch and empty symbols for desorption branch.
Nanomaterials 13 02652 g010
Figure 11. Average fracture toughness (KIc) of pastes with different xGnP concentrations.
Figure 11. Average fracture toughness (KIc) of pastes with different xGnP concentrations.
Nanomaterials 13 02652 g011
Figure 12. Correlation between KIc and 1/ρAC and the fitting curve for cementitious composites with varying xGnP concentrations.
Figure 12. Correlation between KIc and 1/ρAC and the fitting curve for cementitious composites with varying xGnP concentrations.
Nanomaterials 13 02652 g012
Table 1. Mixture compositions used in the present investigation.
Table 1. Mixture compositions used in the present investigation.
Materials
MixtureCEM IWaterxGnPsSP
(g)(g)(g)(g)
M-0.003801140.0-
M-0.053801140.191.52
M-0.103801140.383.04
M-0.203801140.766.08
M-0.403801141.5213.16
Table 2. Exportable data from mercury and nitrogen porosimetry.
Table 2. Exportable data from mercury and nitrogen porosimetry.
MixtureMercury PorosimetryN2 Porosimetry
ε (%)Density (g/cm3)Pore Volume
(cm3/g)
s.s.a 1 (m2/g)Pore Radius
(nm)
Pore
Volume (cm3/g)
Pore
Radius (nm)
s.s.a (BET) (m2/g)
BulkApparent
M-0.0042.142.230.01876.0784–366533.504–20010.08
M-0.05182.132.600.084634.684–265545.504–2007.60
M-0.10182.092.540.086539.184–327046.954–2008.85
M-0.20122.122.440.060627.984–273055.354–20014.04
M-0.40162.092.500.077130.464–269052.304–20011.4
1 Specific surface area.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tziviloglou, E.; Metaxa, Z.S.; Maistros, G.; Kourkoulis, S.K.; Karousos, D.S.; Favvas, E.P.; Alexopoulos, N.D. Electrochemical Impedance as an Assessment Tool for the Investigation of the Physical and Mechanical Properties of Graphene-Based Cementitious Nanocomposites. Nanomaterials 2023, 13, 2652. https://doi.org/10.3390/nano13192652

AMA Style

Tziviloglou E, Metaxa ZS, Maistros G, Kourkoulis SK, Karousos DS, Favvas EP, Alexopoulos ND. Electrochemical Impedance as an Assessment Tool for the Investigation of the Physical and Mechanical Properties of Graphene-Based Cementitious Nanocomposites. Nanomaterials. 2023; 13(19):2652. https://doi.org/10.3390/nano13192652

Chicago/Turabian Style

Tziviloglou, Eirini, Zoi S. Metaxa, George Maistros, Stavros K. Kourkoulis, Dionysios S. Karousos, Evangelos P. Favvas, and Nikolaos D. Alexopoulos. 2023. "Electrochemical Impedance as an Assessment Tool for the Investigation of the Physical and Mechanical Properties of Graphene-Based Cementitious Nanocomposites" Nanomaterials 13, no. 19: 2652. https://doi.org/10.3390/nano13192652

APA Style

Tziviloglou, E., Metaxa, Z. S., Maistros, G., Kourkoulis, S. K., Karousos, D. S., Favvas, E. P., & Alexopoulos, N. D. (2023). Electrochemical Impedance as an Assessment Tool for the Investigation of the Physical and Mechanical Properties of Graphene-Based Cementitious Nanocomposites. Nanomaterials, 13(19), 2652. https://doi.org/10.3390/nano13192652

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