Mobility of Radiogenic Helium in Amphibole

Recently experiments on He extraction from an amphibole by the incremental heating unexpectedly revealed that the He release pattern depends on the heating rate. During slow heating (~4 K·min−1) of the amphibole grains, one smooth peak of the He flux from the mineral was observed; in contrast, during fast heating (~40 K·min−1) an additional sharp peak appeared at a temperature about 750 °C. In order to explain these observations, we developed a model of He diffusion from the amphibole, which allowed the calculated He fluxes from the mineral to be reconciled with those observed. From the modelling we derived: (i) the helium diffusion domain size distribution, and evolution of the distribution in the course of incremental heating; (ii) occurrence of the tensile stresses, operating under enhanced temperatures above 700 °C. The stresses are different in sites with the different local thermal expansion of the crystalline lattice and they increase the He diffusion flux. The model can be applied to other minerals (materials).

In this contribution, we continue investigations of the He sites and mobility in amphiboles. 4 He and 3 He concentrations in amphiboles generally exceed those in the host whole rock samples and in some rocks provide the balance of both He isotopes. Since the early papers of Gerling et al. [22,23] a specific role of amphiboles in the noble gas balance of terrestrial rocks was considered in several contributions [9,[24][25][26][27][28][29][30].
Recently, two new features of the He isotope migration in amphiboles were reported [1].
(1) Even though the light isotope, 3 He, is generally considered as a more movable one [31,32], values of the 4 He retention parameter (RET, which is the He amount retained in a mineral/total amount of radiogenic He produced) in some samples are well below those for 3 He, RET( 3 He) > RET( 4 He). To study He isotope mobility in such samples, the high-resolution step-wise-heating experiments (with a 50 • C heating step) were performed and the diffusion process was applied to model the release fluxes. The modeling showed that the diffusion parameters of 3 He and 4 He atoms, survived in the mineral, can't be applied to those that have been lost [33]. 2 of 13 This result was explained by the different sites of the parent isotopes in the mineral: 6 Li for 3 He and 232 Th, 235 U, 238 U for 4 He. Li atoms are supposed to be uniformly distributed within the amphibole crystalline lattice as components of the structural formula or substituting similar size ions of Mg+, Fe++ and Al+. Therefore, the radiation tracks, yielded by the 6 Li(n, α) 3 H reaction, have an enhanced probability to remain within the amphibole grain and to preserve 3 He atoms, produced by 3 H β-decay. U and Th are the highly incompatible elements; they are often concentrated near the mineral boundaries, in the micro mineral inclusions or within the thin film mineralization. Thus, α-tracks, produced by U and Th decay, often cross the grain surface or fissures. The residence times of He atoms within an open track is as short as ≈1 year (at normal temperature) and 4 He atoms easily migrate from such tracks. Hence, the different sites of the parent isotopes stimulate the preferential release of 4 He atoms (compared to 3 He ones), migrating from a mineral by the non-diffusion mechanisms [33].
(2) Another new feature of He migration from some amphiboles is the dependence of the He release pattern on the heating rate during incremental linear heating [1]. If a slow heating rate, e.g., 6 • C·min −1 , was applied to the amphibole grains (separated from the alkaline granites, the Ponoi massif, the Kola Peninsula), He was released as one smooth peak within a temperature interval from 300 to 900 • C. However, in case of the enhanced rate, e.g., 40 • C·min −1 , a sharp second peak appeared with a maximum at ≈770 • C. These observations were reproduced in two different laboratories (see Figure 4 in [1]). Further studies showed that: (i) such a peculiarity of the release pattern is seen in amphibole samples from several other sites; however, this is not an inherent feature of the amphibole group minerals; (ii) generally, the amphiboles contain micro inclusions of other minerals; (iii) the incremental heating of the crushed samples (a powder particle size ≤ 50 µm) results in the one peak He release pattern, irrespective of the heating rate.
The objectives of this study are to present new data on He release from the amphibole grains and powder by fast and slow incremental heating, and to explain the observed release patterns by an adequate diffusion model that takes into account the diffusion domain size distribution and the processes governing He diffusion in the course of heating of the mineral.

Experimental Procedure
He extraction from a mineral by its heating with a linear incessant increase in the temperature (defined here as incremental heating) was performed in a stainless-steel crucible heated by a resistance furnace up to~1000 • C and controlled by a digital temperature controller with Pt-Pt-Rh thermocouple.
Helium was extracted by the incremental heating from the amphibole grains (which sizes were from 0.25 to 0.63 mm) and powder (from~5 × 10 −3 to~5 × 10 −1 mm).
Gases, extracted from a mineral, were admitted into the inlet-purification system. Chemically active gases and heavy noble gases were removed, using cold traps and Ti-Zr getters. He and Ne were admitted into the mass spectrometer, and the He flux from a mineral was measured. The magnetic sector mass spectrometer (MI 1201) was used for measurements of (i) the flux of He, released from a sample during the incremental heating experiments (in the dynamic mode), and (ii) noble gas isotope compositions and concentrations (in the static mode). The mass resolution of~1000 allowed a complete separation of 3 He + from hydrogen ions H 3 + and HD + . The sensitivity for He was 3 × 10 −5 A torr −1 . Before each run, the mass spectrometer was calibrated by the laboratory standard with the 4 He/ 3 He and 4 He/ 20 Ne ratios being 6.29 × 10 5 and 47, respectively. The 4 He analytical blank was within 4.5 × 10 −14 mol.
The U-Pb age of the granites, obtained from the zircon cores, is 2666 ± 10 Ma. The late Paleoproterozoic stage of the Svekofennian tectono-magmatic activity caused the amphibole facies metamorphism of these rocks. The temperature and pressure of this event approached 600 • C and 6 kbar, respectively. New mineral paragenesis was formed (zircon, titanite, monazite, etc.) and some isotopic chronometers (e.g., K-Ar, Rb-Sr) were reset. U-Th-Pb dating of the zircon shells gives 1802 ± 22 Ma, accepted hereafter as the age of metamorphism [34].

Experimental Results
The He release patterns, measured in the course of the fast and slow incremental heating of the powder and grains of amphibole 23/90, and the heating rates are presented in Figure 1 and The He release patterns reveal three principal features, modelled and discussed below.

1.
In case of the fast heating rate, He is released from the powder and grains under higher temperatures compared with the slow rate.

2.
He is released as a smooth peak within a wider temperature interval compared with that expected from He diffusion from the single size diffusion domain.

3.
During fast isothermal heating of the amphibole grains, a sharp second peak appeared with a maximum at ≈770 • C; the second peak has never been observed if He is extracted from the powder.

He Migration in a Diffusion Domain: The Model
Migration of He atoms in a mineral grain is considered to be governed by a diffusion mechanism operating within a diffusion domain, whose size and shape are unknown. In order to simplify the model, we discuss here three symmetric cases introduced by the shape coefficient γ, which is equal to 0 for a plane, to 1 for a cylinder and to 2 for a sphere. Then, the He concentration C in each site of the domain at time t is described by Equation (1) where ξ is the coordinate, D is the diffusion coefficient, which is assumed to be dependent (within a given mineral) only on the temperature, D ≡ D(T). He is considered to escape readily from the surface of the domain; therefore, the boundary condition is He is considered to be homogeneously distributed within a domain at time t = 0, when the extraction experiment starts, therefore the initial condition is: In expressions (2,3), r is the size of the diffusion domain, i.e., the half-width for the plane, and the radius for the cylinder or sphere. Integrating Equation (1) by the domain volume gives: The left side of (4) expresses the rate of change of the He amount inside the diffusion domain (within a constant factor), and the right side represents the He flux through its surface. Let us introduce the average relative concentration of He within the domain, and estimate the He mass flux through the surface of the domain (see the right part of Equation (4)), replacing the derivative by the finite difference: Substituting the average concentration and the He mass flux from Equations (5) and (6) into (4) gives: Now, we introduce the relaxation time τ, which is the characteristic time of the diffusion, according to the following expression: Substitution of Equation (8) into (7) yields: Equation (9) shows that during the relaxation time τ the amount of He in a diffusion domain decreases by a factor of "e" times (e ∼ = 2.71828). As it follows from Equation (8), the relaxation time τ depends on: the size and shape of the diffusion domain; the diffusion coefficient and therefore on the temperature. Equation (9) is the ordinary differential equation that determines the average He concentration in a diffusion domain at any given time; at t = 0 the concentration C is (in accordance with expression (3)): and the value at any subsequent time is determined by the solution (11) of the problem (9) with the condition (10): Hereafter, we discuss the dimensionless relative He amount in a diffusion domain, i.e., the fraction of He retained at time t, or the fraction lost (see Equation (19)). This comment is also valid for: (i) the case when diffusion takes place in the domains of different size (Equation (12)); (ii) the related parameters, e.g., for the He flux, whose dimension is time −1 . Equation (11) is the simplest version of the model, envisaging an equal size for all diffusion domains, and the respective calculations are not able to reproduce the observed decreasing of the average He amount (expressed as the He flux from the mineral) in the course of incremental heating (see Figure 1 below and related text).

The Model with Variable Sizes of the Diffusion Domains
Now, we consider the multi-size distribution of the diffusion domains and introduce the time-dependent function κ r , representing the density of distribution of the He amount over the size r (radius) of the domains (hereafter, it is termed as the helium domain size distribution function, DSD, m −1 ). Then, the He amount in domains with the radius varying within an interval (r, r + dr) is dC = κ r (t, r)dr, and the total He amount at the time t can be obtained by integration of all domains: As the diffusion processes are independent in the domains with different sizes, then, by analogy with Equations (7) and (9), the following equation is valid: where τ(r) is the relaxation time defined by Equation (8).
Equation (13) is the ordinary differential equation that determines κ r values of the DSD at each time moment. Integrating (13) requires the initial condition: where κ r0 is the DSD at the initial time t = 0. The following normalization takes place: Integrating (13) yields: Equation (16) makes it possible to calculate relationships between the time and: (i) the DSD function κ r ; (ii) the He amount in the sample, expressed by (12). Below, we use the 3-parameter distribution (17) as the initial DSD, where R m is the representative radius of the domains (see Figure 2); α, β are constants (α, β > 0); K is a coefficient, whose value is derived from the normalization (15). Equation (17) is a modification of the Erlang's law of distribution [35], well-known from the theory of probability. The initial distribution κ k0 (r) (Figure 2, see curves corresponding to the temperature 20 • C) shows that for r = 0 the He amount κ k0 (r) = 0. While the radius r increases, κ k0 (r) is increasing as well up to a maximum value for r = R m α is asymmetric. Therefore, the averaged radius of the diffusion domain R at the initial time t = 0, expressed by Equation (18) differs from the representative radius R m . The He amount in a domain is decreasing with time, as it follows from Equation (11), and it falls to zero when the time is long enough. Correspondingly, the function g(t), is the He fraction that has released from a sample at time t; g(t) is increasing with time and tends to 1. The time derivative of g(t) is the He flux from a sample in the course of its incremental heating (Section 3.1): In our experiments, the heating rate is • T ≡ d t T = const and heating temperature varies from an initial room value T 0 ∼ = 20 • C to T max ∼ = 1000 • C. Because of the linear connection between the time and the temperature, it is more convenient to use temperature as the independent variable, i.e., to operate by g(T) and d T g(T) instead of g(t) and • g(t); notice that the following relationship takes place: By analogy with • g(t) we define the derivative d T g as the rate of He release from a sample (also termed here as the He flux, K −1 ), depending on the relationships between the diffusion coefficient and the temperature expressed by the Arrhenius Equation (22) where D 0 , m 2 ·s −1 , is the frequency factor, E, J·mol −1 , is the activation energy and R A = 8.31 J ·mol −1 ·K −1 is the universal gas constant. Equation (22) can also be expressed as the Boltzmann exponent and we use below this expression, where U, J, is the activation energy per 1 molecule, k = 1.38 × 10 −23 J·K −1 is the Boltzmann constant. According to (22) and (23), parameters in the above equations are related as where N A ≈ 6.03 × 10 23 mol −1 is the Avogadro number. The above equations were applied to model He release from the amphibole powder and grains for the slow and fast heating rates. The following fixed values of the model parameters were used for all simulations: the temperature diapason from an initial temperature T 0 = 20 • C to a final T max = 1000 • C; the diffusion constants D 0 = 2 × 10 −5 m 2 ·s −1 , E = 9.94 × 10 4 J·mol −1 . The latter values are similar to those reported in the original papers [30,36] and review [37].

He Release from the Powdered Samples: Comparison of the Experimental and Calculated Results
In this section, we compare the He extraction patterns, obtained by incremental heating of the amphibole powder (Section 3.1), with those derived from the diffusion model (Section 3.2). According to the measurements, the He flux d T g increases with increasing of the temperature, approaches a maximum at the temperature T m , and then smoothly decreases to zero (Figure 1). Two features of the He release patters are important: (i) T m depends on the heating rate of a sample: a higher rate • T shifts the maximum d T g towards a higher temperature; (ii) the He release patterns are almost symmetrical for the slow and fast heating rates.
The above equations were applied to model He release from the amphibole powder and grains for the slow and fast heating rates. The following fixed values of the model parameters were used for all simulations: the temperature diapason from an initial temperature T0 = 20 °C to a final Tmax = 1000 °C; the diffusion constants D0 = 2 × 10 −5 m 2 ·s −1 , E = 9.94 × 10 4 J·mol −1 . The latter values are similar to those reported in the original papers [30,36] and review [37].

He Release from the Powdered Samples: Comparison of the Experimental and Calculated Results
In this section, we compare the He extraction patterns, obtained by incremental heating of the amphibole powder (Section 3.1), with those derived from the diffusion model (Section 3.2). According to the measurements, the He flux dTg increases with increasing of the temperature, approaches a maximum at the temperature Tm, and then smoothly decreases to zero (Figure 1). Two features of the He release patters are important: (i) Tm depends on the heating rate of a sample: a higher rate • T shifts the maximum dTg towards a higher temperature; (ii) the He release patterns are almost symmetrical for the slow and fast heating rates.  Figure 1 shows that the simplest version of the model (Section 3.2), envisaging a single-size diffusion domain with a model-derived radius R = 0.037 mm, reproduces the temperatures of the maximum fluxes (see the calculated curves Powder Slow Model 1, PSM1, and PFM1 for the slow and fast heating rates, respectively, in Figure 1). However, the calculated peaks are narrow and their maxima are significantly higher than those measured.
A good agreement between the experimental and calculated fluxes was achieved using the model 2 with the multi-size diffusion domains (the calculated curves Powder Slow Model 2, PSM2, and PFM2 for the slow and fast heating rates, respectively, in Figure 1). The initial DSD function, κ k0 (r), used for these calculations (Figure 2, see curves labeled by a temperature 20 • C), was derived for α = 0.5, β = 0.7, R m = 0.006 mm (Equation (17)). The above values of α and β are also valid for the DSD calculated for the amphibole grains, but R m , the representative radius of the domains containing the maximum amount of 4 He at the initial time t = 0, approaches 0.064 and 0.037 mm for the slow and fast incremental heating of the grains, respectively (see Section 3.5 below).
The calculated sizes (radii) of the powder diffusion domains vary from~10 −4 to~10 −1 ; these sizes are in a reasonable agreement with those observed for the powder particles, from~10 −3 to~5 × 10 −1 mm. The average radius of the diffusion domain R = 0.033 mm, calculated via Equation (18) for the initial time t = 0 (or the initial temperature, T 0 = 20 • C), is quite similar to the radius R = 0.037 mm, derived from the single-size domain model. The calculated sizes (radii) of the powder diffusion domains vary from ~10 −4 to ~10 −1 ; these sizes are in a reasonable agreement with those observed for the powder particles, from ~10 −3 to ~5 × 10 −1 mm. The average radius of the diffusion domain

He Release from the Grains: Comparison of Experimental and Calculated Results
As in the case with the powdered samples, in this section, we compare the measured He release patterns from the amphibole grains with those calculated for the single-size and multi-size domain diffusion models.  Figure 2 also presents the model-derived evolution of the DSD function κ r (r) in the course of slow and fast heating. As far as the heating proceeds, the distribution curves shift towards the larger domains; the radius corresponding to the maximum value of κ r max is increasing, whereas values of κ r max and areas below the curves κ r (r) are decreasing. A faster heating rate causes increase the temperature of He release from the domains of equal size. For example, when the heating temperature approaches 420 • C, He is still available only in domains with R > 8 × 10 −3 and > 2 × 10 −3 mm for the slow and fast heating rates, respectively ( Figure 2).

He Release from the Grains: Comparison of Experimental and Calculated Results
As in the case with the powdered samples, in this section, we compare the measured He release patterns from the amphibole grains with those calculated for the single-size and multi-size domain diffusion models.
First of all, it should be emphasized that the measured He release fluxes from the powder (Figure 1) and the grains (Figure 3) are quite different, especially for the fast heating rate. In the case of slow heating the temperature interval of He release is wider (approaches ∼ = 600 • C) and is shifted to a higher temperature ( ∼ = 620 • C, Figure 3) compared with the powder release pattern (PSObs in Figure 1). In the case of fast heating (curve GFObs, Figure 3), the He extraction pattern has two peaks: the lower temperature smooth peak with a maximum value at ∼ = 670 • C and the higher temperature sharp peak at ∼ = 770 • C. This difference between the He extraction patterns from a mineral in the course of the slow and fast heating rates was for the first time reported in [1].
To model He fluxes from the grains, we applied: the diffusion model with the temperature range and the diffusion constants used to model He fluxes from the powder (Section 3.3); parameters for the DSD function, presented in Section 3.4.
Comparison of the observed release fluxes (curves GSObs and GFObs) with those calculated for model 1 with the single-size diffusion domains (see curves Grain Slow Model 1, GSM1, and GFM1 for the fast and slow heating rates in Figure 3) shows that we approximately reproduced temperatures of the maximal fluxes, as it was done for the powdered samples ( Figure 1). However, the differences between the widths of the experimental and calculated He mass fluxes are larger, compared to the results for the powdered samples.
For a more complex multi-size-domain model 2, a reasonable agreement between the experimental and calculated fluxes is observed for the slow heating rate (curves GSObs and GSM2 in Figure 3). approximately reproduced temperatures of the maximal fluxes, as it was done for the powdered samples (Figure 1). However, the differences between the widths of the experimental and calculated He mass fluxes are larger, compared to the results for the powdered samples.
For a more complex multi-size-domain model 2, a reasonable agreement between the experimental and calculated fluxes is observed for the slow heating rate (curves GSObs and GSM2 in Figure 3). In the case of the fast heating the experiment-and model-derived release fluxes (curves GFObs and GFM2) are similar up to a temperature of ~700 °C (Figure 3). At the higher temperatures there is a significant mismatch of the fluxes: the measured flux yields a sharp peak in contrast to the calculated flux, which is smoothly decreasing. Such feature of the He release flux indicates a sharp decrease in the He relaxation time τ at temperatures ≥ 700 °C; possible causes of this phenomenon are considered below.

Diffusion and the Tensile Stresses.
Equation (23)  In the case of the fast heating the experiment-and model-derived release fluxes (curves GFObs and GFM2) are similar up to a temperature of~700 • C (Figure 3). At the higher temperatures there is a significant mismatch of the fluxes: the measured flux yields a sharp peak in contrast to the calculated flux, which is smoothly decreasing. Such feature of the He release flux indicates a sharp decrease in the He relaxation time τ at temperatures ≥ 700 • C; possible causes of this phenomenon are considered below.

Diffusion and the Tensile Stresses
Equation (23) determines the diffusion coefficient, D = D 0 e − U kT , via the activation energy U characterizing the atomic bonds of the crystalline lattice. A helium atom starts moving through the lattice only if its energy exceeds the activation energy. However, the moving atom shifts neighboring atoms of the lattice and, thus, reduces locally this potential barrier [38]. Additionally, the potential barrier is decreased in case a crystalline lattice is stretched by the tensile stresses, increasing the inter-atomic intervals. Therefore, in accordance with the thermal-fluctuation theory of destruction [39], the energy U in Equation (23) should be replaced by the difference U − Vσ: where σ is the tensile stress; V is the structural parameter, i.e., the efficient volume activated by the moving atom, whose value exceeds the atomic size by about 2 orders of magnitude. Equation (24) shows that occurrence of the tensile stresses increases the diffusion coefficient: the higher the stresses, the faster the diffusion.

Discussion
Now, we discuss why the tensile stresses appear in a mineral during its heating and estimate their values. It is known that atoms of the crystalline lattice perform periodic vibrations, called the thermal vibrations. Considering an atom as a simple harmonic oscillator, we obtain the following expression for the amplitude A of its thermal vibrations: where ξ and U a are the length and the energy of the chemical bonds, respectively. As it follows from Equation (25), the amplitude of atomic vibrations A depends on the ratio of the atomic thermal energy k T over the atomic bond energy U a : the larger the ratio k T/U a , the higher the amplitude. An increase in the amplitude causes expansion of the crystalline lattice. From the theory of elasticity, it follows that the relative deformation (elongation) de, caused by a change of the temperature dT, is represented as: where γ is the proportionality factor, termed the coefficient of thermal expansion (CTE) of a medium. It appears reasonable to assume that the relative deformation of the medium is proportional to the relative change of the vibration amplitude; namely, de ∼ dA ξ . Then, differentiating Equation (25) with respect to the temperature and comparing the result with Equation (26), we obtain the following expression for γ: indicating the dependence of CTE on the binding energy of atoms: the higher the strength of the interatomic bonds, the lower the thermal expansion of the medium. In Equation (27) [40]. Correspondingly, the CTE value (Equation (27)) at different sites of the amphibole is different, and this initiates the internal stresses. To estimate the difference, we define γ as the CTE of the amphibole as a whole; then, if the temperature is increased by a dT value, the internal stresses should increase approximately by where G is the amphibole elastic modulus. Equation (28) shows that in the course of heating, dT > 0, the positive tensile stresses σ increases in the sites with the local thermal expansion coefficients γ below the average value γ , and, according to Equation (24), the He diffusion is getting faster through these sites. It should be noted that the "equilibrium" temperature, below which the tensile stresses in amphibole are negligible, should be taken into account for correct using of Equation (28). The available data indicate that the process of magmatic crystallization during cooling leads to a gradual formation of the mineral crystalline structure, which is finally set at a certain temperature and does not change upon further cooling. According to thermo-barometric measurements in granites, the last restructuring of an amphibole occurs at a temperature ≈740 • C [41]. Additionally, rocks of the Ponoi massif underwent the amphibolite facies metamorphism, which occurred under temperature slightly below 700 • C [26]. Therefore, we may assume that, at a temperature ∼ = 700 • C, the state of the crystal lattice was at the equilibrium and there were no internal stresses, which are developed under T > 700 • C.
Using reasonable values of the parameters, responsible for occurrence of the tensile stresses (Equations (24) and (28)), and applying the diffusion model with variable sizes of the diffusion domains (Section 3.3), we were able to reconcile the calculated He fluxes with those observed in the course of slow and fast heating of the amphibole grains (compare curves GSObs and GFObs with those obtained by modeling, Grain Slow Model 3, GSM3, and Grain Fast Model 3, GFM3, Figure 3).
The dependence of the diffusion coefficient on the internal thermal stresses (Equations (24) and (28)) should be considered as a qualitative one due to: (i) uncertainty in the values of the structural parameter V, the thermal expansion coefficients γ and γ , and the elastic modulus G; (ii) a simplified presentation of the effects of the internal stresses on the interatomic bond effective energy.
It should be highlighted that the high-temperature intense He flux (curve GFObs in Figure 3) has not been observed in some other amphiboles [1]. Within the framework of the model 3 developed here, two explanations could be proposed.
(1) The tensile stresses originate due to different CTE of the amphibole crystalline lattice and the quartz (other minerals) micro-inclusions in it. In this case Equation (28), the accepted above values of the related parameters and the results of modeling (seen in Figure 3, curve GFM3) could be valid providing a sufficient density of the micro-inclusions; the density is a rather variable "irregular" parameter. Such inclusions are present indeed in amphiboles, separated from the Ponoi alkaline granites (see Figure 8 in [1]). This could explain occurrence of the high-temperature He flux and its variable contributions seen in experiments with the different sample-weights of the same sample (see Figure 4 in [1]). A low density of the micro-inclusions in an amphibole from another locality cause insufficient tensile stresses and absence of the high-temperature He flux.
(2) As it was mentioned above, the tensile stresses occur under enhances temperatures exceeding the equilibrium "formations" temperature. If He has (mainly) released from a sample below this temperature (e.g., because of a small size of the diffusion domains), occurrence of the stresses and the related deformations of the diffusion domains would not affect the He flux. This explains the smooth release pattern of the He flux from the amphibole powder (curve PFObs).

Conclusions
Our results show that gas extraction by the incremental (linear incessant) heating with different rates provide a profound insight on processes governing noble gas mobility in a mineral. Thus, He extraction from some amphibole varieties revealed dependence of the He release pattern on the heating rate [1]. Here, we presented the He release fluxes from the Ponoi amphibole grains and powder in the course of their heating at the slow and fast rates (see Figures 1 and 3). The measured fluxes were compared with those calculated for three different versions of the diffusion model developed in this study.
All set of the experimental results can be explained by the He diffusion model, which includes two new important features: (i) the multi-size distribution of the diffusion domains, (ii) occurrence of the tensile stresses, operating under high temperatures and accelerating He migration from a mineral.
The model allows the diffusion domains size distribution (DSD) to be quantified. This sheds light on the disturbance of the crystalline lattice (e.g., by microfractures, solid and/or fluid micro-inclusions, radiation defects), revealed by the He fluxes from a mineral.
The tensile stresses stretch the crystalline lattice, increasing the inter-atomic intervals due to the different thermal expansion at the different sites of the crystalline lattice. Thus, they increase the diffusion coefficient and accelerate He release from a sample, which process is seen as the high-temperature sharp peak of the He flux.
Even though the observed and model derived fluxes are almost identical, the model solution should be considered as a qualitative one because the parameters involved are poorly known. Further studies are needed to better understand this phenomenon, including experiments on He release from amphiboles (as well as from some other minerals) of different composition, age and evolution history.
The combination of the experiments and modeling, worked out in this contribution, can be applied to other minerals (materials), which undergo destruction in the course of heating. We also would like to emphasize that the incremental heating, applied at different rates, provides a convenient framework for numerical simulations of noble gas migration and thus helps to identify the noble gas sites in minerals and to understand processes responsible for their mobility.
Supplementary Materials: The following are available online at https://www.mdpi.com/2075-163 X/11/1/27/s1, Table S1: He release patterns obtained by slow and fast heating of the powder and grains, alkaline amphibole, The Ponoi massif of alkaline granites, Kola Peninsula.