Use of fine aggregate matrix to analyze the rheological behavior of cold recycled materials

Nowadays, one of the main challenges to a wider application of cold recycling techniques is the lack of reliable information on the mechanical behavior of cold recycled materials (CRM). In this context, measurement and modelling of the complex modulus of CRM mixtures may give an important contribution to the design and analysis of pavements including cold recycled layers. In this study, we analyzed the rheological behavior of CRM mixtures produced using bitumen emulsion and cement through the study of their fine aggregate matrix (FAM). Starting from a fixed CRM mixture composition, we compared different FAM mortars, focusing on the effect of water and air content. Then, we selected a composition as representative of the FAM in the mixture and investigated the evolution of both materials during a fixed curing period. Next, we measured the complex modulus of the CRM mixture and FAM at two curing stages and applied a rheological model to simulate and compare their behavior. Results showed that the properties of CRM mixtures are comparable to those of FAM mortars produced using all the binding agents (bitumen emulsion and cement) and a fraction of the voids contained in the mixture. Despite the huge difference in volumetric compositions, the FAM mortar controlled the curing and the thermo-rheological behavior of the CRM mixture, while the coarse reclaimed asphalt aggregate fraction and the voids mainly affected the asymptotic properties (equilibrium and glassy moduli) and the non-viscous dissipation component.

Abstract Nowadays, one of the main challenges to a wider application of cold recycling techniques is the lack of reliable information on the mechanical behavior of cold recycled materials (CRM). In this context, measurement and modelling of the complex modulus of CRM mixtures may give an important contribution to the design and analysis of pavements including cold recycled layers. In this study, we analyzed the rheological behavior of CRM mixtures produced using bitumen emulsion and cement through the study of their fine aggregate matrix (FAM). Starting from a fixed CRM mixture composition, we compared different FAM mortars, focusing on the effect of water and air content. Then, we selected a composition as representative of the FAM in the mixture and investigated the evolution of both materials during a fixed curing period. Next, we measured the complex modulus of the CRM mixture and FAM at two curing stages and applied a rheological model to simulate and compare their behavior. Results showed that the properties of CRM mixtures are comparable to those of FAM mortars produced using all the binding agents (bitumen emulsion and cement) and a fraction of the voids contained in the mixture. Despite the huge difference in volumetric compositions, the FAM mortar controlled the curing and the thermo-rheological behavior of the CRM mixture, while the coarse reclaimed asphalt aggregate fraction and the voids mainly affected the asymptotic properties (equilibrium and glassy moduli) and the non-viscous dissipation component.

Introduction
Cold recycling of bituminous pavements has gained the interest of an increasing number of road agencies and contractors worldwide. On the one hand, this is due to the increasing emphasis that is placed on sustainability and reuse of end-of-life products like reclaimed asphalt (RA) [1][2][3]. On the other hand, field applications are showing that cold recycled materials (CRM) are an economically and structurally viable solution even for high traffic roadways [4][5][6][7].
In order to promote a wider application of CRM mixtures, a deeper knowledge of their mechanical behavior is critical. First, it must be recognized that CRM were developed and are applied following different composition concepts, leading to diverse mechanical behaviors [8,9]. CRM without any binder or treated only with Portland cement do exists but are not in the scope of this paper. In most cases, bitumen (bitumen emulsion or foamed bitumen) is the main binding agent and it is used in combination with Portland cement or chemical additives, such as lime. If Portland cement is used with dosages greater than about 1% with respect to the mass of aggregates, it has clear effect on the long-term properties of CRM mixtures, hence bitumen and cement shall be considered as co-binders [8,10,11]. In this case, the rheological behavior of CRM mixtures can be characterized by measuring the complex modulus [12,13]. However, compared to HMA, the complex modulus characterization of CRM mixtures must consider additional composition-related issues.
The first is that, due to water evaporation, emulsion breaking and cement hydration, CRM mixtures are subjected to a curing process and evolve from the ''fresh'' to the ''hardened'' state [14]. Therefore, if we are looking for stiffness values to be used in pavement design calculations, measurements must be carried out after a certain curing period. The second issue is that, CRM mixtures contain bitumen from emulsion or foam, as well as RA bitumen. The former is ''fresh'' and provides cohesion to the mixture, thanks to the emulsification of foaming processes. The former is aged, and, without heating, it is not able to blend with the fresh bitumen or to produce cohesion [15]. However, both the fresh and the aged bitumen are time-and temperature dependent materials and therefore they may affect the thermo-rheological behavior of the mixture [12,16,17].
The complex modulus of CRM mixtures can be measured using the same experimental procedures that were developed for HMA mixtures, and CRM mixtures can be considered thermo-rheologically simple materials, i.e. the time-temperature superposition principle can be applied to both the stiffness modulus (absolute value of the complex modulus) and the phase angle [8,13]. In addition, the rheological models developed to simulate complex modulus of HMA, like the Huet-Sayegh model [18] or the 2S2P1D model [19] can also be used for CRM mixtures [12,20,21].
A further improvement to the rheological characterization of CRM mixtures is expected from the application of a multiscale approach and specifically from the study of their fine aggregate matrix (FAM). FAM materials were introduced to study the linear viscoelastic behavior, the fatigue and moisture damage properties of HMA mixtures [22][23][24]. FAM materials are essentially mortars, i.e. mixtures characterized by a nominal maximum aggregate size ranging between 1.18 and 4.76 mm, whose composition should represent the FAM as it exists in the mixture [25][26][27]. FAM mortars have the same gradation of the fine fraction of the original mixture and are normally fabricated using all its bituminous binder, except the small fraction which is part of the thin film adhering to the surface of the coarse aggregate. The air voids content of FAM mortars is generally comprised between 0.5 and 0.7 of the total air voids of the mixture [27].
In the field of cold bituminous materials, which clearly includes CRM mixtures, mortars were considered either as model systems for investigating the properties of bitumen emulsions and their interaction with the cementitious binders [28,29] or, alternatively, to investigate the properties of the CRM mixture itself [30].
The present research follows the second approach and our main objective is to compare the rheological behavior of CRM mixtures to that of mortars representing their FAM. The comparison is based on the measurement and modelling of the complex modulus. In the first part of the experimental study, starting from a fixed CRM mixture composition, we tested different mortar compositions obtained by removing the coarse aggregate and by changing the water and air content. In the second part, we selected a mortar composition to represent the FAM of the CRM mixture and investigated its curing behavior by monitoring the evolution of stiffness and water loss by evaporation. In the third part, we measured the complex modulus of both the CRM mixture and the selected FAM at two curing stages and finally we applied a new rheological model to analyze and compare their behavior.

Materials
The CRM mixtures and FAM mortars tested in this study were produced in the laboratory by mixing RA aggregate, filler, Portland cement, bitumen emulsion and water, at ambient temperature.
The RA aggregate was sampled from a stockpile located in a cold recycling plant in Italy. Its nominal maximum aggregate size was 16 mm (Fig. 1), the particle density was 2482 kg/m 3 , the water absorption was 1.0% and the bitumen content was 4.9% (by dry aggregate mass). In the laboratory, coarse and fine RA aggregate fractions were separated using the 2 mm sieve, their bitumen content was 3.3% and 8.3%, respectively. The FAM mortar was produced using the fine fraction of RA aggregate, passing to the 2 mm sieve (Fig. 1). In this study RA aggregate is considered a ''black rock'', meaning that the aged asphalt is considered a part of the RA particles. The filler was a limestone powder with a particle density of 2650 kg/ m 3 . The Portland cement was a GU type, with compressive strength at 28 days of 43.9 MPa (ASTM C109). The bitumen emulsion was a commercial slowsetting type, designated as C 60 B 10 (EN 13108) specifically produced for cold recycling applications with cement (overstabilized emulsion), with a residual bitumen content of 60%.

Mixing and compaction procedures
For mixing CRM mixtures and mortars, a strict control of the water content is necessary. For this reason, the RA aggregate was dried at 40 ± 2°C until reaching constant mass and then a water amount corresponding to the absorption was added. The wet samples were stored in sealed plastic bags overnight, at room temperature, to allow a uniform moisture distribution and absorption. Mixing was carried out with a mechanical mixer and, during the process, additional water, cement, and bitumen emulsion were gradually added to the aggregate blend [31].
All the specimens were compacted with a gyratory compactor, adopting a constant pressure of 600 kPa, a gyration speed of 30 rpm and an angle of inclination of 1.25°. Molds with a diameter of 100 mm and 150 mm were used for the compaction of mortars and mixtures, respectively.

Composition of the CRM mixture
The dry aggregate blend of the mixture was composed of 94% of RA aggregate and 6% of filler (by mass) and had a grading distribution close to the maximum density curve (Fig. 1). The dosage of bitumen emulsion (residual binder) and Portland cement were 5.0% (3.0%) and 1.5%, respectively (by mass of the dry aggregate blend). Additional water was added to the mixture to help the homogeneous distribution of bitumen droplets and to improve workability and compactability. The total water content of the mixture, Fig. 1 Grading distribution of RA aggregate, CRM mixture and FAM mortar including additional water and emulsion water was 4.0% (by mass of the dry aggregate blend). This can be further subdivided into water absorbed by the RA aggregate (1.0%) and intergranular water (3.0%).
The volumetric composition of the freshly compacted mixture specimens is outlined in Fig. 2 and summarized in Table 1. To simplify the volumetric analysis, we consider the RA aggregate in saturated surface dry (SSD) condition, i.e. its volume is the bulk particle volume (including surface porosity), and its mass includes the absorbed water. We calculate the voids in the mixture V m (volumetric fraction of the non-structural part) as follows: where v is the total volume of the specimen (calculated using the mold diameter and the height measured by the gyratory compactor), v W is the volume of the intergranular water (the volume of the absorbed water is included in the bulk volume of aggregates), and v A is the volume of the air voids. The volume of voids v W þ v A is obtained by subtracting the volume of the solids v S (which includes the bulk volume of aggregates, filler and un-hydrated cement) and the volume of the residual bitumen from emulsion v B , from the total volume of the specimen. Equation 1 is correct only if during the compaction the material loss from the mold is negligible [31,32]. For this reason, the mass of each specimen was measured immediately after the extraction from the mold and compared with the initial mass of the loose mixture cast in the mold. For all specimens, the mass loss during compaction was less than 0.3% and there was no observable leakage of water or fines. In this study, we fixed the voids content of the mixture specimens at V m = 15%. According to the mixture gravimetric composition given above, this comprised voids filled with intergranular water and air voids (V A ¼ v A =v = 9.3%). As an example, Fig. 3 shows the compaction diagrams for the mixture and FAM mortar specimens that were tested in the second part of the experimental study (Sect. 3.2). The slope (k) of the compaction curve (voids Vs Log N) is reported in the plot and can be used as a numerical indicator of the compactability [33]. For the mixture specimens, the selected voids level (V m = 15.0%) was obtained after about 30 gyrations and the compactability (k = 8.66) was in agreement with previous studies [34,35]. For the FAM mortar specimens, the selected voids level (V m = 14.8%, Sect. 2.4) was obtained after about 20 gyrations, and the compactability was lower (k = 6.59) with respect to mixtures.

Composition of the FAM mortar
To obtain the volumetric composition of the FAM mortar, starting from the given mixture composition, first we removed the volume of the coarse aggregate (in SSD condition). We included in the FAM a fraction of the mixture water (at least the emulsion water) and a fraction of the mixture air voids [30]. The remaining fraction of water and air was considered as an independent voids phase. As a result, we considered the mixture as a three-phase composite ( We can consider the FAM mortar as the binding matrix of the mixture, because it contains all the binders (bitumen and Portland cement) and the coarse aggregate as particle inclusions.
As regards the voids, separating the fraction that could be considered part of the FAM from the fraction that could be considered as an independent phase in the mixture required some judgement. For example, we could assume that those voids with characteristic size smaller than about 2 mm could be part of the FAM, whereas larger ones could be considered as an independent phase. However, this posed the additional issue of measuring the size distribution of the voids, which was beyond the scope of the research. We decided to investigate this issue by preparing and testing seven mortars with different voids content V m;FAM , i.e. with different V W;FAM and V A;FAM (Fig. 2). The mortar compositions were obtained applying the principles described above and thus were characterized by the same grading distribution (Fig. 1), gravimetric dosage of emulsion (13.5%) and Portland cement (4.1%).
The volumetric composition of mixture and mortars is summarized in Table 1. The reference (REF) mortar was obtained by removing only the SSD coarse RA aggregate from the mixture and contained all the mixture's voids (water and air). Therefore, it would correspond to a mixture without an independent voids phase (V V = 0). This mortar had a volume of voids 11.2% ? 18.2% = 29.4%, while its calculated gravimetric water content was 9.0%. We mixed this mortar composition in the laboratory, but we were not able to compact a specimen because the uncompacted mortar cast in the gyratory mold had less than 24.9% of voids. In other words, the volume of voids included in the mixture (V m = 15.0%) was too high to fit inside the FAM phase and therefore the multiphase mixture model (Fig. 2) must include the independent voids phase.  The mortar A was prepared with the same gravimetric water dosage of the REF mortar but was compacted until reaching the minimum fraction of air voids of 1.1%. We tried to compact this mortar until reaching zero air voids (saturation condition), but it was impossible because leakage of water started from the bottom and the top of the mold and the compaction had to be stopped to preserve the original composition, as explained in Sect. 2.3 for the mixtures. In other words, a minimum amount of entrapped air voids must be included in the mortar.
For compositions B and C, we reduced the gravimetric water dosage to 8.2%, while for compositions D, E and F the water dosage was reduced to 7.5%. A further reduction was not possible because the mortar became too dry, not homogeneous and almost impossible to compact. Like mortar A, mortars B and D were compacted until reaching the minimum level of entrapped air voids (in this case 1.6%), without leakage of water from the mold. Figure 4 describes the volumetric properties of all the tested mortars: the x-axis is the fraction of air voids in the mortar to the air voids in the mixture (v A;FAM =v A ), while the y-axis is the volumetric fraction of water in the mortar to total water in the mixture (v w;FAM =v W ). The dashed lines drawn in the plot describe the values V m;FAM = constant. The upperright corner of the plot represents the REF composition (v A;FAM =v A ¼ v w;FAM =v W ¼ 1). Since the mortars were produced using all the bitumen emulsion of the mixture, the minimum value on the y-axis corresponds to the water content of the emulsion (for our mixture this was 0.66). Moreover, in Fig. 4 we outlined two areas along the bottom and left margins of the plot, which qualitatively indicate configurations that were impossible to mix or compact in the laboratory, starting from the mixture composition given above.

Testing procedures
In the first part of the experimental study, to evaluate the effect of the voids on the mortar composition, we measured the Indirect Tensile Strength (ITS) of the materials listed in Table 1. Two replicate specimens were tested for each material composition. The test was carried out at 25°C, after 1 day of curing at the same temperature [30]. Following the standard EN 12697-23, we applied a constant rate of deformation (50 ± 2 mm/min) until failure and the ITS was calculated following Eq. 2.
where F MAX is the maximum compressive force sustained by the specimen, H and D are the specimen height and diameter, respectively. In the second part of the study, we selected a specific FAM mortar composition (E) and compared its curing behavior to the curing behavior of the mixture. We produced three replicate specimens for both the mixture and the FAM, having 100 mm diameter and 67.4 mm and 64.7 mm height, respectively. Curing started immediately after compaction, and consisted of a first phase of 14 days at 25°C followed by a second phase of 14 days at 40°C [36]. During the curing phases, we measured the water loss (WL) due to evaporation (Eq. 3) after 6 h, 1, 3, 7, 14, 15, 17, 21 and 28 days.
where M 0 is the specimen mass right after compaction, M i is the specimen mass after i curing days and W TOT is the total mass of water in the specimen.
On the same specimens, starting from the third day of curing, we also measured the Indirect Tensile Stiffness Modulus (ITSM) (Eq. 4) according to the standard EN 12697-26. The test measured the average Fig. 4 Volumetric properties of the tested mortars stiffness at 25°C after the application of 10 ? 5 pulse loads (conditioning ? test pulses) with a rise time of 124 ± 4 ms and a target horizontal strain of 3 9 10 -6 mm/mm.
where F max is the maximum value of the pulse load, d max is the corresponding maximum value of the horizontal deformation, m is the Poisson's ratio (assumed equal to 0.35) and H is the specimen height.
In the third part of the study, we measured the complex modulus E Ã x ð Þ (Eq. 5) of the mixture and the selected FAM mortar composition at the end of the two curing stages, i.e. after 14 days at 25°C and after the additional 14 days at 40°C. Two specimens of each material were prepared with the gyratory compactor, in a 100 mm mold, to the final height of 130 mm and then cored to the diameter 75 mm at the end of curing. The test was carried out using an AMPT PRO system where the axial stress was measured with a load cell, while the axial strain was measured in the middle part of the specimen (measuring base of 70 mm) using three LVDT, placed 120°apart. A haversine compression loading was applied to obtain a target strain amplitude of 30Á10 -6 mm/mm. The testing temperatures were 0, 10, 20, 30 and 40°C, the testing frequencies were 0.1, 0.5, 1, 5 and 10 Hz and 20 loading cycles were applied at each frequency.
where j is the imaginary unit, x = 2pf is the angular frequency (f is the testing frequency in Hz), d x ð Þ is the phase angle and E 0 x ð Þ is the stiffness modulus (Eq. 6): where r 0 and 0 are the steady-state sinusoidal amplitudes of the measured stress and strain signals.

Results and analysis
3.1 Mortar composition Figure 5 shows the ratio of the ITS of the six FAM mortar compositions (ITS FAM ) to the ITS of the mixture (ITS MIX = 0.142 MPa). This ratio is plotted as a function of the air voids fraction in the mortar, v A;FAM =v A , measured after compaction (the same used in Fig. 4). The ITS of the FAM mortars increased when the water dosage decreased and, for each water dosage, it increased when the air content was reduced. Moreover, all the tested FAM mortar compositions had an ITS higher than the mixture. This is because only a small fraction of the air content of the mixture could be included in the mortars (less than 32%), while the remaining fraction should be considered as part of the mixture, within the independent voids fraction (Fig. 2, Phase 3). All the mortar volumetric compositions (A to F) are potentially representative of the FAM properties within the mixture. However, to analyze the effect of curing and to evaluate the rheological behavior, we chose composition E. We considered compositions A, B and D not representative of the FAM within the mixture because they were characterized by a very low air voids content (v A;FAM =v A \ 0.1), whereas the mixture had a relatively high air voids content. We preferred composition E over compositions C and F because it required a higher compaction energy (number of gyrations given in Table 2), closer to the energy required by the mixture. Adopting mortar E to represent the FAM, the volumetric composition of the mixture (Fig. 2) would be as follows: • Phase 1 coarse aggregate, V CA = 49.0%.  Table 1). • Phase 3 independent voids: 8.8% (1.1% water, 7.7% air). Figure 6 shows the relationship between the evolution of WL of the FAM (composition E) and the mixture, from 6 h to 28 days of curing. We observe that WL MIX was always higher than WL FAM and, except the first few hours, the rate of water loss for mixture and FAM was practically identical. The same result was obtained in a previous study on similar materials (see Figure 6 in [30]).

Evolutive behavior of mixture and FAM
Water evaporation from porous media is often described as a two-stage process. In the first stage, the rate of water loss per unit surface is nearly constant and close to the evaporation rate from a free water surface. In the second stage, evaporation is characterized by a falling rate and is controlled by the diffusion of vapor through the pores of the medium [37,38]. The results of this study, and the results presented in [30], clearly indicate that after the first few hours of curing, the evaporation rate decreased for all specimens, thus we were in the diffusion stage. Moreover, since the water loss from the mixture and the FAM mortar was strictly related, the pore structure of the two materials could be similar. However, this hypothesis needs further confirmation that can be obtained only with an experimental characterization of the pore structure of both materials. Figure 7 shows the relationship between the ITSM of mixture and FAM, from 3 to 28 days of curing. At short curing times, the mixture was stiffer than the FAM mortar. This indicated that the reinforcing effect of the CA particles (V CA ) prevailed on the weakening effect of the independent voids (V V ). Afterwards, the difference with FAM reduced and, after about 18 days, the FAM became stiffer than the mixture, highlighting the effect of voids.
To model the relation between the ITSM of mixture and FAM, among the various equations proposed in the literature for particulate composites [39], we chose a slightly modified version of the Hirsch model (Eq. 7), that was originally developed for cement concrete and later applied to HMA [40].  According to the volumetric composition of mortar E, we used V CA ¼ 0:490, V FAM ¼ 0:422, and thus we explicitly considered the independent voids phase in parallel arrangement with the structural part of the material (CA and FAM), as outlined in Fig. 7. Using least squares error minimization, we estimated E CA ¼ 2912 MPa and x ¼ 0:108. The resulting model is plotted in Fig. 7 and gives an excellent fitting of the experimental data. We recall that parallel arrangement implies equality of strain and thus perfect bonding between the FAM and the CA particles, on the other hand, series arrangement implies equality of stress which can be related to poor bonding [39]. The model fitting results suggests that series arrangement prevailed (89.2%), suggesting that in the mixture the bonding between CA particles (RA particles retained on the 2 mm sieve) and FAM was rather poor. This behavior clearly resembles the behavior of HMA, but for CRM mixtures and FAM it is due to the frequency and temperature susceptibility of the two distinct bituminous components: the residual bitumen from the emulsion and the RA bitumen. Moreover, cementitious bonds limited the variability of the stiffness modulus values within about one order of magnitude (700-7000 MPa) and the maximum phase angle to values less 20°. Similar trends can be observed in [12,13,41].

Complex modulus results
In Fig. 8 we also outlined the effect of curing, which can be better observed in Fig. 9, where we directly compared the stiffness modulus and phase angle results (same temperature and same frequency) after 14 and 28 days. Throughout the temperature and frequency range, E 0 increased between 15 and 20%, in agreement with the ITSM measurements shown in Fig. 7. Concurrently, there was a phase angle decrease of about 1°to 4°. Overall, the E Ã variation followed about the same pattern as for temperature decrease and frequency increase [14] which explains why the effect of curing is not clearly visible in Fig. 8. 3.4 Thermo-rheological modeling Figure 8 shows that, for all the specimens and curing conditions, the (E 0 ; d) pairs describe a single curve in the Black diagram. This confirms the validity of the time-temperature superposition principle (variations of temperature and frequency have the same effect on E Ã ) and thus allows the creation of master curves. The master curves of the stiffness modulus (Fig. 10a) were obtained by selecting the reference temperature T ref ¼ 20°C and shifting the data  Fig. 9 Effect of curing on the rheological properties: a stiffness modulus; b phase angle measured at the other temperatures parallel to the logfrequency axis using Eq. 8: where f is the testing frequency, f r is the shifted frequency, normally called reduced frequency and a T ref T ð Þ is the shift factor required to obtain the superposition between the data measured at different temperatures T and at T ref . The shift factors were calculated using the closed form shifting (CFS) algorithm which is based on the minimization of the area between two adjacent isothermal segments and is operator-independent [42].
The master curves of the phase angle (Fig. 10b) were created applying the same shift factors calculated for the stiffness modulus. As can be observed, the superposition is excellent and confirms that the smallstrain response of the tested mixture and FAM can be considered thermo-rheologically simple. Figure 11 shows the calculated shift factors as a function of temperature (values of replicate specimens were averaged). The experimental data are modelled with the Williams-Landel-Ferry (WLF) equation [43]: where C 1 , C 2 are constants obtained by least squares fitting ( Table 2). The complex modulus data were fitted with a modified version of the Huet-Sayegh (HS) rheological model, described by Eq. 10 [12]: where E Ã HS x ð Þ represents the HS model and exp jd 0 ð Þis The HS model is described by Eq. 11: where x ¼ 2pf r is the angular frequency E e and E g are the equilibrium and glassy values of the stiffness modulus, d, h and k are adimensional parameters controlling the shape of the model and s is a characteristic time whose value depends on the selected reference temperature. Table 2 summarizes the estimated model parameters and Fig. 10 shows the fitted model superposed to the experimental data. In Fig. 12 the new model and the HS model are graphically compared using the Cole-Cole and Black diagrams. CRM mixtures show a frequency-and temperaturedependent behavior that can be attributed to the bituminous component, comprising both the emulsion residue and the aged RAP bitumen. This behavior is analogous to that of HMA and can be simulated using Eq. 11. Specifically, the parameters 0\k\h\1 represent the order of derivation of two fractional derivative elements (FDE), also called parabolic dashpots. From a physical point of view, k translates the linear viscoelastic (LVE) behavior in the lowtemperature (high-frequency) range, whereas h translates the LVE behavior in the high-temperature (lowfrequency) range [18]. Their lower limit is 0, representing purely elastic behavior (the FDE becomes a Hookean spring), whereas their upper limit is 1, representing purely viscous behavior (the FDE becomes a Newtonian dashpot).
Along with viscous damping, Eq. 10 also considers an additional frequency-and temperature-independent damping, represented by the constant phase angle d 0 , which was added to the HS model (Eq. 10 and Fig. 12). This is the well-known hysteretic dissipation mechanism [44], that is used to simulate the dynamic response of many materials, including granular (noncohesive) soils [45] and cement concrete [46]. For CRM mixtures, a small but still measurable hysteretic dissipation may be related to the presence of cementitious bonds or to interparticle friction, possibly due to aggregates not fully coated by the bituminous films formed with emulsion setting.
On the one hand, the proposed Eq. 10 is limited because it is applicable only in the frequency domain. On the other hand it is mathematically simple and can be easily extended by replacing the HS model with other LVE models like, for example, the 2S2P1D or the Havriliak-Negami models [47], which are currently used to simulate the behavior of HMA materials.

Comparison between mixture and FAM
First, we highlight that the same values of k (0.122) and h (0.314), along with a slightly variable value of d Fig. 12 Comparison between the new model and the HS model (1.634-2.000), could be used to simulate the behavior of all the mixture and FAM mortar specimens, at both curing conditions (Table 2). This means that the LVE behavior of the mixture was due to its FAM, which is reasonable since the FAM contained most of the bitumen of the mixture. In fact, only the aged bitumen film coating the coarse RA aggregates was not comprised in the FAM and it clearly played a minor role.
For HMA mixtures, the values of k and h normally range between 0.15 and 0.18 and between 0.44 and 0.70, respectively [19,[48][49][50]. The values obtained in the present research, as well as those reported in other studies on CRM mixtures, are clearly lower [12,21]. From the model point of view, this indicates that the two parabolic dashpots of the model are closer to a Hookean spring. In other terms, the rheological behavior of CRM mixtures and FAM mortars was characterized by a less marked viscous component with respect to HMA. This suggests that the aged RA bitumen and the fresh bitumen from the emulsion jointly affected the LVE behavior.
Curing did not affect k and h, while d had only a small variation from 14 to 28 days. However, we observed (Fig. 8) that the effect of curing is similar to that of frequency and temperature and therefore it should be visible mainly in the parameter s, which is a frequency multiplier (Eq. 11). In fact, for all tested specimens, the curing time increase led to an increase in s. This is particularly evident for the two mixture specimens, where the variation is about one order of magnitude (Table 2). From a physical point of view, higher values of s characterize materials which are harder and have lower relaxation ability. The corresponding increase in stiffness modulus and the phase angle reduction has been already highlighted in Fig. 9. An increase in s was also observed when the fraction of RAP-extracted binder was increased in blends of RAP-extracted binder with virgin binder [51,52]. This suggests that the longer curing period (14 days at 40°C) caused some aging of the bituminous binder films that formed upon emulsion breaking and setting.
As regards the phase angle d 0 , describing the frequency-and temperature-independent dissipation, the mixtures showed values from 1.789°to 0.945°, which are similar to those obtained for CRM mixtures tested 8 years after construction [12]. The d 0 values were actually very small, but this is not surprising and simply indicates that, for the given composition, the thermo-viscoelastic dissipation due to the bituminous component certainly prevailed (i.e. the material showed an asphalt-like behavior). For FAM mortars d 0 was even smaller, suggesting that the addition of coarse aggregate and voids was responsible for most of the frequency-and temperature-independent dissipation.
The FAM mortar specimens had higher values of the equilibrium modulus, with respect to mixtures; as regards the glassy modulus the two materials showed similar values. For HMA mixtures, E e and E g are associated to the contribution of the aggregate skeleton and the air voids content [53,54]. From our results, following this interpretation, we may conclude that the reinforcing effect of the coarse RA aggregate was almost perfectly balanced by the presence of the independent voids phase. However, in the hightemperature/low-frequency range (log f r ! À1) we need to be very careful. In that range the bitumen contribution (fresh and aged) to the stiffness becomes negligible, but CRM mixture and mortars also contain cementitious bonds, whose effect probably determine the value of E e .
The rheological behavior of mixture and FAM was also compared using normalized complex modulus curves [53]. We applied the normalization only to the LVE part of the behavior (E Ã exp Àjd 0 ð Þ): Using this definition, we can compare the LVE behavior of mixture and FAM without the influence of E e , E g and d 0 . In Fig. 13 normalized models are superposed to normalized measured data and the superposition between mixture and FAM is very good. The small effect of curing (difference between 14 and 28 days) is highlighted as a small increase in both d and s.

Conclusions
The objective of this study was to compare the rheological behavior of CRM mixtures and FAM mortars. We monitored the curing process by measuring the water loss by evaporation, ITS and ITSM, and performed complex modulus tests at 14 and 28 days of curing. Based on the experimental results and analysis, our conclusions are as follows: • The mixture and the FAM mortar had the same rate of water loss, and their stiffness evolution was strictly linked. In the first curing days the mixture was stiffer, but ultimately the FAM became stiffer.
We used the Hirsch model to simulate this behavior and found that the series arrangement of FAM and coarse RA aggregate prevailed on the parallel arrangement. • The complex modulus results showed that, despite the huge difference in volumetric composition, the thermo-rheological behavior of mixture and FAM mortar was very similar. The effect of curing (increase in stiffness modulus, decrease in phase angle) was similar to that of frequency and temperature. This could be due to aging phenomena that are starting to occur within the bituminous binder films formed upon emulsion setting. • The rheological behavior of mixtures and FAM mortars was characterized by the presence of a small time-and temperature-independent dissipation component. Therefore, the complex modulus results could be successfully simulated with a new model which, besides viscous dissipation, considers a hysteretic dissipation component. The latter may be related to the presence of cementitious bonds or to interparticle friction, possibly due to aggregates not fully coated by the bituminous films formed upon emulsion setting. • The new model showed that, similar to HMA mixtures, the rheological behavior of the CRM mixture was controlled by its FAM. The coarse RA aggregate fraction and the voids mainly affected the asymptotic properties (equilibrium and glassy moduli) and the non-viscous dissipation component.
Author contributions Laboratory testing was performed by SR. Data collection and analysis were performed by SR, CM and AG. The first draft of the manuscript was written by AG and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.

Compliance with ethical standards
Conflict of interest The authors declare that there is no conflict of interest regarding the publication of this article.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.