Influence of Hot Band Annealing on Cold-Rolled Microstructure and Recrystallization in AA 6016

The influence of an intermediate heat treatment at the end of hot rolling and before cold rolling on Cube texture formation during the final solution annealing of AA 6016 is investigated. Three hot bands with different initial grain sizes and textures are considered: the first one without annealing before cold rolling, while the other two hot bands are heat treated at 540 °C for 1 hour in air before being cold rolled. One of the heat-treated hot bands was left to cool down in air and the other inside the furnace. Electron-backscatter diffraction (EBSD) maps of the cold-rolled specimens and crystal plasticity simulations show no difference in the amount of Cube remaining in the microstructure at the end of cold rolling for all three specimens. The initial grain size of the hot band has no influence on the Cube texture fraction left in the microstructure at the end of cold rolling for thickness reductions higher than 65 pct. Nevertheless, the grain size of the hot band affects the shape and distribution of the Cube grains left in the microstructure and the kernel average misorientation in the cold-rolled specimens. Moreover, the heat treatment decreases the intensity of the beta fiber components (Brass, Copper, and S) in the hot band and promotes the formation of a cold-rolled microstructure with a low kernel average misorientation. Both these factors lower the probability of preferential Cube nucleation during solution annealing and keep the Cube volume fraction after recrystallization below 10 pct, while it reaches 25 pct without intermediate annealing.


I. INTRODUCTION
OPTIMIZATION and tailoring of texture is fundamental to reduce anisotropy of aluminum sheets during forming, especially in 6xxx alloys because of their use in automotive components, where surface finish and shape must satisfy strict requirements. Alloys of 6xxx series are of interest in automotive applications because of their good combination of formability, corrosion resistance, and increase of strength during paint baking after forming. [1] Texture evolution depends on casting, [2] hot and cold rolling parameters (strain rate [3] and temperature), deformation history, [4][5][6] alloy composition, [7] and the presence of precipitates. [8,9] After cold rolling, the 6xxx sheets are solution annealed (SA) at a temperature between 500°C and 570°C followed by quenching to dissolve the hardening phases and retain the corresponding alloying elements in solid solution. [1,4] During this solution treatment, recrystallization of the cold-rolled sheet occurs resulting in a modified texture. Cube is considered a critical texture component during the recrystallization treatment because it can show preferential nucleation and growth. [10] Cube preferential growth during recrystallization is found mostly in high-purity aluminum alloys [11] ; however, also 6xxx alloys stabilize a high fraction of Cube during hot rolling which can be difficult to remove during successive cold rolling. [11,12] This fraction of Cube left in 6xxx alloys can show preferential nucleation during the final solution annealing. Recrystallized Cube grains can originate from transition bands, [13] which are thin Cube seeds surrounded by high orientation gradients and that are caused by diverging rotations of unstable orientations according to the Dillamore and Katoh model. [14,15] However, also Cube fragments left in the microstructure at the end of cold rolling and not exactly corresponding to transition bands can act as nucleation sites. [16,17] During in situ annealing experiments and modeling of commercial AlFeSi alloy, Sukhopar and Gottstein [18] have shown that the intensity of Cube texture after recrystallization is due to nucleation at the Cube bands left in the deformed microstructure, while contribution to Cube nucleation from locations outside the Cube bands is negligible. Recrystallization in high purity and highly deformed (98 pct cold rolled) aluminum foil has shown that occurrence of extensive recovery and continuous recrystallization can retard the formation of Cube texture. [19] The recrystallization advantage of Cube grains decreases rapidly when the Cube present in the microstructure is misoriented at least 10 to 15 deg from its exact orientation. [20] Thus, several works were dedicated to increasing the Cube fragmentation during rolling through asymmetric cold rolling [4] or with the adjustment of shear loading. [21,22] Some experimental investigations have suggested that a fraction of Cube texture could be created by S-oriented grains during compression, [23] while elasto-viscoplastic fast Fourier transform (EVP-FFT) simulations have shown that depending on the amount of compression strain, some Cube can originate from random orientations or non-Cube grains misoriented within 10 to 20 deg from ideal Cube, [24] but the fraction of this newly formed Cube during rolling is extremely low and it does not play a key role during recrystallization.
Solutions to reduce the amount of Cube already during casting were suggested through the use of continuous casting (CC) instead of direct casting (DC). For the same aluminum alloy, CC produces ND misoriented Cube microstructures, while DC forms microstructures containing a higher fraction of perfectly oriented Cube. [2] During cold rolling, Brass, Copper, and S develop faster in hot bands from CC than DC; moreover, microstructures from DC maintain a higher fraction of exact Cube orientation at the end of cold rolling. [25] In several aluminum alloys, the amount of Brass texture was increased and the Cube fraction was lowered by tailoring the alloy composition. [7,26,27] Cube recrystallization has been the object of several investigations at microtexture level using bicrystal plane strain compression tests. [28,29] Initially, the attention was focused on the role of grain boundaries and Cube transition bands in the nucleation of new grains. In particular, subgrain coalescence in transition bands helps Cube nuclei in reaching early the critical size to grow. [30] More recently, the influence of different amounts of iron content, deformation strain, and recrystallization temperature on Cube texture in AA 1145 has been presented. [31] The recrystallized texture depends on a competition between Cube, R orientation, and particle stimulated nucleation (PSN). Preferred Cube nucleation is obtained at low strain and low recrystallization temperature because the fast Cube recovery prevails over recrystallization of other texture components; while R-oriented grains nucleate preferably at grain boundaries whose density is favored by high strain deformations. [31] Another microstructure heterogeneity that influences Cube recrystallization is represented by shear bands. In Al-Mg alloys, shear bands possess high-angle grain boundaries that can already migrate at low temperatures (200°C to 350°C) and favor the nucleation of new grains. [32] At these temperatures, the nucleation at shear bands is found to be more rapid than formation of nuclei at other microstructure heterogeneities. [32] Cube recrystallization has been found to be affected by the presence of concurrent precipitation in supersaturated aluminum alloys. [33] In Al-xMn-Fe-Si alloys, concurrent precipitation enhances the presence of P and ND rotated Cube grains, while it suppresses or retards nucleation and growth of exactly oriented Cube. [33] These results were confirmed by Tangen et al. [34] who investigated the influence of concurrent precipitation during recrystallization in a manganese supersaturated aluminum alloy. Tangen et al. [34] concluded that concurrent precipitation leads to an inhomogeneous coarse-grained microstructure which is dominated by P and ND rotated Cube orientations.
The influence of different textures and microtexture for the same alloy system (AA 1050) on Cube recrystallization was investigated by comparing Cube nucleation and growth of plane strained compressed (PSC) microstructure and equal channel angular extruded (ECAE) microstructures. [35,36] Despite the fraction of Cube grains was the same after both deformation routes, a higher intensity of Cube texture was found after recrystallization in plane strain compressed microstructures because of the presence of strong S texture. [35] Recrystallization of microstructures from ECAE deformation did not show Cube growth because of the absence of S orientation in the deformed microstructure. [35] Recently, the influence of pressure on recrystallization of AA 3003 alloy was investigated proving that the presence of pressure during the thermal treatment can refine the recrystallized grains, weaken the pinning effect, and strengthen P and ND rotated Cube. [37] The purpose of this work is to present a solution to limit Cube recrystallization in AA 6016 by introducing an intermediate annealing at the end of hot rolling and before cold rolling. The aim is to investigate how this intermediate annealing changes the microstructure (grain size and texture) during cold rolling and final solution annealing. It is important to assess if there is an influence of hot band grain size and texture on Cube texture evolution in AA 6016 during cold rolling and to experimentally see how the cold-rolled microstructures obtained from different hot bands change after the final solution annealing. Microstructures right after industrial cold rolling are compared with textures from full-field crystal plasticity simulations. This simulation approach allows to investigate the influence of initial hot band texture and grain size on Cube texture evolution during cold rolling. Analysis of distribution and fragmentation of Cube orientation is performed both on experiments and simulations.

A. Experimental Methods
The aluminum AA 6016 used for the hot bands had the composition is shown in Table I RD-ND sections from the three hot bands were prepared by grinding, mechanical polishing with diamond suspensions and oxide polishing suspension (OPS, Struers) with contact forces up to 50 N for at least 18 min. This final polishing was followed by electrochemical polishing using STRUERS A2 electrolyte at 10°C for 4 to 8 seconds and a voltage of 36 V. This procedure has been already used by the authors for the preparation of several AA 6016 specimens. [38] Electron-backscatter diffraction (EBSD) analyses were performed using a scanning electron microscope (SEM; JEOL 7200F FEG-SEM, Tokyo, Japan) equipped with an EBSD-measurement system (Nordlys Nano and Symmetry S2, Oxford Instruments, Abingdon, UK) with a 70 deg tilt specimen at 20 kV. For hot band A, the EBSD scans covered the whole sheet thickness (area of 2000 lm 9 7000 lm with step size of 3 lm). For hot bands B and C, three measurements each of size 820 lm 9 7000 lm with a step size of 4 lm were performed for each sample to obtain a sufficient amount of grains orientations to evaluate texture due to the larger grain size in hot band B and C.
All hot bands with an initial thickness of 7.0 mm were cold rolled down to 1.3 mm (corresponding to a total thickness reduction of % 82 pct). Cold rolling was performed on the industrial plant at AMAG rolling GmbH.
EBSD analyses were performed on RD-ND sections of the cold-rolled specimens to evaluate the Cube texture distribution through thickness and the texture evolution obtained from the three hot bands. Cold-rolled samples for EBSD were prepared through standard grinding with 1000 and 2000 SiC papers before proceeding with 3 lm polishing with diamond suspension. The final polishing was done manually for approximately 1 hour with OPS suspension (Struers) with the addition of 10 pct of hydrogen peroxide (H 2 O 2 ). EBSD maps were obtained using a field emission SEM-JEOL 6500 F operated at a voltage of 20 kV. For each cold-rolled sample, at least three EBSD scans of the full sheet thickness were obtained using a step size of 1.5 lm.
The cold-rolled samples from different hot bands were then solution annealed (SA) in a salt bath for 5 min at a temperature of 520°C followed by water quenching. Fully recrystallized microstructures were obtained in each specimen, and EBSD analyses were performed on two areas of at least 1500 lm 9 800 lm with step size of 3 lm for each sample.
For texture evaluation, unless otherwise stated, orientation distribution functions (ODFs) were calculated using MTEX version 5.2.3 and De la Vallee Poussin kernel with a half-width of 5 deg.

B. Computational Methods
Crystal plasticity simulations were performed to track the influence of initial grain shape of the hot band on Cube texture evolution. For the three hot bands, high-resolution EBSD scans were performed at the center of the aluminum sheet. The microstructure obtained from these EBSD scans was used to generate starting 2D-RVEs to be used for crystal plasticity simulations with the freeware software package DAMASK. [39] The input files were created by writing a pipeline in DREAM.3D software. [40] A first simulation compared an RVE from hot band A with an RVE from hot band B with the same initial Cube fraction. Another simulation was performed on an RVE from hot band A and RVE from hot band C again with the same initial Cube volume fraction. The RVEs were then deformed in plane strain compression applying the compression history given in Table II, while the component of the deformation gradient F TD-TD is equal to 1.0 and the other components F RD-TD , F RD-ND , F TD-ND are all set to zero. During deformation of the polycrystalline RVE, the initial grid of material points becomes distorted because of deformation heterogeneity in each grain. When the distortion becomes too large, the stress equilibrium cannot be found, and convergence is not reached. Thus, to simulate compression higher than 40 to 50 pct of thickness reduction, a remeshing scheme is implemented in DAMASK. During remeshing, a new undistorted grid is created and the state variables from the previous deformation step are mapped onto the newly created grid using a nearest-neighbor mapping algorithm. The next cold rolling step is then performed with the new undistorted grid. During remeshing, the elastic deformation is relaxed. A detailed description of the implemented remeshing scheme is described by Sedighiani et al. [41] and its influence on texture prediction was already assessed, showing that the number of remeshing step does not affect texture evolution. [3] 1. Constitutive model for crystal plasticity simulations The plastic deformation is described through the plastic deformation gradient (F p ) whose evolution is given as: In Eq. [1], the plastic velocity gradient (L p ) can be formulated as sum of the shear rates on all slip systems as: where m a represents the slip direction and n a is the normal to the slip plane of slip system a, while _ c a is the shear rate of the ath slip system and n is the number of slip systems. For cold rolling simulations, the 12 octahedral slip systems (i.e., {111} h110i) are considered in the calculation of the plastic velocity gradient.
For the shear rate, the constitutive behavior of aluminum was described by a pheno-power law, whose parameters are listed in Table III and whose description has already been used for aluminum texture simulations and prediction. [3,42] In Table III, C 11 , C 12 , and C 44 are the elastic constants, _ c 0 is the reference shear strain rate, m is the strain rate sensitivity exponent, s b c is the critical shear stress for slip system b, and s s , a, and h 0 , represent the saturation value of slip resistance, the hardening exponent, and the initial hardening coefficient, respectively. The latent hardening (q ab ) is set to 1.0 for coplanar slip systems and 1.4 for non-coplanar slip systems to make the hardening model anisotropic. In phenomenological constitutive models, the shear rate on slip system a is a function of the critical resolved shear stress (s a c ) and the resolved shear stress for slip system a (s a ) according to: where _ c 0 is the reference shear strain rate and m is the strain rate sensitivity exponent. The evolution of the critical resolved shear stress is expressed as: This equation describes the influence of any set of slip system (b) on the hardening behavior of a (specific) slip system a through the hardening matrix (h ab ). This hardening matrix captures the micromechanical interaction among different slip systems and it is given by: with h 0 , a, and s s representing the self-hardening coefficient, the hardening exponent, and the saturation value of slip resistance, respectively. These parameters are assumed to be identical for all FCC slip systems. The parameter q ab is calculated as: where d ab is the Kronecker delta and q self and q lat are the self and latent hardening, respectively. The hardening matrix contains diagonal terms (q self ) and off-diagonal terms (q lat ). The latent hardening is set to 1.0 for coplanar slip systems and 1.4 for non-coplanar slip systems to make the hardening model anisotropic.

III. RESULTS AND DISCUSSION
A. Hot Rolling Microstructure and Texture The hot band without heat treatment has thin grains elongated along RD, while the annealed hot bands show large grains. From the (111) pole figures, hot band A has a high density of Cube texture with a fully developed beta fiber (Figure 1(d)), while hot bands B (Figure 1(e)) and C (Figure 1(f)) show a high intensity for Cube texture or random orientations, but the beta fiber is suppressed because of the annealing. Figure 2 shows the distribution of Cube grains through the whole thickness for the three hot bands and Figure 3 shows the volume fractions calculated considering up to 15 deg maximum misorientation of the relevant orientations listed in Table IV for the three hot bands. In Figure 3, the amount of Brass, Copper, and S are low in hot bands B and C confirming that the b fiber components have been attenuated by the annealing.
Annealed hot band textures are dominated by Cube, Q and 45 deg ND rotated Cube, while the fractions of beta fiber orientations (Brass, Copper, and S) and R are low, indicating that during annealing before cold rolling, the recrystallization inside the hot-rolled specimen is driven by nucleation at cube bands and shear bands. The contribution of Cube transition bands, strain-induced boundary migration, shear bands, and PSN on texture development during annealing and recrystallization has been investigated by Lu¨cke and Engler [8] and Engler. [43] It was concluded that Cube emerges from nuclei at transition bands that show preferential growth into S and rolling texture components (like Copper), while formation of 45 deg ND rotated Cube and Q orientations is caused by nucleation at shear bands (preferentially shear bands pertaining to S orientations). [43] B. Cold Rolling Microstructure and Texture Figure 4 shows one RD-ND scan and the corresponding Cube distribution map for the cold-rolled samples from each hot band. The other EBSD maps obtained from cold-rolled specimens are included in ''Appendix'' (Figure A-II). In Figure 4, the volume fraction of Cube is comparable for all cold-rolled samples and this suggests that the Cube grain shape and size initially present in the hot bands have no influence on the fraction of Cube texture. However, the Cube morphology and distribution appear different. The cold-rolled sample from hot band A has thin Cube seeds distributed almost through the whole thickness ( Figure 4(b)), while samples from hot bands B and C show thick Cube bands only in some zones across the thickness of the sample (Figures 4(d) and (f)). A similar impact of the hot band grain size on the Cube shape and distribution left in the microstructure at the end of cold rolling was described by Engler. [44] Fine grained hot bands yield to a large number of Cube transition bands parallel to the rolling direction. These transition bands are favorable interfaces for Cube nucleation and, in fine grained samples, they are narrow spaced in the RD-ND section making the Cube nucleation and growth more efficient than PSN or nucleation on shear bands. [44,45] Contrary, in samples from coarse-grained hot bands, few Cube interfaces are present in the RD-ND sections and even after heavy deformation the grains remain thick in the microstructure. As a consequence, Cube bands are too far apart and have too few interfaces to produce many nuclei, and for these coarse microstructures, PSN will be the predominant recrystallization mechanism. [44] Figure 5 shows the fraction of relevant texture components found on each EBSD map (in Figure 4 and in ''Appendix'', Figure A-II) for the cold-rolled samples obtained from the different hot bands. The volume fraction of Cube ranges from 1.5 to 5 pct for all samples suggesting that the larger grain shape in the hot band does not promote more Cube fragmentation during cold rolling than elongated thin grains, but it influences the shape of the remnant Cube grains inside the microstructure. A higher fraction of Brass and S is found for specimens cold rolled from hot band A, while cold-rolled sheets from the annealed hot bands (B and C) exhibit higher fraction of Copper. A higher fraction of Copper at the expenses of S orientation was also found in overaged Al-1.8 pct Cu alloys rolled from hot bands with coarse grain size. [44] The development of rolling texture (i.e., b fiber components) is retarded in coarsegrained materials because a higher level of deformation heterogeneities can occur in such microstructures including deformation bands, shear bands and splitting of  orientations which retard the development of a sharp texture. [44,46] C. Solution Annealing Microstructure and Texture Figure 6 shows one RD-ND section of an EBSD scan for each specimen and the corresponding map of the Cube grain distribution for SA samples. Further EBSD maps used to evaluate texture volume fractions are included in ''Appendix'' (Figure A-III). The Cube distribution maps show more grains pertaining to the exact Cube orientation in the sample from hot band A. However, no preferential growth of Cube occurred as the grain size is the same in all samples, but the higher fraction of Cube in the sample from hot band A seems to be caused by preferential nucleation. Preferential Cube growth with Cube grains much coarser than grains pertaining to other orientations is often found in high-purity aluminum and the strength of Cube texture is known to decrease with alloy content. [47] The recrystallized microstructures in Figure 6 all show similar grain sizes. According to literature, [44] this  is an indication of texture being determined by the superposition of recrystallization occurring at various nucleation sites (Cube bands, shear bands, and PSN). Figure 7 presents the volume fraction of texture components after solution annealing. The cold-rolled sample from hot band A shows the highest Cube fraction among all samples and also preserves a high fraction of S orientation. Samples from hot bands B and C have 8 to 9 pct of Cube fraction, and the volume fractions of the b fiber components (Brass, Copper and S) do not overcome 5 pct. In Figure 7, the higher Cube and S fractions in the sample from hot band A can be caused by a higher fraction of S orientation and b fiber present in the cold-rolled specimen. Cube shows preferential nucleation when surrounded by S orientation. [43] R orientation prevails in the microstructures from hot band A and it can be generated by strain-induced boundary migration (SIBM) but also it is preferentially nucleated in microstructures with high fractions of S texture because of a 40 deg h111i relation with S. [43] Moreover, R is assumed to easily nucleate at grain boundaries. [44] Samples from hot bands B and C have higher fraction of Q and P orientations than the specimen from hot band A. Predominance of P and 45 deg ND rotated Cube in samples from hot bands B and C can be a consequence of PSN. [43,44] However, it was also shown that strong P and Q components can occur at the expenses of Cube when nucleation at shear bands is predominant. [44,48] Microstructure from hot band A has a higher fraction of 22 deg RD rotated Cube which can be attributed to nucleation at Cube transition bands, [44] while 22 deg TD rotated Cube fraction is higher in microstructures from hot bands B and C. This higher fraction of 22 deg TD rotated Cube could be a result of PSN. In AA 3104 alloys, EBSD data have shown that dispersoids with major size lower than 0.6 lm do not affect texture evolution, while significant orientation gradients are found around coarse particles. These orientation gradients are characterized by lattice rotations in TD direction. [49] D.

Crystal Plasticity Simulations Results
For crystal plasticity simulations, two zones were extracted as input from hot band A. The first zone named hot band A (1) which has the same amount of Cube fraction as the initial input microstructure from hot band B. The second zone named hot band A (2) which has the same amount of Cube fraction as the initial input microstructure from hot band C. Figure 8 shows the initial microstructure hot band A (1) and the microstructure from hot band B. Another set of simulations was done taking an initial microstructure hot band A (2) and a microstructure from hot band C (Figure 9), and the evolution of Cube volume fraction was extracted from crystal plasticity simulations. Both Figures 8 and 9 indicate that at 80 pct thickness reduction, the amount of Cube remaining in the microstructure is negligible or even absent according to the simulations.
The evolution of Cube fraction as function of thickness reduction is shown in Figure 10. For a thickness reduction higher than 70 to 80 pct, the Cube fraction does not depend on the initial grain size of the hot band. Moreover, as shown in Figure 10(b), it can also occur that a higher fraction is predicted for the large grain microstructure than for the elongated grains microstructure. This simulation result confirms the experimental fractions shown in Figure 5 where the Cube fraction varied in the same range for all cold-rolled samples. Moreover, in Figure 10(a), for compressions lower than 50 pct, elongated grains show higher Cube fraction than large grains, the opposite holds in Figure 10(b). Thus, no conclusion can be drawn on the effect of initial grain size of hot bands on Cube fraction evolution during cold rolling. A crystal plasticity simulation was repeated by using a 3D RVE with the same initial texture and amount of Cube, but one RVE had elongated grains and the other had equiaxed grains and no difference could be found in Cube fraction evolution for thicknesses equal or higher than 50 pct. Looking at the experimental Cube volume fractions and at the crystal plasticity simulations results, there is no difference in the amount of Cube remaining after cold  rolling that could explain why, after solution annealing, there is a lower Cube volume fraction in samples from hot bands B and C.
The prediction of texture from crystal plasticity simulations was compared with the texture obtained experimentally from cold rolling. The ODFs from experiments were obtained from the full thickness RD-ND EBSD maps, and ODFs were calculated with harmonic series expansion with a half-width of 5 deg. Figure 11 shows the comparison between simulations and cold rolling for specimens from hot band A. The Cube intensity is low in both experiments and simulations, while the intensity is at the maximum for the b fiber components (Brass (Br), Copper (Cu), and S). The same result occurs for microstructures cold rolled from hot bands B and C where b fiber texture components have the highest intensity (Figures 12(a) and (c)). In simulations, there is a tendency to overestimate Goss and underestimate Brass when 2D representative volume elements (RVEs) are used to simulate texture evolution. For a 2D RVE, the component of the deformation gradient F RD-TD is zero globally and also   locally. This tends to underestimate Brass in 2D RVEs, [3] while for a 3D RVE, the F RD-TD is globally zero (as required by a plane strain compression load), but locally can be non-zero. Having F RD-TD both globally and locally equal to zero brings to an overestimation of Goss, [3] which is known to be a stable orientation up to large strains. [28,50] The intensity of ODFs is higher for crystal plasticity simulations than for experiments because of the limited number of orientations available from the simulations. Moreover, the experimental ODFs in Figures 11 and 12 were plot considering the whole RD-ND section, while crystal plasticity simulation is limited to a small area.

E. Cube Neighbors, Dislocation Density, and Particle Size Analyses
From the previous results, it is not a difference in Cube volume fraction remaining in the microstructure at the end of cold rolling that can justify the different volume fractions of Cube in the recrystallized materials.
According to the literature, [31] two theories explain the origin of high Cube fraction at the end of recrystallization: the oriented nucleation and the oriented growth which can be both originated by high mobility grain boundaries where the relation between Cube and the neighbor is 40 deg h111i: This orientation relation would be found between 20 deg ND rotated Cube ( 001 f gh310i) and Copper, [51] while others mention that this orientation relation can be found between the Cube and the S texture component. [24] For aluminum, Hjelen et al. [51] have shown that Cube transition bands promoting Cube recrystallization can be formed inside or at interfaces to Copper orientation. Here, several high-resolution EBSD maps with 80 nm step size were obtained around Cube grains from all three cold-rolled microstructures. The cold-rolled specimen from hot band A (i.e., without heat treatment between hot and cold rolling) has Cube grains which are very often surrounded by S and Copper orientations. EBSD maps obtained from Cube grains in cold-rolled samples from hot bands B and C (Figures 14 and 15 respectively) show Cube grains not directly in contact with S or Copper but instead surrounded by random orientations. Zones extracted from these specimens have also several large precipitates (Al-Fe-Si precipitates up to 10 lm in length) which are close to Cube grains and seem to be fragmented during cold rolling. As discussed by Lu¨cke and Engler, [8] large non-shearable particles (main size > 1 lm) weaken the overall texture creating zones with random orientations due to PSN causing a decrease of Cube fraction in recrystallized microstructures. According to Dons and Nes, [47] AlFeSi precipitates promote PSN and are responsible for texture randomization. The amounts of iron and silicon in the AA 6016 used for the present experiments justify the presence of AlFeSi particles. [47] Even though a pheno-power law was used for the simulations, the shear stresses can be translated into dislocation density using the Taylor equation [s t ð Þ ¼ s 0 þ aGb ffiffiffi q p , with s 0 as the initial shear strength, G the shear modulus of aluminum, b the Burger's vector, a a factor equal to 0.3 and q as dislocation density] and assuming an initial density of the hot bands of the order of 10 8 m À2 . In Figure 16, the dislocation density from crystal plasticity simulations was calculated and compared for all cold-rolled microstructures from hot band A, B, and C. The dislocation density distribution suggests that microstructures from hot band A have slightly higher  dislocation density than microstructures from hot bands B and C, but the difference is not significant and this confirms that at thickness reductions higher than 80 pct, the initial grain size of the hot band does not play a major role in texture and dislocation density evolution. This was confirmed by experiments on polycrystalline Copper where the stored deformation energy of fine and coarse microstructures was compared showing that after 70 pct of reduction, the values of stored energies are the same. [30] Another parameter analyzed to understand if a microstructure is prone to recrystallization is the kernel average misorientation (KAM). The KAM was extracted both from RD-ND sections from crystal plasticity simulations (Figures 17(a) and (b)) and from EBSD maps of the RD-ND section of the whole thickness of the specimen (Figures 17(c) and (d)). The KAM was calculated considering all points in a kernel up to a 3rd order neighbor in scans with the same step size. As shown in Figure 17, both in simulations and experimental EBSD maps, the KAM distribution and maximum peak is at larger angles for the microstructures obtained from cold rolling of hot band A. Thus, cold rolling from hot bands B and C will produce cold-rolled microstructures with lower KAM compared to hot bands without annealing. The lower KAM for specimens from hot bands B and C is another factor decreasing the probability of Cube recrystallization.
The distribution of precipitates dimension was analyzed in hot bands and in the corresponding cold-rolled specimens and their size distribution is plotted as function of the aspect ratio and of the ellipse angle ( Figure 18). The majority of the particles are in the shape of rods and the ellipse angle is 0 deg if the precipitates are aligned with the rolling direction and it increases as they rotate away from this direction. Precipitates pertaining to magnesium and silicon as well as Fe/Mn-containing dispersoids are present; however, their size is too small (< 0.25 lm) and they are not considered in the particle distribution analysis. Only precipitates with an area > 0.22 lm 2 and a major length > 0.6 lm together with a bright appearance in the backscatter electron (BSE) images (originating mainly from Fe) are considered (see ''Appendix'', Figure A-IV). The average particle size in the hot bands is not significantly altered due to annealing (see the averaged sized particle displayed in purple in Figure 18). After cold rolling, the fragmentation is more pronounced in CR-A, while hot bands B and C maintain the same average particle size after cold rolling, based on the mean values of area and major length that are both about 50 pct smaller in CR-A than CR-B and CR-C. The finer and elongated grain size and thus the higher fraction of grain boundaries present in hot band A can be responsible for a more pronounced fragmentation of precipitates due to an increased probability of interaction between dislocations pile up at grain boundaries and precipitates. Moreover, the microstructure of hot band A was not annealed, thus its dislocation density should be higher than the dislocation density of hot band B and C before cold rolling. These two factors could increase the probability of dislocations interacting with precipitates causing more pronounced fragmentation in hot band A.
This precipitates distribution seems to confirm that PSN plays a more important role to texture evolution and randomization during recrystallization of rolled specimens from hot bands B and C than recrystallization of cold-rolled samples stemming from hot band A. Figure 19 shows the 0, 45, and 65 deg sections of the ODF for hot bands A, B, and C and of samples cold rolled and SA. The texture of hot band A has a significant intensity of Cube, but also a fully developed b fiber texture with Brass, Copper, and S components, while hot bands B and C have high intensity for exact Cube and rotated Cube orientations, but the b fiber components are negligible (Figure 19(a)). The intensity of Cube texture after cold rolling and solution annealing is diminished for microstructures cold rolled from hot bands B and C while it increases for the microstructure from hot band A. The maximum ODF intensity for the specimen from hot band A is 28.0 and it is found at the exact Cube orientation, while for both specimens from hot bands B and C, the maximum intensity is lower than 7.0 and it is not found in correspondence to the exact Cube orientation. Moreover, the specimen from hot band A, after solution annealing, has Cube with RD scatter or exact Cube orientation, while samples from hot bands B and C have little Cube with exact orientation and show a scatter of Cube along ND direction. The influence of intermediate annealing during cold rolling was investigated by Engler [52] by applying a heat treatment at 56 pct or at 80 pct thickness reduction on alloy AA 5182. The heat treatment was done at a temperature of 400°C for 4 min in a salt bath. Annealing after a certain amount of cold rolling was effective in reducing the Cube texture [52] ; however, the introduction of the heat treatment before the starting of the cold rolling reduces the Cube intensity even further. Intermediate annealing after more than 50 pct of cold rolling could reduce Cube intensity only of a factor 2, [52] while in the results presented in Figure 19, the maximum ODF intensity has been reduced of a factor 4 and its maximum is not found at the exact Cube orientation. In the work of Engler, [52] the introduction of the intermediate annealing after more than 50 pct of cold rolling should have diminished the rolling texture but also the dislocation density and both factors should have a positive effect in limiting the Cube intensity during the final recrystallization. However, it cannot be excluded that during intermediate annealing after more than 50 pct of cold rolling, the Cube grain size was increased causing a higher fraction of Cube left in the microstructure during the successive 30 pct of thickness reduction. [52] Moreover, PSN in AA 5182 was probably not as strong as in AA 6016.
The work presented here proves that the texture of the initial hot band has an influence on the texture found in the aluminum sheet after cold rolling and annealing. Previously, Dons and Nes [47] tried to investigate the influence of different initial hot band textures on the microstructure and texture obtained after cold rolling and recrystallization in Al-Fe-Si alloys. They started  with a recrystallized hot band having strong Cube texture at the sheet center and random texture at the sheet surface, but after cold rolling and recrystallization, the Cube fraction was the same both at the center and at the surface. [47] It is true however, that their hot band was annealed at a temperature of 285°C for 24 hours, while the recrystallization treatment after cold rolling was done at 400°C for 3 min, so all temperatures are lower than the annealing temperatures applied in the present work. At low recrystallization temperatures, Cube nucleation is favored over other mechanisms (such as strain-induced boundary migration and PSN) and over nucleation of other orientations because of the fast Cube recovery. [44,53]

IV. CONCLUSIONS
The introduction of an intermediate heat treatment at the end of hot rolling and before cold rolling on AA 6016 results in a viable solution to reduce the volume fraction and intensity of Cube texture after solution annealing at the end of cold rolling. The intermediate annealing between hot and cold rolling changes the grain size in the hot band; however, this does not affect the Cube texture fraction during cold rolling for thickness reductions equal or higher than 80 pct. The dislocation density is expected to be similar in all cold-rolled specimens as a consequence of the high thickness reduction. Although the Cube fraction is similar in all cold-rolled samples, strong differences are found in Cube fraction after solution annealing because the initial grain size of the hot bands affects shape, thickness, and distribution of Cube in the cold-rolled microstructure.
For samples cold rolled from annealed hot bands, the lower intensity and fraction of Cube after solution annealing is attributed to the decrease of intensity of b fiber orientations. The removal or decrease of b fiber components lowers the probability of finding the remaining Cube grains surrounded by S or Copper. Moreover, Cube bands and their growth facilitating interfaces are further apart in the microstructures obtained from annealed hot bands. In addition, PSN seems to be favored in cold-rolled specimens from AA 6016 annealed hot bands and it contributes to texture randomization.

FUNDING
Open Access funding enabled and organized by Projekt DEAL.

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://creat ivecommons.org/licenses/by/4.0/.

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://creat ivecommons.org/licenses/by/4.0/.