The authors gratefully acknowledge the financial support from the National Key R&#38;D Program of China (2017YFB0309904), National Natural Science Foundation of China (Grant Nos. 51508090 and 51808188), 973 Program (2015CB655100), State Key Laboratory of High-Performance Civil Engineering Materials (2016CEM005). Also, greatly appreciate Jiangsu Research Institute of Building Science Co., Ltd and the State Key Laboratory of High-Performance Civil Engineering Materials for funding the research project.

1. Preparation of the model concrete with a single ceramic particle
<ol>
	<li>Mold preparation
	<ol>
		<li>Use a brush to clean the mold (25 mm x 25 mm x 25 mm) and ensure that the inner surfaces of the mold are impurity-free.</li>
		<li>Use another brush to uniformly apply diesel oil on the inner surfaces of the mold for easier mold-release. 
		NOTE: Here, we did not use the common mold for mortar or concrete preparation. As the ceramic particle is around 15 mm in diameter, a cubic plastic mold around 30 mm in length is used for sample preparation. Ensure that the size of plastic mold is larger than the ceramic particle.</li>
	</ol>
	</li>
	<li>Molding the model concrete
	<ol>
		<li>Weigh 1,000 g of cement and 350 g of water with an electronic balance (the designed water to cement mass ratio is 0.35).</li>
		<li>Wipe the 5 L mixing pot with a wet towel to moisten it. Add 350 g of water and 1,000 g of cement into the mixing pot sequentially. Place the mixing pot on the mixer and raise it to the stirring position.</li>
		<li>Mix at 65 rpm for 90 s, and let the mixture stand still for 30 s. During this period, scrape off the paste on the inner wall of the pot. Then, mix at 130 rpm for another 60 s.</li>
		<li>Remove the pot from the mixer and put the ceramic particle into the paste, manually mixing it with the cement paste thoroughly.</li>
		<li>Half fill the mold with well mixed fresh cement paste.</li>
		<li>Place the ceramic particle on the top surface of the cement paste and fill the rest of the mold with cement paste. Wipe off excess cement paste with a scraper knife and vibrate the mold on a vibrating table for 1 min at 50 &#177; 3 Hz.</li>
		<li>Seal the mold surface with cling film to prevent moisture evaporation.</li>
	</ol>
	</li>
	<li>Curing
	<ol>
		<li>Cure the specimen in a curing room for 24 h (20 &#177; 1 &#176;C and 95% &#177; 5% relative humidity).</li>
		<li>Remove the specimen from the mold and further cure the specimen for 28 d under the same environmental conditions.</li>
	</ol>
	</li>
</ol>
2. Scanning electron microscopy preparation
<ol>
	<li>Determination of the ceramic particle inside the matrix
	<ol>
		<li>Scan the specimen with X-ray computed tomography to obtain a stack of slices31.</li>
		<li>Roughly choose the slice where the ceramic particle appears to be largest. Fit the boundary of the ceramic particle with a circle and determine the center of the circle as the geometrical center of the ceramic particle. Due to the gray value difference between cement matrix and ceramic particle, a rough boundary of the particle appears on each CT slice (Figure 1).</li>
	</ol>
	</li>
	<li>Cutting
	<ol>
		<li>Cut the cubic specimen into two parts through the geometrical center of the ceramic particle in a cutting machine. Figure 132 is a schematic map showing the cutting direction. 
		NOTE: The ceramic particle was split into two equal parts, while the specimen was not cut into two exactly equal halves. If the ceramic particle is in the exact center of the cubic specimen, the specimen will be split into two equal halves. However, in a real situation, the ceramic particle typically is not in the exact center of the specimen.</li>
	</ol>
	</li>
	<li>Hydration termination
	<ol>
		<li>Immerse the two parts into isopropyl alcohol (&#8805;99.5%) for 3 days at room temperature to remove the unbounded water and terminate the internal hydration process. Replace the isopropyl alcohol solution every 24 h.</li>
		<li>Place the two parts in the vacuum drying even for 7 days to dry the sample at temperature of 40 &#176;C.</li>
	</ol>
	</li>
	<li>Solidifying the microstructure
	<ol>
		<li>Use a finger to smear the inner surface of two cylindrical plastic molds (31 mm in diameter and 25 mm in height) with demolding paste. The molds are all bottom removable.</li>
		<li>Place each piece of the sample into each mold with the surface to be examined facing downward.</li>
		<li>In a paper cup, weigh 50 g of low viscosity epoxy resin and add another 5 g of hardener. Manually stir the mixture with a wooden stick for 2 min.</li>
		<li>Put the mold into the cold mounting machine along with the paper cup with the mixture.</li>
		<li>Start the vacuum on the cold mounting machine and pour the epoxy resin into the mold until it merges with each sample.</li>
		<li>Keep the mold in the cold mounting machine for 24 h until the epoxy resin hardens.</li>
		<li>Remove the bottom of each mold and squeeze out the sample. Store the sample in a vacuum drying oven.</li>
	</ol>
	</li>
	<li>Grinding and polishing
	<ol>
		<li>Grind the sample with SiC paper and alcohol as a lubricant on an automatic polishing machine at the speed of 300 rpm in the following sequence for 3 min each: 180 grit, 300 grit, 600 grit, and 1200 grit.</li>
		<li>Attach the flannelette to the turntable of the automated polishing machine.</li>
		<li>Polish the sample on the flannelette with diamond paste of 3 &#956;m, 1 &#956;m, and 0.25 &#956;m for 15 min at the speed of 150 rpm, each.</li>
		<li>Remove the debris in an ultrasonic cleaner with alcohol as the cleaning solvent after each grinding and polishing step.</li>
		<li>Store each sample in a plastic box of similar size to the sample with each surface to be examined facing up to avoid cause any scratches on the testing surface.</li>
		<li>Keep the boxes containing the samples in a vacuum dry oven32. 
		NOTE: The grinding and polishing process could be completed on an automated polishing machine and at most 6 samples could be polished at the same time. The grinding and polishing time should be carefully chosen to obtain an extremely smooth surface for the SEM without creating height differences between the cement paste and the aggregate. A typical sample is shown in Figure 232.</li>
	</ol>
	</li>
</ol>
3. Backscattered image acquisition and processing
<ol>
	<li>Acquisition
	<ol>
		<li>Spray a thin layer of gold foil on the surface to be examined in a vacuum environment with an automatic sputter coater.</li>
		<li>Place a strip of adhesive tape on the side of the sample to connect the testing surface and opposite surface and place the sample on the test bench with the testing surface facing upward.</li>
		<li>Move the sample to focus the lens on region 1 as labelled in Figure 232.</li>
		<li>Vacuum the SEM and change to backscattered electron mode. Set the magnification at 1,000x and carefully adjust the brightness and contrast before capturing images.</li>
		<li>Move the lens along the direction of the aggregate boundary to another position of the aggregate and take another image. Repeat this moving and imaging process at least 15 times so that enough images can be obtained for statistical analysis.</li>
		<li>Move the lens to region 2 and region 3 and repeat the imaging process. 
		NOTE: Each image should include three phases: the matrix, the aggregate, and the ITZ. Since ITZ is a narrow section existing between another two phases and hard to be distinguished, each image should include both the cement matrix and the aggregate.</li>
	</ol>
	</li>
	<li>Processing
	<ol>
		<li>Pre-treat the image with a best fit and 3 x 3 median filter three times to reduce the noise and enhance the boundary of different phases on ImageJ.</li>
		<li>Manually capture the boundary of the ceramic particle and cut off this part from the original image using ImageJ.</li>
		<li>Roughly determine the upper threshold value of pore phases by setting different threshold values and segmenting the image to compare with the original one.</li>
		<li>Obtain the gray-scale distribution of the remaining part of the image. Choose two approximatively linear parts of the distribution curve just around the roughly determined upper threshold value of pore phases. Fit these two parts with linear curve and the intersection point will be set as the exact upper threshold value of this image (see Figure 3c32).</li>
		<li>Use this value to do the segmentation and compare the binary image with original gray-scale image for final threshold value determination.</li>
		<li>Convert the gray-scale image to a binary image with white (gray value = 255) representing pore phase and black (gray value = 0) representing solid phases. 
		NOTE: The exact determination of threshold value is called the overflow point method33 since the brightness and contrast are kept the same for different images obtained from same sample. Once the upper threshold value is precisely determined, this value could be applied to other images obtained from the same sample.</li>
	</ol>
	</li>
</ol>
4. Data processing
<ol>
	<li>ITZ thickness determination
	<ol>
		<li>Delineate twenty 20 successive strips that are 5 &#956;m in width (use the included strip_delineation.m file), along the captured boundary in the direction of starting from aggregate surface and going into the bulk paste (see Figure 3d32).</li>
		<li>Count the number of pixels with a gray value lower than the threshold in each strip and normalize the values by the number of total voxels contained in each strip. Each normalized value will be viewed as the porosity of each strip.</li>
		<li>Repeat the counting and normalization process for all the images. Average the porosity profiles of the same strip number from different images.</li>
		<li>Draw the porosity distribution graph as a function of distance away from aggregate surface. Determine the inflection point on the curve where porosity become stable as the thickness of the ITZ. 
		NOTE: The number of strips and width of each strip could vary; make sure that the total width of the delineated strips includes all of the ITZ. According to previous research, the ITZ thickness ranges between 20-50 &#956;m13. Even in some model concrete samples with an enlarged ITZ, this value does not exceed 70 &#956;m34,35.</li>
	</ol>
	</li>
	<li>Aggregate surface roughness (SR) characterization
	<ol>
		<li>Save the manually captured boundary as a curve. Fit the irregular boundary with both straight line and circle arc according to the Eq. (1) and Eq. (2) based on least square algorithm. 
		<img align="center" alt="Equation 1" src="/files/ftp_upload/60245/60245eq1.jpg" /> (1) 
		<img align="center" alt="Equation 2" src="/files/ftp_upload/60245/60245eq2.jpg" /> (2) 
		with (a,b) being the center of the fitting circle.</li>
		<li>Define the deviations between original irregular boundary and fitting smooth curve as the surface roughness (SR).</li>
		<li>For straight line, calculate the SRs by averaging the absolute value of the perpendicular distance of the center of each pixel on the boundary to the fitting line: 
		<img align="center" alt="Equation 3" src="/files/ftp_upload/60245/60245eq3.jpg" /> (3) 
		with n being the number of pixels included in each boundary and (xi, yi) being the coordinates of the ith pixel on the boundary.</li>
		<li>For a circle arc, define SRC as: 
		<img align="center" alt="Equation 4" src="/files/ftp_upload/60245/60245eq4.jpg" /> (4)</li>
		<li>Compare the value of SRS and SRC for each boundary and determine the minimum value as the final surface roughness for this curve (use the included surface_roughness_calculation.m file). 
		NOTE: The surface roughness of boundary should be defined against a smooth baseline curve. Both straight line and circle line were used for the following reason. Though the boundary of the spherical ceramic particle appears like a circle in 2D, some local regions appear to be more appeal to a straight line.</li>
	</ol>
	</li>
	<li>K-means clustering
	<ol>
		<li>Define a slope index (SL) to describe the porosity gradient within the interfacial transition zone according to the Eq. (5). 
		<img align="center" alt="Equation 5" src="/files/ftp_upload/60245/60245eq5.jpg" /> (5) 
		where &#966;max is the value of the porosity in the first strip (0 &#956;m to 5 &#956;m) and &#966;min is the value of the porosity in the sixth strip (25 &#956;m to 30 &#956;m).</li>
		<li>Combine the SR and SL of each boundary to be an observation. And for total n boundaries and ITZs, there exist n observations to be saved as a cluster {(SR1,SL1),(SR2,SL2), ... , (SRn,SLn)}.</li>
		<li>Apply a K-means clustering36,37 algorithm (use the included k_means_clustering.m file) to all the observations and subdivide them into 2 clusters: rough and smooth aggregate surface group, respectively.</li>
		<li>Average the porosity distributions of ITZ in rough and smooth cluster, respectively. Compare the average porosity distribution between two clusters. 
		NOTE: Herein, K-means clustering is a method of vector quantization, which is originally used in signal processing and currently widely applied to cluster analysis in data mining. The aim of the method is to subdivide the observations into 2 or more subgroups.</li>
	</ol>
	</li>
</ol>

<ol>
	<li>Scrivener, K. L., Crumbie, A. K., Laugesen, P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+Interfacial+Transition+Zone+(ITZ)+Between+Cement+Paste+and+Aggregate+in+Concrete.">The Interfacial Transition Zone (ITZ) Between Cement Paste and Aggregate in Concrete.</a> Interface Science. 12 (4), 411-421 (2004).</li><li>Scrivener, K. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Backscattered+electron+imaging+of+cementitious+microstructures:+understanding+and+quantification.">Backscattered electron imaging of cementitious microstructures: understanding and quantification.</a> Cement and Concrete Composites. 26 (8), 935-945 (2004).</li><li>Houst, Y. F., Sadouki, H., Wittmann, F. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Influence+of+aggregate+concentration+on+the+diffusion+of+CO2+and+O2.">Influence of aggregate concentration on the diffusion of CO2 and O2.</a> Concrete. , 279-288 (1993).</li><li>Halamickova, P., Detwiler, R. J., Bentz, D. P., Garboczi, E. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Water+permeability+and+chloride+ion+diffusion+in+portland+cement+mortars:+Relationship+to+sand+content+and+critical+pore+diameter.">Water permeability and chloride ion diffusion in portland cement mortars: Relationship to sand content and critical pore diameter.</a> Cement &amp; Concrete Research. 25 (4), 790-802 (1995).</li><li>Yang, Z., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=In-situ+X-ray+computed+tomography+characterisation+of+3D+fracture+evolution+and+image-based+numerical+homogenisation+of+concrete.">In-situ X-ray computed tomography characterisation of 3D fracture evolution and image-based numerical homogenisation of concrete.</a> Cement and Concrete Composites. 75, 74-83 (2017).</li><li>Skar&#380;y&#324;ski, &. #. 3. 2. 1. ;., Nitka, M., Tejchman, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Modelling+of+concrete+fracture+at+aggregate+level+using+FEM+and+DEM+based+on+X-ray+&#181;CT+images+of+internal+structure.">Modelling of concrete fracture at aggregate level using FEM and DEM based on X-ray &#181;CT images of internal structure.</a> Engineering Fracture Mechanics. 147, 13-35 (2015).</li><li>K&#246;nigsberger, M., Pichler, B., Hellmich, C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Micromechanics+of+ITZ-Aggregate+Interaaction+in+Concrete+Part+II:+Stength+Upscaling.">Micromechanics of ITZ-Aggregate Interaaction in Concrete Part II: Stength Upscaling.</a> Journal of the American Ceramic Society. 97 (2), 543-551 (2014).</li><li>Shahbazi, S., Rasoolan, I. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Meso-scale+finite+element+modeling+of+non-homogeneous+three-phase+concrete.">Meso-scale finite element modeling of non-homogeneous three-phase concrete.</a> Case Studies in Construction Materials. 6, 29-42 (2017).</li><li>Ak&#231;ao&#287;lu, T., Tokyay, M., &#199;elik, T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Assessing+the+ITZ+microcracking+via+scanning+electron+microscope+and+its+effect+on+the+failure+behavior+of+concrete.">Assessing the ITZ microcracking via scanning electron microscope and its effect on the failure behavior of concrete.</a> Cement and Concrete Research. 35 (2), 358-363 (2005).</li><li>Chang, H., Feng, P., Lyu, K., Liu, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+novel+method+for+assessing+C-S-H+chloride+adsorption+in+cement+pastes.">A novel method for assessing C-S-H chloride adsorption in cement pastes.</a> Construction &amp; Building Materials. 225, 324-331 (2019).</li><li>Wang, P., Jia, Y., Li, T., Hou, D., Zheng, Q. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Molecular+dynamics+study+on+ions+and+water+confined+in+the+nanometer+channel+of+Friedel's+salt:+structure+dynamics+and+interfacial+interaction.">Molecular dynamics study on ions and water confined in the nanometer channel of Friedel's salt: structure dynamics and interfacial interaction.</a> Physical Chemistry Chemical Physics. 20, 27049-27058 (2018).</li><li>Ma, H., Li, Z. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+Multi-Aggregate+Approach+For+Modeling+The+Interfacial+Transition+Zone+In+Concrete.">A Multi-Aggregate Approach For Modeling The Interfacial Transition Zone In Concrete.</a> ACI Materials Journal. 111 (2), (2014).</li><li>Yun, G., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Characterization+of+ITZ+in+ternary+blended+cementitious+composites:+Experiment+and+simulation.">Characterization of ITZ in ternary blended cementitious composites: Experiment and simulation.</a> Construction &amp; Building Materials. 41 (2), 742-750 (2013).</li><li>Garboczi, E. J., Bentz, D. P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=In+Digital+simulation+of+the+aggregate-cement+paste+interfacial+zone+in+concrete.">In Digital simulation of the aggregate-cement paste interfacial zone in concrete.</a> International Conference on Electric Information and Control Engineering (ICEICE), 2011. , 196-201 (2011).</li><li>Winslow, D. N., Cohen, M. D., Bentz, D. P., Snyder, K. A., Garboczi, E. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Percolation+and+pore+structure+in+mortars+and+concrete.">Percolation and pore structure in mortars and concrete.</a> Cement &amp; Concrete Research. 24 (1), 25-37 (1994).</li><li>Sim&#245;es, T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> Mechanical Characterization of Fiber/Paste and Aggregate/Paste Interfaces (ITZ) in Reinforced Concrete with Fibers. , (2018).</li><li>Xiao, J., Li, W., Sun, Z., Lange, D. A., Shah, S. P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Properties+of+interfacial+transition+zones+in+recycled+aggregate+concrete+tested+by+nanoindentation.">Properties of interfacial transition zones in recycled aggregate concrete tested by nanoindentation.</a> Cement and Concrete Composites. 37, 276-292 (2013).</li><li>Bentz, D. P., Garboczi, E. J., Stutzman, P. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Computer+Modelling+of+the+Interfacial+Transition+Zone+in+Concrete.">Computer Modelling of the Interfacial Transition Zone in Concrete.</a> Interfaces in Cementitious Composites. , 107-116 (1993).</li><li>Kai, L., Wei, S., Changwen, M., Honglei, C., Yue, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Quantitative+characterization+of+pore+morphology+in+hardened+cement+paste+via+SEM-BSE+image+analysis.">Quantitative characterization of pore morphology in hardened cement paste via SEM-BSE image analysis.</a> Construction &amp; Building Materials. 202, 589-602 (2019).</li><li>Ondracek, G., Underwood, E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Quantitative+stereology.">Quantitative stereology.</a> Journal of Nuclear Materials. 42 (2), 237-237 (1972).</li><li>Xu, J., Wang, B., Zuo, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Modification+effects+of+nanosilica+on+the+interfacial+transition+zone+in+concrete:+A+multiscale+approach.">Modification effects of nanosilica on the interfacial transition zone in concrete: A multiscale approach.</a> Cement and Concrete Composite. 81, 1-10 (2017).</li><li>Zhu, Z., Chen, H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> Overestimation of ITZ thickness around regular polygon and ellipse aggregate. , 205-218 (2017).</li><li>Head, M. K., Wong, H. S., Buenfeld, N. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Characterising+aggregate+surface+geometry+in+thin-sections+of+mortar+and+concrete.">Characterising aggregate surface geometry in thin-sections of mortar and concrete.</a> Cement and Concrete Research. 38 (10), 1227-1231 (2008).</li><li>Gao, Y., De Schutter, G., Ye, G., Tan, Z., Wu, K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+ITZ+microstructure,+thickness+and+porosity+in+blended+cementitious+composite:+Effects+of+curing+age,+water+to+binder+ratio+and+aggregate+content.">The ITZ microstructure, thickness and porosity in blended cementitious composite: Effects of curing age, water to binder ratio and aggregate content.</a> Composites Part B: Engineering. 60, 1-13 (2014).</li><li>Erdem, S., Dawson, A. R., Thom, N. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Influence+of+the+micro-+and+nanoscale+local+mechanical+properties+of+the+interfacial+transition+zone+on+impact+behavior+of+concrete+made+with+different+aggregates.">Influence of the micro- and nanoscale local mechanical properties of the interfacial transition zone on impact behavior of concrete made with different aggregates.</a> Cement and Concrete Research. 42 (2), 447-458 (2012).</li><li>Elsharief, A., Cohen, M. D., Olek, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Influence+of+aggregate+size,+water+cement+ratio+and+age+on+the+microstructure+of+the+interfacial+transition+zone.">Influence of aggregate size, water cement ratio and age on the microstructure of the interfacial transition zone.</a> Cement &amp; Concrete Research. 33 (11), 1837-1849 (2003).</li><li>Pan, T., Tutumluer, E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Quantification+of+Coarse+Aggregate+Surface+Texture+Using+Image+Analysis.">Quantification of Coarse Aggregate Surface Texture Using Image Analysis.</a> Journal of Testing &amp; Evaluation. 35 (2), 177-186 (2006).</li><li>Erdogan, S. T., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Three-dimensional+shape+analysis+of+coarse+aggregates:+New+techniques+for+and+preliminary+results+on+several+different+coarse+aggregates+and+reference+rocks.">Three-dimensional shape analysis of coarse aggregates: New techniques for and preliminary results on several different coarse aggregates and reference rocks.</a> Cement &amp; Concrete Research. 36 (9), 1619-1627 (2006).</li><li>Santos, B. O., Valen&#231;a, J., Fowler, D. W., Saleh, H. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Livings+patterns+on+concrete+surfaces+with+biological+stains+using+hyperspectral+images+processing.">Livings patterns on concrete surfaces with biological stains using hyperspectral images processing.</a> Structural Control and Health Monitoring. , (2019).</li><li>Santos, B. O., Valen&#231;a, J., J&#250;lio, E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=In+Classification+of+biological+colonization+on+concrete+surfaces+using+false+colour+HSV+images,+including+near-infrared+information.">In Classification of biological colonization on concrete surfaces using false colour HSV images, including near-infrared information.</a> Optical Sensing and Detection V, International Society for Optics and Photonics. , 106800 (2018).</li><li>Stock, S. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Recent+advances+in+X-ray+microtomography+applied+to+materials.">Recent advances in X-ray microtomography applied to materials.</a> International Materials Reviews. 53 (3), 129-181 (2013).</li><li>Lyu, K., Garboczi, E. J., She, W., Miao, C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+effect+of+rough+vs.+smooth+aggregate+surfaces+on+the+characteristics+of+the+interfacial+transition+zone.">The effect of rough vs. smooth aggregate surfaces on the characteristics of the interfacial transition zone.</a> Cement and Concrete Composites. 99, 49-61 (2019).</li><li>Wong, H. S., Head, M. K., Buenfeld, N. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Pore+segmentation+of+cement-based+materials+from+backscattered+electron+images.">Pore segmentation of cement-based materials from backscattered electron images.</a> Cement &amp; Concrete Research. 36 (6), 1083-1090 (2006).</li><li>Liao, K. -. Y., Chang, P. -. K., Peng, Y. -. N., Yang, C. -. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+study+on+characteristics+of+interfacial+transition+zone+in+concrete.">A study on characteristics of interfacial transition zone in concrete.</a> Cement and Concrete Research. 34 (6), 977-989 (2004).</li><li>Barnes, B. D., Diamond, S., Dolch, W. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+contact+zone+between+portland+cement+paste+and+glass+&#8220;aggregate&#8221;+surfaces.">The contact zone between portland cement paste and glass &#8220;aggregate&#8221; surfaces.</a> Cement &amp; Concrete Research. 8 (2), 233-243 (1978).</li><li>Hamerly, G., Elkan, C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Alternatives+to+the+k-means+algorithm+that+find+better+clusterings.">Alternatives to the k-means algorithm that find better clusterings.</a> Proceedings of the eleventh international conference on Information and knowledge management, ACM. , 600-607 (2002).</li><li>Celebi, M. E., Kingravi, H. A., Vela, P. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> A comparative study of efficient initialization methods for the k-means clustering algorithm. , 200-210 (2013).</li><li>Lu, Y., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Three-dimensional+mortars+using+real-shaped+sand+particles+and+uniform+thickness+interfacial+transition+zones:+Artifacts+seen+in+2D+slices.">Three-dimensional mortars using real-shaped sand particles and uniform thickness interfacial transition zones: Artifacts seen in 2D slices.</a> Cement and Concrete Research. , (2018).</li><li>Gao, Y., De Schutter, G., Ye, G., Huang, H., Tan, Z., Wu, K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Porosity+characterization+of+ITZ+in+cementitious+composites:+Concentric+expansion+and+overflow+criterion.">Porosity characterization of ITZ in cementitious composites: Concentric expansion and overflow criterion.</a> Construction and Building Materials. 38, 1051-1057 (2013).</li><li>Celebi, M. E., Kingravi, H. A., Vela, P. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+comparative+study+of+efficient+initialization+methods+for+the+k-means+clustering+algorithm.">A comparative study of efficient initialization methods for the k-means clustering algorithm.</a> Expert Systems with Applications. 40 (1), 200-210 (2013).</li></ol>

We hereby confirm that this manuscript is our original work and all the authors listed have approved the manuscript and have no interest conflicts on this paper.

The X-CT technique was applied to roughly determine the geometrical center of the ceramic particle to ensure that the analyzed surface is through the equator of the particle. Thus, the overestimation of the ITZ thickness caused by the 2D artifacts could be avoided38. Herein, the accuracy of obtained results is highly dependent on the flatness of the examined surfaces. Generally, a longer grinding and polishing time contributes to an adequately smooth surface for testing. However, due to the varyin...

At the mesoscopic scale, cement-based materials can be regarded as a three-phase composite comprised of the cement paste, the aggregate, and the interfacial transition zone (ITZ) between them1,2. The ITZ is often treated as a weak link since its increased porosity could act as channels for the ingress of aggressive species3,4 or provide easier pathways for crack growth5,6,7,8,9,10,11. Subsequently, it is of great interest to precisely characterize the properties of the ITZ to evaluate and predict the macro performance of the cement-based materials.
To investigate the ITZ, there has been excessive research on its microstructural features, forming mechanisms, and influencing factors12,13,14 using both experimental and numerical methods. Various techniques have been coupled for ITZ characterization including: mechanical tests, transport tests, mercury intrusion porosimetry (MIP) tests15,16 and nano-indentation17. It is widely accepted that the ITZ is mainly caused by the wall effect, as well as water film, micro-bleeding, one side growth, and gel syneresis18.
With the development of digital image processing method (DIP) in the last two decades19, the morphological characteristics of the ITZ (e.g., volume fraction, thickness, and porosity gradient) can be quantitatively determined. Based on examination of the plane sections using scanning electron microscopy (SEM) with a backscattered electron detector (BSE), the three dimensional (3D) features of ITZ can be derived from the 2D results via stereology theory20. Like the SEM-BSE technique, the nano-indentation technique is also based on the examination of polished surfaces, but it more focuses on the elastic modulus of the existing phases21. However, in both SEM-BSE analysis and the nano-indentation test, the ITZ thickness may be overestimated as the examined cross section rarely goes through the normal direction from an aggregate surface22. However, coupling this with fluorescent 3D confocal microscopy, the overestimation of ITZ could be eliminated and a real ITZ porosity and anhydrous cement content could be obtained23.
Previous studies of influencing factors mainly focused on the cement paste, ignoring the role of the aggregate and its surface texture24,25,26. Since the shape and morphological properties of the aggregate have been extensively described based on quantitative analysis of digital slices obtained from SEM or X-ray computed tomography (X-CT)27,28. However, no research focusing on the effect of the aggregate surface texture on the formation of ITZ region has been performed.
Hereby, we present a protocol to investigate the effect of aggregate surface morphology on the ITZ microstructure formation based on quantitative analysis of SEM-BSE images and a K-means clustering algorithm. A model concrete sample was prepared with spherical ceramic particle acting as the coarse aggregate. X-CT was used to roughly determine the relative location of the particle in the opaque cement matrix before halving the sample. After processing to obtained SEM-BSE images, the uneven distribution of ITZ around single aggregate was observed. Also, an index surface roughness (SR) describing the aggregate surface texture at the pixel level was defined. The K-means clustering algorithm, originally used in the area of signal processing and now widely used for image clustering29,30, was introduced to established a relationship between surface roughness (SR) and porosity gradient (SL).

<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<tbody>
		<tr>
			<td>Auto Sputter Coater</td>
			<td>Cressington</td>
			<td>108 Auto/SE</td>
			<td></td>
		</tr>
		<tr>
			<td>Automatic polishing machine</td>
			<td>Buehler</td>
			<td>Phoenix4000</td>
			<td></td>
		</tr>
		<tr>
			<td>Brush</td>
			<td>Huoniu</td>
			<td>3#</td>
			<td></td>
		</tr>
		<tr>
			<td>Cement</td>
			<td>China United Cement Corporation</td>
			<td>P.I. 42.5</td>
			<td></td>
		</tr>
		<tr>
			<td>Cement paste mixer</td>
			<td>Wuxi Construction and Engineering</td>
			<td>NJ160</td>
			<td></td>
		</tr>
		<tr>
			<td>Ceramic particle</td>
			<td>Haoqiang</td>
			<td>&#934;15 mm</td>
			<td></td>
		</tr>
		<tr>
			<td>Cling film</td>
			<td>Miaojie</td>
			<td>65300</td>
			<td></td>
		</tr>
		<tr>
			<td>Cold mounting machine</td>
			<td>Buehler</td>
			<td>Cast N&#39; Vac 1000</td>
			<td></td>
		</tr>
		<tr>
			<td>Conductive tape</td>
			<td>Nissin Corporation</td>
			<td>7311</td>
			<td></td>
		</tr>
		<tr>
			<td>Cup</td>
			<td>Buehler</td>
			<td>20-8177-100</td>
			<td></td>
		</tr>
		<tr>
			<td>Cutting machine</td>
			<td>Buehler</td>
			<td>Isomet 4000</td>
			<td></td>
		</tr>
		<tr>
			<td>Cylindrical plastic mold</td>
			<td>Buehler</td>
			<td>20-8151-100</td>
			<td></td>
		</tr>
		<tr>
			<td>Diamond paste</td>
			<td>Buehler</td>
			<td>00060210, 00060190, 00060170</td>
			<td></td>
		</tr>
		<tr>
			<td>Diesel oil</td>
			<td>China Petroleum</td>
			<td>0#</td>
			<td></td>
		</tr>
		<tr>
			<td>Electronic balance</td>
			<td>Setra</td>
			<td>BL-4100F</td>
			<td></td>
		</tr>
		<tr>
			<td>Epoxy resin</td>
			<td>Buehler</td>
			<td>20-3453-128</td>
			<td></td>
		</tr>
		<tr>
			<td>Hardener</td>
			<td>Buehler</td>
			<td>20-3453-032</td>
			<td></td>
		</tr>
		<tr>
			<td>High precision cutting machine</td>
			<td>Buehler</td>
			<td>2215</td>
			<td></td>
		</tr>
		<tr>
			<td>Image J</td>
			<td>National Institutes of Health</td>
			<td>1.52o</td>
			<td></td>
		</tr>
		<tr>
			<td>Isopropyl alcohol</td>
			<td>Sinopharm</td>
			<td>M0130-241</td>
			<td></td>
		</tr>
		<tr>
			<td>Matlab</td>
			<td>MathWorks</td>
			<td>R2014a</td>
			<td></td>
		</tr>
		<tr>
			<td>Paper</td>
			<td>Deli</td>
			<td>A4</td>
			<td></td>
		</tr>
		<tr>
			<td>Plastic box</td>
			<td>Beichen</td>
			<td>3630</td>
			<td></td>
		</tr>
		<tr>
			<td>Plastic mold</td>
			<td>Youke</td>
			<td>a=b=c=25mm</td>
			<td></td>
		</tr>
		<tr>
			<td>Polished flannelette</td>
			<td>Buehler</td>
			<td>242150, 00242050, 00242100</td>
			<td></td>
		</tr>
		<tr>
			<td>Release agent</td>
			<td>Buehler</td>
			<td>20-8186-30</td>
			<td></td>
		</tr>
		<tr>
			<td>Scanning Electron Microscopy</td>
			<td>FEI</td>
			<td>Quanta 250</td>
			<td></td>
		</tr>
		<tr>
			<td>Scrape knife</td>
			<td>Jinzheng Building Materials</td>
			<td>CD-3</td>
			<td></td>
		</tr>
		<tr>
			<td>SiC paper</td>
			<td>Buehler</td>
			<td>P180, P320, P1200</td>
			<td></td>
		</tr>
		<tr>
			<td>Ultrasonic cleaner</td>
			<td>Zhixin</td>
			<td>DLJ</td>
			<td></td>
		</tr>
		<tr>
			<td>Vacuum box</td>
			<td>Heheng</td>
			<td>DZF-6020</td>
			<td></td>
		</tr>
		<tr>
			<td>Vacuum drying oven</td>
			<td>ZK</td>
			<td>ZK30</td>
			<td></td>
		</tr>
		<tr>
			<td>Vibrating table</td>
			<td>Jianyi</td>
			<td>GZ-75</td>
			<td></td>
		</tr>
		<tr>
			<td>Wooden stick</td>
			<td>Buehler</td>
			<td>20-8175</td>
			<td></td>
		</tr>
		<tr>
			<td>X-ray Computed Tomography</td>
			<td>YXLON</td>
			<td>Y.CT PRECISION S</td>
			<td></td>
		</tr>
	</tbody>
</table>

determination of aggregate surface morphology at the interfacial transition zone (itz)

Here, we present a comprehensive method to illustrate the uneven distribution of the interfacial transition zone (ITZ) around the aggregate and the effect of aggregate surface morphology on the formation of ITZ. First, a model concrete sample is prepared with a spherical ceramic particle in roughly the central part of the cement matrix, acting as a coarse aggregate used in common concrete/mortar. After curing until the designed age, the sample is scanned by X-ray computed tomography to determine the relative location of the ceramic particle inside the cement matrix. Three locations of the ITZ are chosen: above the aggregate, on the side of the aggregate, and below the aggregate. After a series of treatments, the samples are scanned with a SEM-BSE detector. The resultant images were further processed using a digital image processing method (DIP) to obtain quantitative characteristics of the ITZ. The surface morphology is characterized at the pixel level based on the digital image. Thereafter, K-means clustering method is used to illustrate the effect of surface roughness on ITZ formation.

The porosity distribution of ITZ regions above the aggregate, on the side of the aggregate, and below the aggregate are compared and shown in Figure 432. The porosity of the ITZ above the upper surface appears to be smaller than that on the side or above the aggregate, indicating a denser ITZ microstructure, while the ITZ below the aggregate is always the most porous due to micro-bleeding. Figure 432 shows that eve...

Watch this Scientific Journal Video about Determination of Aggregate Surface Morphology at the Interfacial Transition Zone ITZ at JoVE.com

Determination of Aggregate Surface Morphology at the Interfacial Transition Zone ITZ

Hereby, we proposed a protocol to illustrate the effect of aggregate surface morphology on the ITZ microstructure. The SEM-BSE image were quantitatively analyzed to obtain ITZ's porosity gradient via digital image processing and a K-means clustering algorithm was further employed to establish a relationship between porosity gradient and surface roughness.

Determination of Aggregate Surface Morphology at the Interfacial Transition Zone (ITZ)

Here, we present a comprehensive method to illustrate the uneven distribution of the interfacial transition zone (ITZ) around the aggregate and the effect ...

determination-aggregate-surface-morphology-at-interfacial-transition

Research

Education

JoVE Core

JoVE Journal

Civil Engineering

Engineering

Unit 1

Aggregates and Water

SCIENTISTS IN ACTION

7.9K Views.  Southeast University. Hereby, we proposed a protocol to illustrate the effect of aggregate surface morphology on the ITZ microstructure. The SEM-BSE image were quantitatively analyzed to obtain ITZ's porosity gradient via digital image processing and a K-means clustering algorithm was further employed to establish a relationship between porosity gradient and surface roughness.