Review of Several Experimental Methods for Characterization of Micro- and Nano-Scale Pores in Cement-Based Material

Mechanical properties and durability of cement-based materials are largely affected by pore structures. This paper provides an overview of several experimental techniques to characterize pore size distribution and specific surface area, with focus on pores in calcium silicate hydrates. The reviewed experimental techniques are nitrogen and water vapor sorption isotherm, proton nuclear magnetic resonance (1H-NMR) and small-angle scattering (SAS). Different pretreatment methods are compared for sorption measurements. Pore size distribution and specific surface area are analyzed using data from different methods to understand difference and consistency of these methods. It is found that pore size distribution calculated from sorption isotherm is very sensitive to adsorption model. Though specific surface areas from different techniques are quite different from each other, they are all able to detect the microstructural alteration due to long-term drying.


Introduction
Cement paste and concrete are porous material composed from solid skeleton and pores, with the size of pores varying from several micrometers to nanometers. These pores are of great importance for not only mechanical properties but also for durability of cementbased materials. Most of the chemical agents, that are triggers of degradation, invade into the material through pores. Chemical reactions then takes place inside the pores, causing different results based on the size of the pores. For instance, during salt attack in cement paste, crystallization pressure is different in larger pores and smaller pores (Scherer 1999(Scherer , 2004. It was confirmed in Sasano et al. (2018) that drying shrinkage reduced the performance of reinforced concrete structure. This drying shrinkage is caused by three different mechanisms in different pores at different relative humidities (Hansen 1987). In larger pores at higher relative humidity, capillary pressure causes the drying shrinkage (Feldman and Sereda 1970;Powers 1965Powers , 1968. When relative humidity continues to decrease, disjoining pressures of water in smaller pores play main role (Bazant 1972;Wittmann 1973); whereas, at even lower relative humidity in much smaller pores, the shrinkage is caused by surface energy (Feldman and Sereda 1964;Powers 1965). Therefore, characterizing pore size and pores structure is necessary to better estimate the durability of cementitious materials and has been a topic of many studies since decades. This article aims at reviewing such experimental techniques and related theories that are widely used on Ordinary Portland Cement (OPC) pastes. Blended cement pastes will not be considered specifically, but the experimental techniques and theories remain the same as OPC pastes.
Ordinary cement paste is a multi-phase material, composed from unhydrated clinker and hydration products. The hydration products include crystalline phases such as Portlandite, sulfoaluminates and amorphous calcium

Open Access
International Journal of Concrete Structures and Materials   14:55 silicate hydrate (C-S-H). During hydration, calcium oxide and silicate from clinker dissolve into water. Calcium oxides forms layers on which silicates will be paired in the form of dimers (see Fig. 1). Two pairing dimers next to each other can be connected via the bridging position, which can be occupied by silica or aluminum or calcium ions. However, some of the bridging positions are empty or cannot be connected to the next bridging position. Hence, calcium ions can be more concentrated in some places than others, causing more water in these places. As a result, calcium oxide layers will be winded, forming larger spaces between them in some place and smaller spaces in other (Gartner et al. 2017). These larger spaces are called gel pores, while the smaller spaces are known as interlayer space, as shown Fig. 1. Though different from the models in Geng et al. (2017;Richardson 2004) which attribute the gel pores to the imperfect growth of calcium silicate layers, they all suggest that both gel pores and interlayer space are inside of C-S-H, and connected to each other. The connection between gel pores and interlayer spaces is also confirmed by T 2 − T 2 correlation of gel pores and interlayer spaces in 2D proton nuclear magnetic resonance relaxometry (Mcdonald et al. 2007;Valori et al. 2010). Said otherwise, the calcium silicate layers, interlayer water and gel space form together the C-S-H gel. Due to the small number of repeating unit in the calcium silicate layer, experimental techniques for crystallography such as X-ray diffraction cannot be used to study the structure of the C-S-H that appears amorphous. However, it is possible to observe that calcium silica layers are crystal-based structure by transmission electron microscope (TEM), where comes the name of C-S-H nanocrystalline (Richardson 1999). According to Jennings (2000), the size of the interlayer space between calcium silica layers is around 0.1 nm-1 nm; whereas, the size of gel pores are 2 nm-10 nm. On the contrary, Mehta and Monteiro (2006) named all pore spaces in C-S-H gel as interlayer space. Outside of the C-S-H gel, the cement paste contains also capillary pores whose size varies from 10 µm to 10 nm. Capillary pores result from consumption of water during hydration and locate between hydrates and clinkers. Recently, based on relaxation times measured by proton nuclear magnetic resonance relaxometry ( 1 H-NMR), Muller et al. (2013) classified pores into four types: interlayer space, gel pores, interhydrate pores and capillary pores. Among them, interlayer space and gel pores are confined by C-S-H layers. Since most long-term durability issues originate from C-S-H, experimental methods to review in this work are selected with an objective to study gel pores and interlayer spaces of C-S-H.
Many experimental techniques have been used to characterize the pore structure of cementitious materials, including direct methods, in which pore structure is observed by imaging techniques, and indirect methods, in which some probe is used to get information about the pore structure. The range of measurable pore size of each method is limited to only a part of the full pore structure of cement paste. Figure 2 resumes the approximate range of pore size of several different methods.
In direct methods, pore structure is analyzed from 2D or 3D images of samples. The key factor of the imagebased methods is the size of a pixel/voxel, which limits the range of pore size. Smaller pixel/voxel size permits access to smaller pores, but limits, on the other hand, the size of total observation area; hence, representativity of the measurement would be lower. In one of the pioneer works, 2D images of Scanning Electron Microscopy (SEM) with pixel size around 0.5 μm (Lange et al. 1994;Scrivener 1989Scrivener , 2004 were used to estimate total porosity, pore size distribution. Comparing to 2D image of SEM, 3D images provides more topological information on pore network. For instance, Lanzón et al. (2012) used X-ray microtomography to analyze the 3D pore structure  Mater (2020) 14:55 of mortar while 3D microstructure of cement paste was obtained by X-ray microtomography and X-ray nanotomography in Bossa et al. (2015). Laser scanning confocal microscopy (LSCM) can be combined with serial sectioning to analyze influence of 3D pore structure of cement paste (Yio et al. 2015(Yio et al. , 2019. It is also possible to reconstruct 3D image combining the SEM with Focused Ion Beam (FIB) (Holzer et al. 2004). Although the size and shape of the pores can be observed directly with these methods, sophisticated image analysis algorithms are necessary to quantify the pore structure. Among the indirect methods, there are experimental techniques in which a substance is introduced as probe into the sample to measure the pore size and volume, such as Mercury Intrusion Porosimetry (MIP) and gas/ vapor adsorption techniques. MIP is one of the widely used techniques to measure pores structure of cementitious materials (Van Brakel et al. 1981;Muller and Scrivener 2017). MIP was applied on concrete for the first time in 1961 (Edel'man and Sominskii 1961). In MIP, sample is dried prior to the test and then mercury is forced into the pores of the material, with step-by-step increasing pressure. The experimental relation between applied pressure and volume of introduced mercury is then used to infer pore size distribution based on Laplace Equation/ Washburn Equation (see Eq. (5)). MIP is known as a destructive method (Diamond 2000) and recently, using X-ray CT image analysis, Wang et al. (2019) quantified the damage caused by MIP. In addition, pretreatment disturbs the pore structure and measurement results of MIP are affected significantly by pretreatment method (Feldman and Beaudoin 1991). Comparing different pretreatment methods, Konecny and Naqvi (1993) suggested that solvent exchange method keeps better the microstructure for MIP and this was later confirmed by Muller and Scrivener (2017) comparing pore size distribution from MIP with the results from 1 H-NMR. From theoretical point of view, analysis of MIP results from Eq. (5) is also sensitive to the values of surface tension and contact angle between mercury fronts and pore walls of cementitious materials. The value of surface tension between Mercury and cement paste is generally accepted to be equal to 0.48 N/m (e.g., (Powers et al. 1954;Taylor 1997)). However, different values were reported for contact angle, probably because the contact angle is sensitive to many parameters such as the age of sample, roughness of the surface, pretreatment and the purity of Mercury. Shi and Winslow (1985) measured the contact angle between 123 and 135°, and observed that the contact angle may vary depending on cement composition, age of the cement paste and pretreatment. Adolphs et al. (2002) found the contact angle to be dependent on equilibrium relative humidity and wetting-drying cycle of the sample before measurement, and reported values between 142 and 150°. The dependency of contact angle on intrusion-extrusion cycle has been studied in Zeng et al. (2012) by defining a hysteresis factor between contact angle of intrusion cycle and extrusion cycle. Recently, a much lower value of 120° was proposed based on the comparison of MIP and 1 H-NMR (Muller and Scrivener 2017). By comparison between the volume of capillary pores obtained from MIP and that from 1 H-NMR, they also suggested that the so-called ink-bottle effect is minor for the range of pores that can be seen by MIP (Muller and Scrivener 2017), which is not consistent with the point of view in many other articles, such as Diamond (2000). While reviewing MIP technique for cementitious materials, Diamond (2000) argued that the pore size seen in MIP measurement is not the size of pore but the size of pore entry, which is accepted by many others. For instance, Zeng et al. (2019) called the result of MIP as throat size distribution (TSD) and studied the pore size distribution from mercury extrusion data. Seen from Eq. (5), the range of pore size that can be detected by MIP depends on the value of applied pressure. For the maximum pressure applied in common MIP tests, the corresponding minimum pore size that can be detected by MIP is around 4 nm as resumed in Fig. 2. Hence, MIP cannot see internal structure of C-S-H whose size is smaller than this usual range. More suitable method to characterize the C-S-H is sorption isotherms which will be reviewed in detail in Sect. 2.
There are other indirect methods that use water inside the cement paste as probe, such as low-temperature calorimetry (LTC), 1 H-NMR, Small-angle scattering (SAS, includes Small-angle X-ray scattering (SAXS) and smallangle neutron scattering (SANS)). Main advantages of these methods are that the samples can be measured 'as prepared' , so pore structure is not disturbed before measuring. The principle of LTC is to link the pore size to the heat emitted during cooling of pore water, based on Gibbs-Thomson equation (Beddoe and Setzer 1990;Fagerlund 1973;Tombari et al. 2005). In experiment, temperature of sample is decreased gradually from 20 to − 60 °C, while the heat flow is recorded as a function of temperature (Fagerlund 1973). Though the LTC is relatively fast and young samples can be measured, the results are affected by rate of temperature change, alkali concentration in the pore solution (Kjeldsen and Geiker 2008). Furthermore, due to strong adsorption forces, water in interlayer space do not freeze even at − 42 °C (Bager and Sellevold 1986;Brun et al. 1977), at which the last peak of heat emission is observed. Maruyama et al. (2018) quantified that only 20% of the total evaporable water was frozen. Hence, it is difficult to observe pore structure in C-S-H by LTC. Cracks can be caused during ice formation and affect the measurement results. To improve the LTC analysis, Kjeldsen and Geiker (2008) recommended adding nucleation agent and correcting measurement temperature for non-equilibrium between the sample and the reference block. Comparing to LTC, 1 H-NMR and SAS are non-destructive techniques. How to analyze pore structure from 1 H-NMR and SAXS/ SANS will be reviewed in Sects. 3 and 4, respectively.
After reviewing the above-mentioned methods, influence of pretreatment methods on sorption measurement will be compared. Pore size distribution and specific surface area will be analyzed from different experimental methods on white Portland cement paste.

Sorption Isotherm
Sorption isotherm is the oldest and most widely used technique to measure the surface area of cement paste, and probably with the most published results (Thomas et al. 1999). In this section, first, we review different methods that have been used in literature to measure sorption isotherms for cement paste. Then, BET sorption theory will be presented for calculation of specific surface area, followed by estimation of pore size distribution by BJH method and the local density method.

Experimental Methods
Sorption isotherms consist in measuring the amount of condensed liquid of an adsorbate on the surface of an adsorbent, while controlling the relative pressure p/p sat at a constant temperature. Generally, there are gravimetric and volumetric methods to determine the adsorbed amount. In the gravimetric method, samples are put into a chamber with controlled relative humidity. Mass of the sample is weighted until equilibrium is reached between the relative pressure in the chamber and the sample. Then, measurement is moved to the next target relative humidity step by step. In the volumetric methods (also called manometric methods), known amount of adsorbate gas is introduced into sample chamber where the gas pressure is monitored. By adding or removing gas molecules, an equilibrium is achieved at the target relative pressure. Then, amount of adsorption is calculated from ideal gas equation.
Water vapor and nitrogen are the two adsorbates used most widely for the sorption measurement of cement paste. In addition, many other gas have also been used as adsorbate for the sorption measurement of cement paste, for instance, oxygen (Blaine and Valis 1949), argon (Blaine and Valis 1949;Vidick 1987), krypton (Vidick 1987), organic vapors (Mikhail and Selim 1966), ammonia (Odler 2003). Important aspects to consider when choosing an adsorbate are as follows: the heat of adsorption, the size of gas molecules, the polarity of gas molecules and chemical reactivity with cement paste (Thomas et al. 1999). In this review, we focus only on water vapor and nitrogen sorption isotherms.
Generally, the adsorption amount, when expressed in volume of liquid, is less in nitrogen sorption comparing to that in water vapor sorption. Said otherwise, sorption capacity and surface area measured by nitrogen sorption are smaller comparing to that measured by water vapor sorption. An example of comparison is shown in Fig. 3a, on white Portland cement paste, with watercement ratio 0.55 (data from Maruyama et al. (2014)), by the volumetric method. Pretreatment is necessary for the volumetric method since the starting state should be under vacuum in most of the sorption devices. Samples are crushed into powder (sizes of particles are in between 25 μm and 75 μm) prior to the test and dried under vacuum condition at 20 °C. Nitrogen sorption was measured at − 196 °C (BELSORP-miniII MicrotracBEL. Corp) on 75 ± 5 mg sample, while water vapor sorption was measured at 20 °C (VSTAR, Quantachrome instrument) on 25 ± 1 mg sample. As expected, at given relative pressure, volume of adsorption in nitrogen sorption is less than water vapor sorption. Even when nitrogen and Page 5 of 18 Aili and Maruyama Int J Concr Struct Mater (2020) 14:55 water vapor sorption isotherms are compared for same Kelvin radius (defined in Eq. (5)), adsorption in nitrogen is still less than water (see Fig. 3b). Such difference has been compared in many articles (e.g., (Odler 2003;Robens et al. 2004)) and different interpretation exists in the literature. It was suggested that due to larger size of a nitrogen molecule (kinetic diameter of 4.05 Å) compared to a water molecule (kinetic diameter of 3.25 Å), nitrogen sorption cannot access smaller pores, but Feldman and Sereda (1968) argued that this tiny difference in size of molecules cannot explain such big difference of adsorption and suggested that water vapor was structurally incorporated or re-adsorbed into C-S-H interlayer space. It was proposed in Brunauer et al. (1970) that the fact that water molecules possess a distinct dipole, whereas nitrogen molecules do not possess, may play a role. Some others suggested that the difference of adsorption amount between nitrogen and water may come from the fact that nitrogen sorption is usually measured at a temperature much lower than water vapor sorption (Brunauer et al. 1970;Hagymassy et al. 1972;Mikhail et al. 1964). Based on the difference of nitrogen sorption and water vapor sorption Jennings (2000); Tennis and Jennings (2000) suggested the existence of Low-Density (LD) C-S-H, to which nitrogen can access and High-Density (HD) C-S-H, to which nitrogen cannot access. Following their work, though not proved experimentally yet, the most popular view on the difference between nitrogen and water vapor sorption is that nitrogen sorption can access only gel pores; whereas, water vapor can see both gel pores and interlayer space (for instance, Maruyama et al. (2014)). On the contrary, for gravimetric method, measurement can be performed on block samples, not necessarily starting from dry state. When block samples are used, measurement time can be quite long. For instance, up to years for samples with thickness of 5 mm (Baroghel-Bouny 2007), or up to weeks for samples with thickness of 1 mm samples (Maruyama et al. 2015(Maruyama et al. , 2016. As an example, water vapor sorption isotherms measured by different techniques are compared in Fig. 4. The samples are white Portland cement paste, with water-cement ratio 0.55, stored in lime water for at least half year before starting the test. 3 × 13 × 100 mm 3 bar shape sample was dried in different chambers, where relative humidity was controlled by salt-saturated solution, during 1 year, which is noted as the long-term drying method (Maruyama et al. 2014). In the short-term drying method, sorption was measured on sample of size 3 × 3 × 1 mm 3 using RH-TG (denotes for Relative Humidity-Thermo-Gravimetry). Mixture of water vapor and nitrogen was blown constantly into the cell of RH-TG to keep relative humidity at the target value (Maruyama and Rymeš 2019). Each step of relative humidity lasted for 8 h, since more than 90% of mass loss occurred during the first 8 h.
To compare with the long-term and short-term gravimetric methods, the desorption branch of the data from Fig. 3 is also plotted as volumetric method in Fig. 4. The sorption amount was normalized with respect to the vacuum-dried mass at 20 °C. We can see from Fig. 4 that the long-term drying method measures more water loss than the short-term drying, which can be explained by the microstructural change of C-S-H due to the longterm drying proposed in Gartner et al. (2017), Maruyama et al. (2014). This difference between short-term and long-term drying method was also acknowledged in Nguyen et al. (2019) and attributed to slowly surface diffusion in nanopores. As for the volumetric method, the sorption amount at higher relative humidity is lower than the short-term drying method, probably due to two reasons: (i) It is difficult in volumetric method to control relative humidity at high values; hence, large pores cannot be fully filled; (ii) More importantly, due to drying of the sample during the pretreatment of volumetric method, some irreversible changes of microstructure happen (Maruyama et al. 2014;Parrott et al. 1980); Influence of different pretreatment method will be compared further in Sect. 5.1.

Surface Area from Sorption Models
One of the most important information that we get from a sorption isotherm is the surface area of sample. The basic concept of sorption models is to obtain monolayer capacity w m (i.e., amount of adsorption when sample surface is covered fully by one layer of sorptive molecules) and then calculate surface area by: Fig. 4 Water vapor sorption isotherm measured by gravimetric method, i.e., slowly drying (i.e., long-term drying) and RH-TG (i.e., short-term drying) (data from Maruyama et al. (2014) and Maruyama and Rymeš (2019)) and volumetric method (data from Maruyama et al. (2014)).
Page 6 of 18 Aili and Maruyama Int J Concr Struct Mater (2020) 14:55 where M m and s are the molar mass of adsorbate and cross-sectional area of adsorbate molecule, expressed in [kg/mol] and [m 2 ], respectively; N A is the Avogadro's number. The cross-sectional area is measured on the plane of a solid surface, dividing the surface by the number of adsorbed molecules. This cross-sectional area is different from the section calculated from the abovementioned kinetic diameter. The cross-sectional area of water molecule is equal to 0.114 nm 2 and of nitrogen molecule equal to 0.162 nm 2 (Mikhail and Selim 1966). The monolayer capacity is the key parameter to be obtained by sorption models. Brunauer, Emmett and Teller (BET) Theory (Brunauer et al. 1938) is the most widely used adsorption model as it can model quite well the sorption isotherm of cementbased material and also thanks to its simplicity. Like many other sorption models, BET model is based on the kinetic approach in which the number of adsorbed and desorbed gas molecules is related to the temperature and the free energy of adsorption. The main hypothesis of BET model is that the adsorption energy is higher in the first layer of adsorption comparing to all the other layers; and for all layers starting from 2nd layer on, the adsorption energy remains the same and equal to that of bulk liquid phase. At equilibrium, summing up the adsorbed molecules of each layer, total adsorption amount is obtained as a function of relative pressure p/p sat as following: where w is the water content [−]; C BET is the energy constant [−]. For cement pastes, the hypothesis of adsorption of BET theory is generally applicable only in the range of relative humidity 5-35% on the adsorption branch. Hence, the BET Eq. (2) is normally fitted in this range. In practice, the adsorption data are plotted in a linear form, which is known as the BET plot to obtain the two parameters w m and C BET . Then inserting the monolayer capacity w m into Eq. (1), the BET surface S BET area is obtained. As reviewed in Sect. 2.1, water vapor can be adsorbed more than nitrogen, resulting in higher water vapor surface area S H 2 O than nitrogen surface area S N 2 . It is generally assumed, though not proved yet, that the difference between the water vapor surface area S H 2 O and nitrogen surface area S N 2 is the surface area of interlayer space of (Maruyama et al. 2014;.
Although only adsorption branch of sorption isotherms is used to evaluate the BET-specific surface area Maruyama and Igarashi (2011) observed that the desorption branch of water vapor sorption isotherm also follows well the BET theory. This may be due to the fact that, in the desorption branch of sorption experiment, the water adsorbed in the interlayer space of C-S-H layers does not evaporate until very low pressure. The desorption only from gel pores, i.e., outer surface of C-S-H may correspond the hypothesis of the BET theory. Other sorption models have been also used sometimes. Older than BET model, Langmuir model (Langmuir 1918) assumes only monolayer adsorption and suites well for chemisorption. Guggenheim-Andersonde Boer (GAB model) extended the applicability of BET theory on whole range of relative humidity (Anderson 1946;Brunauer et al. 1969;Timmermann 2003). Comparing to BET theory, GAB model introduced one more energy constant to differentiate the adsorption energy of second layer from that of all other subsequent layers. In Anderson (1946) the BET theory was generalized by distinguishing the adsorption energy of each layers. The Langmuir theory, BET theory and GAB model can regarded as one of the specific cases of the generalized BET theory (Anderson 1946;Baroghel-Bouny 1994). Recently, Nguyen et al. (2019Nguyen et al. ( , 2020 modified the BET theory by taking into account the reduction of the adsorption surface with the increase of number of absorbed layers. Despite these progresses, BET surface area remains for the moment the mostly reported values in literature.

Pore Size Distribution by BJH Model
Inferring pore size from sorption isotherm is more laborious than inferring the specific surface area. Though various models have been proposed in literature, in this section, we review basics of BJH method, without going into details of numerical techniques.
During a sorption measurement, condensation of gas is governed by Kelvin's law (Coussy 2011), which was obtained by equating the chemical potential µ L of the condensed liquid phase and µ G of gas phase of the adsorbate. Kelvin's law relates the capillary pressure P c , which is the pressure difference between liquid phase and gas phase, to the relative pressure by: where ρ l is the density of liquid phase, expressed in [kg/ m 3 ]. R is the universal gas constant, equal to 8.314 J/ mol/K and T is the absolute temperature. Due to this capillary pressure P c , the interface between condensed liquid (4) P c = ρ l RT M m ln(p/p sat ), and gas is not flat but forms a curved meniscus. Radius of the capillary meniscus is also related to capillary pressure via Laplace equation (e.g., see Sect. 6.1.7 in (Coussy 2011)) or known also as Washburn's equation (Gardner 1921;Muller and Scrivener 2017): where r k is the mean radius of the capillary meniscus, expressed in [m], θ is the contact angle between adsorbate and solid surface and γ is the surface tension of adsorbate, expressed in [N/m]. For water and nitrogen, the contact angle is considered to be 0. The surface tension depends on temperature. For water, surface tension is equal to 0.072 N/m at 20 °C; for nitrogen, equal to 0.0089 N/m at −196 °C. For a given relative humidity, combining Eqs. (4) and (5), radius r k of meniscus is obtained and named as Kelvin radius. At each step of sorption measurement, when equilibrium is reached, capillary meniscus only with the corresponding Kelvin radius remains. The relation between pore size and Kelvin radius depends on the assumption of the pore geometry. Figure 5 illustrates geometrical relation between Kelvin radius r k and pore diameter d for cylindrical pore and slit pore. In a cylindrical pore, the meniscus has a spherical shape which leads to a pore diameter d = 2r k ; whereas in slit shaped pores, the meniscus is cylindrical and the pore width would be d = r k (Lowell 2005).
The BJH method (Barrett 1951), also called modified Kelvin equation, considers that the condensed liquid (5) r k = − 2γ cosθ P c , has two parts: one is capillary liquid, whose volume depends on relative pressure via Kelvin-Laplace equation; the other is adsorbed layer on pore surface, whose thickness may be considered as a function of relative pressure. When relative pressure changes, both capillary liquid and adsorbed thickness change. Letting sum of the volume change of capillary and adsorbed liquid equal to the experimentally measured adsorbed amount change, an explicit equation of pore size can be obtained. Solving this equation for whole range of adsorption, a pore size distribution can be obtained. Exact expression of BJH method and an approximated expression were derived in Halenda (1951) for cylindrical pore. It should be kept in mind that, considering the geometrical relation illustrated in Fig. 5, the equation of BJH method is different when different pore shape is assumed. For instance, Brunauer et al. (1967) compared volume, surface area and pore size from BJH model for cylindrical, slit and spherical pores. Hence, to summarize, BJH model includes three inputs: pore shape, Kelvin-Laplace equation and a model to estimate adsorption thickness. Different models for adsorption thickness have been used in literature (e.g., see chapter 8 in Lowell (2005)). For instance, both Langmuir sorption model and the BET model in Eq. (2) have been used in Pinson et al. (2015) to estimate thickness of adsorption of water sorption isotherm. Adsorption thickness proposed by Hagymassy et al. (1969) was also widely used in pore structure analysis from water vapor sorption isotherm (Zeng et al. 2014). Again for water sorption isotherm, in calculation where t is the thickness of adsorption in nanometer; parameters K 1 and K 2 are expressed in nanometer. As for nitrogen sorption, Harkins-Jura model (Harkins and Jura 1944), expressed in nanometer, is widely used. When different thickness models are used, BJH model may give totally different pore size distributions. Figure 6 displays an example of pore size distribution from adsorption branch of water vapor sorption displayed in Fig. 3. The parameters of Badmann model were fitted on adsorption branch in the range of relative humidity between 3 and 10%. It can be seen from Fig. 6 that BET model and Badmann model gives very different pore size distribution. Hence, it is important to keep in mind that the adsorption thickness model could be crucial for the evaluation of pore size distribution from BJH method.

Pore Size Distribution by Density Functional Approach
If the BJH method is considered as backward problem, then the corresponding forward problem would be to estimate sorption isotherm from known pore size distribution. Assuming the pore size distribution function as f (φ), 1 the sorption isotherm reads as: 13.99 0.034 − log 10 h r 1/2 , with φ max and φ min are maximum and minimum pore size, respectively; ρ L (φ, h r ) is the quantity of liquid per unit volume of pore with size φ under given relative humidity. ρ L (φ, h r ) is called the density function, depends on pore size and should not be confused with the density of bulk liquid ρ l , which is an intrinsic property of the liquid itself. If the density function ρ L (φ, h r ) can be estimated for each pore size a priori, then assuming some type of pore size distribution function, by back-fitting Eq. (8) on the measured sorption isotherm, we can obtain the pore size distribution. The main task of this type of density functional approach is to estimate the density function ρ L (φ, h r ).
The most common way in literature of cementitious material is to calculate ρ L (φ, h r ) as the sum of condensed liquid (capillary part) and adsorbed liquid, same way as BJH (e.g., (De Burgh et al. 2016;Jiang et al. 2019;Wang et al. 2018)). Other more sophisticated ways within the frame of statistical physics were also used in (Lastoskie et al. 1993a, b) on sorption isotherm of carbon.
From thermodynamic point of view, the ensemble of liquid molecules in a pore at a given relative humidity and temperature is a Grand Canonical ensemble. Said otherwise, temperature T , volume V and chemical potential µ of the system are given. Thermodynamic equilibrium state of such a grand canonical ensemble can be simulated by minimizing the grand potential Ω(T , V , µ) (a measurement of energy, similar to Gibbs free energy or enthalpy) of the system in molecular dynamics or by Grand Canonical Monte Carlo (GCMC) simulation. Lastoskie et al. (1993a, b) proposed the non-local density functional theory (NLDFT) by considering interaction between liquid-liquid and liquid-solid molecules. They compared the density function ρ L (φ, h r ) as well as pore size distribution f (φ) obtained from NLDFT with those obtained by Kelvin-Laplace equation, BJH model and GCMC. It was acknowledged in Lowell (2005) that NLDFT and GCMC are the most accurate methods for pore size analysis of very narrow pores. As for cementitious materials, Wenzel et al. (2017) estimated pore size distribution of C-S-H from argon sorption isotherm using NLDFT and compared with TEM observation results. GCMC were used to estimate amount of adsorption in virtual 3D microstructure of C-S-H, for noble gases in Pellenq and Levitz (2002) and for water in Ioannidou et al. (2016). However, these methods are not widely used in pore size analysis of cement paste, probably due to the difficulty in identification of interaction forces between C-S-H layers and water molecules.

Proton Nuclear Magnetic Resonance Relaxometry
Proton nuclear magnetic resonance relaxometry is a nondestructive technique (Scrivener et al. 2016) to investigate the confined state of water in pores of cement paste. The 1 H-NMR was applied on cement paste for the first time in Japan in Kawachi et al. (1955) to study the ratio between chemically bounded water and free water as a function of hydration. Later, the relaxation time was measured to study evaporable water of cement paste in Blinc et al. (1978). Since then, many authors (for instance (Apih et al. 2001;Barbic et al. 1982;Blinc et al. 1988;Halperin et al. 1994;Jehng et al. 1996;McDonald et al. 2005;Schreiner et al. 1985;Valori et al. 2010)) contributed to the progress of the study of evaporable water state in cement paste by relaxation measurement of 1 H-NMR. General review of application of 1 H-NMR in cement paste can be found in the review papers by Valori et al. (2013). Now, experimental protocol of relaxation measurement and data analysis of 1 H-NMR experiment have been established and well documented (see chapter 7 in Scrivener et al. (2016)). In natural state, the proton ( 1 H) like all other atomic nuclei has a spin-spin rotation moment. The directions of the spin of the proton are randomly distributed in the water inside of the sample. In 1 H-NMR, external magnetic pulses are applied to align spin directions inducing magnetic field. Once the external magnetic field is removed, it takes certain time for the spin directions to gradually go back to the original random state, which can be measured by the decay of the signal intensity of the induced magnetic field. From the decay of signal, T 2 relaxation (spin-spin relaxation) time can be calculated. This T 2 relaxation time depends on the confining state of water, from which we can infer information about pore size.
In 1 H-NMR experiment, no pretreatment is required and sample can be used as prepared. For cement paste, generally, two series of magnetic fields are applied to measure the relaxation times: Quadrature Echo (QE) sequence (or, called as solid-echo (SE)) for chemically bounded water and CPMG pulse sequence for evaporable water. For the details, readers can refer to chapter 7 in Scrivener et al. (2016). This section focuses on the evaluation pore size and specific surface area from the typical results of 1 H-NMR experiment.
In CPMG measurement, the decay of signal intensity of induced magnetic field is deconvoluted and four peaks are observed at different T 2 relaxation times. Figure 7 (top) shows T 2 relaxation time of a matured white Portland cement with water-cement ratio 0.55 (data from ). Four peaks observed in Fig. 7 correspond to (from left to right) interlayer water, gel water, interhydrate water and capillary water (Muller et al. 2013). Size of the peak corresponds to the amount of water. Knowing the amount of chemically bounded water from QE sequence, volume of water in each type of pores can be computed. Two models are proposed in literature for estimation of the pore size: relaxivity method and amplitude method.
The relaxivity method is known as fast exchange model as it is based on "fast exchange" hypothesis: to relax the spin alignment, a water molecule in the bulk of a pore has to exchange first its position with a water molecule on surface of pore. Supposing that this exchange process is much faster than the relaxation of spin, relaxation rate in a given pore should be proportional to the ratio S/V of surface over volume (D'Orazio et al. 1990). This hypothesis reads as following: where T bulk 2 and T surf 2 are the T 2 relaxation times of bulk water and surface adsorbed water, respectively; ε is the thickness of monolayer adsorbed water, taken to be the height of one water molecule, 0.28 nm. If the pore shape is assumed to be slit, the size (width) of pore can be computed as: Then, combining the pore size with the volume of corresponding type of water, surface area can be computed.
Key factor for applying the fast exchange model is the value of relaxation time T surf 2 of surface adsorbed water, which can be measured by 1 H-NMR experiment on material dried to mono-layer coverage. The difficulty lies in the definition of mono-layer covered dried state. In  Mater (2020) 14:55 addition, the measured relaxation time T surf 2 depends on cement type, mixture proportion and the device of 1 H-NMR. The density of elements that can be easily magnetized has also a large impact on the relaxation time T surf 2 . For example, Halperin et al. (1994) measured the surface relaxation time on a sample with water-cement ratio 0.43 which was dried to lower than 20% of saturation degree and obtained 0.04 ms. Bohris et al. (1998) reported a value of 0.01 ms on cement paste with water-cement ratio 0.3 ms and 0.5 ms; whereas, 0.085 ms was measured on mortar sample by (Valckenborg et al. (2001). For white Portland cement with water-cement ratio of 0.4, Muller et al. (2013)  The amplitude model, proposed in McDonald et al. (2010), is based on the amplitude of NMR signals rather than the relaxation time. In this model C-S-H is considered as layered structure. When pores are saturated with water, it is assumed that all water molecules are in fast exchange and only one relaxation time is observed. When the sample is dried, two relaxation times are supposed: one from water on the pore surface, which has lower relaxation time but increasing amplitude with increase of drying; the other from bulk water molecules, which has same relaxation time as saturated state but decreasing amplitude with increase of drying. Therefore, 1 H-NMR experiments should be performed at different stages of drying. McDonald et al. (2010) showed that the pore sizes from amplitude model were in a good agreement with the fast exchange model. Given the fact that more assumption and experiments are needed in amplitude model comparing to fast exchange model, only the fast exchange model will be used in discussion section.

Small-Angle Scattering
Small-angle X-ray scattering was used to determine specific surface area of cement paste for the first time by (Winslow and Diamond 1974). Since then, both SAXS and SANS were widely applied to characterize complex microstructure of cementitious materials (Allen et al. 2007;Allen and Livingston 1998;Allen and Thomas 2007;Bogner et al. 2020;Maruyama et al. 2017b;Thomas et al. 1998Thomas et al. , 1999). An incident beam of either X-ray (in SAXS) or neutron (in SANS) is passed through the material. Due to the interaction between microstructure of the material and the incident beam, small component is scattered out from the direction of the incident beam. The intensity of X-ray or neutron is measured as a function of scattering angle. In SAXS, X-rays interact with outer shell electrons; whereas, in SANS, the neutrons interacts with atomic nuclei. For porous materials such as cement paste, X-ray probes the solid-pore interface where the scattering contrast is maximum. Measured X-ray scattering pattern I(q) (i.e., intensity pattern) is a form of Fourier transform of the solid-pore interface (Brisard et al. 2019;Thomas et al. 1999). Allen et al. (Allen 1991) concluded that the C-S-H gel is the primary source of the scattering. One of the main advantages of SAXS and SANS is that no pretreatment is required and samples can be measured as prepared. Generally, we have to use block samples but not powder, because surface roughness and particle size will add some scattering signal on the data. The range of microstructure that can be measured by SAXS and SANS is wide, up to 3 orders of magnitude (Maruyama 2017b). However, the main difficulty remains on data analysis and exploration of the results. In the following, methods of surface area calculation are presented only for SAXS.
SAXS results are expressed as the relation between scattering intensity I (expressed usually in m −1 ) and scattering vector q (expressed in nm −1 ), which is obtained by q = 4π sin (ϕ/2)/ , where ϕ is half of the scattering angle and is the wave length of X-ray. The relation I(q) follows a power-law in a limited range of q for many materials. In Porod regime, i.e., for instance such as 0.1 nm −1 < q < 1 nm −1 (Porod et al. 1982), the scattering intensity is proportional to the volumetric specific surface area S v (expressed in m 2 /m 3 ). In practice, a linear relationship can be fit between Iq 4 and q 4 , to obtain the volumetric specific surface area: where �ρ e is the scattering length contrast between solid and pore water (expressed in [m −2 ]), i.e., between C-S-H and water in the case of cement paste. There are two major issues in applying the Eq. (11): the range of Porod regime and the value of scattering length contrast.
The range of Porod regime is affected by the typical size of particles of interest that scatter the X-ray. Originally, it was proposed that the range of for Porod regime is in between 0.1 nm −1 and 1 nm −1 (Porod et al. 1982). Later, Thomas et al. (1998), by comparing the coefficient of correlation for different range of q on SANS data measured on cement paste, suggested that Porod regime is in between 1.4 nm −1 and 2 nm −1 for C-S-H.
Regarding the scattering length contrast, its value can be calculated from the density and chemical composition of the particles. Noting that number of electrons in one molecule as n e , density as ρ , electron density would be ρn e N A /M w (expressed as number of electron per volume, m −3 ). Multiplying this electron density by Thomson (11) Iq 4 = 2π |�ρ e | 2 S v + Cq 4 , Page 11 of 18 Aili and Maruyama Int J Concr Struct Mater (2020) 14:55 radius r T = 2.82 × 10 −15 m , X-ray scattering length reads as ρ e = ρn e N A r T /M w (expressed in m −2 ). Since the exact chemical composition of C-S-H is not fully defined, calculating the X-ray scattering length causes problems. Moreover, the impact of the scattering length contrast in Eq. (8) is significant. For instance, scattering length of C 1.7 SH 1.5 is 12% higher compared to that of C 1.7 SH 2.5 . In the above Porod regime method, microstructure is obtained from scattering. The inverse problem, i.e., from microstructure to scattering is less ambiguous and analytical fractal models with specified shape for building blocks of C-S-H were proposed in literature. For instance, supposing spheroid building block, Thomas et al.(Pearson et al. 1983) reproduced appropriately the SANS intensity pattern in the range of q lower than Porod regime. Based on the colloidal model of C-S-H (CM-II) proposed by Jennings et al. (2008), Chiang et al. (2012) assumed disk-shape building blocks with layered structure inside and developed fractal disk model to evaluate globule size in the CM-II model. Knowing the average diameter φ d and thickness t d of disk from the fractal disk model, surface area of C-S-H building block per volume of amorphous C-S-H gel can be estimated. Neglecting surface area of other crystalline hydrates and unreacted clinkers, and if the volume fraction of amorphous C-S-H gel χ can be estimated from other techniques, the volumetric surface area can be calculated as: It is worth noting that the surface area calculated from Eq. (8) differs from that of Porod regime, due to (1) the range of q is lower meaning that the scattering particle size are larger; (2) only surface of building block is taken into account and, hence, interlayer space is neglected. That is to say, Eq. (8) reflects more the surface area of gel pores; whereas Eq. (11) is more corresponding to the interlayer space.
An example of SAXS intensity profile is displayed in Fig. 8. The experimental data are taken from Maruyama et al. (2017b) where the authors measured SAXS on mature white Portland cement paste with water-cement ratio 0.5, kept under water for 1.5 years. Equation (11) is fitted for 3.5 nm −1 < q < 5 nm −1 which corresponds to size of 0.2 nm − 0.3 nm . The illustrated region for fractal disk model is also taken from Maruyama et al. (2017b).

Discussion
This section is dedicated to compare different methods based on various experimental results previously published by the authors (e.g., (Maruyama et al. 2014, (12) Maruyama et al. 2017b;Sugimoto et al. 2017)). The samples were prepared from white Portland cement, with water-cement ratio 0.55. Except from NMR samples and sorption measurement on 3-day-old samples, all other samples were from the same batch. All samples were stored in lime-saturated water until age of 6 months, and then stored in different relative humidity chambers for at least 1 year before test.

Comparison of Different Pretreatment on Sorption Isotherm Test
Like MIP, sorption isotherm results are affected by the pretreatment method. Different pretreatment methods had been compared for MIP, for instance, Galle (2001). In Muller and Scrivener (2017), comparing the MIP pore size distribution with 1 H-NMR results, solvent exchange was suggested as the best method to keep the original microstructure of the sample. Similarly, Bogner et al. (2020) also concluded that solvent exchange and then drying at 105 °C is the best way to keep the original microstructure of the sample for MIP and nitrogen sorption measurement. However, the influence of pretreatment on the sorption isotherm might be different, since the pores that are seen in the sorption measurement are much smaller than those seen by MIP. In Sugimoto et al. (2017), different pretreatment methods were applied on both nitrogen sorption isotherm and water vapor sorption isotherm: drying under vacuum at 20 °C during 3 h or at 105 °C during 30 min; solvent exchange with isopropanol or acetone, and then vacuum dry at 20 °C during 3 h or at 105 °C during 30 min. Freeze drying for one week was applied also on nitrogen sorption samples. In each method, the vacuum drying time, either at 20 °C or at 105 °C, should depend on the sample amount. For instance, in this series of test (Sugimoto Fig. 8 Typical SAXS profile measured on white Portland cement, with indication of Porod regime and Fractal disk region (experimental data from ).
Page 12 of 18 Maruyama Int J Concr Struct Mater (2020) 14:55 et al. 2017), drying time that kept best the original microstructure of the sample was found to be 1 h for 25 ± 1 mg of sample, 3 h for 75 ± 5 mg of sample. Considering microstructural change with age of specimen, sorption isotherms were measured on both early age (3-day-old) and mature (5-year-old) cement pastes. Using Eq. (3), BET surface area was calculated. The results are compared in Fig. 9 for nitrogen sorption and in Fig. 10 for water vapor sorption. From Figs. 9 and 10, it can be seen that vacuum drying at 105 °C provided the highest surface are for both nitrogen and water vapor BET surface area. However, it is known that ettringite and monocarbonate lose water when dried at temperatures even lower than 100 °C (see chapter 5.3 in (Scrivener et al. 2016)). Vacuum drying at 105 °C may reduce number of bounded water molecules per mole of ettringite from 32 to approximately 10-13 (Zhou and Glasser 2001). Rymeš and Maruyama (2018) compared water vapor sorption isotherm on vacuumdried sample at 20 °C and 105 °C. They showed that the total sorption amount and desorption branch of sorption isotherm changed due to decomposition of ettringite. Therefore, stability of hydrates should also be taken into consideration when comparing vacuum drying at 20 °C and 105 °C.
As for solvent exchange method, another drawback is that these solvents are strongly adsorbed on cement paste (Beaudoin et al. 1998;Thomas 1989) and could cause some change in phase composition of cement paste. For instance, it was observed by thermogravimetric (TG) analysis that isopropanol has reacted with portlandite (Beaudoin et al. 1998(Beaudoin et al. , 2000Mikhail and Selim 1966;Zhang and Scherer 2011) As for acetone, similar conclusion was achieved also by TG in Beaudoin et al. (1998Beaudoin et al. ( , 2000, Taylor and Turner (1987), Zhang and Scherer (2011).

Comparison of Pore Size Distribution from Different Methods
Nitrogen sorption and water vapor sorption see different parts of the pore structure. As seen in Sect. 2.3, different theories can be applied to analyze pore size form a given sorption isotherm. Besides, the interpretation of the pore size from sorption isotherm and NMR is based on different hypothesis. This section aims at comparing pore size from different experiments and different analysis methods. In Fig. 11, pore size distributions of mature white Portland cement pastes from various methods are compared. All samples are stored in lime-saturated water before test. From nitrogen sorption data (see Fig. 3), we calculated the pore size distribution using three different methods, all with slit pore assumption: (1) BJH model, with Harkins-Jura model (Eq. (7)) as adsorption thickness model;   Fig. 11 Comparison of pore size distribution from different techniques. Sorption isotherm data from Maruyama et al. (2014), NMR data from . The volume of pores obtained from nitrogen sorption should be read from left Y-axis, while pore volume from water vapor sorption and NMR should be read from right Y-axis.
Page 13 of 18 Aili and Maruyama Int J Concr Struct Mater (2020) 14:55 (2) GCMC simulation, (3) NLDFT, both using BELMas-ter7 software (MicrotracBEL. Corp) and properties of carbon for solid-liquid interaction parameters. Pore size distribution is also computed from water vapor sorption isotherm (see Fig. 3) using BJH model with slit pore shape and BET model (Eq. (2)) as the adsorption thickness model. These results are compared with NMR pore size distribution, calculated from fast exchange model (Eq. (10), surface relaxation time taken as 0.176 ms). It can be seen from Fig. 11 that pore size distributions obtained from the same experimental results (here, nitrogen sorption isotherm) could be different. BJH method suggests no peaks in its range; while, both GCMC and NLDFT simulation suggest three peaks. The three peaks of GCMC are at around 1.2 nm, 3 nm and 10 nm, respectively; whereas those of NLDFT are at around 1.2 nm, 3 nm and 5 nm, respectively. The range of pore size obtained from water vapor sorption is smaller comparing to those from nitrogen sorption and NMR, which might be due to the difficulty in establishing adsorption thickness model adequately. In spite of the difference, pore size distributions from different methods show some consistency in Fig. 11. The pore size 7 nm of interhydrate water from NMR is in between the largest size from GCMC and that from NLDFT. The pore size of 3 nm of gel water from NMR corresponds well to the middle peak of both GCMC and NLDFT. In addition, the size of interlayer space from NMR corresponds approximately to the smallest peak position from GCMC and NLDFT of nitrogen sorption isotherm, as well as the first peak of the water vapor sorption isotherm. However, this consistency of interlayer size should be interpreted carefully, since it was believed that nitrogen molecules cannot enter interlayer space and therefore, information about interlayer space size cannot not be obtained from the nitrogen sorption isotherm. As a conclusion, the calculated pore size distribution of cement paste is affected not only by the experimental technique itself but also by the interpretation method of experimental results. Combining multiple techniques would be useful to get a more comprehensive understanding of the pore structure.

Comparison of Surface Area from Different Methods
Similar to pore size distribution, surface areas measured by different techniques differ from each other, depending on the experimental technique and analysis method. Objective of this section is to compare surface area obtained from different methods. Maruyama et al. (2014) concluded that the microstructure of C-S-H alternates due to the long-term slowly drying. Since the surface area of cement paste originates mainly from C-S-H, the equilibrium relative humidity of the sample (i.e., relative humidity of the storage chamber) is chosen as a comparison index for the surface area. Figure 12 displays the surface area of white Portland cement paste with water-cement ratio 0.55, obtained from different techniques. Nitrogen surface area S N 2 and water vapor surface area S H 2 O are computed using BET theory from the experimental results in Fig. 3. With the experimental protocol and used devices of the laboratory, the relative error of nitrogen surface area S N 2 and water vapor surface area S H 2 O were around ± 2.5%. Combining the simulation results of SAXS-Fractal Disk model of Maruyama et al. (2017b) and volumetric composition of cement paste obtained from XRD-Rietveld analysis, inserting them into Eq. (12), surface area per volume is calculated. To convert this surface area per volume to surface area per mass of dried hardened cement paste, dry density of cement paste is assumed to be equal to 2.1 g/cm 3 for all equilibrium relative humidity. It is worth noting that the scattering profile variation was within 3% and that the errors related to fitting of disk diameter and thickness, given in Fig. 6 of Maruyama et al. (2017b), were less than 3%. Also using the SAXS data of Maruyam et al. (2017b), we calculated surface area per volume with Porod law Eq. (11). The range of scattering vector is chosen to be in between 3 nm −1 and 5 nm −1 to correspond the size of C-S-H layers. Scattering length contrast is computed from chemical composition of C-S-H supposing that number of water molecules (number of H in C-S-H) as a linear function of relative humidity, reducing from 4 in saturated state to 1.5 in 11% of equilibrium relative humidity. As for 1 H-NMR, surface area is calculated for each type of pores based on the fast exchange model with slit shape of pores. Only saturated sample is analyzed, since when sample is dried, 1 H-NMR cannot Fig. 12 Change of the surface area as a function of equilibrium relative humidity, obtained from different techniques: sorption isotherm data from Maruyama et al. (2014), SAXS data from Maruyama et al. (2017b), NMR data from .
Page 14 of 18 Aili and Maruyama Int J Concr Struct Mater (2020) 14:55 see those dried pores. With the repetition number of solid-echo and CPMG signal measurement in the test, the cumulated error of 1 H-NMR surface area related to quadrature echo and CPMG measurement and that to relaxation time T 2 obtained by deconvolution algorithm ILT (Inverse Laplace Transform) was around ± 5%. It can be seen from Fig. 12 that the surface areas from different method differ from each other significantly. As expected, nitrogen surface area is lower than water vapor surface area, probably due to the fact that nitrogen molecules cannot enter interlayer space. The surface area obtained from SAXS, both fractal disk model and Porod regime, are significantly higher, which is consistent with conclusion in Thomas et al. (1999). 1 H-NMR data show surface area higher than sorption data. For the sample stored in lime water, surface area of interlayer space (204 m 2 /g) alone from 1 H-NMR is almost equal the water vapor surface area (196 m 2 /g), while surface area of gel pore (118 m 2 /g) is similar to nitrogen surface area (113 m 2 /g). If we assume that the difference of batch is neglectable, the fact that interlayer surface area (204 m 2 /g) from 1 H-NMR is much higher than the difference between water vapor sorption isotherm and nitrogen sorption isotherm ((196-113 = 83 m 2 /g)) suggests that interlayer space may not be the only difference of water vapor sorption and nitrogen sorption isotherms. Though there is discrepancy between surface area measured by different techniques, all of them are able to detect the decrease of surface area due to the microstructural change of C-S-H structure during long-term drying. Hence, it would be useful to combine various methods to characterize microstructure of cement paste.

Conclusion
This article reviewed three experimental methods that are used for analysis of pore size distribution and specific surface area of cement based materials: sorption isotherm, proton nuclear magnetic resonance relaxometry and small-angle scattering. Different pretreatment methods are compared for sorption isotherm, for both early age and mature cement paste. We also compared pore size distribution and specific surface area of the same cement paste, using data of different experimental method from literature and analyzed difference and consistency between different methods. Following conclusions are drawn: • For both early age and mature cement paste, drying under vacuum at 105 °C gives the highest specific surface area in nitrogen sorption and water vapor sorption. However, decomposition of ettringite when dried at 105 °C should be paid attention since the sorption capacity is affected. • Pore size distribution computed using BJH model from sorption isotherm is sensitive to adsorption thickness model. • Pore size distribution computed using GCMC or NLDFT from nitrogen sorption isotherm showed similar peak position as the pore size of gel and interhydrate computed from 1 H-NMR. • Values of specific surface area from sorption isotherm, 1 H-NMR and SAXS differ from each other significantly. However, all of these methods detected reduction of specific surface area due to alternation of C-S-H microstructure under longterm drying. • To be able to compare and/or combine the reviewed experimental techniques in view of characterization of microstructure of cement pastes, our findings suggest that experiments be performed on samples from same casting and curing condition, at the same age with least possible disturbance on the original microstructure. It would be helpful to design full set of testing campaign for various parameters that cause microstructural changes, such as age, temperature, drying, etc.