1 Introduction

Bituminous binders undergo changes in properties with time, even in the course of annealing, storage and/or use at room temperature—a phenomenon long known from practice and often called ageing. The underlying microstructural changes may be connected to two different processes occurring in parallel: (1) chemical reactions resulting in irreversible processes e.g. UV-irradiation assisted oxidation; leading to an increase and/or decrease in molecular size of single components of the bitumen and (2) physical aging or reversible processes such as re-arrangements of molecular aggregates, aggregation of molecules and crystallization of components within the mixture of compounds within the bitumen [1]. Both processes lead to an increase in hardness or to stiffening in time, the first one in a rather long time-scale (years). The second process is usually much faster (minutes up to several months).

The study of the process of physical hardening during extended storage at isothermal conditions dates back to 1936 when a report by Traxler and Schweyer [2] was published. Different asphalt binders were investigated and a substantial increase in viscosity with time of storage at room temperature was observed. The phenomenon was called in this time “age hardening”. In a subsequent study, the process of increasing viscosity during isothermal storage was explained by the structural sol–gel transition [3] according to the micellar model of bitumen structuring [4].

Increase of the stiffness of several asphalt binders after isothermal storage was found in a later investigation by Brown et al. [5]. This process was called “steric hardening”. The effect was explained by the collapse of the molecular free volume at low temperature, an isothermal sol–gel transition and wax crystallization. The term “physical hardening” of asphalt binders was first used by Blokker and Van Hoorn [6] who explained the process by wax crystallization and asphaltene aggregation.

In the current view, physical hardening which may be also termed as “spontaneous hardening” is connected with the establishment of (new) interactions between polar and apolar components leading to the formation of discrete crystalline domains as observed earlier. This process takes days to weeks if bitumen is stored at room temperature [5, 7, 8] respectively at low temperatures, but still above T g [912] and is basically a thermally activated process. Hardening effects during the cooling of a bitumen melt are also discussed by Masson et al. [8]. In this study, a four stage ordering within bitumen upon cooling from the melt was suggested. In the first stage, all fractions order rapidly into a weakly organized phase. In the second stage, low molecular weight saturated segments crystallize.

In the third stage, high molecular weight saturated segments crystallize. In the fourth stage, resins and asphaltenes order into a mesophase. The third and fourth stages are responsible for the room-temperature (steric) hardening of bitumen. The formation of microstructures during cooling down from an isotropic liquid as initial condition at high temperatures as well as the gradually melting during heating from a two phase structure show good agreement with the thermal signals recorded by differential scanning calorimetry (DSC) in a recent study [13].

It is certainly of interest to quantify the kinetics of the formation of the micro-phases during cooling of bitumen as it occurs in the applications of bitumen, e.g. a road pavement, to better predict performance of asphalt in models and to understand the functional properties and their development.

1.1 Experimental part

The bituminous binders used in this study are labelled Q8 70–100 and Total 70–100 and were supplied by Kuweit Petroleum Research & Development (Rotterdam, Netherlands) and Latexfalt B. V (Koudekerk, The Netherlands). The first bitumen was derived from Kuweit Crude Export. The penetration (PEN) value was 82, the softening point 45.8 °C and the Fraass breaking point −15 °C. The content of SARA fractions is as follows: Saturates: 3.8 %, Aromatics: 59.6 %, Resins: 22.2 %, Asphaltenes: 14.6 %. The other bitumen had similar properties with a content of SARA fractions is as follows: Saturates: 5.9 %, Aromatics: 29.5 %, Resins: 47.0 %, Asphaltenes: 17.6 %. Bitumen films were prepared by the application of a bead of bitumen onto a steel disk, followed by a heating stage at 120 °C for 1–2 min in an oven to ensure a low viscosity of the binder. The resulting film thickness was about 300 µm. Prior to the AFM experiments, the sample was quenched with a rate of 30 K/min to 60 °C. Atomic Force Microscopy observations were performed on a NanoSurf Easyscan 2 apparatus operated in Tapping® mode combined with a Linkam TP 94 microscopy hot stage for heating/cooling and annealing experiments. The Nanosensor silicon SPM probes sensors (type XYNCHR) had tip radii less than 10 nm, the nominal spring constant and resonance frequency of the cantilever were respectively, 42 N/m and 330 kHz. The height and phase contrast images were recorded with a set-point amplitude of 60–65 % in order to minimize the interaction of the tip with the soft and sticky sample surface (soft tapping). A typical scan size was 12.5 × 12.5 µm; the recording time per frame was approx 2 min. For the isothermal calorimetric experiment, a TA instruments DSC Discovery Series was used. The sample was heated to 120 °C, held for 10 min and quenched to 60 °C and isothermally hold for 60 min.

2 Results

Kinetics of polymer crystallization and or demixing of polymers in a blend are “classically” investigated by means of isothermal DSC [14, 15].

The performed isothermal calorimetric experiments however did not provide suitable data, a treatment of the data to extract crystallization kinetics as usually applied for polymer crystallization [14, 15] is not feasible due to the extremely small heat effects connected with the steric hardening and the different and widely spread time constants (see Fig. 1) for the occurring processes (minutes to weeks).

Fig. 1
figure 1

Heat flow recorded during an isothermal DSC experiment of bitumen Q8 70–100 at 60 °C after quenching with ca. 30 K/min from 120 °C

However, the development of the microstructures during an isothermal experiment can be followed by AFM. For AFM-experiments in a fast changing environment, quick scanning is essential to prevent creeping phenomena distorting the scanning during recording a snapshot of the crystallization process. The fastest scanning rate proved to provide acceptable images, which are applicable for studying crystallisation kinetics. The pictures obtained are subsequently linked to the isothermal DSC scan (see Fig. 2) and analyzed with respect to the distribution of surface area belonging to the catana/peri phase (i.e. the crystalline ‘particle’ phase) and to the perpetua phase (i.e. the more viscous ‘matrix’ phase).

Fig. 2
figure 2

Heat flow recorded during an isothermal DSC experiment of bitumen Q8 70–100 at 60 °C after quenching with ca. 30 K/min from 120 °C combined with the by AFM recorded phase contrast pictures

The different stages of the process of physical ageing; nucleation and growth of the disperse phase (catana/peri-phase) as well as the diffusion and growth of its particles can be regarded as being similar to those processes occurring during crystallization in solutions which are supersaturated by cooling. Consequently, the kinetics of the formation of the equilibrium structures in bitumen can be treated within the framework of (heterogeneous) nucleation and crystallization theory.

The kinetics of isothermal crystallization of melts of low molecular compounds as well as of polymers can be described satisfactorily by the Avrami Eq. (1):

$$x_{\tau } = x_{\infty } \left( {1 - {\text{e}}^{{ - k\tau^{n} }} } \right)$$
(1)

With x τ as the degree of crystallization at the time τ; x as the degree of crystallization at the end of the studied crystallization process, k as the crystallization rate constant; n is a coefficient depending on the growth of the crystals (the Avrami exponent) and hence describing the mechanism of formation.

The experimentally determined value of the Avrami exponent allows assumptions concerning the mechanism of crystallization.

Following Lorenzo et al. [10], this exponent can be considered to be composed by two terms (2):

$$n = n_{\text{d}} + n_{\text{n}}$$
(2)

with n d representing the dimensionality of the growing crystals (can only have integer values 1–3) and n n as the time dependency of nucleation (in case of homogeneous nucleation 1 and of heterogeneous nucleation 0).

However, since in some cases there is a mixed nucleation also non-integer values of the Avrami exponent may be obtained. The larger the value of n is, the more perfectly the crystals are formed. For processes with n equal to 3 or 4 it is characteristic to find the growth of three-dimensional spherolites from instantaneously forming or gradually developing nuclei.

Previous investigations of crystallization kinetics of bitumen binders in bulk covering a time scale of days and using dilatometry resulted in a value of the Avrami parameter equal to 1, independent on the composition or hardness (PEN-grade) of the bitumen samples [16]. It was observed, that the change in bitumen density occurred most rapidly at 50 °C i.e., at the temperature which is also characteristic for the maximum change in shear stress (softening). The growth of the dense structures was described to be controlled by diffusion processes in the system as characteristic of a crystallization process. However, the degree of order and the perfection of the particles of the new phase which was formed were not high. Complete crystallization was absent, the observed change in density was small (about 0.55 %) and the time constant was about 6 h [16].

The degree of crystallization x τ can be determined from the ratio of peri-to perpetua-phase, while analyzing changes in dimensions of the forming and growing of crystals of the peri-phase.

Looking at the data obtained from AFM phase contrast pictures of the samples during the course of isothermal crystallization, firstly the fast development of flat, sheet-like crystals in time can be observed. These will be also the very first signs of crystallization since the temperature drop at the air interface will be the largest and the crystals formed by saturates display a lower density than the bulk bitumen (Fig. 4). However, since the overall degree of crystallization is rather small, there will be enough area for the crystals to form at the surface without interfering with each other.

The course of crystallization can be therefore described by using the development of relative area as registered by the phase contrast signal during scanning using the AFM instead of determining the volume of the growing crystalline phase in time, Eq. (3) can be applied:

$$x_{\tau v} = \frac{{V_{\tau } - V_{\infty } }}{{V - V_{\infty } }}\sim x_{\tau s} = \frac{{A_{\tau } - A_{\infty } }}{{A - A_{\infty } }}$$
(3)

with x τv as as the degree of crystallization in bulk (volume) and x τs as as the degree of crystallization as observable on the surface.

The shortest time to record a complete image was approx. 2 min; pictures were taken immediately after cooling down (quenching) from 120 °C with a speed of ≥30 K/min and reaching 60 °C. The development of crystallinity (area of crystalline catana/peri phase) in time is plotted in Fig. 3.

Fig. 3
figure 3

Development of surface crystallinity in time and fit according to Eq. (1)

The data obtained from the phase contrast images can be satisfactorily fitted using Eq. (1). Additionally, the kinetic constant can be calculated from the fit, a value of k = 5 min. Also the half time of initial crystallization τ 1/2 ~3 min. for the bitumen Q8 and 2 min 25 s for the Total bitumen was found.

The variation in n-value (Avrami parameter during the isothermal crystallization can now be computed according to Eq. (4).

$$n = \frac{{ - \tau \left( {\frac{{{\text{d}}x_{\tau } }}{{{\text{d}}\tau }}} \right)}}{{\left[ {\left( {x_{\infty } - x_{\tau } } \right)\;\ln \;\left( {1 - \frac{{x_{\tau } }}{{x_{\infty } }}} \right){\kern 1pt} } \right]}}$$
(4)

and is plotted in Fig. 4 together with schematical pictures for the evolution of crystallisation in time.

Fig. 4
figure 4

Variation in n value during formation of the catana/peri phase at 60 °C and schematical pictures showing the evolution of crystallisation and of the crystal morphology/location of crystalline phase in time

3 Conclusions

The variation in n value with increasing relative crystallinity provides insights into the mechanism of crystallisation and the result with respect to crystal structure and perfection. It can be seen in Fig. 4 that the Avrami exponent changed progressively from 4 to 1.8 during the initial development of crystallinity with an average value of about 2.6.

Crystallisation starts almost immediately after or instantaneously during cooling. An earlier study by Masson et al. [17]. proposed liquid liquid demixing within bitumen upon cooling resulting finally in steric hardening. The beginning of the crystallization is homogeneously nucleated and three-dimensional, which changes in the course of crystallization to 2-dimensional crystal growth and mixed nucleation, both homogeneous and heterogeneous. The primary crystallization process is very fast and followed by subsequent hardening events as reported by Pechenyi and Kusnetzov [16] and by Masson et al. [8, 17]. This especially the case in the Total bitumen sample, which crystallizes very fast and reaches also an initially higher degree of crystallinity compared to the Q8 sample. In this case, the originally homogenously 3 dimensionally nucleated crystals changes to 2-dimensional crystals, which grow in a very short time.

The analogous nature of the kinetics of formation of the equilibrium structures in bitumina and of the crystallization process makes it possible to regard both these processes as first order transitions. The equilibrium structures of the peri-phases should consequently be classified as semi-crystalline. The formation of the peri-phase is connected with an increase in stiffness/hardness and in density. Consequently crystallization causes deformation creep and internal stresses, which promotes cracking.