Property Improvement of Additively Manufactured Ti64 by Heat Treatment Characterized by In Situ High Temperature EBSD and Neutron Diffraction

Among various off-equilibrium microstructures of additively manufactured Ti-6Al-4V alloy, electron beam powder bed fusion, in which three dimensional metallic objects are fabricated by melting the ingredient powder materials layer by layer on a pre-heated bed, results in a specimen that is nearly free of the preferred orientation of the α-Ti phase as well as a low beta phase fraction of ∼1 wt%. However, when further heat treatment of up to 1050 ∘C was applied to the material in our previous study, a strong texture aligning the hexagonal basal plane of α phase with the build direction and about 6% β phase appeared at room temperature. In this study, to understand the mechanism of this heat treatment, the grain level microstructure of the additively manufactured Ti-6Al-4V was investigated using in situ high temperature EBSD up to 1000 ∘C, which allows the tracking of individual grains during a heat cycle. As a result, we found a random texture originating from the fine grains in the initial material and observed a significant suppression of α phase nucleation in the slow cooling after heating to 950 ∘C within the α and β dual phase regime but close to the the β-transus temperature at ∼980 ∘C, which led to a coarse microstructure. Furthermore, the texture resulting from phase transformation of the additively manufactured Ti-6Al-4V assuming nucleation at the grain boundaries was modeled, using the double Burgers orientation relationship for the first time. The model successfully reproduced the measured texture, suggesting that the texture enhancement of the α phase by the additional heat treatment derives also from the variant selection during the phase transformation and nucleation on grain boundaries.


Introduction
Titanium alloy Ti-6Al-4V (wt.%) offers excellent formability, fatigue and creep strength, originating from the balanced α and β-Ti crystallographic phases and therefore is widely used in the aerospace industry [1,2]. Recently, with the advent of additive manufacturing (AM) technology [3], a lot of research regarding fabrication and characterizations of this light-weighted titanium alloy was reported [4][5][6][7]. Powder bed fusion is one of the major AM processes, in which three-dimensional metallic objects are fabricated by melting the ingredient powder materials layer by layer on a platform of a pre-heated powder bed. Our preliminary bulk characterization, using pulsed neutron diffraction, revealed that an additively manufactured sample processed by electron beam powder bed fusion (EB-PBF) method showed little preferential orientation of the α phase, whereas that of laser powder bed fusion method showed significant preference towards the hexagonal basal plane, becoming oriented along the build direction [8]. Furthermore, the weight percentage of the β phase of the sample processed with EB-PBF method was found to be ∼1 wt.%, significantly lower than conventional hot-rolled samples, resulting in weight fractions of 5-7 wt.% for the β phase [9,10]. Figure 1 shows the phase fractions of α and β phase in Ti-6Al-4V predicted using the Thermo-Calc [11]. Based on the equilibrium calculation of crystallographic phase fraction as a function of temperature, the phase fraction of the α phase is nearly 100% at room temperature but becomes zero at temperatures higher than 940 • C. The cooling rate during the AM process is estimated to be fast enough [12] for the microstructure to be off equilibrium, which typically leads to a remaining high temperature β phase different from the experimental data of an as-built material processed with the EB-PBF method [13]. Since the phase balance is critical to the mechanical properties of the Ti-6Al-4V, the mechanism of the β phase suppression needs to be clarified to control the material properties. With this motivation, our previous neutron diffraction study [13] was expanded to an in situ high temperature environment up to 1050 • C using a heating chamber, where the microstructure was characterized as a function of temperature, including the α to β to α transformation. A peculiar finding was that starting from a preferred orientation of hexagonal basal planes aligned with the build direction with a maximum pole density of ∼2.5 mrd, after the cycling through the phase transformation, the maximum pole density of the (002) poles aligned parallel to the build direction doubled to ∼5 mrd. Furthermore, the β phase fraction retained at room temperature increased from ∼1 wt.% to ∼6 wt.% before and after the heat treatment, respectively. The measured texture failed to be reproduced by a conventional texture modeling without considering variant selection [13]. There are many additional studies [14][15][16][17][18][19] on the influence of heat treatment on microstructure and mechanical properties of the Ti-6Al-4V. Katzarov [20] successfully modeled that a high cooling rate results in a large amount of α phase and a suppressed rate of nucleation of the α phase depending on the diffusion of vanadium in the material.
The texture enhancement was also reported in conventional Ti-6Al-4V samples [10,21]. Several papers have discussed possible variant selection in β to α transformation [22][23][24][25]. Lonardelli et al. suggested that a preferable grain growth of existing α phase grains is also responsible for variant selection. Obasi et al. suggested that only the α variant nucleated in the early stage can survive when the β phase grain is large, leading to a strong α texture [10,21]. Although these various studies indicate that the texture evolution relates to phase transformation, the difficulty in direct observation of the transformation hinders deep understanding of the mechanism [26]. In this manuscript, two new experimental and modeling approaches are applied to further understand the microstructure evolution of the Ti-6Al-4V processed by the EB-PBF method: the grain level analysis using in situ heating EBSD method, and the phase transformation texture modeling taking into account two neighboring grain pairs simultaneously, which satisfy the Burgers orientation relationship (BOR).

Experimental Methods
The Ti-6Al-4V alloy investigated in this work was fabricated to the cylinder shaped parts as shown in Figure 2 from used grade 5 powder by an EB-PBF method, using an Arcam A2X machine (Arcam AB, Mölndal, Sweden) at the Institute of Metals, Technion, Israel. The chemical compositions of the powder used are shown in Table 1. The original powder compositions reported by the vendors are within the criteria which meet the ASTM requirements listed in reference [27]. The used powder and the as-built materials had the same high oxygen values (oxygen content exceeded ASTM F2924 standard requirements [27]). One layer of powder (∼70 µm) was spread on a platform of a pre-heated plate. Each powder layer was pre-heated in order to attain a temperature between 550 • C and 700 • C (∼0.5 T m , and T m is the temperature of the melting point in Kelvin) by using a de-focused electron beam. Secondary pre-heating up to 0.8 T m selectively heated only the designated surface areas that were melted eventually by electron beam. Hence, the temperature of the entire build part during the AM process continuously varied within seconds between about 700 • C at the current printed layer and about 400 • C at the lowest (first) layer (note that this temperatures are estimates, as these temperatures were never measured to the best of our knowledge). By repeating the process, the full height of the product was fabricated. The acceleration voltage of the electron beam remained constant at 60 kV. Vacuum conditions were 10 −5 mbar of the initial vacuum with a needle valve regulating a constant 10 −3 mbar helium environment. The electron beam path and the parameters (e.g., beam speed and current) were guided by an algorithm (proprietary) formulated by Arcam. The purpose of this algorithm was to maintain a constant heat deposition (J/mm 2 ) in each area in a single layer; additionally, all layers should have had a similar heat deposition. In the current work, the geometries were simple, and therefore, equal heat deposition was achievable in a cylindrical part. There is no texture difference as a function of distance from powder bed [13]. On the other hand, a microstructure difference upon location (build placement) was observed [27]. Mechanical properties of the as-built material as determined by tensile testing are shown in Table 2. Even with the large amount of oxygen, the EB-PBF material used in this study showed high elongation, possibly due to a fine microstructure characteristic of the EB-PBF method. More information was reported in our previous manuscripts [13,27].

Neutron Diffraction
One of the rod pillar samples cut from the product was used for a bulk analysis, using HIPPO (high pressure preferred orientation time-of-fight diffractometer) at the pulsed LANSCE (Los Alamos Neutron Science Center) neutron source [28,29]. Diffraction histograms were collected at different temperatures for a sample prepared from the middle part of a rod by exposing the sample to a neutron beam with a circular cross section of diameter of 10 mm for 20 min (Figure 3).The incident beam was parallel to the build direction of the sample. Diffraction histograms were analyzed by means of Rietveld analysis [30], using the MAUD software [31] to evaluate the phase fraction in weight percent, as well as the crystallographic texture following the procedures described in [32]. For the texture analysis, the orientation distribution function (ODF) was described by the E-WIMV algorithm in MAUD [32], using a resolution of 7.5 degrees. Diffraction histograms were obtained as a function of temperature starting at room temperature and heating up to 1050 • C, using a vanadium sample holder, vanadium heating elements and heat shields in the furnace. Vanadium has a negligible coherent neutron cross section and thus, does not contribute reflections to the measured diffraction histograms [13,29]. Vacuum conditions were 10 −6 torr. The sample was held at each temperature for about 60 min. Presumably due to an unstable microstructure during temperature equilibration, the histograms obtained at the early stage of temperature holding resulted in a poor Rietveld analysis. Therefore, only histograms from the last 15 min-5 min each at the three sample rotation angles of 0 • , 67.5 • , and 90 • -were used ( Figure 3b). The detailed procedure can be found in [13].

EBSD
Direct observation of microstructure evolution was carried out by the in situ heating EBSD measurement [33,34] with an FE-SEM (JSM-7001F, JEOL, Japan) operating at 15 kV of the accelerating voltage. A heating stage system (HSEA-1000, TSL Solutions Co., Ltd., Japan) was used, and the heating stage was mounted on a SEM sample stage.A sheet sample was cut from the middle part of an as-built rod sample by using wire electric discharge machining to the size of 5 mm × 7 mm × 1 mm (length × width × thickness). The thickness direction corresponds to the build direction. One surface of the sample was polished by a standard method with a final step, using 0.04 µm colloidal silica particles.
The temperature was measured by thermocouples attached to the surface of the sample as shown in Figure 4a.
Heat patterns used in the experiment are shown in Figure 4b. The heating rate was about 2 • C/s, which is about 8 times faster than the HIPPO experiment. The temperature changed in steps and the microstructure was scanned at each temperature. Once the temperature reached each target value, the EBSD pattern was obtained. Due to the thermal expansion, the sample moves inside the heating stage, especially when the temperature is changing. Therefore, before obtaining EBSD patterns, the sample position was adjusted, using a SEM sample stage, ensuring the observation area to be the same during the entire measurement.
The EBSD datasets were subjected to a clean-up procedure, using neighbor confidence index (CI) correlation with a minimum CI of 0.05, followed by a single iteration of grain dilation with a minimum grain size of 5 pixels for the in situ measurements and grain tolerance angle of 5 • . After the clean-up process, misindexed data were replaced with the neighbor data that had a higher CI, and every data point was recognized as being part of a grain with a minimum grain size of the designated pixel numbers. The fraction of points changed in each data set was less than 2%. The EBSD data analysis was performed, using an OIM analysis software version 5.3.1 (EDAX Inc., Markham, ON, Canada). For the ex situ EBSD measurement before heating, an area of 2000 µm× 2000 µm was scanned with a step size of 1.0 µm, and in the in situ measurement, an area of ∼400 µm × ∼400 µm was scanned with a step size of 0.4 µm. From the EBSD analysis, pole figures were constructed by a series expansion of generalized spherical harmonics, using smoothing parameters with Gaussian half-width of 10 • , and plotted with an equal area projection method.

Model for Transformation Texture Formation
The texture transformation upon α → β → α phase transformation was modeled based on so-called double Burgers orientation relation (DBOR) as proposed in the literature [35]. In DBOR, the variants that satisfy BOR:{0001} α //{011} β , < 1120 > α //< 111 > β [36], in which there are 6 and 12 variants in α → β and β → α transformations, respectively, or relationship near to it (up to about 10 degrees of deviation allowed) preferentially nucleate and grow. Two types of growth modes were considered. One is that DBOR variants nucleated on grain boundaries grow into the parent grain with which the variants hold BOR, and the growth is confined within the parent grain similar to the martensite transformation (see Figure 5a). The second mode is more extended in which the DBOR variants grow not only by the migration of the coherent interfaces, but also by semi-coherent or incoherent interfaces into the grains with which they hold only near-BOR or even no orientation relation (EDBOR) as shown in Figure 5b. The DBOR variants can appear on five and four kinds of special boundaries in α and β parent grain structures, respectively, as listed in Table 3 [35]. Therefore, the probability for a particular variant to be chosen according to DBOR, ρ(g), is related to the density of those special boundaries in the parent phase, where g specifies the orientation of a parent grain as well as the variant to be chosen. The ρ(g) is formulated for the former growth mode as follows:

Phase Transformation Type Rotation between Adjacent Grains Number Ratio Angle
Axis Deviation from <1210> Here, f (g) is ODF of the parent phase,∆g is the crystal rotation due to BOR, g c k and g p i are rotational operators for crystal symmetry for child and parent crystals, N k and N i are the numbers of the operators (12 and 24 for α and β phases), respectively. The first term in (1) is for the density of the special boundaries, and ζ k is the weight factor for each type of special boundaries, which is assumed to be unity in this study. The second term in (1), (2), is the term for the growth mode. The parameter ω is the strength parameter of the variant selection by DBOR, which is discussed below. Obviously when ω is zero, ρ(g) is unity and the variant selection is absent. When ω is unity, ρ c (g) is close to zero and the nucleation and growth are all governed by the proposed DBOR mechanism. For the latter growth mode, EDBOR, the second term in (1), ρ c (g) is replaced by a constant δ, which is determined to maintain total volume of material during transformation [35]. In this study, the EDBOR calculation was used. The values of ζ k were chosen to be unity for α → β for simplicity. Those for β → α were chosen to be 0, 3, 1 and 1 for the I, II, III and IV types of the special boundaries that were determined for Ti in the previous study [35], respectively. The values of ζ k determined for α → β in the previous study were not used since they were thought to be influenced by macro-zones caused by rolling processes [35]. The ρ(g) in Equations (1) and (2) can be readily expanded by spherical harmonics [35], and the method for transformation texture calculation based on the harmonic expansion of ODF [37] can be used for the texture prediction. Therefore, the observed experimental textures measured on HIPPO were expanded by harmonics, using recalculated pole figures by MAUD, and the transformation textures were calculated using the DBOR model. In all calculations, the orthorhombic sample symmetry was applied in which the symmetry axes are along the c-axis direction and the normal direction in the initial pole figures. The harmonic expansion of the observed textures was performed, using a method described in the literature [38] for α phase and the software called "standard ODF" [39] for β phase, which was truncated at the 28th and 22nd order for α and β, respectively. The transformation texture calculation and the expansion of ρ(g) were truncated at the 20th and 28th orders, respectively.

Microstructure Evolution
Selected regions of neutron diffraction histograms of Ti-6Al-4V measured at different temperature are shown in Figure 6a. The β phase reflection, {110}, is not readily visible at 200 • C in the heating process. With increasing temperature, this reflection clearly appears, shifts due to thermal expansion and the repartitioning of alloying elements, and grows in intensity due to an increasing weight fraction of the β phase. Contrary to that, the α phase related reflections decrease with increasing temperature, until they vanish entirely at 1050 • C. Using the Rietveld analysis, we identified a hysteresis behavior of the β phase fraction and texture during the measurement heat cycle [13] as indicated by a small, but significant above the background, intensity of the β phase reflection {110} at 200 • C in the cooling process. Although it is almost zero before heat treatment, which is reasonable in light of thermodynamic equilibrium (Figure 1), the β-phase remains about 6% at room temperature after the heating and cooling cycle.
Texture analysis revealed an obvious enhancement of the hexagonal basal plane of α phase {0001} parallel to the build direction after the heat cycle. Figure 6b shows changes in the maximum {0002} α and {011} β pole density as a function of temperature. The pole density slightly decreases during the heating up to 900 • C at which the α phase fraction is approximately 70% [13]. However, after it is retransformed from the full β phase in the cooling process, the α phase exhibits a jump in the pole density. Besides that, the preferential orientation of {0001} along the build direction remains. This jump cannot be explained by a texture simulation where the variants are randomly selected, suggesting the presence of a variant selection mechanism as is well studied in the other materials [40]. Figure 7a shows inverse pole figure maps and phase maps during heat treatment analyzed by the EBSD method. The initial microstructure includes fine grains and partly shows the Widmanstätten structure. There is no obvious change in the microstructure during heating up to 800 • C, but the β phase increases at 900 • C and higher. Thousands of fine α phase grains observed at room temperature transform to only several coarse β phase grains at 900 • C. Nucleation of the transformation is generally considered to occur near grain boundaries [41]. The high cooling rate of the AM process from the melt to the temperature of the remaining part (∼700 • C) would presumably increase the number of nucleations at the grain boundaries, explaining the observed fine α microstructure. When the temperature reached 1000 • C, the microstructure consisted entirely of β phase. Furthermore, the grains at this stage grew excessively, resulting in only a single grain in the observed area as shown in Figure 7. During the cooling process, the α phase fraction gradually increased with decreasing temperature, but the texture changed drastically after heat treatment (Figure 7b,c), unlike the neutron diffraction data [13]. This is due to the local characterization of the EBSD method, resulting for large-grained materials in poor grain statistics for texture characterization, much different from the bulk and averaging characterization by neutron diffraction, providing much better grain statistics, even for large-grained materials as is the case here. The significant grain growth of the β phase, together with the fact that the observed grain has a free surface where the EBSD is applied, precludes observation of β grain boundaries in the observed area such that the texture memory (of the previously observed α) did not occur in the field of view of the EBSD instrument. Thus, the texture memory could not be observed in the EBSD analysis. It should be noted that the bulk texture of β phase will not change significantly by the normal grain growth and, therefore, the bulk texture of α phase did not change drastically (but slightly changed) as shown in the neutron diffraction data [13].  Figure 8 shows the changes in the phase fraction in volume as a function of temperature as measured by EBSD. A small amount of β phase ∼0.5% was detected before heat treatment, which was also reported in reference [42]. The irreversible response of crystallographic phase appears as shown at 950 • C. After cooling to room temperature, ∼5% of the β phase grains were retained, which is consistent with the neutron diffraction data.
The first EBSD measurement revealed that the prior β grain boundary did not change significantly up to 950 • C. Therefore, to investigate the possible texture memory of the α phase during heat treatment, the second heat pattern with the maximum temperature of 950 • C (Figure 4b) was applied to the material, thus avoiding complete transformation to the β phase. Figure 9a shows the inverse pole figure maps and phase maps obtained in the second heat treatment. There is no indication of movement of β grain boundary at 950 • C.
Due to texture memory following the BOR, the original β grain formed at high temperature in the AM process is assumed to be reconstructed. The microstructure after heat treatment in this second heating experiment below the α/β transus is similar to that before heat treatment but has a slight (but an important) difference as discussed below. About the same amount of β phase ∼5% remains after heat treatment as the material exposed to the first heat treatment crosses the α/β transus, which is equivalent to the value investigated in the previous experiment [13].  Figure 9b shows the {1010} pole figures of the α phase. Texture does not change drastically after heat treatment, suggesting an effect of texture memory. However, the α phase texture at 950 • C shows a slightly higher intensity than before heating, which remains, even after cooling. Figure 10 shows elements distributions before and after heat treatment without crossing the α/β transus analyzed by EPMA. As reported before [43,44], the concentration of elements V and Al is obviously enhanced in some regions after heat treatment in addition to the observation that some regions are depleted of Ti. Regions belonging to β Ti with b.c.c. crystal structure show lower Ti density than α Ti with h.c.p. crystal structure. In regions where V is concentrated and both Al and Ti concentrations are decreased, the β phase is stabilized and hard to be transformed to the α phase during cooling process. The increased β phase fraction after the heat treatment compared to the as-built material is therefore explained by elemental repartitioning. As studied in a forged material [45], exposure to high temperature around the α/β transus followed by slow cooling allows the V to stabilize the β phase and leads to the hysteresis behavior of the phase fraction during heat treatment. In the EB-PBF process, on the other hand, the as-built material is pre-heated to ∼700 • C only for seconds, which is insufficient for the V to diffuse and stabilize the β phase. Our previous data [13] revealed that the lattice parameter of the β starts to deviate from the low temperature behavior at ∼600 • C, indicating that in addition to thermal expansion, lattice parameter changes due to vanadium redistribution occur. Our findings suggest, therefore, that vanadium redisitribution, enabling the stabilization of the β phase, only occurs at temperatures higher than at least 600 • C.
Previous studies reported that the cooling rate influences the resulting microstructure during the β to α phase transformation [12,[45][46][47]. Reference [19] reported that the martensitic phase appears with the cooling rate 410 • C/s and more, while another work indicated that the martensitic phase decomposes into α +β in the EB-PBF process [48], even with an estimate that the cooling rate in the EB-PBF process is much faster than 410 • C/s [12]. Wang et al. [49] showed that the scanning speed of the electron beam during the build affects the microstructure and that the martensitic phase α' exhibits fine microstructure ∼3 µm. Figure 11 shows the build direction (BD) inverse pole figure map of an AM fabricated sample. The microstructure shows the basket weave morphology. The average minor axis size of grain is 1.30 ± 1.56 µm, which is approximately the same as an AM fabricated sample reported in reference [47]. However, the electron beam AM material used in this study exhibits a mixture of large α colony phase and small grains as shown in the upper area of the figure. We assume that the fine microstructure corresponds to the martensitic α' phase.   The cooling rate in the heating experiment in this study is about 4 • C/s (from 750 • C to 300 • C) and, therefore, is much lower than in the EB-PBF process, where the martensitic phase does not occur. A careful inspection of the inverse pole figure map (Figure 9a) revealed that fine α grains found in the as-built sample disappear after the transformation from β phase, leading to a more homogeneous and larger grain size of the α phase as shown in Figure 12a.   Figure 13d that there is little orientation deviation inside each grain of 20 • C, suggesting that the α grains observed at 950 • C are a part of α grains of 20 • C and that the increase in the α phase fraction during cooling derives from the grain growth of existing α phase rather than the nucleation of new α grains.
The maximum {0002} pole figure density in Figure 12 suggests that the texture during the grain growth of the cooling process from 800 • C to 20 • C does not change significantly. The slight increase in the {0002} pole density of the α phase likely derives from grain growth during the heating process or β to α transformation at the higher temperature. To analyze the texture in detail, the microstructure of an AM fabricated sample was divided into two groups: an area with large grains (>5 µm) and the other area with small grains ( 5 µm). The texture of each microstructure was investigated, using OIM analysis. The area fraction of the grains with dimension 5 µm and less, which we consider the martensitic phase, was 9.2%. Figure 14 shows the inverse pole figure maps and corresponding pole figures. The microstructure with large grains shows much higher {0002} pole density along the build direction than that with small grains. Given that the small grains in the AM processed sample are the martensitic phase, they exhibit a random texture due to the large number of α' variants from the β phase. The small grains shrink easily in the process of grain growth to reduce the boundary energy in the system. Therefore, the α' grains are consumed by grain growth in the process of the heating, which leads to texture enhancement as shown in Figure 12b. During the cooling process, new α' grains are not generated, due to the low cooling rate as discussed above and, thus, the strong texture remains during the cooling process.
However, the mechanism of the texture enhancement of α phase transformed from 100% β phase as measured in our neutron diffraction study is still not understood. In fact, the degree of the texture enhancement of α phase after heat treatment is much larger when the microstructure is exposed to full β (Figure 6b) rather than α and β dual phase temperature regions during heat treatment ( Figure 12). Unlike the EBSD measurement with the second heat treatment not crossing the α/β transus above, the nucleation of the α phase after full transformation to β occurs inevitably in the cooling process from full β phase to increase the phase fraction. In what follows, we apply the new transformation model with the double Burgers orientation relationship (DBOR) to address the texture evolution in the α → β → α phase transformation to explain the observed strengthening of the texture. The results of transformation texture prediction by using the observed parent textures and EDBOR are shown as pole figures and cross-sections of ODF in Figures 15 and 16, respectively. Note that the predicted β textures were used to calculate the final α textures. The values of ω in Equation (1) were determined to be 0.15 and 1.5 for α → β and β → α transformations by minimization of root mean square deviation between ODFs of predicted and observed child textures. The agreement between the predicted and observed textures of high temperature β and final α is excellent. This suggests that (1) the phase transformations for not only β → α, but also α → β should be caused by nucleation and growth of child phases (not by growth of residual child phases), (2) the nucleation takes place on grain boundaries of parent phases and the variants are at least partly selected to reduce the interfacial energy, having BOR or near-BOR with multiple parent grains, and (3) whereas the variant selection for α → β is weak, that for β → α is rather strong in the AM material studied here, which is different from the behavior observed in rolled Ti [35]. The smaller value of ω, 0.15 for α → β, suggests that the nucleation should occur without variant selection on general grain boundaries that are not the special boundaries for DBOR. As described in the previous research [35], the probability for the DBOR variants to be available on the plane grain boundaries in texture-less materials is only about 13% for α → β transformation, when angular deviation up to 10 • is allowed for BOR on one side of the grain boundary. Although textures can alter the probability, the observed texture of α is relatively weak and, more importantly, the main texture component is of a c-axis fiber, in which only the type I special boundary is allowed; the other types of special boundaries have rotation axes near 1210 (see Table 3). Nucleation on triple junctions or grain corners can increase this probability by twice or more. However, there are still many parent grains in which no DBOR variant is available on surrounding boundaries; then, non-DBOR variants nucleate instead. Therefore, the value of ω should be much smaller than unity as predicted. For β → α, the above probability for DBOR variants in textureless materials is larger and 19% on plane boundaries. Furthermore, whereas the observed texture of β is not strong either, it is of a 110 fiber in that the types I, II and IV of special boundaries for DBOR are allowed. The probability thus significantly increases so that an ideal DBOR variant selection should be statistically possible. Hence, the texture memory upon α → β → α phase transformation in the present AM material can be well explained by the DBOR mechanism.

Conclusions
The microstructure evolution of an additively manufactured Ti-6Al-4V during heat treatment was explored by means of in situ EBSD, complementing data from high temperature in situ neutron diffraction. An enhancement of the α phase texture and an increase in the phase fraction of β phase after heat treatment were observed in the EBSD experiment, although the degree of enhancement of the α texture was smaller than the bulk analysis, using neutron diffraction. However, this deviation can be explained by the poor grain statistics of the EBSD technique. Due to the advantage in grain level analysis, the EBSD analysis plus EPMA revealed the following new outcome.
(i) The rapid solidification phenomena, which is a typical AM powder bed outcome, dictates a unique fine microstructure in as-built material. A heating cycle, analyzed utilizing EBSD, reveals that the smaller α grains, which exhibit a weak texture, disappear in the early stage of heating, while the relatively large α grains remain at high temperature. Even after cooling from 950 • C, most of the α grains remain and coarsen by grain growth, while fine α grains as shown in the as-built material are not frequently observed possibly due to a slow cooling rate. Thus, the α phase texture develops slightly after a heat treatment with the maximum temperature of 950 • C, where several percentages of α phase remain.
(ii) Vanadium concentrates in the β phase after heat treatment, which is not observed in the as-built material. The slow cooling in the heat cycle is considered to allow for the vanadium atoms to diffuse and stabilize the β phase, leading to a suppression of α nucleation and a few weight percentage points of the retained β phase at room temperature.
Furthermore, the texture simulation that follows the double Burgers orientation relationship successfully reproduced the experimentally observed texture enhancement of the α phase after transformation from the full β phase at 1050 • C as analyzed from the neutron diffraction experiment. This suggests that the variant selection at the grain boundary as well as the transformation from the 100% β phase is a key to the preferential α texture formation. Funding: The authors are grateful to Pazy Foundation for its partial support of the present study (grant 2020-ID147). This work was performed under the auspices of the National Nuclear Security Administration-Israel Atomic Energy Commission Memorandum of Understanding for science cooperation. And The APC was funded by S.C.V. (LANL).