Next Article in Journal
An Investigation of the Microstructure and Fatigue Behavior of Additively Manufactured AISI 316L Stainless Steel with Regard to the Influence of Heat Treatment
Next Article in Special Issue
Modeling of Phase Equilibria in Ni-H: Bridging the Atomistic with the Continuum Scale
Previous Article in Journal
Surface Analysis of Uncoated and PVD Coated Punch at the Hole-Flanging Process
Previous Article in Special Issue
On the Mn–C Short-Range Ordering in a High-Strength High-Ductility Steel: Small Angle Neutron Scattering and Ab Initio Investigation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Linking Ab Initio Data on Hydrogen and Carbon in Steel to Statistical and Continuum Descriptions

1
Forschungszentrum Jülich GmbH, Institute for Energy and Climate Research, IEK-2, 52428 Jülich, Germany
2
Jülich-Aachen Research Alliance (JARA Energy), RWTH Aachen University, 52056 Aachen, Germany
*
Author to whom correspondence should be addressed.
Metals 2018, 8(4), 219; https://doi.org/10.3390/met8040219
Submission received: 28 February 2018 / Revised: 23 March 2018 / Accepted: 23 March 2018 / Published: 27 March 2018
(This article belongs to the Special Issue First-Principles Approaches to Metals, Alloys, and Metallic Compounds)

Abstract

:
We present a selection of scale transfer approaches from the electronic to the continuum regime for topics relevant to hydrogen embrittlement. With a focus on grain boundary related hydrogen embrittlement, we discuss the scale transfer for the dependence of the carbon solution behavior in steel on elastic effects and the hydrogen solution in austenitic bulk regions depending on Al content. We introduce an approximative scheme to estimate grain boundary energies for varying carbon and hydrogen population. We employ this approach for a discussion of the suppressing influence of Al on the substitution of carbon with hydrogen at grain boundaries, which is an assumed mechanism for grain boundary hydrogen embrittlement. Finally, we discuss the dependence of hydride formation on the grain boundary stiffness.

1. Introduction

One of the challenges for the modeling of materials, like metals and steel, is that a description on a single length scale is not sufficient in many cases. The bridging of scales, however, is not a routine task and requires to accurately link expressions of the different modeling levels. Nowadays, ab initio methods are able to deliver a large variety of material parameters and process routes with high accuracy, but are confined to systems with about ∼100 atoms at most. At the other extreme, continuum methods are able to capture e.g., long ranged elastic effects, but they require material parameters from the lower scales. In between, statistical and thermodynamic perspectives can help to bridge the scales, and to link also equilibrium to out-of-equilibrium regimes. Often, the parameters to be transferred across the scales involve some aspects of coarse graining, as the level of detailed (atomistic) information is lost on the higher scales.
Here, we demonstrate how selected scale transitions between the atomic and continuum scale can be performed, also in a quantitative manner. All these aspects are related to failure mechanisms in high and medium manganese steel, with a special focus on hydrogen. Especially during the last 20 years, while isolated scale bridging investigations on this topic have been conducted, the general need for a comprising scale bridging approach became increasingly clear when light could be shed onto the complex interplay of effects on multiple length and time scales [1,2,3]. Nevertheless, we keep some aspects more generic and discuss e.g., also carbon as interstitial element. The following four features are discussed within this article:
  • We demonstrate that volume dependent solution enthalpies of carbon in iron are fully consistent with the theory of linear elasticity.
  • We show how solution enthalpies of hydrogen can be used for statistical considerations, and how aluminum affects the occupation of octahedral sites in austenitic steel.
  • We introduce a simple mean field scheme to estimate energies of a grain boundary populated with hydrogen and carbon and use that to discuss hydrogen enrichment at carbon saturated grain boundaries in ferrite and the role of aluminum as protective agent.
  • We propose a simple model for grain boundary stiffening or weakening, which affects the probability for hydride formation near grain boundaries in comparison to the bulk.
These applications cover different aspects of importance to hydrogen embrittlement and illustrate how to establish a quantitative link between ab initio simulations, statistical descriptions and continuum theories.

2. Scale Bridging for Strain Dependent Solution Enthalpies in Fe-C

In this section, we demonstrate how data, which is obtained from ab initio or empirical potential calculations, can be transferred to the continuum level. Usually, this transfer is not straightforward, and we use here a particularly simple case to demonstrate this.
The starting point is published data on the solution enthalpy of carbon in bcc iron [4]. This paper contains a rich collection of strain dependent energies for different interatomic potential descriptions, and it is useful to compare the predictions to a continuum perspective, which allows using the computed data also for larger scale modeling.
First, we note that C in bcc iron occupies octahedral positions in equilibrium, and induces a tetragonal distortion. This effect is considered in the careful ab initio based computations of Hristova et al. [4]. They additionally considered the situation where they neglect this anisotropy and assumed a cubic symmetry. It turns out that the results do not change significantly, if the full relaxation is considered. For the present purpose, it has the advantage that we can use the simpler cubic symmetry for the comparison, which leads to more compact analytical expressions on the continuum scale. In the end, we briefly discuss the transfer to the more general situation. In addition, we note that hydrogen in austenite, which is discussed in the other parts of this manuscript from the perspective of a scale bridging, does not exhibit such tetragonal distortions, and therefore the expressions can directly be transferred to this case.
For the calculation of hydrogen solution enthalpies, the authors of [4] use a 3 × 3 × 3 supercell consisting of 54 iron atoms. For the numerical parameters related to the ab initio calculations, we refer to the description in [4]. One of the octahedral positions is occupied by a carbon atom. From the arising energy differences, the solution enthalpy is defined as
H e x c ( V ) = E F e 54 C ( V ) E F e 54 ( V ) E C .
Here, all energies are T = 0 K bare energies from the underlying simulations. E F e 54 C ( V ) is the volume dependent energy of the system with an additional carbon atom and E F e 54 ( V ) the same without the carbon. Notice that both expressions are taken at the same volume. We note that carbon wants to expand the lattice, therefore the first situation has a higher compressive strain as the latter. For a handshaking to continuum methods, it is important to take into account such effects, as we will see below. Finally, E C is a reference potential from a “carbon reservoir”. For the data related to density functional theory calculations (DFT), Hristova et al. use diamond as carbon reservoir [4]. As we are mainly interested here in the elastic rather than the chemical contributions, this reference term turns out to be of minor importance, but in general it is important to consider it e.g., for thermodynamic considerations.
On the continuum level, a natural approach for the total energy of the system is an additive decomposition into chemical and elastic contributions,
E ( ϵ i j , c ) = E e l ( ϵ i j , c ) + E c ( c ) ,
where c is the carbon concentration, and ϵ i j the strain tensor components, as defined below. A central difference between the discrete and continuum energies is that for large scale applications it is often more convenient to normalize the energy per volume rather than per particle. The reason is that, in more complex situations, where inhomogeneous strain and concentration distributions occur in non-equilibrium states, e.g., phase field approaches can be a versatile approach for modeling the time dependent microstructural evolution. For that, a generating energy functional is written as space-integral over an energy density, which is therefore defined with dimension energy per volume.
For the isotropic expansion (under the constraints mentioned above), the elastic energy can be expressed as
E e l ( ϵ i j , c ) = 9 2 B V 0 ( ϵ i j ϵ i j 0 ( c ) ) 2 ,
using the linear theory of elasticity and an isotropic volume expansion, as explained in detail below. Here, V 0 is the undeformed volume in the Lagrangian reference state, B the bulk modulus, and ϵ i j are the linear elastic strain tensor components
ϵ i j = 1 2 u i x j + u j x i ,
expressed through the displacement field components u i . Here, the strain is diagonal with ϵ i j = ϵ δ i j and δ i j = 1 for i = j and 0 otherwise. The eigenstrain ϵ i j 0 ( c ) expresses the stress-free expansion of the lattice through the presence of a finite carbon concentration c. As long as the concentration is not too high, it is appropriate to assume Vegard’s law, hence the eigenstrain scales linearly with c. In addition, if we neglect the tetragonal distortion, the expansion is uniform, hence the eigenstrain takes the tensorial form
ϵ i j 0 = χ c δ i j = ϵ 0 δ i j ,
with the Vegard coefficient χ . For the definition of the concentration, we can use c being the ratio of populated and total number of octahedral positions, hence c is in the range from 0 to 1.
In general, the bulk modulus is also concentration dependent. However, for low concentrations, this effect can typically be neglected. This is supported by the DFT data in [4]. Both with and without one carbon atom the bulk modulus is reported to be B = 189 GPa . Conceptually, a concentration dependence of the elastic constants leads to a higher order correction of the leading terms.
According to [4], the single carbon atom in the 54 atom host matrix leads to a volume increase of 1.7%, hence the strain in each direction is 1 / 3 of this value, ϵ 0 = 0.0047 . The volume of the reference cell is 11.42 Angstrom 3 per iron atom.
The corresponding expression to the solution enthalpy (1) is in continuum representation
H e x c = E ( ϵ , c ) E ( ϵ , 0 ) E C E e l ( ϵ , c ) c c = 0 c + E c ( 0 ) c E C
after Taylor expansion. This gives
H e x c = 9 B V 0 ϵ χ c + E c ( 0 ) c E C = 9 B V 0 ϵ ϵ 0 + E c ( 0 ) c E C .
Under isotropic volume change the strain becomes ϵ = 1 3 ( V / V 0 1 ) , and, according to the above expression, the solution enthalpy should be a linear function of the strain or volume and decay with increasing volume. Interestingly, this behavior is qualitatively satisfied by the DFT data in [4], but not at all by several of the used empirical potentials (see Figure 1). The authors argue that these potentials are therefore unsuitable for describing the strain effects of hydrogen in iron. Here, we confirm this statement from a continuum perspective, as apparently these potentials disagree with predictions from elasticity theory. We can furthermore perform a quantitative comparison using the slope H e x c ( V ) , which in the continuum picture becomes
d E e x c d V = 3 B ϵ 0 .
Inserting the material parameters given above gives the prediction H e x c ( V ) = 0.0201 eV / Angstrom 3 , which is in perfect agreement with the DFT value (identical to the given accuracy, see Table 3 in [4]), and deviates from the empirical potentials, as expected. The continuum result for the strain dependence of the solution enthalpy is depicted in Figure 1. Only for large strains, where linear theory of elasticity is no longer applicable [5], deviations occur with the DFT results.
Finally, we mention that consideration of tetragonality is a rather straightforward extension of the preceding discussion. For that, one uses the more general expression
E e l = 1 2 V C i j k l ( ϵ i j ϵ i j 0 ) ( ϵ k l ϵ k l 0 )
for the elastic energy, and the eigenstrain depends then on the carbon concentrations in the different sublattices. In practise, ϵ i j 0 becomes a superposition of the tetragonal expansions in the different directions.
Overall, this example shows how the transfer from discrete atomistic descriptions to continuum pictures involving elastic effects can be done. From this, one can obtain a closed analytical expression with all material parameters, which can subsequently be used for higher scale modeling.

3. Hydrogen in Fe-Mn-Al Alloys

In this section, we discuss certain thermodynamic aspects of hydrogen in austenitic Fe-Mn-Al alloys. The starting point is the work by von Appen et al. on hydrogen solubility in Fe-Mn alloys [6]. Figure 2 shows a collection of the solution enthalpy of hydrogen. Similarly to the discussion above, it is defined as
H = E Fe 1 x Mn x H E Fe 1 x Mn x μ H ,
where the reference chemical potential μ H is half of the energy of an isolated H 2 molecule. Notice that the calculations were obtained for fully relaxed configurations, hence in particular the pressure vanishes, and consequently the volumes of the two energy terms with and without hydrogen differ. However, from the considerations above, we can conclude that the Legendre transformation between fixed volume and fixed pressure can be performed using theory of elasticity, without computationally expensive additional ab initio simulations. Obviously, the energy essentially depends linearly on the number of manganese neighbors around the H occupied octahedral site (see Figure 2). We note that the same linear energy trend is even more pronounced if the solution enthalpy is plotted versus the Voronoi volume. We can therefore conclude that the leading contribution is due to elastic effects.
We can use this information to obtain some insights into the equilibrium occupation of the different octahedral sites by hydrogen. For an alloy with the composition Fe 1 x Mn x , we assume a random distribution of Fe and Mn atoms, which ignores potential near ordering effects. Therefore, the probability for an octahedral site to have exactly n manganese neighbors is given by
p ( n ) = 6 n ( 1 x ) 6 n x n .
The occupation of a single octahedral site follows a Fermi-Dirac distribution,
c ( E ) = 1 1 + exp ( β [ E ( n ) μ ] ) ,
where 0 < c ( E ) < 1 is the probability of the site to be populated with hydrogen, β = 1 / k T is the inverse temperature, μ the equilibrium chemical potential of hydrogen (not to be confused with μ H above), and E is the solution enthalpy for the given atomic environment. Therefore, the equilibrium hydrogen concentration in the Fe-Mn system becomes
c = n = 0 6 6 n ( 1 x ) 6 n x n 1 1 + exp ( β [ E ( n ) μ ] ) .
As the concentrations are typically low, the additional consideration of more than one hydrogen atom in the supercell should not be necessary. With the linear approximation E ( n ) = a b n (see Figure 2, using a = 66 meV and b = 29 meV ), we can relate the average concentration c to the chemical potential. If we consider (here, for example, for x = 0.2 ) c = 0.5 ppm at room temperature, which is a concentration that may already trigger hydrogen embrittlement, the equilibrium chemical potential is around μ = 351 meV . An octahedron with 8 Mn neighbors is rare within such an alloy, but energetically most favorable. Its average population with hydrogen is around c ( n = 6 ) 60 ppm , which is two orders of magnitude higher than the average concentration c. Similarly, the energetically most unfavorable configuration, where hydrogen is surrounded by Fe only, has an expected H occupation of c ( n = 0 ) 0.06 ppm , which is one order of magnitude lower than the average value. The octahedral positions with a large number of Mn atoms may therefore be considered as hydrogen traps, in particular as their energy is about 130 meV lower than the average site. However, in equilibrium, even these sites have only a minute H occupation, and therefore these sites cannot trap a large amount of hydrogen. This is a general feature of hydrogen with such low concentrations. The chemical potential is very low, and hence all site occupations are in the tail E μ k T of the Fermi-Dirac distribution. Typical trap depths are not sufficient to provide a significant increase of the hydrogen concentration beyond a dilute limit, which would require reaching the regime E μ k T .
Aluminum changes the energy landscape rather strongly, as worked out for specific atomic arrangements in [7,8]. Both author groups use different concentrations of their alloys, which is captured by a 32 atom fcc supercell. Song et al. use Fe 26 Mn 5 Al 1 , Hüter et al. Fe 23 Mn 8 Al 1 , and both groups come to similar conclusions despite the different stoichiometry, which suggests that the situation is more general. We follow the notation in [8], where x-y-z denotes an octahedral environment with x iron, y manganese and z aluminum atoms; this notation, however, cannot distinguish between different atomic arrangements, but we stay with it for convenience. In the following discussion, we focus on the Fe 23 Mn 8 Al supercell. This system allows for distinguishing 6-0-0, 5-0-1, 4-2-0 and 3-2-1 neighbourhood configurations for the hydrogen atom.
The starting point is the trans (high energy) 4-2-0 configuration in Figure 2 (notice that the simple notation 4-2-0 does not reflect whether Al is present in a next-to-nearest neighbor position, therefore this information needs to be provided additionally). In contrast to the Fe-Mn alloys discussed above, the energy is here not only determined by the nearest neighbor configuration. A 4-2-0 neighborhood, where an Al atom occupies a next-to-nearest neighbor position, diminishes the hydrogen solution enthalpy by about 60 meV [7]. It is therefore speculated that this octahedral configuration may act as an H trap; in view of our discussion above, we stress that this effect is however rather weak.
More pronounced, however, is the effect that a direct neighborhood of Al to an H atom is energetically highly unfavorable. Configurations like 3-2-1 or 5-0-1 are about 130–140 meV less favorable than the potential 4-2-0 “trap” site. These octahedra are populated with a concentration that is about two orders of magnitude lower than the average concentration at room temperature; we may therefore effectively say that such sites are blocked.
Compared to the Fe-Mn system, Al therefore leads to strong alternations of the entire energy landscape. Changes of atomic neighborhoods by replacing one Fe atom by Al raises the energy in several cases by around 80 meV. In contrast, only one extreme case in the Fe-Mn system, for a transition between the trans 4-2-0 to the mere 3-3-0 environment, leads to a comparable shift of about −80 meV. In all other cases, the energetic shifts are much milder in the Fe-Mn system. Altogether, we can therefore conclude that the entire solution enthalpy landscape becomes significantly rougher in the Fe-Mn-Al case by the addition of rather small amounts of Al. Already for about 3 at% of Al almost 20% of all octahedral cages are neighboring an aluminum atom and are therefore energetically unfavorable for hydrogen.

4. Ferritic Grain Boundary Hydrogen Embrittlement: Estimates From Ab Initio Datasets

In this section, we discuss a possible phase separation into hydrogen enriched regions and hydrogen depleted regions at ferritic grain boundaries, and the function of aluminum as protective agent against grain boundary related hydrogen embrittlement.
An important contribution to the macroscopic interpretation of hydrogen embrittlement is the phase formation based on available ab initio data. Especially in the case of hydrogen in steel, this approach becomes complex, as at least carbon, hydrogen and iron need to be considered, and the required amount of ab initio data quickly becomes large. Especially for grain boundaries, which require rather expensive electronic scale simulations, the data set is sparse. Therefore, the development of possible approximative schemes can be helpful in this context.
Our first approach employs a mean field approximation for grain boundary energies, based on data published in [4,9,10] for varying carbon and hydrogen populations. The most investigated grain boundary in these publications is Σ 5 ( 310 ) [ 001 ] in bcc iron, and the authors present—amongst other results—segregation and grain boundary energies for carbon and hydrogen at grain boundaries, where the most attractive type of grain boundary site is populated. In the considered supercell, four of these sites exist per grain boundary, i.e., any combination of altogether four hydrogen or carbon atoms per grain boundary corresponds to a fully saturated grain boundary with respect to this site. For details about the determination of the most attractive sites in a grain boundary state of minimal energy, we refer to [4,9,10,11]. The grain boundary energies reported in [9] are defined as
γ G B n X = E G B F e + 2 n X E b u l k F e 2 n μ X 2 A ,
where γ G B n X is the energy of a grain boundary with interfacial area A and n atoms of type X in it. E G B F e + 2 n X is the absolute energy of the used simulation cell, which includes two grain boundaries of the simulated type due to periodicity constraints. By E b u l k F e , the absolute energy of a simulation cell with the same number of iron atoms is defined. Note that all simulations were relaxed using zero stress boundary conditions for the supercell, i.e., the only elastic contributions of the resulting energies are those confined to atomistic length scales around the grain boundary. μ X is the reference chemical potential, for hydrogen that corresponds to H 2 molecules, for carbon to diamond, as discussed also in the preceding sections.
The scheme to estimate grain boundary energies for various hydrogen and carbon populations of the reported Σ 5 ( 310 ) [ 001 ] grain boundaries is a summation of the averaged energetic contributions per atom relative to the fully saturated grain boundary for each pure population, which can be understood as a mean field approach. This can be expressed as
γ G B n H + m C = n 4 γ G B 4 H + m 4 γ G B 4 C + 4 n m 4 γ G B ,
where the required grain boundary energies reported in [9] are γ G B = 1.55 Jm 2 , γ G B 4 H = 1.27 Jm 2 , and γ G B 4 C = 0.58 Jm 2 . Besides the reported grain boundary energies for zero to four sites occupied by carbon or hydrogen, also a co-segregated grain boundary energy is provided, for two hydrogen atoms and one carbon atom. This co-segregated grain boundary energy is an interesting test data point for the approximative scheme, which would predict the grain boundary energy for the reported population as 1/4 of a carbon saturated grain boundary, 2/4 of a hydrogen saturated grain boundary and 1/4 of an empty grain boundary. In Table 1 below, we see the comparison of predicted mean field grain boundary energies and ab initio calculated grain boundary energies for the testable data points in [9].
Apparently, the mean field description works excellently for hydrogen, while the discrepancies for carbon can reach up to 0.09 J/m 2 . This has to be put in context of the variance in reported ab initio results for grain boundary energies that correspond to identical definitions and therefore should give identical numbers. We provide here the grain boundary energies reported for two grain boundaries. The reporting references include [9,10,11,12,13,14,15], and the reported grain boundary energies range from 1.49 Jm 2 to 1.63 Jm 2 ( Σ 5(310)) and from 0.43 Jm 2 to 0.47 Jm 2 ( Σ 3(112)).
While the approximative scheme can be tested for the Σ 5 ( 310 ) [ 001 ] grain boundary in the especially interesting regime of co-segregated grain boundaries, we are currently not aware of grain boundary energies for co-segregated sites for any other grain boundary. Consequently, we continue with a focus on the Σ 5 ( 310 ) [ 001 ] grain boundary.
As the next step, we obtain the energies required for the determination of grain boundary phase separation. From the findings reported in [9], we know that the ground state of the grain boundary is carbon saturated with respect to the most attractive sites. These sites are also most attractive to hydrogen, thus the phase equilibrium in the ternary Fe-C-H system at the grain boundary reduces to the coexistence of regions in the grain boundary with varying hydrogen and carbon fractions, where all attractive sites are occupied by one of those species.
When the mixture of carbon saturated (hydrogen free) phase ( α ) and hydrogen saturated phase ( β ) forms, the averaged Gibbs energy per host atom is given as
g ( x H , T ) = ( 1 v ) g α ( x α , T ) + v g β ( x β , T ) ,
where g α is the averaged Gibbs energy per host atom in the hydrogen-free phase, g β is the averaged Gibbs energy per host atom in the hydride phase, x H , x α , x β are the averaged hydrogen concentration, hydrogen concentration in the carbon-saturated phase and hydrogen concentration in the hydride phase. We assume that we have an ideal statistical distribution of the hydrogen atoms along the available sites, and restrict the description to the configurational entropy contribution. The Gibbs energies of pure iron ( α ), x α = 0 and a pure hydride phase ( β , fully saturated, x β = 1 ) are
g α , β = g α , β 0 ( x α , β ) + g c ( x α , β , T ) ,
where, in the low concentration regime,
g c k T x H ln ( x H / x 0 )
is the dominant contribution. We require equal chemical potentials for the heterogenous mixture of hydrogen saturated and carbon saturated grain boundary phases and the α -phase, so we obtain for low temperatures
g α ( 0 ) + g β ( x β 0 ) x x β g α ( 0 ) x α 0 + k T ln ( x H / x 0 ) ,
where we used v x H / x β 0 , with x α , β 0 denoting the respective phase equilibrium concentration. For the solulibility limit, we get
x H x 0 exp Δ G k T ,
with Δ G = g α ( 0 ) + g β ( x β ) ) / x β 0 + g α / x x = 0 . The derivative of the Gibbs energy at zero concentration of hydrogen, g α / x x = 0 , is obtained as finite difference approximation from the available ab initio data. The expression for the enthalpy difference Δ G becomes more familiar when we rewrite it in terms of classical defect energies E α , β f , and we identify the relevant contributions for Δ G ,
Δ G = E α f E β f , E α f = E G B ; N M t o t ; 3 C : 1 H E G B ; N M t o t ; 4 C μ H r e f + μ C r e f , E β f = E G B ; 1 / x 0 β t o t ; 1 H E G B t o t ; 1 C 1 μ H r e f + 1 μ C r e f , E β f = 1 4 E G B ; N M t o t ; 4 H E G B t o t ; 4 C 4 μ H r e f + 4 μ C r e f .
In the last line, we took advantage of the scaling behaviour of the energy with respect to the particle numbers. By E G B ; N M t o t ; n C : m H , we define the total energy of the grain boundary simulation cell including n carbon atoms and m hydrogen atoms at the attractive sites at each of the two grain boundaries, and N M host atoms. To calculate the differences in energies per atom from the differences in grain boundary energies, Δ E = 3.207 eV m 2 J 1 Δ γ G B + Δ μ r e f holds for the given grain boundary area of A = 51.371 × 10 20 m 2 . The values obtained here for a separation in the considered Σ 5 ( 310 ) [ 001 ] grain boundary into fully carbon saturated and fully hydrogen saturated phases lead to a hydrogen solubility at the interface governed by the enthalpy difference ( Δ G 0.06 eV ) , where the reference chemical potentials correspond to diamond and H 2 .
When we relate the enthalpy difference to the bulk solution energies in the preferred octahedral site for carbon, at 0.74 eV [12], and the preferred tetrahedral site for hydrogen, at 0.25 eV [11], we can estimate that the grain boundary will be co-segregated by hydrogen and carbon without separation into hydrogen-free and carbon-free regions. The reason is a relative shift in the chemical potentials of about 0.49 eV in favor of a partial substitution of carbon with hydrogen at the grain boundary.
From this perspective, it is interesting to consider the influence of aluminum, which is known as protective agent against hydrogen embrittlement. The ab initio data we use to estimate the influence of aluminum has been published in [16]. In this publication, there were no grain boundary or bulk solution energies for hydrogen in Fe-Al reported, but the diffusion barriers of hydrogen in carbide-free bainite have been calculated. Since Fe-(Al)-H was chosen as a representative system, we use the reported results as the upper limit for the repulsive interaction of Al with hydrogen in the direct neighborhood. Depending on the Al fraction in the supercell, the shift in the hydrogen solution energy at the transition point along the tetrahedral-tetrahedral path varies and reaches up to 350 meV.
We choose an intermediate value of 150 meV that reflects two aspects: the transition point shift due to Al alloying—at least in fcc Fe—is around twice as high as the shift of the low energy occupied sites solution energies due the Al alloying [8], and even with Al segregating to ferritic grain boundaries [17], we do not assume a local peaking to 50 atomic percent of Al. However, even with 150 meV shift per hydrogen atom, the preference of hydrogen to segregate to the grain boundaries is suppressed.
Besides this thermodynamic picture, we briefly consider possible kinetic scenarios for the substitution of carbon with hydrogen in grain boundaries. While the carbon atoms act as a strengthening element in grain boundaries, as reported in detail in [4], they also increase the local electronic density. Specifically, as reported also in [9], the carbon atom in the grain boundary shows overlapping bands with its adjacent Fe atoms and a spherically symmetric local increase of the charge density distribution (see Figure 3 in [9]). It is possible that the hydrogen atoms in the metallic matrix, which are partially negatively polarised, lead to a local charge redistribution at these spots of increased electronic density, providing access to a substitution of the carbon by hydrogen. We note that this mechanism will require a more detailed quantum mechanical study of the dynamic bonding behavior at this site.

5. Hydride Precipitation near Grain Boundaries

As described in Section 2, the application of a tensile strain can lower the solution enthalpy of interstitial elements and therefore favor segregation. If the material is homogeneously deformed, also the shift of the solution enthalpy is the same everywhere. For a grand canonical ensemble, the hydrogen concentration raises uniformly. In contrast, for a canonical ensemble with a fixed total amount of hydrogen, the concentration distribution is unaffected. This is different for non-uniform deformation, e.g., localized tensile strains near crack tips or dislocations, where the equilibrium hydrogen concentration can strongly increase and eventually lead to localized hydride formation.
As we have worked out recently, solubility limits can also be affected by the proximity of surfaces and interfaces even in the absence of external stresses [18]. At the solubility limit, a secondary phase (hydride) forms, which typically has an elastic mismatch with the parent phase. This induces long-ranged coherency stresses, which tend to suppress the phase separation. If free surfaces are close to the potentially forming precipitate, the stresses can partially be released, and therefore the elastic energy penalty for two-phase situations is reduced. Consequently, there is a preference for phase separation near surfaces, or, in other words, the solubility limit for hydrogen is reduced in these regions.
We have worked out a scale bridging description of these processes, where all relevant material parameters are derived from density functional theory, and the long-ranged elastic distortions are treated on the continuum level. Both descriptions are linked by a statistical mechanics and phase field perspective for the prediction of stress and surface affected solubility limits. A particular finding is that, in the low concentration and low temperature regime, which is most relevant for hydrogen in metals and steel, the solubility limits between bulk and surface differ by a factor
s = exp γ Δ G e l bulk k T .
This expression is independent of chemical bonding effects and only depends on elastic parameters. It contains the bulk elastic energy
Δ G e l bulk = E 1 ν χ 2 Ω N 0 c β 0 ,
with Young’s modulus E, Poisson ratio ν , Vegard coefficient χ and low temperature hydrogen concentration c β 0 in the precipitate phase. In Equation (22), γ is a dimensionless geometrical factor of order unity, which describes the near surface relaxation (for 0 < γ < 1 ), which can be obtained e.g., by finite element simulations [19]. We note that this prediction is based in particular on a strain dependent energy including hydrogen, analogous to carbon in Fe as analyzed in Section 2. For a detailed discussion, we refer to [18].
We can use these results to discuss the possibility of precipitation, in particular hydride formation, near grain boundaries. For that, we assume that the grain boundary appears on mesoscopic scales as a layer with a locally different mechanical strength. The inclusion is assumed to have a diagonal eigenstrain relative to the matrix phase, and the elastic constants are taken as equal. A sketch of the geometry is shown in Figure 3.
A spherical precipitate is placed at a certain distance from either a free surface or a grain boundary. If the precipitate is neither close to the free surface nor to the grain boundary, the elastic energy reaches the bulk value, which can be calculated analytically under the given conditions (see [18] for details). For different distances of the precipitate to a free surface or the grain boundary and different elastic properties of the grain boundary, the elastic energy of the system is calculated with the finite element method. Close to the free surface, the elastic energy drops with the distance of the precipitate’s center of mass to the surface h (see Figure 4).
When the precipitate approaches the grain boundary, the elastic energy can either increase or decrease, depending on whether the elastic constant E g b in the grain boundary layer is effectively higher or lower than in the bulk. In the limit that E g b / E tends to zero, the grain boundary acts effectively as a free surface, and we recover the same elastic energy relaxation. On the contrary, for E g b E , the grain boundary becomes rigid, and we receive the same situation as for nucleation near a fully confined surface, as discussed in detail in [18].
We note that, due to the computational expense, limited datasets of usable mechanical data for grain boundaries are available. In [9], the dependence of the theoretical strength σ t h of the Σ 5 ( 310 ) [ 001 ] grain boundary on the hydrogen and carbon population is reported. By varying the carbon and hydrogen content of this grain boundary, σ t h varies by 10–50 percent when the empty grain boundary is populated with hydrogen or carbon. In particular, an energetic preference is reported in [9] for a segregation of hydrogen to the carbon saturated grain boundary to replace carbon, and an associated decrease of the theoretical strength by about 20 percent is determined. We note that this value is subject to variations arising due to orientational changes and elastic anisotropy in real samples. We assume that these effects lead effectively to a local alternation of the ratio E g b / E in the spirit of our continuum picture.
The energy modification 1 γ = G e l / G e l bulk can then be inserted into Equation (22) to obtain the increase or decrease of the solubility limit. Near a stiff grain boundary, hydride formation is less likely.
In conjunction with the role of Al as a protective agent against grain boundary embrittlement discussed in the previous section, the findings in this section suggest the possibility of an increased weakness against hydride formation close to grain boundaries when the segregation of Al to grain boundaries prevents the substitution of carbon by hydrogen. However, we mention that a detailed understanding of the role of Al on carbon population of attractive grain boundary sites is presently lacking to give a more complete estimate here. In particular, the possible configurations of Al that range from a disperse distribution to highly localised concentrations require further investigations on the quantum mechanical scale.
We remark that one can interpret the Boltzmann factor in Equation (22) as a kinetic modification of the energy barrier for a hydride nucleation process. Apart from phase separation triggered by nucleation events, for higher concentrations, also barrierless transformations driven by spinodal decomposition are conceivable. The mixing enthalpy (per unit volume) is given by
g m i x ( c ) = α N 0 Ω c ( 1 c )
as a function of the hydrogen concentration c, using the number of metal atoms N 0 per unit cell of volume Ω and the interaction strength parameter α that can be obtained by ab initio simulations. The dominant configurational contribution to the lattice gas entropy is in the low concentration regime
μ c = k T ln c .
According to the seminal work by Cahn and Hilliard [20,21], the chemical spinodal is determined by the condition μ ( c ) = 0 . Asymptotically for T 0 and for the low concentration branch,
c k T 2 α .
In contrast, the coherent spinodal, which is comparable to the bulk elastic solubility limit, is determined by the condition [20]
μ ( c ) + 2 E χ 2 Ω ( 1 ν ) N 0 = 0 .
In the low temperature limit for the low concentration branch, we get
c k T 2 α 2 E χ 2 Ω ( 1 ν ) N 0 1 .
Surface spinodal modes have been discussed [22], and they can appear in between the coherent and the surface spinodal in the vicinity of a free surface. From these expressions, we can conclude that, for hydrogen related problems in steel, the required concentrations for spinodal decomposition are typically too high, and phase separation should therefore be initiated by nucleation, as discussed above.
Within this section, we have focused on the influence of mechanical characteristics of grain boundaries on the hydride formation in the vicinity of a grain boundary. Analogously to the thermodynamic description of the influence of elastic relaxation close to free surfaces in [18], the stiffness of the grain boundary can lead to a local decrease or increase of the hydrogen solubility limit, hence favoring or suppressing hydride precipitation from the host matrix.

6. Discussion

We have presented several interlinked scenarios for the continuum level interpretation of ab initio data on hydrogen in steel in the context of hydrogen embrittlement, also in competition with carbon. The connecting motif is the high complexity of the hydrogen embrittlement process at grain boundaries, which requires understanding of all the presented phenomena in a scale-bridging fashion to develop a quantitative connection to macroscale transport and mechanical phenomena.
The connection of ab initio data on the dependence of the solution enthalpy of carbon on the strain state to a continuum scale description supports the discussion of the influence of grain boundaries on hydride formation (or precipitation in the vicinity of grain boundaries in general), also depending on the mechanical properties of the grain boundary. Here, we recognise that, depending on the hydrogen segregation at the grain boundary, a change of the stiffness by about 20 percent can be expected—and this implies a possible competitive embrittlement mechanism either in the vicinity of the grain boundary or in the grain boundary by hydride formation or enrichment. A more conclusive prediction will require more ab initio data, especially since energies can significantly depend on the type of grain boundary, misorientation, etc.
Furthermore, we have described the statistical connection of the ab initio data on the influence of Al on hydrogen solution energies to its local concentration in steel. This also contributes to the discussion of the protective influence of Al against hydrogen embrittlement at grain boundaries. We have introduced an approximative scheme to estimate the energies of grain boundaries for varying levels of co-segregation by hydrogen and carbon, with an accuracy that is comparable to the variation within published ab initio data. We note that, in a carbon saturated ground state of the considered grain boundary, the hydrogen enrichment is driven by a segregation energy. It gives a relative preference to hydrogen instead of carbon by about 100 meV, and we estimate the repulsive effect of Al on H in the grain boundary to be about 150 meV at sufficiently high Al concentrations. From this, a suppressive effect of Al on the hydrogen accumulation at the grain boundary can be expected. For higher Al concentrations, competing phases may form that can also affect the intergranular fracture behavior.

7. Materials and Methods

All numerical approaches for obtaining the simulation data have been obtained using approaches that are explained in detail in the main cited references. In particular, the Vienna ab initio simulation package (VASP, version 5.4.1, Computational Materials Physics, University of Vienna, Vienna, Austria, 2015) has been used for density functional theory calculations [23,24,25]. Finite element method calculations have been performed with the open source software FreeFEM++ (FreeFEM++, version 3.59) [19].

8. Conclusions

To conclude, we mention that the possibilities to connect available ab initio data sets to continuum descriptions require more data from the ab initio side. Especially for grain boundary hydrogen embrittlement, a comprising study would benefit from data on the dependence of carbon and hydrogen solution energies in at least two different grain boundaries on the presence of aluminum and the strain state. However, these kinds of calculations are computationally expensive, though certain synergy effects are accessible—the theoretical strength can be considered as a byproduct of such simulations. Still, the configurational complexity of the corresponding simulation cells due to the distribution of Al as a substituting element and carbon and hydrogen as interstitials creates a vast set of configurations to be covered. For sure, such data will help to shed light on the complex interplay of different mechanisms, which may contribute to the understanding and prevention of hydrogen embrittlement in the future.

Acknowledgments

This work has been supported by the Collaborative Research Center 761 “Stahl ab initio” of the German Research Foundation. The authors gratefully acknowledge the computing time on the supercomputer JURECA at Forschungszentrum Jülich. The authors thank Tilmann Hickel for many valuable discussions.

Author Contributions

The results in this paper have been worked out together by all authors. The manuscript has been written together by all authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Venezuela, J.; Liu, Q.; Zhang, M.; Zhou, Q.; Andrej, A. A review of hydrogen embrittlement of martensitic advanced high-strength steel. Corros. Rev. 2016, 34, 153–186. [Google Scholar] [CrossRef]
  2. Jemblie, L.; Olden, V.; Akselsen, O.M. A review of cohesive zone modelling as an approach for numerically assessing hydrogen embrittlement of steel structures. Philos. Trans. R. Soc. A 2017, 375, 20160411. [Google Scholar] [CrossRef] [PubMed]
  3. Robertson, I.M.; Sofronis, P.; Nagao, A.; Martin, M.L.; Wang, S.; Gross, D.W.; Nygren, K.E. Hydrogen embrittlement understood. Metall. Mater. Trans. A 2015, 46, 2323–2341. [Google Scholar] [CrossRef]
  4. Hristova, E.; Janisch, R.; Drautz, R.; Hartmaier, A. Solubility of carbon in a-iron under volumetric strain and close to the Σ5(310)[001] grain boundary: Comparison of DFT and empirical potential methods. Comput. Mater. Sci. 2011, 50, 1088–1096. [Google Scholar] [CrossRef]
  5. Hüter, C.; Friák, M.; Weikamp, M.; Neugebauer, J.; Goldenfeld, N.; Svendsen, B.; Spatschek, R. Nonlinear elastic effects in phase field crystal and amplitude equations: Comparison to ab initio simulations of bcc metals and graphene. Phys. Rev. B 2016, 93, 214105. [Google Scholar] [CrossRef]
  6. Von Appen, J.; Dronskowski, R.; Chakrabarty, A.; Hickel, T.; Spatschek, R.; Neugebauer, J. Impact of Mn on the Solution Enthalpy of Hydrogen in Austenitic Fe-Mn Alloys: A First-Principles Study. J. Comput. Chem. 2014, 35, 2239–2244. [Google Scholar] [CrossRef] [PubMed]
  7. Song, E.J.; Bhadeshia, H.K.D.H.; Suh, D.-W. Interaction of aluminium with hydrogen in twinning-induced plasticity steel. Scr. Mater. 2014, 87, 9–12. [Google Scholar] [CrossRef]
  8. Hüter, C.; Dang, S.; Zhang, X.; Glensk, A.; Spatschek, R. Effects of Aluminum on Hydrogen Solubility and Diffusion in Deformed Fe-Mn Alloys Adv. Mater. Sci Eng. 2016. [Google Scholar] [CrossRef]
  9. Tahir, A.M.; Janisch, R.; Hartmaier, A. Hydrogen embrittlement of a carbon segregated Sigma 5 (310) [001] symmetrical tilt grain boundary in α Fe. Mater. Sci. Eng. A 2014, 612, 462–467. [Google Scholar] [CrossRef]
  10. Jiang, D.E.; Carter, E.A. Diffusion of interstitial hydrogen into and through bcc Fe from first principles. Phys. Rev. B 2004, 70, 064102. [Google Scholar] [CrossRef]
  11. Du, Y.A.; Ismer, L.; Rogal, J.; Hickel, T.; Neugebauer, J.; Drautz, R. First-principles study on the interaction of H interstitials with grain boundaries in α and γ-Fe. Phys. Rev. B 2011, 84, 144121. [Google Scholar] [CrossRef]
  12. Wang, J.; Janisch, R.; Madsen, G.K.H.; Drautz, R. First-principles study of carbon segregation in bcc iron symmetrical tilt grain boundaries. Acta Mater. 2016, 115, 259–268. [Google Scholar] [CrossRef]
  13. Cak, M.; Sob, M.; Hafner, J. First-principles study of magnetism at grain boundaries in iron and nickel. Phys. Rev. B 2008, 78, 054418. [Google Scholar] [CrossRef]
  14. Bhattacharya, S.K.; Tanaka, S.; Shiihara, Y.; Kohyama, M. Ab initio perspective of the [110] symmetrical tilt grain boundaries in bcc Fe: application of local energy and local stress. J. Mater. Sci. 2014, 49, 3980–3995. [Google Scholar] [CrossRef]
  15. Gao, N.; Fu, C.-C.; Samaras, M.; Schäublin, R.; Victoria, M.; Hoffelner, W. Multiscale modelling of bi-crystal grain boundaries in bcc iron. J. Nucl. Mater. 2009, 385, 262–267. [Google Scholar] [CrossRef]
  16. Li, Y.; Chen, C.; Zhang, F. Al and Si Influences on Hydrogen Embrittlement of Carbide-Free Bainitic Steel. Adv. Mater. Sci. Eng. 2013, 382060. [Google Scholar]
  17. Lejcek, P.; Hofmann, S. Prediction of enthalpy and entropy of grain boundary segregation. Surf. Interface Anal. 2002, 33, 203–210. [Google Scholar] [CrossRef]
  18. Spatschek, R.; Gobbi, G.; Hüter, C.; Chakrabarty, A.; Aydin, U.; Brinckmann, S.; Neugebauer, J. Scale bridging description of coherent phase equilibria in the presence of surfaces and interfaces Phys. Rev. B 2016, 94, 134106. [Google Scholar] [CrossRef]
  19. Hecht, F. New development in FreeFem++. J. Numer. Math. 2012, 20, 251–266. [Google Scholar] [CrossRef]
  20. Cahn, J.W. Coherent Fluctuations and Nucleation in Isotropic Solids. Acta Metall. 1962, 10, 907–913. [Google Scholar] [CrossRef]
  21. Cahn, J.W.; Hilliard, J.E. Free Energy of a Nonuniform System. I. Interfacial Free Energy. J. Chem. Phys. 1958, 28, 258–267. [Google Scholar] [CrossRef]
  22. Tang, M.; Karma, A. Surface Modes of Coherent Spinodal Decomposition. Phys. Rev. Lett. 2012, 108, 265701. [Google Scholar] [CrossRef] [PubMed]
  23. Kresse, G.; Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. Rev. B 1993, 47, 558. [Google Scholar] [CrossRef]
  24. Kresse, G.; Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 1996, 6, 15–50. [Google Scholar] [CrossRef]
  25. Kresse, G.; Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B. 1996, 54, 11169. [Google Scholar] [CrossRef]
Figure 1. Solution (excess) enthalpy of carbon in bcc iron as function of volume. The continuum theory prediction is compared to DFT and empirical potential data. See [4] for details of the atomic scale simulations and the different empirical potentials.
Figure 1. Solution (excess) enthalpy of carbon in bcc iron as function of volume. The continuum theory prediction is compared to DFT and empirical potential data. See [4] for details of the atomic scale simulations and the different empirical potentials.
Metals 08 00219 g001
Figure 2. Hydrogen solution enthalpy as function of the number of Mn neighbors. For some cases, there are different energy values, which are due to different atomic arrangements, e.g., trans or cis configurations (see [6] for a detailed discussion). The dashed curve is a linear fit. The black arrow points to the trans configuration with two neighboring Mn atoms at opposite ends of the octahedron surrounding the H atom. This configuration is the reference for the Fe-Mn-Al-H system.
Figure 2. Hydrogen solution enthalpy as function of the number of Mn neighbors. For some cases, there are different energy values, which are due to different atomic arrangements, e.g., trans or cis configurations (see [6] for a detailed discussion). The dashed curve is a linear fit. The black arrow points to the trans configuration with two neighboring Mn atoms at opposite ends of the octahedron surrounding the H atom. This configuration is the reference for the Fe-Mn-Al-H system.
Metals 08 00219 g002
Figure 3. Spherical precipitate in the vicinity of a free surface or grain boundary. The Young’s modulus of the grain boundary is effectively either higher or lower than the value in the bulk. The sketch shows a two-dimensional cut through the three-dimensional system.
Figure 3. Spherical precipitate in the vicinity of a free surface or grain boundary. The Young’s modulus of the grain boundary is effectively either higher or lower than the value in the bulk. The sketch shows a two-dimensional cut through the three-dimensional system.
Metals 08 00219 g003
Figure 4. Elastic energy of the two-phase system with a precipitate near a free surface, in the bulk, or close to a hard or soft grain boundary. The energy is shown relative to the bulk energy value. The location of the grain boundary is at h / R = 6 , and we use a grain boundary thickness d / R = 0.8 .
Figure 4. Elastic energy of the two-phase system with a precipitate near a free surface, in the bulk, or close to a hard or soft grain boundary. The energy is shown relative to the bulk energy value. The location of the grain boundary is at h / R = 6 , and we use a grain boundary thickness d / R = 0.8 .
Metals 08 00219 g004
Table 1. Comparison of calculated and estimated grain boundary energies.
Table 1. Comparison of calculated and estimated grain boundary energies.
Grain BoundaryCalculated Value/(Jm 2 )Estimated Value/(Jm 2 )Error/(Jm 2 )
γ G B 1 H 1.491.480.01
γ G B 2 H 1.421.410.01
γ G B 3 H 1.341.34<0.01
γ G B 1 C 1.251.310.06
γ G B 2 C 0.900.930.03
γ G B 3 C 0.730.820.09
γ G B 2 H 1 C 1.081.170.09
The calculated grain boundary energies reported in [9] and their estimated values based on the mean field approximation. We do not compare the energies for four carbon and four hydrogen atoms as they coincide by definition of the approximation.

Share and Cite

MDPI and ACS Style

Weikamp, M.; Hüter, C.; Spatschek, R. Linking Ab Initio Data on Hydrogen and Carbon in Steel to Statistical and Continuum Descriptions. Metals 2018, 8, 219. https://doi.org/10.3390/met8040219

AMA Style

Weikamp M, Hüter C, Spatschek R. Linking Ab Initio Data on Hydrogen and Carbon in Steel to Statistical and Continuum Descriptions. Metals. 2018; 8(4):219. https://doi.org/10.3390/met8040219

Chicago/Turabian Style

Weikamp, Marc, Claas Hüter, and Robert Spatschek. 2018. "Linking Ab Initio Data on Hydrogen and Carbon in Steel to Statistical and Continuum Descriptions" Metals 8, no. 4: 219. https://doi.org/10.3390/met8040219

APA Style

Weikamp, M., Hüter, C., & Spatschek, R. (2018). Linking Ab Initio Data on Hydrogen and Carbon in Steel to Statistical and Continuum Descriptions. Metals, 8(4), 219. https://doi.org/10.3390/met8040219

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop