covid
Buscar en
Journal of Applied Research and Technology. JART
Toda la web
Inicio Journal of Applied Research and Technology. JART A first principles calculation of Ni-16Cr and Ni-16Mo alloys
Información de la revista
Vol. 15. Núm. 1.
Páginas 78-82 (febrero 2017)
Compartir
Compartir
Descargar PDF
Más opciones de artículo
Visitas
2686
Vol. 15. Núm. 1.
Páginas 78-82 (febrero 2017)
Short communication
Open Access
A first principles calculation of Ni-16Cr and Ni-16Mo alloys
Visitas
2686
Ashish Pathaka, Kumud Kant Mehtab, Ashok Kumar Singha,
Autor para correspondencia
Singh_ashok3@rediffmail.com

Corresponding author.
a Defence Metallurgical Research Laboratory, Hyderabad 500058, India
b RDAQA (M), c/o DRDO, Hyderabad 500058, India
Este artículo ha recibido

Under a Creative Commons license
Información del artículo
Resumen
Texto completo
Bibliografía
Descargar PDF
Estadísticas
Figuras (1)
Tablas (3)
Table 1. Density, equilibrium lattice constants and formation energy/atom for Ni-16Cr and Ni-16Mo alloys.
Table 2. Elastic constants and anisotropy factor (A) for Ni-16Cr and Ni-16Mo alloys.
Table 3. Bulk modulus (B), Shear modulus (G) and Young's modulus (E) and Poisson's ratio (υ) for Ni-16Cr and Ni-16Mo alloys.
Mostrar másMostrar menos
Abstract

The present paper describes the structural stability and mechanical properties of the γ phase for Ni-16Cr and Ni-16Mo alloys using first principles density functional theory (DFT) within generalized gradient approximation (GGA). The equilibrium lattice constant values of these alloys are in good agreement with the experimental data obtained by X-ray diffraction technique. The formation energy per atom values suggests that the γ phase of both alloys is stable. The Ni-16Mo alloy is energetically more stable than Ni-16Cr. Both alloys satisfy the Born stability criteria in terms of elastic constants and are associated with ductile behaviour based on shear to bulk modulus ratios. Both alloys show an anisotropic behaviour. The Ni-16Mo alloy exhibits higher anisotropy than the Ni-16Cr one.

Keywords:
Alloys
First-principles
Crystal structure
Physical properties
Texto completo
1Introduction

Ni-based solid solution alloys such as Inconel 617, Inconel 625, Inconel 686, Hastelloy S, Hastelloy C, Hastelloy C-4, Hastelloy C-22, Hastelloy C-276, Nicrofer 5923, among others, consist of Cr and Mo as major alloying elements (Davis, 2000). The industrial applications of these alloys are largely in wrought forms such as plates, sheets, strips, extruded tubes, wires and flanges due to their excellent hot and/or cold deforming capability with moderate to high work-hardening rate and weldability. The wrought forms of these alloys are therefore produced by thermo mechanical processing of the cast ingot. The microstructures and mechanical properties of Ni-16Cr and Ni-16Mowt. % (Ni-18Cr and Ni-8Moatom %) alloys have been investigated recently (Mehta, Mukhopadhyay, Mandal, & Singh, 2015). The properties of both alloys are quite different in hot rolled and annealed condition. The contents of Cr and Mo lie in the range of 10–20 (wt. %) in all the alloys mentioned above apart from other alloying elements. The limits of solubility of Cr and Mo in Ni are 35–40wt. % and about 20wt. %, respectively (Takeuchi & Inoue, 2005). In addition, the size difference between the Ni-Mo is fairly large in comparison to that of the Ni-Cr in view of the fact that the atomic radii (Å) of Ni, Cr and Mo are 1.25, 1.28 and 1.36, respectively (Delehouzee & Deruyttere, 1967). A linear relation between flow stress and lattice parameter change has been observed for any single solute element (Cr, Mo, W, Fe, Co and Cu) in nickel (Tawancy & Al-Hadhrami, 2012). The Ni-based solid solution alloys having around 16wt. % of both the Cr and Mo display marked improvements in thermal stability, fatigue life, high temperature strength, low expansion characteristics, oxidation and corrosion resistance (Tawancy & Al-Hadhrami, 2012). In addition, most of modelling work of solid solution Ni-based alloys has been carried out around this composition (Mishima, Ochiai, Hamao, Yodogawa, & Suzuki, 1986; Roth, Davis, & Thomson, 1997).

An understanding of these binary alloys in terms of their stability, lattice constants and mechanical properties using first principles calculation is therefore quite important. This can in fact be utilized to tailor the compositions and properties of multi-component solid solution alloys with optimized properties. The present paper is thus concerned with a first principles study of the Ni-16Cr and Ni-16Mo alloys.

2Experimental details

The alloy pancakes of 600g with nominal compositions Ni-16Cr and Ni-16Mo (wt. %) were prepared using non-consumable vacuum arc melting. The melting was repeated four times to ensure chemical homogeneity. The average analyzed chemical compositions of both alloys obtained by wet chemical and electron probe microanalyzer (EPMA-SX 100, CAMECA, France) techniques are Ni-16.4Cr and Ni-16.3Mo (wt. %). The density of both alloys has been obtained by Archimedes principle. Pancakes were unidirectionally hot rolled (HR) up to 70% reduction at 1150°C and then air cooled (AC). The materials were deformed 10% in each pass and reduction direction was kept strictly unidirectional. The samples were placed back into furnace for 15min after each pass to re-attain the rolling temperature. Hot rolled materials were then annealed at 1000°C and 1100°C for Ni-18Cr and Ni-8Mo alloys, respectively for 5min in vacuum furnace and subsequently furnace cooled in vacuum.

Microstructural characterization of both alloys in hot rolled (3mm sheet) and annealed condition was done using optical microscope (OM). The specimens were prepared following standard metallographic technique used for Ni-based alloys and etched with Kallings reagent (5gm CuCl2 in 40ml HCl and 60ml methanol solution). X-ray diffraction (XRD) study of hot rolled and annealed materials was performed using a Philips 3020 diffractometer with CuKα radiation equipped with graphite monochrometer.

3Computational details

A first principles ultrasoft (US) pseudo-potential method in the frame work of the density functional theory (DFT) has been utilized with Quantum Espresso code in present calculation (Giannozzi et al., 2009; Hohenberg & Kohn, 1964; Kohn & Sham, 1965; Vanderbilt, 1990). Perdew–Burke–Ernzerhof formulation within generalized gradient approximation (GGA) has been considered for exchange-correlation effects (Perdew, Burke, & Ernzerhof, 1996). Both the plane wave cut-off energy and cut-off density of 40Ry and 400Ry, respectively have been chosen in all calculations after the verification of the convergence for plane wave. The wave functions between the real and reciprocal lattices have been converted using Fast Fourier transform algorithm (Goedecker, 1997). Consequently, conjugate gradient algorithm has been employed within the frame work of self-consistency (Gonze, 1996; Payne, Teter, Allan, Arias, & Joannopoulos, 1992). The Monkhorst–Pack scheme has been utilized for the integration over the Brillouin zone (BZ) (Monkhorst & Pack, 1976). The optimum value used for the present calculations is 10×10×10 after ensuring the convergence. The differences between energies or forces in two consecutive steps are less than 1.36×10−4meV or 0.36eV/Å, respectively have been presumed for convergence.

Both the formation energy per atom and the lattice constant ‘a’ of Ni-16Cr and Ni-16Mo have been computed using equivalent atom % (Ni-18Cr and Ni-8Mo) and given in Table 1. The equilibrium lattice constant has been obtained by varying lattice constant and calculating the energy corresponding to each lattice constant to get the E vs. a curve. The lattice constant corresponding to minimum energy in the curve is regarded as an equilibrium lattice constant. The energy of formation per atom has been defined as

where ENiyM is the total system energy of the Ni-16Cr and Ni-16Mo alloys having ‘x’ Ni-atoms, ‘y’ M (Cr/Mo), ENig and EMg are the total energy per atom in their ground states for Ni and M (Cr/Mo) atoms, respectively. The (x+y) denotes the total number of atoms considered in the unit cell that are 16 and 12 for the Ni matrix of Ni-18Cr and Ni-8Mo alloys, respectively.

Table 1.

Density, equilibrium lattice constants and formation energy/atom for Ni-16Cr and Ni-16Mo alloys.

Alloys  Density (gm/cc)  Lattice constants (Å)Formation energy/atom (eV) 
  Experimental  Experimental  Present study   
Ni-16Cr  8.620  3.5366  3.5040  −0.0060 
Ni-16Mo  9.098  3.5690  3.5580  −0.0396 

The ElaStic code has been utilized to calculate the elastic constant (Golesorkhtabar, Pavone, Spitaler, Puschnig, & Draxl, 2013). This code is able to calculate the full second-order elastic stiffness tensors for any crystal structure from ab initio total-energy and/or stress calculations.

4Results and discussion

The density of both alloys is given in Table 1. This reflects the density values of Ni-16Cr and Ni-16Mo are smaller and larger than that of pure Ni (8.9g/cc), respectively. It is to be noted that the density values of pure Cr and Mo are 7.19 and 10.28g/cc, respectively. It reflects that the density of the present alloys depends mainly on the density of individual alloying elements and their respective weight fractions.

The optical microstructures of Ni-16Cr and Ni-16Mo are shown in Figure 1. The microstructure of the Ni-16Cr alloy in hot rolled and annealed condition reveals equiaxed grains with large variation in grain size (Fig. 1a). On the other hand, the Ni-16Mo alloy display nearly uniform microstructure with annealing twins within equiaxed grains (Fig. 1b). XRD patterns show the presence of γ phase (Ni matrix, Fm 3¯ m) only in both alloys (Fig. 1c). The lattice parameters of the γ phase of both alloys obtained from XRD patterns are given in Table 1. The lattice parameter of the Ni-16Mo alloy is larger than that of the Ni-16Cr one. This observation is not surprising since the atomic radii of Ni, Cr and Mo are 1.25, 1.28 and 1.36Å, respectively.

Fig. 1.

Optical microstructures of hot rolled and annealed alloys: (a) Ni-16Cr, (b) Ni-16Mo and (c) corresponding XRD patterns.

(0.33MB).

The lattice parameter values of the γ phase of both alloys have also been calculated using first principles calculation and given in Table 1. These are in good agreement and follows similar trend with lattice parameter values obtained by XRD technique.

The formation energy per atom values of Ni-16Cr and Ni-16Mo alloys obtained by first principle calculation are given in Table 1. The formation energy per atom of both alloys are negative indicating that the γ phase of the present alloys is stable. This also reflects that the Ni-16Mo alloy is more stable than the Ni-16Cr one since the formation energy per atom of former is lower than the later. Incidentally, the stacking fault energy (SFE) of Ni-16Mo is less than two third of the Ni-16Cr (Mehta et al., 2015). This follows similar trend to those of formation energy per atom of both alloys.

The elastic properties of single crystal of the γ phase of both alloys are described by C11, C12 and C44 and given in Table 2. The elastic constant values of Ni-16Cr are higher than those of the Ni-16Mo alloy. A comparison of the elastic constant values of pure Ni (C11=261, C12=151 and C44=132GPa) and present alloys indicates that the equal amount of Cr addition is more effective than Mo to increase the compliances. It also appears that this effect is opposite to the stability of the γ phase (in terms of formation energy per atom) in these alloys and extent of solid solution of both elements in Ni.

Table 2.

Elastic constants and anisotropy factor (A) for Ni-16Cr and Ni-16Mo alloys.

Alloys  C11 (GPa)  C12 (GPa)  C44 (GPa)  A 
Ni-16Cr  310  178  141  6.994 
Ni-16Mo  270  175  116  9.522 

The nature of metallic bonding of both alloys can be predicted based on Cauchy pressures (Olijnyk & Jephcoat, 2000). This is defined as C12–C44<0. The negative and positive Cauchy pressures are sign of more directional and metallic bonding, respectively. The calculated values of Cauchy pressures for the γ phase of both experimental alloys are positive. This clearly points towards the presence of predominant metallic bonds in these alloys.

The values of elastic constants can be used to predict the stability of phases. The Born stability criterions are given in Eqs. (2)–(4) (Born, 1940; Fedorov, 1968). The corresponding calculated values are given in parenthesis.

C44>0 (Ni-16Cr: 141; Ni-16Mo: 116)
C11–¿C12¿>0 (Ni-16Cr:132; Ni-16Mo:95)and
C11+2C12>0 (Ni-16Cr: 666; Ni-16Mo:620)

The values given in parenthesis of Eqs. (2)–(4) suggest that the γ phase of these alloys is mechanically stable. This is also in agreement with available experimental binary phase diagrams wherein the γ phase is stable at Cr and Mo contents of present alloys (Massalski, Okamoto, Subramanian, & Kacprzak, 2001).

The anisotropy factor (A) for the γ phase of these alloys is defined below and corresponding calculated values are given in Table 2.

For cubic materials

The values of A should be 1 for isotropic materials and away from 1 point towards the extent of anisotropy in elastic constants. The values of A obtained in present study clearly indicate that the single crystal of both alloys exhibit considerable amount of anisotropy. However, the extent of anisotropy present in alloy Ni-16Mo is higher than that of the Ni-16Cr.

The effective elastic modulus of isotropic polycrystalline cubic materials can be evaluated from the elastic constants by following two approximations namely, the Voigt and Reuss (Reuss & Angnew, 1929; Voigt, 1928). These approximations provide information about the upper and lower limits that are

where B is the bulk modulus while GV and GR are shear modulus values obtained by Voigt and Reuss approximations, respectively.

The average value of these two estimates, as mentioned above, is given by Hill approximation (Hill, 1952). The Voigt–Reuss–Hill (VRH) average values are given by

where BH, GH, E and ν are the bulk modulus, shear modulus, Young's modulus and Poisson's ratio, respectively.

The calculated values of B, G and E for both alloys are given in Table 3. The Ni-16Cr alloy exhibits higher values of modulus than the Ni-16Mo one. A comparison of modulus values of pure Ni (B=188, G=101 and E=257GPa) with the present alloys indicates that the addition of Cr increases these values. On the other hand, an addition of Mo increases B and decreases G and E values. Mehta et al. have recently investigated the mechanical properties of these alloys along three sample directions, i.e., longitudinal (L or 0°), 45° (specimen axis at 45° to the rolling direction) and transverse (T or 90°) directions (Mehta et al., 2015). The Young's modulus values of these alloys have been calculated from engineering stress-strain curves along all the three directions. The average values of the Young's modulus of the Ni-16Cr and Ni-16Mo alloys are 242 and 216GPa, respectively. These are in good agreement with Young's modulus values obtained in the present study by first principles calculation. The Poisson's ratio of alloy Ni-16Mo is larger than that of the Ni-16Cr. These values are also larger than the Poisson's ratio of pure Ni (0.27).

Table 3.

Bulk modulus (B), Shear modulus (G) and Young's modulus (E) and Poisson's ratio (υ) for Ni-16Cr and Ni-16Mo alloys.

Alloys  BV  BR  B  GV  GR  G  E  G/B  υ 
  GPa   
Ni-16Cr  223  223  223  111  96  104  270  0.466  0.30 
Ni-16Mo  202  202  202  90  74  82  217  0.406  0.32 

The values of shear and bulk modulus can also be utilized to predict the brittle and ductile behaviour of materials (Pugh, 1954). This can be envisaged by taking the ratio of G and B. The ratio (G/B)>0.57 is associated with brittle otherwise related with ductile behaviour. The G/B ratios of both alloys are reasonably lower than the 0.57. This indicates that these alloys are ductile. These are in agreement with the experimental results of Mehta et al. (2015). They have recently investigated the orientation dependent flow behaviour of these two binary alloys and reported high ductility in tensile tested samples.

5Conclusions

The main conclusions are:

  • 1.

    Based on formation energy per atom values and Born stability criteria in terms of elastic constants, the γ phase of both alloys is stable.

  • 2.

    The equilibrium lattice constant values of these alloys are in good agreement with the experimental data.

  • 3.

    The Ni-16Mo alloy is energetically more stable than the Ni-16Cr one while the latter one is mechanically more stable than the former one.

  • 4.

    Both alloys show anisotropic and ductile behaviour. The Ni-16Mo alloy exhibits higher anisotropy than the Ni-16Cr one.

Conflict of interest

The authors have no conflicts of interest to declare.

Acknowledgements

The authors are grateful to the Ministry of Defence of Government of India for financial support. The authors are indebted to Director DMRL Hyderabad for his encouragement. They also thank Director ANURAG, Hyderabad for the provision of computational facilities and Dr. R. Sankarasubramanian and Shri A. Mondal for their kind support.

References
[Born, 1940]
M. Born.
On the stability of crystal lattices. I.
Mathematical Proceedings of the Cambridge Philosophical Society (Vol. 36, No. 02), pp. 160-172
[Davis, 2000]
J.R. Davis.
ASM Specialty handbook: Nickel, cobalt, and their alloys.
ASM International, Member/Customer Service Center, (2000), pp. 442
[Delehouzee and Deruyttere, 1967]
L. Delehouzee, A. Deruyttere.
The stacking fault density in solid solutions based on copper, silver, nickel, aluminium and lead.
Acta Metallurgica, 15 (1967), pp. 727-734
[Fedorov, 1968]
F.I. Fedorov.
Theory of elastic waves in crystals.
Plenum, (1968),
[Giannozzi et al., 2009]
P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, C. Cavazzoni, …, A. Dal Corso.
QUANTUM ESPRESSO: A modular and open-source software project for quantum simulations of materials.
Journal of Physics: Condensed Matter, 21 (2009), pp. 395502
[Goedecker, 1997]
S. Goedecker.
Fast radix 2, 3, 4, and 5 Kernels for fast Fourier transformations on computers with overlapping multiply—Add instructions.
SIAM Journal on Scientific Computing, 18 (1997), pp. 1605-1611
[Golesorkhtabar et al., 2013]
R. Golesorkhtabar, P. Pavone, J. Spitaler, P. Puschnig, C. Draxl.
ElaStic: A tool for calculating second-order elastic constants from first principles.
Computer Physics Communications, 184 (2013), pp. 1861-1873
[Gonze, 1996]
X. Gonze.
Towards a potential-based conjugate gradient algorithm for order-N self-consistent total energy calculations.
Physical Review B, 54 (1996), pp. 4383
[Hill, 1952]
R. Hill.
The elastic behaviour of a crystalline aggregate.
Proceedings of the Physical Society. Section A, 65 (1952), pp. 349
[Hohenberg and Kohn, 1964]
P. Hohenberg, W. Kohn.
Inhomogeneous electron gas.
Physical Review, 136 (1964), pp. B864
[Kohn and Sham, 1965]
W. Kohn, L.J. Sham.
Self-consistent equations including exchange and correlation effects.
Physical Review, 140 (1965), pp. A1133
[Massalski et al., 2001]
T.B. Massalski, H. Okamoto, P. Subramanian, L. Kacprzak.
ASM International, (2001),
[Mehta et al., 2015]
K.K. Mehta, P. Mukhopadhyay, R.K. Mandal, A.K. Singh.
Microstructure, texture, and orientation-dependent flow behavior of binary Ni-16Cr and Ni-16Mo solid solution alloys.
Metallurgical and Materials Transactions A, 46 (2015), pp. 3656-3669
[Mishima et al., 1986]
Y. Mishima, S. Ochiai, N. Hamao, M. Yodogawa, T. Suzuki.
Solid solution hardening of nickel—Role of transition metal and B-subgroup solutes.
Transactions of the Japan Institute of Metals, 27 (1986), pp. 656-664
[Monkhorst and Pack, 1976]
H.J. Monkhorst, J.D. Pack.
Special points for Brillouin-zone integrations.
Physical Review B, 13 (1976), pp. 5188-5192
[Olijnyk and Jephcoat, 2000]
H. Olijnyk, A.P. Jephcoat.
The E2g phonon and the elastic constant C44 in hexagonal van der Waals bonded solids.
Journal of Physics: Condensed Matter, 12 (2000), pp. 10423
[Payne et al., 1992]
M.C. Payne, M.P. Teter, D.C. Allan, T.A. Arias, J.D. Joannopoulos.
Iterative minimization techniques for ab initio total-energy calculations: Molecular dynamics and conjugate gradients.
Reviews of Modern Physics, 64 (1992), pp. 1045
[Perdew et al., 1996]
J.P. Perdew, K. Burke, M. Ernzerhof.
Generalized gradient approximation made simple.
Physical Review Letters, 77 (1996), pp. 3865
[Pugh, 1954]
S.F. Pugh.
XCII. Relations between the elastic moduli and the plastic properties of polycrystalline pure metals.
The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 45 (1954), pp. 823-843
[Reuss and Angnew, 1929]
A. Reuss, Z. Angnew.
A calculation of the bulk modulus of polycrystalline materials.
Math Methods, 9 (1929), pp. 55
[Roth et al., 1997]
H.A. Roth, C.L. Davis, R.C. Thomson.
Modeling solid solution strengthening in nickel alloys.
Metallurgical and Materials Transactions A, 28 (1997), pp. 1329-1335
[Takeuchi and Inoue, 2005]
A. Takeuchi, A. Inoue.
Classification of bulk metallic glasses by atomic size difference, heat of mixing and period of constituent elements and its application to characterization of the main alloying element.
Materials Transactions, 46 (2005), pp. 2817-2829
[Tawancy and Al-Hadhrami, 2012]
H.M. Tawancy, L.M. Al-Hadhrami.
A nanocrystalline Ni2 (Cr, Mo) intermetallic with potentially useful combination of properties for gas turbine seal ring applications.
Journal of Materials Engineering and Performance, 21 (2012), pp. 1374-1379
[Vanderbilt, 1990]
D. Vanderbilt.
Soft self-consistent pseudopotentials in a generalized Eigenvalue formalism.
Physical Review B, 41 (1990), pp. 7892
[Voigt, 1928]
W. Voigt.
Lehrbuch der Kristallphysik (mit Ausschluss der Kristalloptik).
1st ed., B.G. Teubner, (1928),

Peer Review under the responsibility of Universidad Nacional Autónoma de México.

Copyright © 2016. Universidad Nacional Autónoma de México, Centro de Ciencias Aplicadas y Desarrollo Tecnológico
Descargar PDF
Opciones de artículo