Strength and Electrical Properties of Cementitious Composite with Integrated Carbon Nanotubes

The main objective of this work was to study the effects of carbon nanotubes (CNTs) on the strength and electrical properties of cement mortar. Molecular dynamic simulations (MDSs) were carried out to determine the mechanical and electrical properties of a cementitious composite and its associated mechanisms. To model the atomic structure of a calcium silicate hydrate (C-S-H) gel, tobermorite 11Å was chosen. Single-walled carbon nanotubes (SWCNTs) embedded in a tobermorite structure were tested numerically. In particular, it was concluded that a piezoelectric effect can be effectively simulated by varying the concentration levels of carbon nanotubes. The deformation characteristics were analyzed by subjecting a sample to an electrical field of 250 MV/m in the z-direction in a simulation box. The results indicated a progressively stronger converse piezoelectric response with an increasing proportion of carbon nanotubes. Additionally, it was observed that the piezoelectric constant in the z-direction, denoted by d33, also increased correspondingly, thereby validating the potential for generating an electrical current during sample deformation. An innovative experiment was developed for the electrical characterization of a cementitious composite of carbon nanotubes. The results showed that the electrostatic current measurements exhibited a higher electric sensitivity for samples with a higher concentration of CNTs.


Introduction
Concrete, a very old material, has been widely used since the 19th century. It is the most commonly used material in the construction of civil engineering structures. Since climatic warming is a hard fact, its use is questioned because of its significant carbon footprint. Today, researchers are working to improve its physical and mechanical characteristics, as well as its durability, by mixing cement paste with different additive powders or nanoparticles. These nanoadded substances are generally used with superplasticizers. One of the most useful nanoadditives these days is carbon nanotubes (CNTs). The exceptional mechanical, electrical, and thermal properties of carbon nanotubes (CNTs) [1][2][3] make them beneficial for applications in various areas of civil engineering. The use of CNTs improves the mechanical properties of concretes, as reported in the literature [4,5]. However, objective of recent research is not only to focus on improving the mechanical properties but also to view CNTs as an adaptive and evolving material. Inducing the piezoresistive effect in cementitious materials to monitor structures can be achieved via the use of conductive fillers such as carbon nanofibers, carbon nanotubes, and graphene [6][7][8][9][10]. Azhari, F. et al. [8] considered cement-based sensors, including carbon nanofibers or carbon nanotubes, to change the resistivity, with good results for sensors with both nanotubes and nanofibers. Zhao, P. et al. [9] showed an improvement in the piezoelectric properties when carbon nanotubes were dispersed in cement-based composites in the range of 0 to 0.9 vol% (volume percentage). Gong et al. 2011 [10] created composites using Portland cement and PZT powders with the inclusion of modified CNTs ranging from 0 to 1.3 vol%. The best outcome was achieved when the concentration of the modified CNTs was 0.6 vol%. The results showed that the highest constant d 33 (see Equation (2) below) was equal to 62 pC/N for a 0.3 vol% CNT.
Electrical conductivity is one of the properties that the microstructure of cement cannot naturally develop. The resistivity of concrete varies from 105 Ω/mm to a range of 1012 Ω/mm between wet and dry concrete, respectively, changing from a semiconductor to an insulator in these two phases, thereby allowing evaporable water to be considered as having an important role in electrical conductivity. The ion concentration in water induces an ionic association that leads to the formation of C-S-H gel, consequently decreasing the mobility of the ions due to the electrical insulation layer in the cement grains. After modeling a concrete sample, the conductivity of the aqueous phase, facilitated by the dissolved ions, was significant. So, the resistivity was very low; however, this resistivity increased with the hydration of the concrete, indicating a proportionality between the evolution of resistivity and the evolution of strength. This variation in resistivity could help predict the performance and thus indicate the age factor of structures [1]. At this point, the modification of concrete to obtain electrical properties is very promising for structural health monitoring applications, for example. In this study, we combined experimental and numerical research to investigate the effects on mechanical and electrical properties caused by the modification of cement with carbon nanotubes. The first part of our study focuses on numerical simulations. In these simulations, the dynamics occur through discrete interactions; in our case, this occurs with atoms, which are called particles. The number of particles used in the simulations in this study varied depending on the number of carbon nanotubes incorporated into the tobermorite structure. These simulations were made possible by the ongoing development of the LAMMPS program for molecular modeling, which is a computational tool that continuously improves its ability to accurately represent the physical, chemical, and biological characteristics of systems composed of interacting molecules and atoms [11].
In addition, innovative experimental tests were conducted to measure the electrical currents in prismatic mortar samples of 40 × 40 × 80 mm with carbon nanotubes in weight percentages (wt.%) by weights of cement of 0.00 wt.%, 0.006 wt.%, and 0.018 wt.%. The effectiveness of CNTs in modifying the properties of concrete is largely determined by their level of dispersion [12]. There is a general consensus that standard mixing or manual mixing methods are insufficient for achieving the uniform dispersion of CNTs [13]. Consequently, various dispersion techniques, such as sonication and CNT functionalization, have been devised to achieve stable and uniform aqueous CNT dispersions [14]. It should be noted, however, that the CNT dispersion obtained through sonication tends to re-agglomerate over time, resulting in a decrease in the dispersion quality. For the experimental test, we utilized a Masterbatch CW2-45, which was produced by Arkema Group Co. To disperse the premix, we employed a high-speed bead mill mixer with 55% carboxymethylcellulose and 45% MWCNT. This Masterbatch can be used directly in the water needed to prepare the cement paste, which helps to avoid the coagulation of carbon nanotubes in the dispersion. This technology has been patented for use in cement systems [15].
Therefore, the aim of this paper was to focus on the modification of the cement paste structure via the addition of CNTs. To achieve this, we aimed to determine the atomic structure of the C-S-H/CNT composite and simulate it.

Molecular Dynamics Simulation
To elucidate the initiation of electrical conductivity and strength improvement in the cementitious composite with CNTs, an understanding of the interaction processes between nanoadditives and cement paste is needed. The electrical and mechanical properties of the cement paste containing CNTs were determined using a molecular dynamics study. The piezoelectric effect can be obtained using applied strains. Under compression or tension, the atoms change position in the crystal structure and generate a voltage. The polarization along the r-axis can be obtained via the following process: where q i is the charge of the ith atom, r i is the r-coordinate of the ith atom, N is the number of atoms, and V is volume.
When an electric field is applied to generate deformations, the piezoelectric effect is simulated. It can be represented by matrices [11] to obtain the strain resulting from the elastic, thermal, electric, and magnetic phenomena. The tensors for variables such as the stress σ kl , electric field E k , magnetic field H l , and temperature differential ∆T can be represented as follows: where S ijkl is the elastic compliance, d kij is piezo-electric, q lij is piezo-magnetic, and α ij is the thermal expansion tensor. Apart from the aforementioned equation, we can represent the effects of interest in this study using another equation that eliminates the magnetic effect. This can be expressed as: where D i is the electric displacement field, d ijk is the piezo-electric tensor, ik is the material electric permittivity, and p i is the pyro-electric tensor. For this study, we selected tobermorite as the main constituent of cement after hydration, as it has a similar structure to that of the C-S-H gel. Tobermorite is a crystalline and natural mineral with the chemical formula Ca 5 Si 6 0 16 (OH) 2 4H 2 O. There are different subspecies of tobermorite categorized by their interlayer spacing, with tobermorite 11Å being the most commonly used model to describe cement paste in MD simulations [16][17][18]. Its initial structure is monoclinic and has vectors a = 6.735Å, b = 7.385Å, and c = 22.487Å, and angles between vectors α = 90 • , β = 90 • , and γ = 123.25 • , as determined by Merlino in 1999 [19].
To create the initial configuration system, a Python code was utilized. The triclinic unit cell basis can be expressed in a matrix format using the following formula: The vector format was used to include the fractional positions and the positions were extended in the three cartesian dimensions using three loops based on the system size to create a supercell, as shown in Figure 1. We employed tobermorite 11Å supercells of size 5 × 5 × 3(Å), containing a total of 9683 to 9899 atoms in this project.
The vector format was used to include the fractional positions and the positions were extended in the three cartesian dimensions using three loops based on the system size to create a supercell, as shown in Figure 1. We employed tobermorite 11Å supercells of size 5 × 5 × 3(Å), containing a total of 9683 to 9899 atoms in this project. We opted for a single-wall carbon nanotube, specifically the armchair (3,3)-type with a diameter of 4.073Å and length of around 70Å, due to its ideal diameter-to-insertion ratio within the tobermorite structure. The nanotubes were inserted into the tobermorite structure along the z-axis (Figure 1) after establishing a vacuum in the structure. In the context of this investigation, the word "vacuum" denotes the empty space generated to accommodate the nanotube insertion into the tobermorite structure while preserving the system's neutral charge. This was accomplished by refraining from allocating unit cells to a specific location in the xy plane, resulting in the absence of atoms within the simulation box. In Figure 2, to create a void in the xy plane and along the z-direction for the 5 × 5 × 3 supercell to insert one CNT, unit cells U331, U332, and U333 were not considered in the code. As a result, the code created unit cells around the nanotube. When the unit cell's geometry, particularly its cross-sectional dimensions, closely matches the targeted hole diameter, the ensuing vacuum within the monoclinic-shaped tobermorite 11Å exhibits a cylindrical appearance after the minimization procedure, as illustrated in Figure 2. Therefore, it could be advantageous to consider the dimensions of the unit cell during the creation of the hole to obtain the desired shape of the void. There are alternative approaches to incorporating CNTs into tobermorite, apart from the method mentioned previously. For example, in a study by Eftekhari M. et al. [20], a cylindrical diameter of a hole, which was initially zero, was incrementally increased using the "fix ident cylinder" command of LAMMPS [21]. Generating a vacuum is also possible by removing the less stable atoms while keeping the molecular neutral following the verification of the charge [22]. Plassard C. et al. [23] noted the absence of silica tetrahedral in tobermorite like C-S-H and accordingly identified the silicon atoms with lower stability within the silicon chains. By eliminating the less stable silicon atoms, the system's charge can remain neutral while altering the geometry of the silica tetrahedron to a silica trigonal pyramidal structure. We opted for a single-wall carbon nanotube, specifically the armchair (3,3)-type with a diameter of 4.073Å and length of around 70Å, due to its ideal diameter-to-insertion ratio within the tobermorite structure. The nanotubes were inserted into the tobermorite structure along the z-axis ( Figure 1) after establishing a vacuum in the structure. In the context of this investigation, the word "vacuum" denotes the empty space generated to accommodate the nanotube insertion into the tobermorite structure while preserving the system's neutral charge. This was accomplished by refraining from allocating unit cells to a specific location in the xy plane, resulting in the absence of atoms within the simulation box. In Figure 2, to create a void in the xy plane and along the z-direction for the 5 × 5 × 3 supercell to insert one CNT, unit cells U331, U332, and U333 were not considered in the code. As a result, the code created unit cells around the nanotube. When the unit cell's geometry, particularly its cross-sectional dimensions, closely matches the targeted hole diameter, the ensuing vacuum within the monoclinic-shaped tobermorite 11Å exhibits a cylindrical appearance after the minimization procedure, as illustrated in Figure 2. Therefore, it could be advantageous to consider the dimensions of the unit cell during the creation of the hole to obtain the desired shape of the void. There are alternative approaches to incorporating CNTs into tobermorite, apart from the method mentioned previously. For example, in a study by Eftekhari M. et al. [20], a cylindrical diameter of a hole, which was initially zero, was incrementally increased using the "fix ident cylinder" command of LAMMPS [21]. Generating a vacuum is also possible by removing the less stable atoms while keeping the molecular neutral following the verification of the charge [22]. Plassard C. et al. [23] noted the absence of silica tetrahedral in tobermorite like C-S-H and accordingly identified the silicon atoms with lower stability within the silicon chains. By eliminating the less stable silicon atoms, the system's charge can remain neutral while altering the geometry of the silica tetrahedron to a silica trigonal pyramidal structure. Materials 2023, 16, x FOR PEER REVIEW 5 of 14 To govern the interaction between the carbon atoms of CNT, we adopted the Tersoff potential [24]. The cumulative energy associated with the interatomic forces between all the atoms within the tobermorite structure is derived from the summation of all the various types of atomic interactions. Such interactions between the atoms of tobermorite and CNT are described in numerous studies, such as in [25,26].
To examine the incidence of the concentration of CNT in the cement paste, systems without СNT and one, two, and three CNTs that yield 0.00%, 1.22%, 2.55%, and 4.00% of CNTs by weight (wt.%) of tobermorite were prepared ( Figure 3). To determine the concentrations, one can calculate the product of the total number of atoms present in the system and their respective atomic masses. The simulation employed periodic boundary conditions to simulate a larger system and prevent the boundary effect. After minimization, an NPT ensemble was utilized with a time step of 1 femtosecond at a temperature of 300 K and zero pressure to regulate the temperature and pressure of the simulation box.

Carbon Nanotubes
Multi-wall carbon nanotubes (MWCNT) from Arkema company were used in this study. The product used for this experiment was Graphistrength ® CW 2-45, which is a masterbatch with a high concentration of 45% by weight of pre-dispersed carbon nanotubes Graphistrength ® C100 (Figure 4) with the typical properties presented in Table  1, perfectly dispersed in carboxyl methyl cellulose (CMC). To govern the interaction between the carbon atoms of CNT, we adopted the Tersoff potential [24]. The cumulative energy associated with the interatomic forces between all the atoms within the tobermorite structure is derived from the summation of all the various types of atomic interactions. Such interactions between the atoms of tobermorite and CNT are described in numerous studies, such as in [25,26].
To examine the incidence of the concentration of CNT in the cement paste, systems without CNT and one, two, and three CNTs that yield 0.00%, 1.22%, 2.55%, and 4.00% of CNTs by weight (wt.%) of tobermorite were prepared ( Figure 3). To determine the concentrations, one can calculate the product of the total number of atoms present in the system and their respective atomic masses. To govern the interaction between the carbon atoms of CNT, we adopted the Tersoff potential [24]. The cumulative energy associated with the interatomic forces between all the atoms within the tobermorite structure is derived from the summation of all the various types of atomic interactions. Such interactions between the atoms of tobermorite and CNT are described in numerous studies, such as in [25,26].
To examine the incidence of the concentration of CNT in the cement paste, systems without СNT and one, two, and three CNTs that yield 0.00%, 1.22%, 2.55%, and 4.00% of CNTs by weight (wt.%) of tobermorite were prepared ( Figure 3). To determine the concentrations, one can calculate the product of the total number of atoms present in the system and their respective atomic masses. The simulation employed periodic boundary conditions to simulate a larger system and prevent the boundary effect. After minimization, an NPT ensemble was utilized with a time step of 1 femtosecond at a temperature of 300 K and zero pressure to regulate the temperature and pressure of the simulation box.

Carbon Nanotubes
Multi-wall carbon nanotubes (MWCNT) from Arkema company were used in this study. The product used for this experiment was Graphistrength ® CW 2-45, which is a masterbatch with a high concentration of 45% by weight of pre-dispersed carbon nanotubes Graphistrength ® C100 (Figure 4) with the typical properties presented in Table  1, perfectly dispersed in carboxyl methyl cellulose (CMC). The simulation employed periodic boundary conditions to simulate a larger system and prevent the boundary effect. After minimization, an NPT ensemble was utilized with a time step of 1 femtosecond at a temperature of 300 K and zero pressure to regulate the temperature and pressure of the simulation box.

Carbon Nanotubes
Multi-wall carbon nanotubes (MWCNT) from Arkema company were used in this study. The product used for this experiment was Graphistrength ® CW 2-45, which is a masterbatch with a high concentration of 45% by weight of pre-dispersed carbon nanotubes Graphistrength ® C100 (Figure 4) with the typical properties presented in Table 1, perfectly dispersed in carboxyl methyl cellulose (CMC).     Graphistrength ® CW 2-45 is presented in the form of pellets with the key characteristics given in Table 2.  Table 3 summarises the composition of clinker and some physical properties of the cement used in this work, provided by the Vicat group, France. Standard sand (ISO standard sand) was certified in accordance with EN 196-1, which is siliceous, and its finest fractions with the grading given in Table 4 were used. Mortar mix with 0.006 wt.% of MWCNTs was chosen as a principle mix based on the results of [28] that show the best concentration of CNT in cement concrete using the same carbon nanotubes from Arkema. We hypothesized that this represents 1.22% of the CNTs by weight of tobermorite. To obtain the proportionality of the nanotube concentration compared to the simulation and obtain additional results, we prepared a mix with 0.018 wt.% of MWCNTs that represent 4.00% of CNTs by weight of tobermorite.
We used 40 × 40 × 160 cm 3 prismatic molds to prepare the samples with two 2D inox meshes 2 × 2 cm placed in the normal axis 2 cm apart from each other to simulate the electrodes of the capacitor (see Figures 5 and 6) and homogenizing via a vibrating table. Table 5 shows the mixing proportions. After curing for 28 days, the samples were cut into two parts to obtain 6 of each concentration of MWCNTs of size 40 × 40 × 80 cm 3 (see Figure 5).

Testing Procedure
The samples were heated by an infra-red lamp to simulate the sun, and the corresponding electrostatic current was measured with a Keithley model 6482 dualchannel picoammeter. The temperatures were measured using a thermal camera FLIR Systems AB model 60b x 1.0 ( Figure 5). As shown in Figure 5, the samples with and without MWCNTs look identical. Only a scanning electron microscope (SEM) analysis can be used to differentiate them.
The Keithley model 6482 dual-channel picoammeter was connected to the two electrodes of the capacitor ( Figure 6). The leakage current was measured when the sample was subjected to a temperature field.

Dynamic Molecular Study
The Young's modulus of the built model was evaluated and compared with previous research using Tables 6 and 7. The system was equilibrated before undergoing deformation via a uniaxial tension by fixing one side of the structure and stretching the other side at a strain rate of 0.00001/fs. The initial length of the simulation box in the zdirection, as shown in Figure 3, was saved before this process.

Testing Procedure
The samples were heated by an infra-red lamp to simulate the sun, and the corresponding electrostatic current was measured with a Keithley model 6482 dualchannel picoammeter. The temperatures were measured using a thermal camera FLIR Systems AB model 60b x 1.0 ( Figure 5). As shown in Figure 5, the samples with and without MWCNTs look identical. Only a scanning electron microscope (SEM) analysis can be used to differentiate them.
The Keithley model 6482 dual-channel picoammeter was connected to the two electrodes of the capacitor ( Figure 6). The leakage current was measured when the sample was subjected to a temperature field.

Dynamic Molecular Study
The Young's modulus of the built model was evaluated and compared with previous research using Tables 6 and 7. The system was equilibrated before undergoing deformation via a uniaxial tension by fixing one side of the structure and stretching the other side at a strain rate of 0.00001/fs. The initial length of the simulation box in the zdirection, as shown in Figure 3, was saved before this process.

Testing Procedure
The samples were heated by an infra-red lamp to simulate the sun, and the corresponding electrostatic current was measured with a Keithley model 6482 dual-channel picoammeter. The temperatures were measured using a thermal camera FLIR Systems AB model 60b × 1.0 ( Figure 5). As shown in Figure 5, the samples with and without MWCNTs look identical. Only a scanning electron microscope (SEM) analysis can be used to differentiate them.
The Keithley model 6482 dual-channel picoammeter was connected to the two electrodes of the capacitor (Figure 6). The leakage current was measured when the sample was subjected to a temperature field.

Dynamic Molecular Study
The Young's modulus of the built model was evaluated and compared with previous research using Tables 6 and 7. The system was equilibrated before undergoing deformation via a uniaxial tension by fixing one side of the structure and stretching the other side at a strain rate of 0.00001/fs. The initial length of the simulation box in the z-direction, as shown in Figure 3, was saved before this process. After carrying out separate tensile tests on tobermorite 11Å and CNT (3,3), the combination of tobermorite 11Å + CNT was subjected to the same test to determine the effect of CNT concentration on the tensile resistance of tobermorite 11Å. Table 8 and Figure 7 present the results for Young's modulus, tensile strength, and stress evolution due to applied strain for the tobermorite-CNT models. The initial stress peak at 0.1 strain (Figure 7) indicates bond breaking in the tobermorite atoms (Figure 8b), followed by an increase in stress as only the CNTs were strained (Figure 8c) until the CNTs ruptured at around 0.25 strain (Figure 8d).   [29] 103.25 Molecular Dynamics Lushnikova [30] 78.39 Molecular Dynamics R. Shahsavari et al. [31] 82.82 Ab initio  [33] 0.510 Molecular Dynamics X. Lu, Z. Hu [34] 0.989-1.058 Finite Element Analysis P Subba Rao [35] 1.

22-1.28 Finite Element Analysis
After carrying out separate tensile tests on tobermorite 11Å and CNT (3,3), the combination of tobermorite 11Å + CNT was subjected to the same test to determine the effect of CNT concentration on the tensile resistance of tobermorite 11Å. Table 8 and Figure 7 present the results for Young's modulus, tensile strength, and stress evolution due to applied strain for the tobermorite-CNT models. The initial stress peak at 0.1 strain (Figure 7) indicates bond breaking in the tobermorite atoms (Figure 8b), followed by an increase in stress as only the CNTs were strained (Figure 8c) until the CNTs ruptured at around 0.25 strain (Figure 8d).    Subsequently, the piezoelectric effect was simulated. The system was first equilibrated, and then the lengths of the simulation box were monitored while an external electric field of 250 MV/m was applied in the z-direction (as shown in Figure 3). The evolution of the box length in the z-direction for the four systems is compared and presented in Figure 9. The strain tensor was determined between simulation steps 4000 and 10,000, resulting from an electric field of 250MV/m applied in the z-direction (refer to Figure 3). As there was no external stress or temperature change, Equation (2) for the converse piezoelectric effect in a triclinic structure could be expressed as follows: = This can be solved for , yielding: The piezoelectric constant calculated for each of the samples is shown in Table 9. Subsequently, the piezoelectric effect was simulated. The system was first equilibrated, and then the lengths of the simulation box were monitored while an external electric field of 250 MV/m was applied in the z-direction (as shown in Figure 3). The evolution of the box length in the z-direction for the four systems is compared and presented in Figure 9. Subsequently, the piezoelectric effect was simulated. The system was first equilibrated, and then the lengths of the simulation box were monitored while an external electric field of 250 MV/m was applied in the z-direction (as shown in Figure 3). The evolution of the box length in the z-direction for the four systems is compared and presented in Figure 9. The strain tensor was determined between simulation steps 4000 and 10,000, resulting from an electric field of 250MV/m applied in the z-direction (refer to Figure 3). As there was no external stress or temperature change, Equation (2) for the converse piezoelectric effect in a triclinic structure could be expressed as follows: This can be solved for , yielding: The piezoelectric constant calculated for each of the samples is shown in Table 9. The strain tensor ε ij was determined between simulation steps 4000 and 10,000, resulting from an electric field of 250 MV/m applied in the z-direction (refer to Figure 3). As there was no external stress or temperature change, Equation (2) for the converse piezoelectric effect in a triclinic structure could be expressed as follows: This can be solved for d 33 , yielding: The piezoelectric constant calculated for each of the samples is shown in Table 9. The tobermorite 11Å structure alone displays a higher ultimate tensile strength than the models Tob-CNT (1.22 wt.%) and Tob-CNT (2.55 wt.%). However, the Tob-CNT model (4.00 wt.%) has the highest tensile strength of 14.27 GPa and the highest Young's modulus of 165.46 GPa, as shown in Table 8. It is observed that the models Tob-CNT (1.22 wt.%) and Tob-CNT (2.55 wt.%) exhibit a decreased maximum tensile strength. In all cases where CNTs were incorporated into tobermorite, a higher Young's modulus was observed. The initial stress peak, which is associated with the collapse of the tobermorite structure, occurs at a strain of 0.1. Figure 9 shows multiple stress peaks, which indicate that collapse occurs for strains ranging from 0.25 to 0.3 when CNTs are present. Despite a discontinuous collapse in the modeling when translating reality, cracking should not be separated. This numerical result can be interpreted as a discontinuity in the overall strain. Under the application of an electric field, the samples follow displacements, and the sample oriented in the same direction (z) of the electric field is analyzed. The piezoelectric constant is aligned in this direction, and the highest strain of 0.002155 was observed in Tob-CNT (2.55 wt.%). The samples with 1.22 wt.% and 4.00 wt.% displayed similar results. Without CNTs, the response to the applied electrical field is weak. By taking into account these numerical considerations, it should be experimentally possible to cause the deformation of a mortar sample containing CNT by a strain field (temperature field, for example), and thus generate an electric field. That is what we have tested in the following section.

Experimental Study
A total of 137 measurements of current were carried out on 6 samples of size 40 × 40 × 80 cm 3 (see Figure 5) with 0.00 wt.% of MWCNTs. For 0.006 wt.% and 0.018 wt.% of MWCNTs, 6 samples were prepared and 109 measurements of current were carried out on each mixture. A solar lamp was used to simulate these materials. It can produce UVA or UVA-UVB. This infrared lamp was chosen to simulate sun exposure. The temperatures were measured using a thermal camera ( Figure 5). The values of the current obtained from the experimental test in η A with respect to the temperature for all the samples are displayed in Figure 10. Only the measurements where the temperature was stabilized at 20 • C (ambient temperature), 30 • C, 35 • C, 45 • C, and 50 • C were retained in the range of ±2 • C. The results show great dispersion (see Figure 10) because electrical measurements at such low levels are very sensitive to the environment. An electrically neutral enclosure, controlled by temperature and humidity, should have been used. This is planned as the next step of this research. Nevertheless, the tendencies (see Figure 10) show that the current increases for the same temperature when the percentage of wt.% of MWCNTs increases in the mixture, which was expected, given the dynamic molecular study.
As shown in Equation (2), the strain (ε) of a material can be expressed as the change in temperature (∆T) multiplied by the coefficient of thermal expansion (α). For mortar, the coefficient of thermal expansion is in the order of 10 × 10 −6 / • C [36]. Table 10 uses α∆T to express the increment of the measured current concerning the strain by averaging the samples that have the same CNT concentration. The average electrostatic current increment for a strain of 10 × 10 −5 and ∆T= 10 • C increases between 1.28 ( η A) and 1.51 ( η A) when the percentage increases between 0.000 wt.% and 0.018 wt.%.  Dynamic molecular and experimental studies demonstrated that the piezoelectric effect can be obtained in a mortar containing MWCNTs. As in certain ceramics, piezoelectricity is the electric charge that accumulates in solid materials in response to applied stresses and reciprocally. This electricity, resulting from the displacement caused by an increase in temperature, can be capped. The piezoelectric phenomenon arises due to the interdependent relationship between the mechanical and electrical properties of materials. By including MWCNTs in a mortar, we can achieve a piezoelectric effect. This piezoelectric effect can be used as a sensor to analyze the field of stress detection caused by loads or temperature.

Conclusions
This study examined how the inclusion of carbon nanotubes affects the electrical properties of cement composites. To investigate this, a molecular dynamics simulation was performed using tobermorite 11Å as a representative of calcium silicate hydrate gel, combined with armchair (3,3)-type single-wall carbon nanotubes. Thus, numerical simulations of tobermorite mixtures with different percentages of carbon nanotubes were prepared. First, the molecular dynamics study showed that the piezoelectric effect can be simulated for four different concentrations. The deformation in the simulation box was then analyzed by subjecting it to an electrical field of 250 MV/m in the z-direction, The results showed a stronger converse piezoelectric response when the percentage of carbon nanotubes was increased. The deducted piezoelectric constant in the z-direction d33 increased in the same manner, proving in theory that current can be generated when samples are subjected to strains.
The experimental investigation was conducted using samples of mortar reinforced with different amounts (0.00 wt.%, 0.006 wt.%, 0.018 wt.%) of tangled multi-walled carbon nanotubes from Arkema. The samples were heated by an infra-red lamp to simulate the  Dynamic molecular and experimental studies demonstrated that the piezoelectric effect can be obtained in a mortar containing MWCNTs. As in certain ceramics, piezoelectricity is the electric charge that accumulates in solid materials in response to applied stresses and reciprocally. This electricity, resulting from the displacement caused by an increase in temperature, can be capped. The piezoelectric phenomenon arises due to the interdependent relationship between the mechanical and electrical properties of materials. By including MWCNTs in a mortar, we can achieve a piezoelectric effect. This piezoelectric effect can be used as a sensor to analyze the field of stress detection caused by loads or temperature.

Conclusions
This study examined how the inclusion of carbon nanotubes affects the electrical properties of cement composites. To investigate this, a molecular dynamics simulation was performed using tobermorite 11Å as a representative of calcium silicate hydrate gel, combined with armchair (3,3)-type single-wall carbon nanotubes. Thus, numerical simulations of tobermorite mixtures with different percentages of carbon nanotubes were prepared. First, the molecular dynamics study showed that the piezoelectric effect can be simulated for four different concentrations. The deformation in the simulation box was then analyzed by subjecting it to an electrical field of 250 MV/m in the z-direction, The results showed a stronger converse piezoelectric response when the percentage of carbon nanotubes was increased. The deducted piezoelectric constant in the z-direction d 33 increased in the same manner, proving in theory that current can be generated when samples are subjected to strains.
The experimental investigation was conducted using samples of mortar reinforced with different amounts (0.00 wt.%, 0.006 wt.%, 0.018 wt.%) of tangled multi-walled carbon nanotubes from Arkema. The samples were heated by an infra-red lamp to simulate the sun, and the corresponding electrostatic current was measured with a Keithley model 6482 dualchannel picoammeter. In parallel, the temperatures were measured using a thermal camera with a precision of ±2 • C. The results showed that electrostatic current measurements exhibited a higher electric sensitivity for samples with a higher concentration of CNTs. These experimental results are very promising when compared with molecular dynamics studies, and such effects can be used for analyses of stress detection caused by loads or temperatures in concrete structural elements.