Next Article in Journal
Effect of Ball Milling Parameters on the Refinement of Tungsten Powder
Next Article in Special Issue
Ab Initio Guided Low Temperature Synthesis Strategy for Smooth Face–Centred Cubic FeMn Thin Films
Previous Article in Journal
Failure Analysis of PHILOS Plate Construct Used for Pantalar Arthrodesis Paper II—Screws and FEM Simulations
Previous Article in Special Issue
Linking Ab Initio Data on Hydrogen and Carbon in Steel to Statistical and Continuum Descriptions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling of Phase Equilibria in Ni-H: Bridging the Atomistic with the Continuum Scale

by
Dominique Korbmacher
1,
Johann Von Pezold
1,
Steffen Brinckmann
1,
Jörg Neugebauer
1,
Claas Hüter
2,3 and
Robert Spatschek
2,3,*
1
Max-Planck-Institut für Eisenforschung GmbH, D-40237 Düsseldorf, Germany
2
Institute for Energy and Climate Research, Forschungszentrum Jülich GmbH, D-52428 Jülich, Germany
3
Jülich Aachen Research Alliance (JARA), RWTH Aachen University, D-52056 Aachen, Germany
*
Author to whom correspondence should be addressed.
Metals 2018, 8(4), 280; https://doi.org/10.3390/met8040280
Submission received: 29 March 2018 / Revised: 13 April 2018 / Accepted: 16 April 2018 / Published: 18 April 2018
(This article belongs to the Special Issue First-Principles Approaches to Metals, Alloys, and Metallic Compounds)

Abstract

:
In this paper, we present a model which allows bridging the atomistic description of two-phase systems to the continuum level, using Ni-H as a model system. Considering configurational entropy, an attractive hydrogen–hydrogen interaction, mechanical deformations and interfacial effects, we obtained a fully quantitative agreement in the chemical potential, without the need for any additional adjustable parameter. We find that nonlinear elastic effects are crucial for a complete understanding of constant volume phase coexistence, and predict the phase diagram with and without elastic effects.

1. Introduction

The understanding of phase equilibria and transitions is essential for the comprehension of many processes in nature as well for a theory guided search for novel materials with superior properties. The modeling of two-phase coexistence in binary alloys is based on the common tangent construction: In equilibrium, the chemical potentials of each species in the two phases of interest have to be equal. As shown in pioneering work by several authors [1,2,3], elastic effects can significantly modify the phase coexistence behavior, especially in solid phases, e.g., due to density differences.
For a true multi-scale modeling of phase equilibria and transitions in complex materials, ranging from macroscopic dimensions down to the nanoscale, an efficient and accurate matching between the atomistic simulations and formal thermodynamic and continuum concepts is critical. However, despite significant progress, there is still a substantial gap between these two levels. The purpose of the present article is therefore to seamlessly connect the atomistic and continuum scale, illustrated for the Ni-H system. This system is used as a prototype in atomistic simulations for the understanding of hydrogen embrittlement phenomena [4] and rechargeable batteries [5,6,7]. Mechanical deformations play a crucial role, and thus a careful thermodynamic inspection is essential for a thorough understanding. In previous work [8], we connected ab initio modeling of coherent phase equilibria in the presence of lattice strains and interfacial proximity effects to continuum descriptions. Here, in contrast, the focus is on a matching between Monte Carlo (MC) simulations using empirical potentials and bulk effects on the continuum level. Additionally, we consider interfacial energy contributions, which were neglected in [8].
Experimentally, it is known that metallic nickel absorbs very limited amounts of hydrogen under ordinary conditions, but is known to form a nearly stoichiometric hydride NiH at high hydrogen fugacities. Fukai et al. performed in situ measurements of the lattice parameter of Ni-H alloys over a wide range of hydrogen partial pressures. This way they constructed a pressure–composition–temperature diagram, encompassing the region of two-phase coexistence [9,10].
One of the central outcomes of the present study is the importance on nonlinear elastic effects in constant volume simulations, as large stresses can arise [11]. Moreover, our approach also provides important fundamental insights into the theory of phase equilibria in coherent solid-state systems, as it elucidates quantitatively the role of different energetic and entropic contributions.
The present work is complementary to atomistic modeling [12] and cluster expansion approaches [13,14], which aim at a description on the atomic level. Cluster expansion approaches are based on an energy description of an Ising-like model of alloys and are powerful methods for predicting phase diagrams and ordering phenomena; long-ranged elastic effects, however, are difficult to capture in such approaches. Therefore, in the present work, we focused on the bridging between a hybrid molecular static and Monte Carlo approach for finding the equilibrium configurations on the one hand and long-wavelength continuum theories on the other hand. This link across the scales is needed to capture in a closed description both microscopic effects related to H–H interactions as well as long-ranged elastic effects due to phase separation. The derived scale bridging energy functional can then serve as a direct input e.g., to phase field simulations. The resulting quantitative scale-bridging is the major outcome of this paper. Step by step, we introduce the various contributions on the continuum scale to match the atomistically determined chemical potential of hydrogen, which is a key quantity in the description of phase coexistence.

2. Methods

The starting point of our study is the determination of the equilibrium spatial distribution of the interstitial H atoms in the metallic matrix, employing grand canonical Monte-Carlo simulations and molecular statics. The simulations are performed within a periodic 10 × 10 × 10 fcc supercell containing 4000 Ni atoms at constant volume (corresponding to the equilibrium volume of pure Ni) and temperature. H atoms are introduced on the octahedral interstitial sites according to the Metropolis algorithm [15]. The energies of the trial steps are determined by relaxing the ionic degrees of freedom—thus fully capturing elastic interactions—using the LAMMPS simulation package [16,17] in conjunction with a modified version [4] of the Ni-H EAM potential by Angelo et al. [18,19]. The effects of temperature are thus only considered to account for the configurational entropy contributions to the free energy, while vibrational contributions are disregarded. This description has previously been used to model hydrogen mediated dislocation–dislocation interaction [4].
The continuum modeling involves the construction of a free energy functional, which is derived step by step in the following section. This partly requires solving continuum mechanics problems, for which we used in particular finite element modeling using ABAQUS.

3. Results

3.1. Monte Carlo Modeling

In the grand canonical simulations, the chemical potential is prescribed, and the number of hydrogen atoms can vary. Depending on the H chemical potential dilute or condensed H distributions are obtained. The condensed hydride precipitates remain coherent and adopt characteristic shapes depending on the bulk H concentration, as shown in Figure 1.
A peculiarity, which motives the present scale bridging investigations, is the functional form of the hydrogen chemical potential μ H as function of the average H concentration, see dots in Figure 2, which will be discussed in detail below. Classically, we would expect (in a strain free situation) the chemical potential to be constant in the entire two-phase region as a consequence of the common tangent construction. Obviously, this is not the case here, and it is one of the goals of this paper to shed light on this effect.

3.2. Free Energy Formulation

Our starting point for transferring the atomic scale behavior to the continuum level is the free energy for the single phase material,
F = μ 0 N H + F c + F H-H + F el
with the number of hydrogen atoms N H and the elastic free energy F el , the configurational free energy F c and the H–H interaction F H H . μ 0 is the solvation energy needed to insert an isolated hydrogen atom into the (empty) matrix, in contrast to the aforementioned chemical potential μ H , which also includes mutual interactions, elastic effects, etc.

3.3. Solvation Energy and H–H Interaction

For the formation of a hydride phase, it is essential to include a description of the lattice mediated interaction between the hydrogen atoms [4]. We use the Margules model [5], F H-H / N Ni = α c 2 / 2 + β c 3 / 3 , where N Ni is the number of nickel atoms and c = N H / N Ni the (homogeneous) hydrogen concentration; it is normalized to 1 for a crystal where all interstitial octahedral sites are occupied. The parameters are determined from the Monte Carlo calculations by averaging over several random mixtures of fully relaxed homogeneous solutions. We point out that, for this matching, we only consider homogeneous “lattice gas” configurations and not phase separated states, which involve additionally (coherency) elastic and interfacial effects, and which will be parametrized below. These reference calculations are performed under constant pressure conditions, P = 0 , to suppress (external) stress effects. The resulting curve F ( c ) is shown in Figure 3. It is fitted by a polynomial F f i t / N Ni = a 0 + a 1 c + a 2 c 2 + a 3 c 3 . Identification with the above parameters therefore uniquely gives μ 0 = a 1 , α = 2 a 2 and β = 3 a 3 ( a 0 is an irrelevant zero-point energy). Hence, we obtain μ 0 = 2 . 148 eV relative to monatomic hydrogen in vacuum, α = 0.751 eV and β = 0.355 eV . These numbers differ slightly from those given in Ref. [5] due to a different pair potential cutoff to stabilize the hydride, as this leads to a positive value of the elastic constant C 44 , see Ref. [4] for details. The positive sign of α favors an increase of the local hydrogen concentration, thus leading to the hydride formation.
The chemical potential is defined in the usual way,
μ H = F N H T , V , N Ni ,
and is split into additive contributions analogous to the free energy in Equation (1). We therefore get for the chemical potential contribution by the H–H interaction μ H H = α c + β c 2 . It is shown in Figure 4 for both the present fitting parameters and the ones given in Ref. [5], exhibiting only small differences between the predictions.

3.4. Configurational Entropy

The configurational free energy F c stems from the different possibilities to occupy the interstitial sites with hydrogen (or the hydride with vacancies), and is therefore given by
F c = k B T N Ni [ ( 1 c ) ln ( 1 c ) + c ln c ]
under the assumption of equal occupation probability for all octahedral sites [20]. In the low concentration regime, it gives the dominant nontrivial contribution, as elastic and interaction effects are negligible there, see Figure 5.
We note that the configurational term is the only temperature dependent one in the present description. To verify this dependence we also computed in the Monte Carlo simulations cases with different temperatures, and the results are shown in Figure 6, showing very good agreement.

3.5. Maxwell Construction

Thus far, the system is described for a single phase state only, and elastic effects are suppressed by a free volume expansion. Due to homogeneity, elastic stresses as well as interfacial contributions do not appear, and the description is therefore complete on this level.
According to the usual thermodynamic picture, a single phase state is unstable in regions where the slope of the chemical potential is negative, μ H ( c ) < 0 , and phase separation should occur there. As can be seen in Figure 6, this is the case for a wide concentration regime for low temperatures, whereas for high temperatures the hydrogen solubility limit is significantly larger.
Phase coexistence on this level is described by Maxwell’s equal area rule, which states that phase separation sets in at the intersection points of a horizontal line with the S-shaped van der Waals loop, cutting it into two equal areas above and below this Maxwell line. This is shown in Figure 2 for T = 300 K . From this, we get the constant chemical potential in the two-phase region, μ M = 2 . 405 eV for T = 300 K . In equilibrium, the system therefore enters the two-phase region at the first intersection of the horizontal Maxwell line with the S-shaped single phase curve. Apparently, this happens already for very low concentrations in the ppm range. Nevertheless, on the following ascending part of the curve μ H ( c ) , the system is still metastable. We point out that, according to the principles of the two-phase construction, the dual phase hydrogen chemical potential should be constant and equal to μ M until the phase transformation is completed. Obviously, this is not the case for the atomistic data with the fixed volume constraint.

3.6. Elastic Effects

Before entering into a detailed analysis of elastic effects, let us briefly discuss them in relation to the H–H interaction term. The latter is sometimes also introduced as long-range elastic interaction, and we want to stress that we do not double count effects here. The H–H interaction term expresses that an H atom locally deforms the lattice, and therefore makes the placement of another H atom in the vicinity energetically favorable. In a mean field context, this is expressed by the term ∼ c 2 in F H H . If the site occupancy gets higher and the octahedral sites more and more filled, the energy for placing more H atoms into the lattice increases again, as expressed by the + c 3 term. Although this consideration makes reference to elastic deformations, this is distinct from the elastic effects considered in the following. To make this point more explicit, we consider the following two situations, for simplicity both with homogeneous hydrogen distributions, i.e., spatially constant concentration c:
First, for fixed pressure P = 0 , we increase the hydrogen concentration. This leads to a widening of the lattice, as discussed in the following and shown in Figure 7. However, since the system remains stress free on a mesoscopic or macroscopic level, σ i j = 0 , the elastic energy (see Equation (5)) in the continuum mechanics sense remains zero. Nevertheless, the total energy of the system is changed due to the H–H interaction, and this is expressed through F H-H , which is independent of the elastic stress state.
Second, for fixed concentration (and total number of hydrogen atoms), we change the external pressure of the system. Then, the H–H interaction term does not change, but the mesoscopic elastic energy does.
The difference between the H–H interaction term and the explicit elastic term becomes most prominent in two-phase situations. Then, the hydride has a larger equilibrium lattice constant and therefore distorts the surrounding, coherently connected matrix phase. This leads to long-ranged elastic deformations, and the range of these interactions of the order of the precipitate size, which is captured by the mesoscopic elastic energy term introduced below, but not by the H–H interaction term. For further discussion of this separation of microscopic and mesoscopic elastic contributions, we refer to [8]. There, it was shown that the elastic energy due to microscopic deformations around individual H atoms and long ranged strains resulting from the formation of mesoscopic precipitates decompose additively into microscopic and mesoscopic contributions without the appearance of a cross term.
The large deviations between the atomistic data and the continuum model clearly indicate that elastic and interfacial effects are critical and cannot be neglected. We remind that the basis for the Maxwell construction—or here equivalently the common tangent construction—is based on the assumption that the two-phase energy in a macroscopic system decomposes additively into the contributions of the two phases, which do not influence each other. This means, that they are described by free energy densities, which depend only on the local concentration. This condition, however, is violated in the presence of elastic effects. Here, e.g., a density variation in one phase also influences the other, as the global pressure in the system changes.
In the following, we discuss the modeling mainly from the continuum perspective. The technical aspects of the parameter matching between the scales are given in the appendices.
We focus first on the low-temperature regime, where the hydrogen solubility limit is low (see Figure 13), and higher temperatures will be discussed later. In the present regime, the phases are almost stoichiometric, and the chemical potential of the hydrogen in the two-phase region is modified. From a minimization of the total free energy in the two-phase region, we obtain the modified chemical potential (see Appendix A)
μ H ( c ) = μ M + μ el ( c ) + μ s ( c ) .
In contrast to the Maxwell term μ M , the elastic and interfacial contributions, μ el and μ s , respectively, do depend on the concentration, and therefore the chemical potential is no longer constant in the two-phase region, as observed in the Monte Carlo simulations (see Figure 2).

3.6.1. Linear Elastic Effects

The pure nickel and the hydride exhibit a substantial lattice mismatch, and therefore, in the two-phase region, elastic stresses arise. The linear elastic free energy density per unit volume is given by
f el = 1 2 C i j k l ( ϵ i j ϵ i j 0 ) ( ϵ k l ϵ k l 0 ) ,
where ϵ i j = ( i u j + j u i ) / 2 is the strain derived from the displacements u i , and the eigenstrain ϵ i j 0 = ε 0 δ i j represents the isotropic volume expansion of the hydride [21]. It is concentration dependent, and we extracted it again from T = 0 K simulations of homogeneous states with P = 0 ; it is well described by ε 0 ( c ) = b 1 c + b 2 c 2 + b 3 c 3 with fitting parameters b 1 = 0.1027 , b 2 = 0.0598 and b 3 = 0.0189 , going beyond a linear dependence in the spirit of Vegard’s law. The eigenstrain is directly related to the concentration dependent lattice constants with a Ni = 3.520 and a Ni H = 3.738 via ε 0 ( c ) = [ a ( c ) a Ni ] / a Ni , since we use here the relaxed, hydrogen-free nickel as reference configuration (see Figure 7).
The elastic constants for the pure phases are listed in Table 1. We note that all values are given with respect to Ni as reference configuration; see Appendix B for a discussion of this issue. The determination of the elastic constants from the Monte Carlo simulations is described in detail in Appendix C.
To better understand the role of the elastic effects, it is instructive to inspect an analytical isotropic linear elastic model, where we assume for simplicity the elastic constants to be equal in both phases. We use a spherical sample of radius R consisting of nickel, with a hydride inclusion of radius R H in its center. The radial displacement component is continuous at the coherent interface and vanishes at the outer boundary due to the volume constraint. At low temperatures, the concentration is in good approximation c = 0 in Ni and c = 1 in the hydride. This gives (see Appendix D)
μ lin . el = a Ni 3 ( 3 λ + 2 G ) [ 2 G + c ( 3 λ + 2 G ) ] ( Δ ε 0 ) 2 ( λ + 2 G ) N 0 ,
where N 0 = 4 is the number of octahedral sites per unit cell, Δ ε 0 = ε 0 ( 1 ) ε 0 ( 0 ) , and λ and G are Lamé coefficient and shear modulus respectively. In this approximation, the elastic contribution in Equation (6) is a linear function of the concentration, and it depends quadratically on the eigenstrain difference Δ ε 0 . Since the material is relaxed for c = 0 , the elastic effects therefore give only an additive constant for c = 0 , which is related to the elastic hysteresis due to the internal stresses that arise upon precipitate formation. For the isotropic elastic constants we use the values for nickel and λ = C 12 and G = ( C 11 C 12 ) / 2 for the solid curve μ M + μ lin . el in Figure 8.
The calculation can also be done using the elastic constants of the individual phases, leading to the long-dashed curve shown in the same graph (see Appendix D). A numerical finite element solution with ABAQUS for spherical precipitates in a cubic box with fixed volume, taking into account the full cubic elasticity and periodic boundary conditions, leads to the short-dashed curve shown in Figure 8. It is very close to the analytical results. Since we still observe a significant discrepancy to the Monte Carlo data, we conclude that the consideration of linear elasticity is not sufficient to explain the slope of the chemical potential in the two-phase region.

3.6.2. Nonlinear Elastic Effects

The reason for the discrepancy between the atomistic data and the continuum modeling involving linear elasticity is the appearance of large compressive stresses for higher hydrogen concentrations (at fixed volume), and therefore the linear elastic approximation breaks down. Instead, we have to take into account nonlinear elastic effects; they include nonlinearities in the stress-strain relationship, leading effectively to expressions C i j k l ( { ϵ m n } ) in the elastic free energy density. Additionally, geometrical nonlinearities appear due to the volume deformation and the nonlinear strain tensor [11],
ϵ i j = 1 2 u j x i + u i x j + u k x i u k x j .
We have estimated that these effects are small in comparison to the change in the constitutive law, and therefore effectively captured them in the modified elastic constants. However, it turns out to be important to carefully use consistently the same reference configuration for the Lagrangian formulation of elasticity for both phases. We note that the chemical potential contains now also a contribution due to the strain dependence of the elastic constants. We mention in passing that we do not observe plastic deformations in the atomistic simulations. From a careful analysis of the Monte Carlo simulations we find that the nonlinear corrections depend only on the volume change and not on shear effects, thus the elastic constants depend only on the trace of the (local) strain tensor, which simplifies the expressions. We therefore write them as series expansion
C = C 0 1 + d 1 C tr ( ϵ ϵ 0 ) + d 2 C [ tr ( ϵ ϵ 0 ) ] 2 +
for all elastic constants or combinations of them. From the atomistic simulations we find the coefficients of the nonlinearities for C 11 C 12 , C 44 and the bulk modulus K, see Table 2 and Appendix C. From them, all elastic constants can be expressed using the relation 3 K = C 11 + 2 C 12 . We note that for the present geometry bulk compression is the dominant effect, and the material becomes stiffer in this regime of negative strain (compression). Including the nonlinear contributions the continuum chemical potential shows now a very satisfactory agreement with the Monte Carlo data, see the dotted line in Figure 8.

3.7. Interfacial Effects

In a final step, we take into account interfacial effects. As already visible from the agreement of the Monte Carlo data and the nonlinear elastic curve, their contribution is obviously small. In isotropic approximation, the surface energy for a spherical inclusion is given by F s = 4 π γ R H 2 with the surface energy γ . From planar interface calculations, we obtain γ = 0.11 Jm 2 for (100) interfaces and γ = 0.17 Jm 2 for (111) interfaces. From the fitting of two-phase data, we independently get an excellent agreement using γ = 0.17 Jm 2 for the spherical inclusion.
The magnitude of interfacial terms should be compared to bulk contributions to see their influence and to get an impression on the dimension of critical nuclei. For an order of magnitude estimate, we compare elastic bulk contributions in the approximation in Equation (6) to interfacial contributions for spherical precipitates for small concentrations in Equation (A42), where the influence of interfacial terms is largest. For c 1 the elastic and interfacial contributions are comparable for
N H 4 γ 3 ( 2 G + λ ) 3 n 0 π 3 a Ni 3 ( Δ ε 0 ) 6 G 3 ( 2 G + 3 λ ) 3 ,
which is the case for about 23 H atoms. This indicates that interfacial contributions are indeed small for the present scenario in comparison to the elastic bulk terms.
For a linear isotropic material with equal elastic constants in both phases and a purely dilatational eigenstrain, the elastic energy does not depend on the shape of the inclusion but only on its volume fraction [22]. Although the conditions for this rigorous statement are not exactly fulfilled here, we find only a small dependence of the elastic energy on the precipitate shape in the Monte Carlo data, as exemplarily shown for spherical hydride inclusions in a Nickel matrix and vice versa in Figure 9. Consequently, the shape of the precipitate is determined by the interfacial energy terms alone, see Appendix E for explicit calculations. This means, that for each given average concentration c, which is related to a precipitate size, the equilibrium shape with the minimum interfacial energy appears. For isotropic surface energy we expect spherical inclusions for c < 4 π / 81 0.16 , followed by tubes up to c = 1 / π 0.32 , and slabs thereafter. The situation is symmetric for c > 0.5 with then the hydride being the matrix phase. Thus we expect nickel tubes in the range 1 1 / π < c < 1 4 π / 81 and spheres beyond this concentration. In our Monte Carlo calculations, we indeed find all these structures in the correct ordering (see Figure 1 and Figure 10), but the transition points differ slightly due to anisotropy and nonlinear elastic effects. Based on the observed structures and the calculated interfacial energy we added this contribution to the nonlinear elastic chemical potential in Figure 10. The transitions between different precipitate shapes explain the small discontinuities in the chemical potential, as computed from the Monte Carlo simulations. Apparently, with all the aforementioned effects taken together, we get an excellent agreement with the Monte Carlo data in the entire two-phase region.
Finally, we briefly comment on the role of surface stress. In isotropic approximation, the surface energy for a spherical inclusion is given by F s = 4 π γ 0 R H 2 with the surface energy γ 0 . The surface stress contribution is F β = 4 π R H 2 β ( ϵ θ θ + ϵ ϕ ϕ ) in spherical coordinates, with the scalar surface stress coefficient β [23,24]. We obtain—in the same approximation as for Equation (6)—due to the additional pressure difference between the precipitate and the matrix Δ σ r r = 2 β / R H also a contribution from the elastic bulk energy Δ F el [25]
F β + Δ F el = 8 3 π 3 λ + 2 μ λ + 2 μ β Δ ε 0 R H 2 ,
where we have neglected a small destabilizing term proportional to β 2 for | β | 3 K Δ ε 0 R H (we have estimated from the atomistic data β to be of the order 1 Jm 2 , thus this condition is fulfilled), and also assumed c 1 , since spherical precipitates appear for small concentrations only, as discussed below. Consequently, these terms have the same radius scaling as the surface energy term ∼ R H 2 , and despite their partial origin from a bulk energy, they appear effectively as surface term; therefore, we treat them via a renormalized interfacial energy density γ .
With this, the parametrization of the continuum model is completed, and it is used in the following for predicting hydrogen partial pressures, volume relaxation effects, and for the prediction of phase diagrams.

3.8. High Concentration Limit

We can also correctly describe the limit of high hydrogen concentrations, where almost all octahedral sites are filled with H. For 1 c 1 , the system is again in a single phase state, with a dilute distribution of vacancies, thus the steep divergency of the chemical potential is due to the configurational degrees of freedom. Nevertheless, the central difference to the dilute limit c 1 is that the material is here severely under (nonlinear) stress (as we focus on fixed volume situations, where pure Ni is stress free, hence stresses are maximized for c = 1 ) , thus elastic effects play a major role here. In contrast to the phase separated region we need precise knowledge of the concentration dependence of the elastic constants, their nonlinearities and the eigenstrain. Consideration of all these effects leads to the good agreement between the Monte Carlo data and the continuum picture, as shown in Figure 10.

3.9. Conversion to Partial Pressures

To relate the chemical potentials to experimentally accessible quantities, it is useful to translate them to hydrogen (H 2 ) partial pressures. In the following, we use a subscript H 2 to discriminate between quantities related to the gaseous hydrogen molecules and the monatomic hydrogen in solid solution. From the ideal gas equation and the Maxwell relation
μ H 2 p H 2 T , N H 2 = V N H 2 T , p H 2 = k B T p H 2
we obtain by integration the chemical potential of H 2 up to an additive term μ H 2 0 ,
μ H 2 ( T , p H 2 ) = μ H 2 0 ( T ) + k B T ln p H 2 p ref
with an arbitrary reference pressure p ref . From the equilibrium between H 2 and dissociated hydrogen follows
μ H 2 = 2 μ H .
This leads to
p H 2 = p ref exp 2 μ H μ H 2 0 ( T ) k B T .
We therefore need an expression for μ H 2 0 ( T ) to link the chemical potential to a hydrogen partial pressure.
In the dilute limit c 1 , the chemical potential is (see also Figure 5)
μ H = μ 0 + k B T ln c .
In combination with Equations (12) and (13), we therefore get Sievert’s law
c = p H 2 p ref 1 / 2 exp μ 0 1 2 μ H 2 0 ( T ) k B T .
This is in agreement with experimental findings [26] for the solubility of H 2 in Ni
c = 2 c 0 p H 2 p 0 1 / 2 exp 1485 . 167 K T + 3 . 9881 ,
expressed in terms of the quantities used here, with c 0 = 5.8694 · 10 7 and p 0 = 133.3224 Pa (the prefactor 2 stems from the dissociation of a H 2 molecule in two H atoms). This allows to identify an expression for the offset potential μ H 2 0 ( T )
μ H 2 0 ( T ) = 2 μ 0 k B · 2970.334 K + 7.9762 k B T = 4.552 eV + T · 6.873 · 10 4 eV K 1 ,
where we have chosen p ref = p 0 / ( 4 c 0 2 ) . With these identifications, Equation (14) is the desired relation, which is shown in Figure 10. Apparently, the H 2 partial pressures needed for the higher concentrations are extremely high and not reachable in experiments under normal conditions, making computer simulations particularly useful there. In addition, we note that for such high pressures the used ideal gas description of H 2 breaks down. On the other hand, the H 2 partial pressures needed for formation of a hydride are much lower for free volume situations, which will be discussed below (see Figure 11). We can therefore conclude that volumetric constraints, or—more generally—compressive stresses can suppress the formation of hydrides, which can be responsible for the embrittlement of the metal.

3.10. Volume Relaxation

Instead of a fixed volume constraint, one may also consider a situation with free expansion, P = 0 . This does not mean that elastic effects disappear, because in the two-phase region internal coherency stresses still arise. Thus, the same elastic barrier appears for the first nucleation of the precipitate, but thereafter the elastic contribution to the chemical potential decreases with concentration, see Figure 11. Since the stress effects are much lower, elastic nonlinearities are negligible here, and the curve is calculated in linear isotropic approximation. The qualitatively different behavior shows the important role of external boundary conditions. In contrast to scenarios with free expansion, the volume constraint stabilizes two-phase equilibria. This behavior can intuitively be understood as follows: When a precipitate is inserted into the matrix, coherency stresses arise, which increase the elastic energy. If the precipitate grows and finally fills the whole sample, the freely expanding material is homogeneous and therefore stress free. This implies that the elastic energy decays with the hydride volume fraction, and therefore the elastic contribution to the chemical potential appears with negative sign. Therefore, we get a hydrogen chemical potential with a negative slope in the two-phase region, which corresponds to an unstable situation. We have verified this prediction using Monte Carlo simulations for fixed and free volume at T = 300 K , see Figure 12.
We note that volume constraints can therefore stabilize precipitates. In a simulation, this is useful for probing situations with phase coexistence, which would otherwise not appear, as either the system is homogeneously in the dilute state or fully saturated with H for given chemical potential. This is expressed here by a negative slope of the chemical potential of H. The present physically motivated approach is therefore an alternative to bias potentials [12] to probe configurations with phase coexistence.

3.11. Phase Diagram

Based on the temperature dependence of F c , we can predict the entire bulk phase diagram without elastic effects via the Maxwell construction, see Figure 13. We note that this phase diagram is derived solely from T = 0 K data for the H–H interaction. To verify it also for higher temperatures, we have compared the analytical expression for the chemical potential to the single-phase data from grand-canonical Monte Carlo simulations with free volume relaxation, where elastic and interfacial effects do not appear, and find an excellent agreement (see Figure 6). We note that the purpose of the present work is to match a continuum description of Ni-H to Monte Carlo simulations based on an empirical interatomic potential. The agreement between the two theoretical approaches is excellent. Compared to experimental values, however, the description is less suitable for high temperatures, and therefore it is not accurate in this regime. The reason is—apart from potential limitations of the interatomic potential—the use of the (lattice) Monte Carlo molecular statics description, which neglects in particular phonon excitations.
An application that we can readily extract from the above results is the phase diagram taking into account elastic effects. In particular, we consider the situation of fixed volume and determine the two-phase region by minimization of the total free energy. For the sake of simplicity, we assume a slab geometry with the phase normal in [ 100 ] direction (see Figure 14), which is the most relevant pattern in the medium concentration regime, as discussed above in Figure 10. In contrast to other geometries, homogeneity implies the constancy of concentrations in each phase. We neglect the elastic nonlinearities, which allows for a straightforward solution of the underlying one-dimensional elastic problem. As the elastic constants were determined mainly for c = 0 and c = 1 , we assume additionally that they interpolate linearly for arbitrary concentrations 0 < c < 1 . This is mainly relevant for higher temperatures, where the appearing phases have larger solubilities for hydrogen or vacancies, thus they appear with intermediate concentrations. As will be discussed below, it turns out that the concentration dependence of the elastic constants is not important, and therefore this approximation is legitimate.
We inspect the case of a constant total volume which is stress free for pure nickel, i.e., the lattice constant for a homogeneous system is a Ni . Consequently, for the two-phase system with cubic symmetry, the only non-vanishing displacement component is u x , which is a linear function of x in each phase, u x ( i ) = a i x + b (the phases are enumerated by i). Hence, the only nontrivial strain component is ϵ x x ( i ) = a i ; stresses follow from Hooke’s law with cubic symmetry,
σ x x ( i ) = C 11 ( c i ) ( a i ε 0 ( c i ) ) 2 C 12 ( c i ) ε 0 ( c i ) ,
σ y y ( i ) = σ z z ( i ) = C 12 ( c i ) [ a i 2 ε 0 ( c i ) ] C 11 ( c i ) ε 0 ( c i ) .
Because the concentration is homogeneous in each phase for this effectively one-dimensional situation, the linear bulk elastic equations are automatically fulfilled by the above ansatz. The coefficients a 1 , b 1 , a 2 , b 2 are determined by the boundary and matching conditions, u x ( 2 ) ( x = 0 ) = 0 , u x ( 1 ) ( x = L ) = 0 , u x ( 1 ) ( x 0 ) = u x ( 2 ) ( x 0 ) and σ x x ( 1 ) ( x 0 ) = σ x x ( 2 ) ( x 0 ) . The elastic energy density is in each phase
f el ( i ) = 2 C 12 ( c i ) ε 0 ( c i ) [ 2 a i + 3 ε 0 ( c i ) ] + C 11 ( c i ) [ a i 2 2 a i ε 0 ( c i ) + 3 ε 0 ( c i ) 2 ] .
From that, we can calculate the elastic energy per unit area,
F el = f el ( 1 ) ( L x 0 ) + f el ( 2 ) x 0 .
The total free energy is then F = μ 0 N H + F c + F H H + F el . Additionally, the concentrations c 1 and c 2 in the two phases are related to the (fixed) average concentration c via the lever rule,
L c = ( L x 0 ) c 1 + x 0 c 2 .
We can therefore investigate F / L as function of c 1 , c 2 for given values of T , c , and the phase fractions are determined by the lever rule. The concentration domain can be restricted to c 1 < c < c 2 , and the free energy is minimized with respect to c 1 and c 2 in equilibrium. Whenever the minimum is located inside this domain, the system is in a two-phase state, whereas the minimum is located on the border c 1 = c or c 2 = c for a single phase state. By performing this minimization in the entire c , T plane one obtains the elastic phase diagram.
Already at this level we see the tremendous influence of the elastic effects in Figure 13, with an enlargement of the two-phase region towards higher temperatures and the hydrogen rich side. This unusual behavior—one would intuitively expect the suppression of phase separation since coherency stresses are energetically unfavorable—is due to the fact that the eigenstrain is a concave function of the concentration, ε 0 ( c ) < 0 (see Figure 7), thus volumetric deviations from Vegard’s law play the dominant role here.
To investigate this peculiarity in more detail, we look at the elastic energy depending on the phase concentrations. We start with a situation which differs from the true Ni-H problem. Instead, we assume a linear dependence of the eigenstrain on concentration (see Figure 7) (Vegard’s law). Here, we find that single-phase states have the lowest elastic energy in a wide region of the phase diagram and are therefore energetically favorable. Thus, phase separation is suppressed by elastic effects in this fictitious situation, leading to an increased H solubility limit, and this is shown in Figure 15. In contrast, for the true concentration dependence of the lattice constant on the H concentration in Figure 7, we obtain a reduced solubility limit of hydrogen in the nickel matrix. As mentioned above, this is a counterintuitive result, as the appearance of phase separation should usually increase the elastic energy and therefore make this state thermodynamically unfavorable, in contrast to our observation. Here, the situation is opposite: Single-phase states are more expensive from point of view of elastic energy, and the hydrogen solubility limit is therefore reduced. The explanation is that the single-phase case has a higher average eigenstrain, and therefore the mechanical volume work is higher than for the phase separated case for the anticipated fixed volume ensemble.
As a result, we find that the ϵ 0 ( c ) dependence has a tremendous influence on the phase diagram. This effect is much more pronounced that the concentration dependence of the elastic constant. If we assume, instead of the linear interpolation between the values for Ni and Ni-H, just concentration independent constants using the values of Ni, we find only a small shift of the binodal (see Figure 15). We can therefore conclude that the concentration dependence of the lattice constant is the most relevant quantity for the influence of elastic effects on hydrogen solubility limits.

4. Discussion and Conclusions

Understanding the thermodynamics of phase separation across the scales is a key for deriving macroscopic material properties from fundamental microscopic descriptions. Often, such a transfer is not straightforward, but, with careful descriptions, it is possible to obtain not only macroscopic parameters but also describe energy functionals in agreement with the underlying behavior on atomic scales. We demonstrated such an approach for the Ni-H system with a special focus not only on classical thermodynamics but also strong mechanical and interfacial effects.
In detail, we have seamlessly matched the results of hybrid molecular static and Monte-Carlo simulations in a binary system with full consideration of elastic effects to a continuum description for large scale simulations. The accurate identification of configurational, interaction, nonlinear elastic and surface effects allows for a thorough understanding, which will be equally useful and applicable to other mechanisms or material systems. Presently, the description is most useful for the low temperature regime below the Debye temperature. The consideration of vibrational degrees of freedom and defect formation, as well as the consideration of more complex phase diagrams will be the subject of future investigations. It should be pointed out that the accuracy of the predictions is limited by the atomic scale description, to which the continuum energy functional is tailored. The inclusion of the effects mentioned before will presumably require the consideration of additional energy terms and variables such as defect concentration.
The determination of phase diagrams demonstrates the opportunities arising from the transfer of atomistic data to the mesoscale: The quantitative determination of the free energy allows for example to construct the entire phase diagram in Figure 13 computationally much more efficiently than using atomistic simulations which require large system sizes and long relaxation times, especially in the room temperature regime. Moreover, the description serves as an accurate basis for mesoscale simulation methods such as phase field [27].

Acknowledgments

This work was supported by the DFG Collaborative Research Center 761 Steel ab initio. The authors gratefully acknowledge the computing time on the supercomputer JURECA at Forschungszentrum Jülich.

Author Contributions

All authors contributed equally to this work.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviation is used in this manuscript:
MCMonte Carlo

Appendix A. Thermodynamic Framework

To understand the influence of elastic effects on the phase diagram and the chemical potential, we start with considering the free energy of the system. We use the following notations: N Ni is the total number of nickel atoms, N Ni d is the number of Ni atoms in the dilute phase, N Ni h = N Ni N Ni d is the number of Ni atoms in the hydride, and n 0 is the number of nickel atoms per unit cell (4 for fcc). Similarly, N H is the total number of hydrogen atoms, with N H d of them being in the dilute phase and N H h = N H N H d in the hydride. In the fully saturated hydride, we have N 0 hydrogen atoms per unit cell (here, N 0 = 4 ), which corresponds to the normalization c = 1 . Thus, we define the hydrogen concentration in the dilute phase c d = N H d n 0 / ( N Ni d N 0 ) , and in the hydride, c h = N H h n 0 / ( N Ni h N 0 ) . The average concentration is c = N H n 0 / ( N Ni N 0 ) .
The free energy of the system is
F = N Ni d f d ( c d ) + N Ni h f h ( c h ) + F el ( N Ni d , N Ni h , c d , c h ) .
Here, f d is the free energy density (per host atom) of the dilute phase, and f h is the same for the hydride. The elastic energy depends on the volume fractions and the concentrations in each phase. Here, we already suppressed the degree of freedom of the shape of the hydride inclusion and assume that the free energy is already minimized with respect to it. Notice that we consider here only situations with fixed global strain (i.e., fixed volume), therefore the free energy is the appropriate thermodynamic functional.
As a further approximation of the elastic energy for low temperatures, we take into account that the dilute phase is almost free of hydrogen ( N H d 0 ), whereas the hydride is almost fully saturated ( N H h N Ni h N 0 / n 0 ) . This implies N H h N H . Therefore, we obtain the simplification F el N Ni f el ( c ) with the intensive elastic energy density f el .
For the determination of the two-phase equilibrium, we proceed in the usual way: The free energy has to be minimized with respect to the phase fractions and the partitioning. Therefore, we get
F N H d N H , N Ni , N Ni d = μ d ( c d ) μ h ( c h ) = 0 ,
where we defined the chemical potentials
μ d / h = n 0 N 0 d d c d / h f d / h ( c d / h ) .
The equality of grand potentials becomes
F N Ni d N H , N Ni , N H d = f d ( c d ) c d μ d f h ( c h ) + c h μ h = 0 .
Notice that the elastic terms do not appear in the common tangent construction in the framework of the above approximations, since the concentration is assumed to be fixed in this approximation, which holds here for low temperatures.
Next, we calculate the chemical potential of the hydrogen in the two-phase region. It is given by
μ H = d F d N H N Ni = F N H N Ni , N Ni d , N H d + F N Ni d N Ni , N H , N H d d N Ni d d N H N Ni + F N H d N Ni , N H , N Ni d d N H d d N H N Ni = F N H N Ni , N Ni d , N H d ,
since the last two terms vanish by the equilibrium conditions in Equations (A2) and (A4). Here, we obtain
μ H = μ h + n 0 N 0 f el c ,
which corresponds to Equation (4) in the main text. Interfacial effects can be considered in a similar way.

Appendix B. Reference State

For the elastic problem, we describe the material in the Lagrangian reference frame of the stress-free nickel. The eigenstrain of the hydride is therefore ϵ i j 0 = [ a Ni H a Ni ] / a Ni δ i j ; notice that also for the hydride the (eigen-)strain is defined relative to the lattice constant of the relaxed Ni. To make this point more transparent, we denote all quantities with Ni-H as reference state by a tilde, e.g., a dilatational strain by ϵ ˜ i j = [ a a Ni H ] / a Ni H , thus apart from an additive eigenstrain contribution ϵ i j = ϵ ˜ i j × a Ni H / a Ni . Since stresses are defined as force per area in the reference configuration, we have the transformation rule σ i j = σ ˜ i j × ( a Ni H / a Ni ) 2 . From Hooke’s law, σ i j = C i j k l ϵ k l and σ ˜ i j = C ˜ i j k l ϵ ˜ k l , we therefore conclude C i j k l = C ˜ i j k l × a Ni H / a Ni . This means that that the elastic constants of the hydride relative to nickel as reference state appear by a factor a Ni H / a Ni higher than their nominal values for Ni-H reference configuration.

Appendix C. Elastic Constants

Here, we briefly explain the extraction of the elastic constants from the Monte Carlo simulations. Notice that in the absence of vibrational contributions the elastic constants are temperature independent and are therefore calculated at T = 0 K . In all cases, the relaxed equilibrium configuration of pure Ni is taken as reference state, see Appendix B.

Appendix C.1. Bulk Modulus

To get the bulk modulus K for the phases, we strain the system isotropically with ϵ x x = ϵ y y = ϵ z z = ϵ . The resulting energies are represented in Figure A1. In linear elasticity, the energy increases quadratically with the strain,
f el = 9 2 K ϵ 2 ,
and we determine the bulk modulus K from the curvature of the measured energy–volume curves. Hence, we get for pure Nickel K Ni = 180.5 GPa and for the fully saturated hydride K Ni H = 230.0 GPa .
Figure A1. Energy density as function of the external hydrostatic strain ϵ . With the relation in Equation (A7), we get the bulk modulus K in quadratic approximation. The atomistic data are taken from simulations at T = 0 K . For high strains deviations from the linear elastic response can be seen.
Figure A1. Energy density as function of the external hydrostatic strain ϵ . With the relation in Equation (A7), we get the bulk modulus K in quadratic approximation. The atomistic data are taken from simulations at T = 0 K . For high strains deviations from the linear elastic response can be seen.
Metals 08 00280 g0a1

Appendix C.2. Monoclinic Distortion

The monoclinic distortion describes shear effects to get the elastic constant C 44 both for Ni and the hydride. One gets the translational movement of an atom from the original location r = r ( x , y , z ) to the new position r = r ( x , y , z ) ,
r = x + α y / 2 α x / 2 + y 1 + α 2 / ( 1 α 2 ) z
with a deformation parameter α . With the displacement u = r r and the linear strain tensor ϵ i j = ( i u j + j u i ) / 2 , we obtain
ϵ x y = α 2 and ϵ z z = α 2 1 α 2 = O ( α 2 ) ,
whereas all other strains are zero and the quadratic term in α can be neglected in a linear approximation. Thus, the elastic energy is
f el = 1 2 σ i j ϵ i j = 1 2 C 44 α 2 .
Figure A2 shows the energy density of both phases depending on the strain ϵ x y for the monoclinic distortion, and the elastic constant C 44 is obtained from the curvature at the origin. We obtain a value of C 44 = 134 . 3 GPa for Nickel and C 44 = 33 . 5 GPa for the hydride. Notice that for pure shear the quadratic energy–deformation dependence holds up to several percent of strain.
Figure A2. The results for the monoclinic distortion are shown. The atomistic data is fitted via a quadratic ansatz.
Figure A2. The results for the monoclinic distortion are shown. The atomistic data is fitted via a quadratic ansatz.
Metals 08 00280 g0a2

Appendix C.3. Orthorhombic Distortion

The orthorhombic distortion delivers the result for the term C 11 C 12 for both phases. It is described by
r = ( 1 + α ) x ( 1 α ) y 1 + α 2 / ( 1 α 2 ) z ,
hence
ϵ x x = α , ϵ y y = α and ϵ z z = α 2 1 α 2 ,
where ϵ z z can be neglected again with the same argument as above. Consequently, the elastic energy is
f = 1 2 σ i j ϵ i j = ( C 11 C 12 ) α 2 .
The elastic constants C 11 and C 12 are
C 11 = K + 2 3 ( C 11 C 12 ) ,
C 12 = K 1 3 ( C 11 C 12 ) .
From the quadratic fit in Figure A3, we get for nickel C 11 = 250.7 GPa , C 12 = 145.5 GPa , and for the hydride C 11 = 295.3 GPa , C 12 = 197.3 GPa .
Figure A3. The results for the orthorhombic distortion are shown. The atomistic data are fitted via a quadratic ansatz.
Figure A3. The results for the orthorhombic distortion are shown. The atomistic data are fitted via a quadratic ansatz.
Metals 08 00280 g0a3

Appendix C.4. Nonlinear Elastic Coefficients

To demonstrate the extraction of the elastic nonlinearities, we inspect the mechanical response of pure Ni under bulk compression or expansion, as shown in Figure A4.
For larger strains, we see deviations from the linear response, i.e., a quadratic energy–strain relationship. For tensile stresses the material appears effectively softer, whereas it is stiffer under compression. In the strain regime shown in the plot, the elastic response is well described by
f el = 9 2 K ϵ 2 1 + d 1 K tr ϵ + d 2 K ( tr ϵ ) 2 + d 3 K ( tr ϵ ) 3 + ,
in the spirit of Equation (8). The coefficients d i K are given in Table 2.
For pure shear during a monoclinic transformation, the energy is well described by a quadratic strain dependence, as shown before in Figure A2, and the same holds for the volume preserving orthorhombic transformation (see Figure A3). This demonstrates that the nonlinearities depend only on the trace of the strain tensor. We have performed monoclinic and orthorhombic transformations at different lattice units (i.e., bulk compression or expansion), to extract the coefficients in the spirit of Equation (8). This is exemplarily shown in Figure A5 for the monoclinic deformation of pure Ni. This allows to extract the elastic constants, here e.g., C 44 as function of the trace of the strain tensor. This dependence is shown in Figure A6. The resulting data for both phases are then described by Equation (8) with the coefficients in Table 2. The corresponding plot for the orthorhombic deformation is shown in Figure A7.
Figure A4. The free energy density as a function of the external strain ϵ for the pure Nickel. The quadratic ansatz is only true for small strains, whereas for higher strain we have to use a nonlinear function to fit the Monte Carlo data.
Figure A4. The free energy density as a function of the external strain ϵ for the pure Nickel. The quadratic ansatz is only true for small strains, whereas for higher strain we have to use a nonlinear function to fit the Monte Carlo data.
Metals 08 00280 g0a4
Figure A5. Monoclinic deformation of pure nickel for different bulk lattice constants. The fits (curves) are parabolic, and the data (points) are well described without further nonlinear corrections, in agreement with the statement that corrections depend only on tr ϵ . The latter is contained in the lattice constant dependence of the curvatures, leading to strain-dependent elastic constants C 44 ( tr ϵ ) . The energy of the bulk compression is subtracted in the plots.
Figure A5. Monoclinic deformation of pure nickel for different bulk lattice constants. The fits (curves) are parabolic, and the data (points) are well described without further nonlinear corrections, in agreement with the statement that corrections depend only on tr ϵ . The latter is contained in the lattice constant dependence of the curvatures, leading to strain-dependent elastic constants C 44 ( tr ϵ ) . The energy of the bulk compression is subtracted in the plots.
Metals 08 00280 g0a5
Figure A6. The elastic constant C 44 as function of the trace of the strain tr ϵ for Nickel and the hydride is shown. For tr ϵ = 0 , we obtain the values: C 44 , N i = 134 . 3 GPa and C 44 , H = 33 . 5 GPa .
Figure A6. The elastic constant C 44 as function of the trace of the strain tr ϵ for Nickel and the hydride is shown. For tr ϵ = 0 , we obtain the values: C 44 , N i = 134 . 3 GPa and C 44 , H = 33 . 5 GPa .
Metals 08 00280 g0a6
Figure A7. The difference C 11 C 12 depending on the trace of strain tr ϵ for Nickel and the hydride is shown.
Figure A7. The difference C 11 C 12 depending on the trace of strain tr ϵ for Nickel and the hydride is shown.
Metals 08 00280 g0a7

Appendix D. The Isotropic Model

We assume that the material is isotropic, and that a spherical nucleus of the hydride with radius R H forms within the dilute matrix. To simplify the geometry, we assume that the system is also spherical with radius R, thus the entire problem is rotational invariant. Hence, the displacement field depends only on the radius coordinate and has only a radial component. We note that we treat here only the linear elastic case, where corrections due to higher order terms in the strain as introduced in Equation (8) are not considered. Nevertheless, we note that the nonlinear case can also be treated in a similar way. However, for going beyond the one-dimensional case for an inclusion inside a cubic cell, we use finite element simulations instead.
In the outer dilute phase, the displacement is given by
u r d = a d r + b d / r 2 ,
with constants a d and b d . This ansatz satisfies the linear elastic bulk equations. In the hydride (inner phase), we use similarly
u r h = a h r .
The strains are given in spherical coordinates by
ϵ r r ( i ) = u r ( i ) , ϵ θ θ ( i ) = ϵ ϕ ϕ ( i ) = u r ( i ) / r ,
and all non-diagonal terms vanish. The stresses in the dilute phase are
σ i j d = λ d ϵ k k d δ i j + 2 G d ϵ i j d ,
and in the hydride
σ i j h = λ h ( ϵ k k h 3 Δ ε 0 ) δ i j + 2 G h ( ϵ i j h Δ ε 0 δ i j )
with the eigenstrain Δ ε 0 δ i j .
At the interface coherency and stress balance demands
u r d ( R H ) = u r h ( R H ) , σ r r d ( R H ) = σ r r h ( R H ) .
At the outer boundary, the fixed volume condition implies
u r d ( R ) = 0 .
These are altogether three conditions to determine the constants a d , b d and a h . The elastic energy densities (per unit volume) are
f el d = 1 2 λ d ϵ k k d 2 + G d ϵ i j d 2
and
f el h = 1 2 λ h ( ϵ k k h 3 Δ ε 0 ) 2 + G h ( ϵ i j h Δ ε 0 δ i j ) 2 .
The total free energy is
F el = 4 π 0 R H r 2 f el h d r + 4 π R H R r 2 f el d d r ,
which is a lengthy expression. Using the conversions
4 3 π R 3 = N Ni n 0 a 3 , c = R H R 3
with the lattice unit a Ni , we can calculate the elastic contribution to the chemical potential
μ el = n 0 N 0 f el c .
In the general case, it becomes
μ el = ( 3 a Ni 3 ( Δ ε 0 ) 2 ( 3 λ h + 2 G h ) ( c 2 ( 3 λ d + 2 G d ) × ( 3 λ d 3 λ h + 2 G d 2 G h ) + 4 G d ( 3 λ h + 4 G d + 2 G h ) + 2 c ( 3 λ d + 2 G d ) ( 3 λ h + 4 G d + 2 G h ) ) ) / [ 2 ( 3 λ h + 4 G d + c ( 3 λ d 3 λ h + 2 G d 2 G h ) + 2 G h ) 2 N 0 ] .
In particular for equal elastic constants ( λ d = λ h = λ , G h = G d = G )
μ el = a Ni 3 ( Δ ε 0 ) 2 ( 3 λ + 2 G ) ( 2 G + c ( 3 λ + 2 G ) ) ( λ + 2 G ) N 0 .

Appendix E. Interfacial Effects

The shape of the precipitate forming in the two-phase region is mainly determined by interfacial energy. For isotropic surface energy, a spherical inclusion minimizes the interfacial energy for small radii R H L , with L being the edge length of the enclosing box with periodic boundary conditions. However, for R H L , the formation of other interface contours is energetically favorable, since we use periodic boundary conditions in the atomistic and FEM simulations.
To find the different ranges and the limiting concentrations at which the forms change, we have to compute the minimum surface energy. The surface energy is defined as
F s = γ S
where S defines the size of the interfaces and γ the surface energy coefficient.
First, we look for the spherical inclusion of the hydride with the radius R H . In this case, the surface energy is
F s = γ S = γ 4 π R H 2 .
The concentration (volume fraction) is given as
c = 4 π R H 3 3 L 3 ,
so the surface energy is
F s ( c ) = γ L 2 ( 36 c 2 π ) 1 / 3 .
The surface energy for a hydride tube with radius r and length L is
F s = γ 2 π r L .
with the respective concentration
c = π r 2 L L 3
we obtain
F s ( c ) = 2 γ L 2 ( c π ) 1 / 2 .
In the case of slabs, the surface energy is independent of the concentration c:
F s = 2 γ L 2 .
The cases with c > 1 / 2 are analogous, with the role of the matrix and the hydride being exchanged. Figure A8 shows the surface energy depending on the concentration for the respective case. Up to c = 4 π / 81 0.16 , it should be a spherical inclusion of the hydride in the nickel matrix. For 4 π / 81 c 1 / π , the inclusion of the hydride should be a tube followed by plates up to c = 1 1 / π . The analytical results predict a symmetric behavior for concentrations c > 0 . 5 , so for concentrations 1 1 / π c 1 4 π / 81 there should be a tube and for c > 1 4 π / 81 we expect a nickel sphere.
Figure A8. Surface energy depending on the concentration for different geometric cases. The thick line segments indicate the lowest energy configurations.
Figure A8. Surface energy depending on the concentration for different geometric cases. The thick line segments indicate the lowest energy configurations.
Metals 08 00280 g0a8
With the calculated surface energy F s , we compute now the chemical potential with μ s = d F s / d N H . Then, we have to rewrite the surface energy in an expression which depends on the number of hydrogen atoms N H . For the first case, the hydride sphere inclusion, the radius can be expressed as
R H = a Ni 3 N H N 0 4 π 1 / 3 ,
so
F s ( N H ) = a Ni 2 γ π 1 / 3 6 N H N 0 2 / 3 .
The chemical potential then becomes
μ s ( N ) = 2 · 2 2 / 3 a Ni 2 γ π / 3 1 / 3 N 0 N H / N 0 1 / 3 ,
and finally with c = N H n 0 / ( N Ni N 0 )
μ s ( c ) = 2 · 2 2 / 3 a Ni 2 γ π / 3 1 / 3 N 0 c N Ni / n 0 1 / 3 .
Similarly, for tubes, we get
μ S = a Ni 2 π γ c N Ni / n 0 2 / 3 N 0 ,
for slabs
μ S = 0 ,
for inverse tubes
μ S = a Ni 2 π γ N 0 ( 1 c ) N Ni / n 0 2 / 3 ,
and for inverse spheres
μ S = 4 a Ni 2 π 1 / 3 γ N 0 6 ( 1 c ) N Ni / n 0 1 / 3 .

References

  1. Williams, R.O. The calculation of coherent phase equilibria. Calphad 1984, 8, 1–14. [Google Scholar] [CrossRef]
  2. Cahn, J.W.; Larché, F. A simple model for coherent equilibrium. Acta Metall. 1984, 32, 1915–1923. [Google Scholar] [CrossRef]
  3. Johnson, W.C.; Vorhees, P.W. Phase equilibrium in two-phase coherent solids. Metall. Trans. A 1987, 18, 1213–1228. [Google Scholar] [CrossRef]
  4. Von Pezold, J.; Lymperakis, L.; Neugebauer, J. Hydrogen-enhanced local plasticity at dilute bulk H concentrations: The role of H–H interactions and the formation of local hydrides. Acta Mater. 2011, 59, 2969–2980. [Google Scholar] [CrossRef]
  5. Haftbaradaran, H.; Song, J.; Curtin, W.A.; Gao, H. Continuum and atomistic models of strongly coupled diffusion, stress, and solute concentration. J. Power Sour. 2011, 196, 361–370. [Google Scholar] [CrossRef]
  6. Hüter, C.; Fu, S.; Finsterbusch, M.; Figgemeier, E.; Wells, L.; Spatschek, R. Electrode-Electrolyte Interface Stability in Solid State Electrolyte Systems: Influence of Coating Thickness under Varying Residual Stresses. AIMS Mater. Sci. 2016, 4, 867. [Google Scholar] [CrossRef]
  7. Markus, I.M.; Prill, M.; Dang, S.; Markus, T.; Spatschek, R.; Singheiser, L. High temperature investigation of electrochemical lithium insertion into Li4Ti5O12. Phys. Chem. Chem. Phys. 2016, 18, 31640. [Google Scholar] [CrossRef] [PubMed]
  8. 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]
  9. Shizuku, Y.; Yamamoto, S.; Fukai, Y. Phase diagram of the Ni-H system at high hydrogen pressures. J. Alloys Comp. 2002, 336, 159. [Google Scholar] [CrossRef]
  10. Fukai, Y.; Yamatomo, S.; Harada, S.; Kanazawa, M. The phase diagram of the Ni-H system revisited. J. Alloys Comp. 2004, 372, L4–L5. [Google Scholar] [CrossRef]
  11. Hüter, C.; Friák, M.; Weikamp, W.; 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]
  12. Sadigh, B.; Erhart, P. Calculation of excess free energies of precipitates via direct thermodynamic integration across phase boundaries. Phys. Rev. B 2012, 86, 134204. [Google Scholar] [CrossRef]
  13. Laks, D.B.; Ferreira, L.G.; Froyen, S.; Zunger, A. Efficient cluster expansion for substitutional systems. Phys. Rev. B 1992, 46, 12587. [Google Scholar] [CrossRef]
  14. Kurta, R.P.; Bugaev, V.N.; Ortiz, A.D. Long-Wavelength Elastic Interactions in Complex Crystals. Phys. Rev. Lett. 2010, 104, 085502. [Google Scholar] [CrossRef] [PubMed]
  15. Metropolis, N.; Rosenbluth, A.W.; Rosenbluth, M.N.; Teller, A.H. Equation of State Calculations by Fast Computing Machines. J. Chem. Phys 1953, 21, 1087–1092. [Google Scholar] [CrossRef]
  16. Plimpton, S.J. Fast Parallel Algorithms for Short-Range Molecular Dynamics. J. Comp. Phys. 1995, 117, 1–19. [Google Scholar] [CrossRef]
  17. LAMMPS Molecular Dynamics Simulator. Available online: http://lammps.sandia.gov (accessed on 17 April 2018).
  18. Angelo, J.E.; Moody, N.R.; Baskes, M.I. Trapping of hydrogen to lattice-defects in nickel. Model. Simul. Mater. Sci. Engr. 1995, 3, 289. [Google Scholar] [CrossRef]
  19. Baskes, M.I.; Sha, X.W.; Angelo, J.E.; Moody, N.R. Trapping of hydrogen to lattice defects in nickel. Model. Simul. Mater. Sci. Eng. 1997, 5, 651–652. [Google Scholar] [CrossRef]
  20. Haasen, P. Physical Metallurgy; Cambridge University Press: Cambridge, UK, 2003. [Google Scholar]
  21. Landau, L.D.; Lifshitz, E.M. Theory of Elasticity; Pergamon Press: Oxford, UK, 1987. [Google Scholar]
  22. Mura, T. Micromechanics of Defects in Solids; Springer: Berlin/Heidelberg, Germany, 1987. [Google Scholar]
  23. Marchenko, V.I. Theory of the equilibrium shape of crystals. Zh. Eksp. Teor. Fiz. (USSR) 1981, 81, 1141. [Google Scholar]
  24. Fischer, F.D.; Waitz, T.; Vollath, D.; Simha, N.K. On the role of surface energy and surface stress in phase-transforming nanoparticles. Progress Mater. Sci. 2008, 53, 481–527. [Google Scholar] [CrossRef]
  25. Cahn, J.W.; Larché, F. Surface Stress and Chemical Equilibrium of Small Crystals. II. Solid. Particles Embedded in a Solid Matrix. Acta Metall. 1982, 30, 51–56. [Google Scholar] [CrossRef]
  26. Armbruster, M.H. The Solubility of Hydrogen at Low Pressure in Iron, Nickel and Certain Steels at 400 to 600. J. Am. Chem. Soc. 1943, 65, 1043–1054. [Google Scholar] [CrossRef]
  27. Brener, E.A.; Marchenko, V.I.; Spatschek, R. Influence of strain on the kinetics of phase transitions in solids. Phys. Rev. E 2007, 75, 041604. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Selected atomistic simulation results for different hydrogen concentrations with fixed volume at T = 300 K . The Ni matrix is not shown, and the hydrogen atoms are visualized by the dots. The equilibrium shape of the precipitates in a periodic system depends on the H concentrations. First, hydride spheres appear at low concentrations, here for c = 0.22 (a); followed by a hydride tube (b) with c = 0.23 . For higher saturations (here c = 0.78 ), slabs are observed (c). Beyond c = 0.5 , the situation is inverted and Ni precipitates form inside the hydride matrix. Snapshot (d) shows a Ni tube for c = 0.85 . The last possible shape—a Ni sphere inside a hydride matrix—is not shown.
Figure 1. Selected atomistic simulation results for different hydrogen concentrations with fixed volume at T = 300 K . The Ni matrix is not shown, and the hydrogen atoms are visualized by the dots. The equilibrium shape of the precipitates in a periodic system depends on the H concentrations. First, hydride spheres appear at low concentrations, here for c = 0.22 (a); followed by a hydride tube (b) with c = 0.23 . For higher saturations (here c = 0.78 ), slabs are observed (c). Beyond c = 0.5 , the situation is inverted and Ni precipitates form inside the hydride matrix. Snapshot (d) shows a Ni tube for c = 0.85 . The last possible shape—a Ni sphere inside a hydride matrix—is not shown.
Metals 08 00280 g001
Figure 2. Chemical potential of hydrogen as function of H concentration at T = 300 K . The S-shaped curve is the single phase prediction, the horizontal dashed line the Maxwell construction. For them, elastic effects are not taken into account, in contrast to the Monte Carlo data, which is for fixed volume (lattice constant of pure Ni). For high concentrations, compressive elastic effects become dominant. Apart from the dilute regions c 1 and 1 c 1 , the Monte Carlo data correspond to two-phase configurations.
Figure 2. Chemical potential of hydrogen as function of H concentration at T = 300 K . The S-shaped curve is the single phase prediction, the horizontal dashed line the Maxwell construction. For them, elastic effects are not taken into account, in contrast to the Monte Carlo data, which is for fixed volume (lattice constant of pure Ni). For high concentrations, compressive elastic effects become dominant. Apart from the dilute regions c 1 and 1 c 1 , the Monte Carlo data correspond to two-phase configurations.
Metals 08 00280 g002
Figure 3. Energy as function of the H concentration for a homogeneous and pressure free system, to determine the solvation energy μ 0 and the H–H interaction. Notice that neither entropic effects are present at T = 0 , nor elastic or interfacial effects. The Monte Carlo data is averaged over several configurations.
Figure 3. Energy as function of the H concentration for a homogeneous and pressure free system, to determine the solvation energy μ 0 and the H–H interaction. Notice that neither entropic effects are present at T = 0 , nor elastic or interfacial effects. The Monte Carlo data is averaged over several configurations.
Metals 08 00280 g003
Figure 4. The chemical potential contribution μ H-H as function of the hydrogen concentration in a homogeneous system.
Figure 4. The chemical potential contribution μ H-H as function of the hydrogen concentration in a homogeneous system.
Metals 08 00280 g004
Figure 5. Dilute limit of the chemical potential, which is dominated by the configurational entropy. The data are for T = 300 K and for fixed volume with the equilibrium lattice constant of pure Ni. The chemical potential is (apart from the offset μ 0 ) dominated by the configurational contribution μ c , whereas elastic and H–H interaction terms, μ el and μ H-H are negligible there. The accuracy of the continuum description in comparison to the Monte Carlo data are finally about 20 meV per atom in the entire concentration range 0 < c < 1 .
Figure 5. Dilute limit of the chemical potential, which is dominated by the configurational entropy. The data are for T = 300 K and for fixed volume with the equilibrium lattice constant of pure Ni. The chemical potential is (apart from the offset μ 0 ) dominated by the configurational contribution μ c , whereas elastic and H–H interaction terms, μ el and μ H-H are negligible there. The accuracy of the continuum description in comparison to the Monte Carlo data are finally about 20 meV per atom in the entire concentration range 0 < c < 1 .
Metals 08 00280 g005
Figure 6. Chemical potential in the single phase state for different temperatures. The curves are the analytical descriptions based on the expression μ 0 + μ H-H + μ c using the parameters of the H–H interaction and the solvation energy. The points are the data from the grand canonical Monte Carlo simulations of homogeneous, single-phase states. At low temperatures, the solubility is very low, and therefore the concentration range of the single phase equilibrium states limited to the dilute regions.
Figure 6. Chemical potential in the single phase state for different temperatures. The curves are the analytical descriptions based on the expression μ 0 + μ H-H + μ c using the parameters of the H–H interaction and the solvation energy. The points are the data from the grand canonical Monte Carlo simulations of homogeneous, single-phase states. At low temperatures, the solubility is very low, and therefore the concentration range of the single phase equilibrium states limited to the dilute regions.
Metals 08 00280 g006
Figure 7. Lattice constant as function of the hydrogen concentration. The Monte Carlo results are obtained from homogeneous free volume simulations, averaged over several configurations. The fit describes the nonlinear concentration dependence, in contrast to a linear interpolation (Vegard’s law).
Figure 7. Lattice constant as function of the hydrogen concentration. The Monte Carlo results are obtained from homogeneous free volume simulations, averaged over several configurations. The fit describes the nonlinear concentration dependence, in contrast to a linear interpolation (Vegard’s law).
Metals 08 00280 g007
Figure 8. Elastic contribution to the chemical potential in various approximations: The solid curve expresses the isotropic prediction by Equation (6), using equal elastic constants for both phases (we use the values of pure Ni). The long-dashed curve is the same isotropic model, but this time using the different elastic constants of both phases. This prediction is almost identical to a numerically calculated expression based on cubic elasticity, which assumes a spherical inclusion in a cubic box with periodic boundary conditions (short-dashed curve). In contrast to all these linear elastic approximations, the dotted line also considers elastic nonlinearities. Its prediction is very close to the Monte Carlo data.
Figure 8. Elastic contribution to the chemical potential in various approximations: The solid curve expresses the isotropic prediction by Equation (6), using equal elastic constants for both phases (we use the values of pure Ni). The long-dashed curve is the same isotropic model, but this time using the different elastic constants of both phases. This prediction is almost identical to a numerically calculated expression based on cubic elasticity, which assumes a spherical inclusion in a cubic box with periodic boundary conditions (short-dashed curve). In contrast to all these linear elastic approximations, the dotted line also considers elastic nonlinearities. Its prediction is very close to the Monte Carlo data.
Metals 08 00280 g008
Figure 9. Energy per Ni atom as function of the hydrogen concentration. The atomistic data are calculated at T = 0 K without Monte Carlo steps but atomic relaxation for fixed volume. A spherical precipitate is placed in the center to the simulation cell (with periodic boundary conditions). The irrelevant integration constant E 0 for pure Ni is subtracted, as well as the trivial term μ 0 , thus only the H–H interaction, elasticity and surface contributions remain. The shape of the precipitate is not important for the elastic energy; we used both Ni-H and Ni inclusions in the atomistic simulation, visualized by open and filled spheres.
Figure 9. Energy per Ni atom as function of the hydrogen concentration. The atomistic data are calculated at T = 0 K without Monte Carlo steps but atomic relaxation for fixed volume. A spherical precipitate is placed in the center to the simulation cell (with periodic boundary conditions). The irrelevant integration constant E 0 for pure Ni is subtracted, as well as the trivial term μ 0 , thus only the H–H interaction, elasticity and surface contributions remain. The shape of the precipitate is not important for the elastic energy; we used both Ni-H and Ni inclusions in the atomistic simulation, visualized by open and filled spheres.
Metals 08 00280 g009
Figure 10. The chemical potential of hydrogen at T = 300 K . The dots are the data from the Monte Carlo simulations. The chemical potential including nonlinear elastic effects is very close to the Monte Carlo data, and, together with interfacial effects, the agreement is excellent. The contribution from the interfacial terms is largest for low concentrations, as there the precipitates are small and therefore have the largest surface-to-volume ratio. The dotted curves near c = 0 and c = 1 are the analytical predictions for the single phase dilute limits. The H 2 partial pressure is based on an ideal gas model, which is not realistic for higher pressures and only serves for illustrational purposes.
Figure 10. The chemical potential of hydrogen at T = 300 K . The dots are the data from the Monte Carlo simulations. The chemical potential including nonlinear elastic effects is very close to the Monte Carlo data, and, together with interfacial effects, the agreement is excellent. The contribution from the interfacial terms is largest for low concentrations, as there the precipitates are small and therefore have the largest surface-to-volume ratio. The dotted curves near c = 0 and c = 1 are the analytical predictions for the single phase dilute limits. The H 2 partial pressure is based on an ideal gas model, which is not realistic for higher pressures and only serves for illustrational purposes.
Metals 08 00280 g010
Figure 11. Comparison of the chemical potential for fixed volume and in a stress-free situation, P = 0 , T = 300 K . The volume constraint leads to a positive slope of the chemical potential in the two-phase region, thus stabilizing phase separation of given chemical potential. In contrast, two phase states are unstable for given μ H , as shown here using the continuum prediction.
Figure 11. Comparison of the chemical potential for fixed volume and in a stress-free situation, P = 0 , T = 300 K . The volume constraint leads to a positive slope of the chemical potential in the two-phase region, thus stabilizing phase separation of given chemical potential. In contrast, two phase states are unstable for given μ H , as shown here using the continuum prediction.
Metals 08 00280 g011
Figure 12. Evolution of the hydrogen concentration as function of the number of Monte Carlo steps for T = 300 K . Grand-canconical simulations with given chemical potential μ H are used. (Top) Fixed volume using the equilibrium lattice constant of Ni. Due to the positive slope of the chemical potential in the two-phase region (solid curve and data points in Figure 11) stable two phase states with arbitrary average concentrations can be obtained in thermodynamic equilibrium. (Bottom) For free volume, the miscibility gap is inaccessible for the given chemical potential in the grand canonical simulations, and instead the equilibrium state is almost free of hydrogen or almost fully saturated, as the solubilities are low at low temperatures (see phase diagram Figure 13). The reason for this behavior (which is accompanied by a hysteresis) is the negative slope of the chemical potential in the two-phase region (see Figure 11), which implies that these configurations are unstable. All simulations are started as a hydrogen free system, and formation of the hydride therefore demands to overcome a nucleation barrier.
Figure 12. Evolution of the hydrogen concentration as function of the number of Monte Carlo steps for T = 300 K . Grand-canconical simulations with given chemical potential μ H are used. (Top) Fixed volume using the equilibrium lattice constant of Ni. Due to the positive slope of the chemical potential in the two-phase region (solid curve and data points in Figure 11) stable two phase states with arbitrary average concentrations can be obtained in thermodynamic equilibrium. (Bottom) For free volume, the miscibility gap is inaccessible for the given chemical potential in the grand canonical simulations, and instead the equilibrium state is almost free of hydrogen or almost fully saturated, as the solubilities are low at low temperatures (see phase diagram Figure 13). The reason for this behavior (which is accompanied by a hysteresis) is the negative slope of the chemical potential in the two-phase region (see Figure 11), which implies that these configurations are unstable. All simulations are started as a hydrogen free system, and formation of the hydride therefore demands to overcome a nucleation barrier.
Metals 08 00280 g012
Figure 13. Phase diagram of the Ni-H systems, based on parameters extracted from Monte Carlo simulations. The solid line is the binodal without consideration of elastic effects, whereas the dash-dotted line is the same with elastic effects for fixed volume.
Figure 13. Phase diagram of the Ni-H systems, based on parameters extracted from Monte Carlo simulations. The solid line is the binodal without consideration of elastic effects, whereas the dash-dotted line is the same with elastic effects for fixed volume.
Metals 08 00280 g013
Figure 14. Sketch of the hydride phase (light grey) forming as a slab in the nickel matrix with fixed volume. Due to translational invariance, we can assume the hydride phase to be located in the left part of the cube, starting at x = 0 .
Figure 14. Sketch of the hydride phase (light grey) forming as a slab in the nickel matrix with fixed volume. Due to translational invariance, we can assume the hydride phase to be located in the left part of the cube, starting at x = 0 .
Metals 08 00280 g014
Figure 15. The phase diagram depending on the eigenstrain is shown. The binodal and elastic binodal curves are the same as in Figure 13. The other two curves show the tremendous influence of the concentration dependence of the eigenstrain or lattice constant. If, instead of the concave dependence ε 0 ( c ) < 0 , a linear relationship is assumed (Vegard’s law), the solubility of hydrogen is drastically enhanced. In contrast, the concentration dependence of the elastic constants plays only a minor role: The case of different elastic constants anticipates a linear interpolation of the elastic constants between the values for Ni and Ni-H, whereas, for equal constants, both phases are assumed to have the same elastic constants (using values for Ni).
Figure 15. The phase diagram depending on the eigenstrain is shown. The binodal and elastic binodal curves are the same as in Figure 13. The other two curves show the tremendous influence of the concentration dependence of the eigenstrain or lattice constant. If, instead of the concave dependence ε 0 ( c ) < 0 , a linear relationship is assumed (Vegard’s law), the solubility of hydrogen is drastically enhanced. In contrast, the concentration dependence of the elastic constants plays only a minor role: The case of different elastic constants anticipates a linear interpolation of the elastic constants between the values for Ni and Ni-H, whereas, for equal constants, both phases are assumed to have the same elastic constants (using values for Ni).
Metals 08 00280 g015
Table 1. Elastic constants of pure nickel and the fully saturated hydride, as obtained from the EAM potential (see Appendix C for details). We note that the relaxed configuration of Ni is used as reference state.
Table 1. Elastic constants of pure nickel and the fully saturated hydride, as obtained from the EAM potential (see Appendix C for details). We note that the relaxed configuration of Ni is used as reference state.
C ij NiNi-H
C 11 250.7 GPa295.3 GPa
C 12 145.5 GPa197.3 GPa
C 44 134.3 GPa33.5 GPa
Table 2. Coefficients for the nonlinear elastic corrections. If higher order coefficients are missing, the expansion is truncated already at lower order.
Table 2. Coefficients for the nonlinear elastic corrections. If higher order coefficients are missing, the expansion is truncated already at lower order.
Elastic ConstantMaterial d 1 d 2 d 3 d 4 d 5
C 44 Ni−3.25−7.04
C 11 C 12 Ni−2.01−12.51
KNi−1.000.53−4.492.5750.03
C 44 Ni-H−22.80−47.20
C 11 C 12 Ni-H−9.77−14.69
KNi−2.253.3711.513.55−51.00

Share and Cite

MDPI and ACS Style

Korbmacher, D.; Von Pezold, J.; Brinckmann, S.; Neugebauer, J.; Hüter, C.; Spatschek, R. Modeling of Phase Equilibria in Ni-H: Bridging the Atomistic with the Continuum Scale. Metals 2018, 8, 280. https://doi.org/10.3390/met8040280

AMA Style

Korbmacher D, Von Pezold J, Brinckmann S, Neugebauer J, Hüter C, Spatschek R. Modeling of Phase Equilibria in Ni-H: Bridging the Atomistic with the Continuum Scale. Metals. 2018; 8(4):280. https://doi.org/10.3390/met8040280

Chicago/Turabian Style

Korbmacher, Dominique, Johann Von Pezold, Steffen Brinckmann, Jörg Neugebauer, Claas Hüter, and Robert Spatschek. 2018. "Modeling of Phase Equilibria in Ni-H: Bridging the Atomistic with the Continuum Scale" Metals 8, no. 4: 280. https://doi.org/10.3390/met8040280

APA Style

Korbmacher, D., Von Pezold, J., Brinckmann, S., Neugebauer, J., Hüter, C., & Spatschek, R. (2018). Modeling of Phase Equilibria in Ni-H: Bridging the Atomistic with the Continuum Scale. Metals, 8(4), 280. https://doi.org/10.3390/met8040280

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