Recently synthesized (Ti1-x Mo x )2AlC (0 ≤ x ≤ 0.20) solid solutions: deciphering the structural, electronic, mechanical and thermodynamic properties via ab initio simulations.

The structural, electronic, mechanical and thermodynamic properties of (Ti1-x Mo x )2AlC (0 ≤ x ≤ 0.20) were explored using density functional theory. The obtained lattice constants agree well with the experimental values. The electronic band structure confirms the metallic nature. Strengthening of covalent bonds due to Mo substitution is confirmed from the study of band structure, electronic density of states and charge density mapping. The elastic constants satisfy the mechanical stability criteria. Strengthening of covalent bonds leads to enhanced mechanical properties. (Ti1-x Mo x )2AlC compounds are found to exhibit brittle behavior. The anisotropic nature of (Ti1-x Mo x )2AlC is revealed from the direction dependent Young's modulus, compressibility, shear modulus and Poisson's ratio as well as the shear anisotropic constants and the universal anisotropic factor. The Debye temperature, minimum thermal conductivity, Grüneisen parameter and melting temperature of (Ti1-x Mo x )2AlC have been calculated for different Mo contents. Our calculated values are compared with reported values, where available. This journal is © The Royal Society of Chemistry.

Introduction

Intrinsically nanolaminated MAX phase compounds have drawn significant attention from the materials science community because of their unusual combination of physical properties such as damage tolerance, high stiffness, resistance to corrosion and thermal shock like ceramics as well as fairly high electrical and thermal conductivity, and machinability like metallic compounds.[1,2] This outstanding combination of these properties makes MAX compounds suitable for extensive use in a number of applications mentioned in the literature.[1-6] In the 1960s, the powder forms of these carbides/nitrides were first discovered,[7-10] but interest was reinvigorated after the studies of Barsoum and coworkers in the 1990s[6,11] which disclosed their interesting properties in detail for the first time. The physical properties of these systems originate mainly due to the layered structure of MAX phases that have A-layers (Sn, Al, Ge, etc.) sandwiched in between MX sheets (e.g., Ti2C) with the stacking of atomic layers along the c-axis.[12,13] To date, a number of compounds (around 80) belonging to MAX phases have already been synthesized,[1,14] theoretically around 155 are known[15] out of a possible 665 viable MAX phases (MAX, n = 1–4)[16] and the number of reports performed theoretically as well as experimentally are increasing tremendously.[3,17-38] The study of the existing and potentially viable MAX phase nanolaminates has established itself as an important subfield of materials science research.

The Ti2AlC compound belongs to the Ti–Al–C system representing a large family of solids with overall advantages of MAX phases and thus becomes one of widely studied compounds from the 211 MAX phase materials.[39] The M–Al–X phases are most attractive for high-temperature applications, such as Ti2AlC and Cr2AlC because of their high oxidation resistant (up to 1400 °C) resulted from the slow growing protective alumina (Al2O3) layers.[40,41] The coefficient of thermal expansion (8.2 × 10−6 K−1) of Ti2AlC[42] is almost equal to that of α-Al2O3 (9.3 × 10−6 K−1),[43] gives it another advantage by making it resistant to thermal stresses upon thermal cycling.[44] However, low hardness and strength limit its potential use as a structural component like other binary carbide ceramics.[45] Therefore, an increase of hardness and bonding strength along with melting temperature of Ti2AlC will be highly advantageous to extend its domain of applications.

Scientific community are always keen on the reporting of new MAX phase solid solutions, one of the efficient ways of altering the properties of existing MAX phase materials, opening new wing of research and widening their field of applications. Enhancement of many important properties such as oxidation,[46-50] fracture toughness,[51,52] strength[48,51] and self-healing properties[47,53] are reported in the solid solutions. Attempts have also been made to increase the hardness and strength of Ti2AlC both theoretically and experimentally.[42,45,54-57] A considerable increase in elastic modulus is reported by Nb substitution in Ti2AlC.[42] Strengthening of (M and M′ = Ti, V, and Cr) solid solutions has been reported by Wang et al.[56] Meng et al.[57] have investigated Ti2AlC and (Ti0.8V0.2)2AlC where enhancement of Vickers hardness, flexural strength, and shear strength were reported due to V substitution. A considerable improvement in compressive strength has already been reported for Ti2AlC1−N.[58] Recently, Pan et al.[59] have successfully synthesized the Mo substituted Ti2AlC solid solution where a significant increase of Vickers hardness, flexural strength, and fracture toughness were reported. The atom Mo has the possibility to form (Ti1−Mo)2AlC solids solutions because of its smaller atomic radius compared to Ti. The earlier reports also proved the ability of Mo to modify the microstructural and mechanical properties of MAX phase materials.[45,60] Moreover, we have investigated the effect of M mixing on the physical properties of (Zr,Ti)2AlC in our previous study.[27] Thus, the mechanical and thermodynamic properties are expected to be upgraded by substituting Mo into Ti2AlC.

Therefore, we aimed to study the structural, electronic, mechanical and thermodynamic properties of (Ti1−Mo)2AlC (0 ≤ x ≤ 0.20) solids solutions theoretically for the first time where the modification of mechanical and thermodynamic properties due to alloying with Mo are explored in details. The value of x is limited up to 0.20 and only the experimentally synthesized compositions are chosen to keep away from the controversy of the thermodynamical stability. The rest of the paper has been structured as follows. Section 2 describes the computational methodology in brief. Results are presented, analyzed and discussed in Section 3. Finally, important conclusions from this study are summarized in Section 4.

Calculation methods

The CAmbridge Sequential Total Energy Package (CASTEP) code[61] was used to calculate the physical properties of (Ti1−Mo)2AlC solids solutions via the density functional theory (DFT)[62,63] based on the plane–wave pseudopotential method. The generalized gradient approximation (GGA) of Perdew–Burke–Ernzerhof (PBE) is used as an exchange and correlation functional.[64] The pseudo atomic calculation was performed for C – 2s22p2, Al – 3s23p1, Ti – 3s23p63d24s2 and Mo – 4s24p64d55s1 electronic orbitals. The cutoff energy of 400 eV and a k-point mesh of 12 × 12 × 3 (ref. 65) were set for good convergence using the Broyden Fletcher Goldfarb Shanno (BFGS) geometry optimization[66] technique. The self-consistent convergence of the total energy was 5 × 10−6 eV per atom, and the maximum force on the atom was 0.01 eV Å−1. The maximum ionic displacement was set to 5 × 10−4 Å, and a maximum stress of 0.02 GPa was used. The Virtual Crystal Approximation (VCA) approach was employed to model the solid solutions. In the VCA the crystal preserves the primitive periodicity with mixture of atoms. VCA allows the study of crystals with mixed atoms by preserving unit cell of the unsubstituted parent material. The disorder due to alloying is modeled by defining mixture atoms or hybrid atoms where one atomic site is occupied by one of the two different types of atoms. The relative concentrations can be set for any number of atoms involved, where the total concentration of the substituted and substituent atoms (for instance, Ti and Mo, respectively, in (Ti1−Mo)2AlC solid solutions) is always 100%. In practice, this approximation works well when the level of alloying substitution is not very high and the Vegard's law is reasonably held. These conditions are satisfied for the alloys considered here. Moreover, a number of works have been published recently using the VCA[27,28,67-70] which resulted in reliable estimates of thermo-physical and electronic band structure properties of MAX phase compounds.

Results and discussion

Structural properties

The MAX phases are crystallized in hexagonal structure with the space group P63/mmc. There are 8 atoms (for 211 system) in the unit cell of which 4M, 2A and 2X atoms are contributed to the total 8 atoms. In case of Ti2AlC, the position of carbon atoms is 2a, Ti atoms is 4f and Al atoms is 2c Wyckoff positions,[6] as shown schematically in Fig. 1(a). In order to calculate the ground state physical properties, we have to obtain the equilibrium structure by optimizing the total energy of the solids under consideration. Optimized lattice constants are presented in Table 1 together with values obtained by others.[71-75]

Fig. 1

(a) Crystal structure (unit cell) of the Ti2AlC compound and (b) lattice constants a (Å) and c (Å) for different Mo contents (the linear fits to the lattice constants are also shown).

The obtained lattice parameters (a and c), c/a ratio, and cell volumes (V) for the (Ti1−Mo)2AlC (0 ≤ x ≤ 0.20) solid solutions, together with values obtained by previous studies

x	a (Å)	c (Å)	c/a	V (Å3)
0.00	3.0670	13.737	4.478	111.94
	3.058a	13.642a	4.461a	106.62a
	3.066b	13.694b	4.466b	111.49b
	3.04c	13.60c
	3.053d	13.640d	4.468d
	3.06e	13.67e	4.46e
0.05	3.0620	13.677	4.466	111.06
0.10	3.0552	13.636	4.463	110.23
0.15	3.0461	13.618	4.470	109.43
0.20	3.0334	13.606	4.485	108.66

Expt.[71]

GGA.[72]

LDA.[73]

GGA.[74]

GGA.[75]

The calculated ground state lattice constants of Ti2AlC: a and c are just 0.29% and 0.69%, respectively, higher than the experimentally obtained values. This ensures very high level of consistency of present first-principles calculations. Fig. 1(b) illustrates the variation of a and c with Mo contents from which it is clear that the lattice constants a and c decrease with Mo content. The decrement of lattice constants can be largely attributed to the difference in ionic radii of Mo (0.69 Å) and Ti (0.74 Å). Since Mo is smaller than the Ti, its substitution causes shrinkage of the unit cell. The decrease of lattice constant c with Mo content is reported in experiment while the values of a have been found to vary by little.[59] However, the compositions considered here follows the Vegard's law reasonably well.[76]

Electronic properties

Fig. 2(a) shows the electronic band structure of Ti2AlC in which no band gap is observed between valence and conduction bands and the density of states (DOS) at Fermi level is 3.04 states per eV [Fig. 2(c)]. Accordingly, Ti2AlC is found to be metallic in nature which is able to conduct electrically as well as thermally, which is in good accord with the experimental results of Barsoum et al.[77] The electrical conductivity exhibited by Ti2AlC is direction dependent. The paths Γ–A, H–K and M–L are along the c-direction which exhibited much less degree of energy dispersion. On the other hand, the paths A–H, K–Γ, Γ–M and L–H are in the basal planes, exhibiting much greater degree of energy dispersion.[73] This shows that the electronic effective mass tensor is higher for c-axis conduction compared to that in the ab-plane. Therefore, anisotropic electrical conductivity is predicted for Ti2AlC i.e., the electrical conductivity is higher in the basal plane than that along the c-direction.

Fig. 2

The electronic band structure of (a) Ti2AlC and (b) (Ti0.95Mo0.05)2AlC. Total and partial DOS of (c) Ti2AlC and (d) the total DOS for different Mo contents.

Fig. 2(b) shows the band structure of (Ti0·95Mo0.05)2AlC exhibiting almost similar characteristics except that the position of some of the dispersion curves. These curves are observed to be shifted slightly to lower energy. The band structures (not shown) for other compositions (x = 0.10, 0.15 and 0.20) also show similar nature with shifting of the curves towards lower energy which increases with increasing Mo contents as can be observed from Fig. 2(d).

Fig. 2(c) shows total and partial DOS of Ti2AlC which gives better understanding of bonding characteristics in the system. The partial DOS is used to disclose the contribution from different electronic states. The dominating contribution comes from Ti-3d with a very small contribution from Al-3p orbital at the Fermi level. The hybridization among different states is observed in different energy ranges. The peak in the energy range −9 eV to −11.5 eV comes from hybridization of C-2s and Ti-4s states with a prime contribution from C-2s states. The peak in the range −1.8 eV to −3.6 eV comes from the hybridization of C-2p and Ti-3d states. This bonding is primarily covalent in nature resulting from strong mixing between C-2p and Ti-3d states which dominates the higher energy range. The peak just below the Fermi level is attributed from the hybridization of Ti-3d and Al-3p states that results in a comparatively weaker covalent bonding than the previous one.[78] The energy range 0 to 5 eV is dominating by the metal-to-metal dd interactions and anti-bonding states. The covalent bond among the unlike atoms determines the overall bonding strength and consequently the mechanical properties of the MAX phases. The total and partial DOS of Mo substituted compositions are also calculated (not shown) and revealed the same qualitative nature of bonding except the positions of the peaks (energy range of hybridization) that comes from the contribution of Mo-4d sates. The shifting of these peaks can be visualized in Fig. 2(d). Fig. 2(d) shows the total DOS for different Mo contents in which two vertical red solid lines are drawn to refer the position of the peaks for the x = 0 compound. These two reference lines clearly demonstrate the shifting of the peaks towards lower energy with increasing Mo content. The hybridized energy states always give an indication of the bonding strength. In general, lower the energy states stronger the bonding strength. Therefore, lowering of the position of the peaks indicates the better hybridization among different states and stronger covalent bonds between them. The strengthening of the covalent bonds can also be observed from Fig. 3. Fig. 3 shows the charge density mapping (CDM) (in the units of e Å−3) in the (101) crystallographic plane. CDM is an important tool to describe the bonding nature of solids. The occupation positions of atoms have been labeled in the figure. A very clear change is observed in the charge density at the Ti/Mo positions with varying level of Mo substitution. It is quite clear from Fig. 3 that the accumulation of charges at the Ti sites increased to a considerable extent owing to Mo substitution and it increases gradually with increasing Mo contents. Therefore, the strength of covalent bond is expected to increase because of the better mixing of d and p electronic states. Hence, the dominating covalent bond between Ti/Mo-d and C-p states as well as comparatively weaker covalent bond between Ti/Mo-d and Al-p states is assumed to be strengthened. Hence, it is reasonable to expect that the overall bonding strength within (Ti1−Mo)2AlC (0 ≤ x ≤ 0.20) solid solutions will increase owing to Mo substitution and this should be reflected in the enhanced mechanical properties presented in the following section.

Fig. 3

The electronic charge density mapping for different Mo contents.

Mechanical properties

The knowledge of elastic properties e.g. stiffness constant, elastic moduli, elastic anisotropy factors and ductile or brittle characteristic is of scientific interest for practical use of engineering materials. The elastic constants studied here are obtained from linear finite strain–stress method as implemented in the CASTEP.[79] The stiffness constants (C) of (Ti1−Mo)2AlC (0 ≤ x ≤ 0.20) solid solutions are calculated and presented in Table 2. The confirmation of mechanical stability is essential for studying the mechanical properties of a solid and we have checked the mechanical stability conditions of (Ti1−Mo)2AlC (0 ≤ x ≤ 0.20) solid solutions using the Born criteria[80] for hexagonal system: C11 > 0, C11 − C12 > 0, C44 > 0, (C11 + C12)C33 − 2C13 > 0. To judge the accuracy of our calculations, the reported elastic constant values of Ti2AlC are also presented herein. Comparison with the calculated values presented here ensures the reliability of our theoretical estimations.

The elastic constants C (GPa), bulk modulus, B (GPa), shear modulus, G (GPa), Young's modulus, Y (GPa), Poisson ratio, υ, Pugh ratio, G/B and Cauchy pressure of (Ti1−Mo)2AlC (0 ≤ x ≤ 0.20) solid solutions

x	C 11	C 12	C 13	C 33	C 44	B	G	Y	υ	G/B	Cauchy pressure	Ref.
0.00	302	59	58	269	107	135	113	265	0.17	0.84	−48	This study
	308	55	60	270	111	137	118	275a	0.17a	0.86a		Ref. 56 and 74
	309	72	68	284	112	146	115	273	0.19	0.78a		Ref. 72
	301	69	62	266	106	139						Ref. 81
	306	59	59	272	109	138						Ref. 82
	302	68	64	268	107	141	111	264	0.18	0.79a		Ref. 83
0.05	304	68	66	277	116	143	116	274	0.18	0.81	−49
0.10	307	55	73	280	122	145	119	280	0.18	0.82	−67
0.15	304	54	84	272	128	148	120	283	0.18	0.81	−74
0.20	323	62	83	278	131	154	124	293	0.18	0.81	−69

Calculated using published data.

Additionally, the uniaxial directional elastic behaviors can also be understood from the elastic constants. Like the usual trends in MAX phases, the values of C11 are observed to be higher than C33 for (Ti1−Mo)2AlC solid solutions. The C11 and C33 measure the resistance required to compress the compound along[100] (along a-axis) direction and along [001] (along c-axis) direction, respectively. The higher values of C11 and C33 are found in Mo substituted compositions revealing the increased bonding strength both along a and c directions. Another important stiffness constant C44 also known as shear stiffness, by itself corresponds to a change in shape keeping the volume unchanged and provides with the information regarding hardness directly.[84] The C44 is also found to be increased by 8.4%, 14.0%, 19.6% and 22.4% for x = 0.05, 0.10, 0.15 and 0.20, respectively. The obtained values of C13 and C12 are close to each other as usually observed in MAX phases.[21]

The mechanical parameters such as bulk moduli (B), shear moduli (G), Young's moduli (Y), Poisson's ratio (υ) are calculated using the formalism that can be found elsewhere.[85] The calculated elastic moduli are found in good agreement with the reported values as shown in Table 2. The two main mechanical parameters of polycrystalline solids are bulk and shear moduli. The ability to resist volume deformation and plastic deformation is represented by B and G which are found to be increased by 5.9%, 7.4%, 10.0%, 14.1% and 2.6%, 5.3%, 6.2%, 9.7% for x = 0.05, 0.10, 0.15 and 0.20, respectively. Another important mechanical parameter for polycrystalline solids is the Young's modulus, by itself representing the stiffness of solids where the higher value indicates the stiffer nature and vice versa. The value of Y is found to be increased by 3.3%, 5.6%, 6.8%, 10.6% for x = 0.05, 0.10, 0.15 and 0.20.

As we have already discussed that the direct information regarding hardness can be known from C44 which is found to be increased by 8.4%, 14.0%, 19.6% and 22.4% for x = 0.05, 0.10, 0.15 and 0.20, respectively. The enhancements of mechanical parameters B, G and Y are considerable due to alloying with Mo and are illustrated in Fig. 4. Although, these parameters cannot predict directly the actual hardness of solids, but all these elastic moduli are usually higher for comparatively harder solids.[86] It should be noted that Pan et al.[59] have measured the Vickers hardness, flexural strength and fracture toughness and found those to be increased with Mo contents with a maximum increase of 44, 34 and 136%, respectively for (Ti0.80Mo0.20)2AlC.

Fig. 4

The variation of C44, B, G and Y with Mo contents in (Ti1−Mo)2AlC solid solutions.

The brittleness of (Ti1−Mo)2AlC (0 ≤ x ≤ 0.20)

The brittleness or ductility of the studied materials has been assessed by calculating the Pugh ratio (G/B), Poisson's ratio (υ) and Cauchy pressure as presented in Table 2. The ratio G/B is named after the name of Pugh, a widely used parameter of brittleness or ductility of solids. A G/B value of 0.571 indicates the border line between the brittleness and ductility.[87] Another widely used parameter to predict brittleness or ductility is the Poisson's ratio where the value 0.26 indicates the border line between the brittleness or ductility as proposed by Frantsevich et al.[88] A value of G/B greater than 0.571 (or υ lower than 0.26) indicates the brittle character of solids while a value of G/B lower than 0.571 (of υ greater than 0.26) exhibits the ductile character of the solids. The values presented in Table 2 revealed the brittle character of the studied solid solutions. The difference between two stiffness constants (C12 − C44) is known as the Cauchy pressure (CP). This parameter is assumed to define the bonding nature in solids depending on the value of CP (either positive or negative).[89] For ductile metallic material the values of CP is positive while a negative value points toward the covalent bonding with brittle nature with directional bonding character. Furthermore, the magnitudes of negative values exhibit the level of dominance of covalent bonding in the overall makeup of the solid. The large negative values CP for (Ti1−Mo)2AlC solid solutions revealed the dominant role of directional covalent bonding in association with brittle nature within the studied compounds. The CP increases with elevated level of Mo implying that Mo substitution makes covalent contribution stronger in (Ti1−Mo)2AlC.

The elastic anisotropy

The compositions considered here are mechanically anisotropic. The study of anisotropy has been considered because of its involvement with the many physical processes of engineering importance, i.e., the creation of plastic deformation in solids, enhanced charged defect mobility, dynamics of micro-scale cracking in ceramics etc.[90]

The values of Y, B, G and υ are direction depend that can be presented by three-dimensional (3D) plots. The plotting of 3D contour graphs is an attractive way of presenting the mechanical anisotropy where a perfect sphere designates perfect isotropy of solids and departure from it spherical shape measures the level of anisotropy.[91] The 3D plots of Y, B (in terms of compressibility), G and υ of (Ti1−Mo)2AlC solid solutions obtained using the ELATE program[92] are shown in Fig. 5–8. These plots were generated from the calculated values of C of (Ti1−Mo)2AlC. The anisotropy level of the studied materials can be understood from the figures where the most anisotropic composition is (Ti0.85Mo0.15)2AlC. Various anisotropy indices are also presented in Table 3 where the value of A is 1 for isotropic solids which represents a perfect sphere in the ELATE generated profile.

The minimum and the maximum values of the Young's modulus, 1/bulk modulus, shear modulus, Poisson's ratio and anisotropy indices of (Ti1−Mo)2AlC solid solutions

x	Y min (GPa)	Y max (GPa)	A	(1/B)min (TPa−1)	(1/B)max (TPa−1)	A	G min (GPa)	G max (GPa)	A	υ min (GPa)	υ max (GPa)	A
0.00	251	282	1.12	2.3	2.7	1.15	108	121	1.12	0.15	0.19	1.25
0.05	254	279	1.10	2.2	2.5	1.13	112	118	1.05	0.16	0.19	1.16
0.10	250	289	1.15	2.3	2.4	1.04	109	125	1.14	0.13	0.22	1.78
0.15	233	296	1.27	2.2	2.2	1.01	100	128	1.27	0.10	0.27	2.69
0.20	242	306	1.26	2.0	2.3	1.12	107	131	1.22	0.12	0.26	2.06

The three shear anisotropic factors A1, A2 and A3 correspond to three independent elastic shear constants for hexagonal crystals and are calculated using the formalism that can be found elsewhere.[28] The values of Ai's are not 1 (for isotropic solids, Ai = 1), revealed the anisotropic nature of the (Ti1−Mo)2AlC solid solutions. The universal anisotropic index of (Ti1−Mo)2AlC solid solutions are also calculated and found to have non-zero value, revealing the anisotropic nature of the (Ti1−Mo)2AlC solid solutions[93] (Table 4).

Anisotropic factors, A1, A2, A3, and universal anisotropic index AU of (Ti1−Mo)2AlC solid solutions

x	A 1	A 2	A 3	A U
0.00	1.04	0.88	0.91	0.22
0.05	0.95	0.98	0.94	0.23
0.10	0.86	0.97	0.83	0.19
0.15	0.74	1.02	0.75	0.19
0.20	0.77	1.00	0.78	0.15

Thermodynamic properties

The important thermodynamic properties such as Debye temperature, minimum thermal conductivity, Grüneisen parameter and melting temperature of the (Ti1−Mo)2AlC solid solutions are studied in this section.

One of the fundamental characteristics parameters of solids is the Debye temperature (ΘD) which promotes understanding of a number of important physical parameters by correlating the mechanical properties with thermodynamic properties such as phonons, lattice vibrations, enthalpy, thermal conductivity, melting point, specific heat etc.[94] A simple but widely used method (details can be found elsewhere[95-97]) proposed by Anderson[98] has been used to calculate ΘD using the average sound velocity in (Ti1−Mo)2AlC solid solutions. The sound velocities were determined from the elastic constants and crystal density of the compounds under consideration.

Table 5 displays the obtained values of sound velocities and ΘD for different Mo content along with the previously reported values for Ti2AlC to assess the reliability of our calculations. The values of transverse sound velocity is much lower than the longitudinal sound velocity;[99] more energy is required to vibrate neighboring atoms during the propagation of transverse wave that results in the loss of wave energy and hence reduced the wave velocity. The calculated values of ΘD presented in Table 5 are observed to vary in opposing trend of mechanical properties. As we know, in general, ΘD are higher for harder solids[100,101] and the hardness is assumed to be increased with increased elastic constants and moduli but the ΘD is observed to be decreased with Mo contents. This apparent contradiction vanishes when we take into consideration of the much larger atomic weight of Mo (95.95) compared to Ti (46.86). In fact, within the harmonic approximation, the Debye frequency (equivalently the Debye temperature) varies inversely with the square-root of the atomic weight. Therefore, the slight reduction in the ΘD with increasing Mo content does not imply that the lattice is becoming less stiff, it rather implies that the lattice is becoming relatively stiffer with Mo substitution when the huge difference in mass between Mo and Ti is properly taken into account. This follows from the harmonic approximation where the frequency (ν) of the phonon mode can be expressed roughly as, ν ∼ (k/M)1/2. Here, k is the stiffness constant or force constant and M is the mass of the atomic species.

Calculated density (ρ), longitudinal, transverse and average sound velocities (vl, vt, and vm, respectively), Debye temperature (ΘD), minimum thermal conductivity (Kmin), Grüneisen parameter (γ) and melting temperature (Tm) of (Ti1−Mo)2AlC solid solutions

x	ρ (g cm−3)	v l (m s−1)	v t (m s−1)	v m (m s−1)	Θ D (K)	K min (W m−1 K−1)	γ	T m (K)
0.00	3.997	8453	5316	5853	723	8.47	1.18	1663
	4.013a	8640a	5350a	5900a	729a
0.05	4.172	8446	5272	5808	719	8.45	1.21	1681
0.10	4.349	8356	5230	5761	715	8.43	1.20	1695
0.15	4.526	8248	5148	5672	706	8.34	1.21	1674
0.20	4.705	8238	5133	5656	705	8.35	1.22	1740

Ref. 72.

The MAX phase materials are promising candidates for high temperatures applications. Therefore, study of high temperature physical properties is of scientific interest. Like, ΘD, another fundamental physical parameter is the lattice thermal conductivity (Kph), resulting from role of phonon–phonon scattering on thermal transport. At high enough temperatures, the lattice thermal conductivity Kph assumes a minimum constant value (Kmin). The phonons become completely uncoupled and the thermal energy is distributed evenly among the atoms within the solid.[102] We have calculated the Kmin using the following equation:[103], where kB, vm, NA and ρ are Boltzmann constant, average phonon velocity, Avogadro's number and density of crystal, respectively and presented those in Table 5. The values Kmin are observed to decrease very slowly with increasing Mo contents. This is largely due to the increase in the crystal density with increasing level of Mo substitution. The Grüneisen parameter (γ) provides the information regarding anharmonic effect in crystalline solids attributed from lattice dynamics. The Grüneisen parameter is calculated using Poisson's ratio: .[104] The calculated values lie in the range of 0.85 to 3.53 for Poisson's ratio within the range 0.05–0.46 for diverse class of polycrystalline materials.[105] Relatively low and almost constant values of γ of (Ti1−Mo)2AlC imply that anharmonic effects are low in all these alloys. The knowledge of melting temperature (Tm) of MAX phases is very important as they are potential candidates for high temperature applications. We have calculated the melting temperature of (Ti1−Mo)2AlC solid solutions via the empirical formula:[106]. The increased values of Tm owing to Mo substitution indicates that compared to Ti2AlC, Mo substituted alloys are more suitable for high temperature applications.

Conclusions

In this paper Mo-substituted MAX phase solid solutions (Ti1−Mo)2AlC (0 ≤ x ≤ 0.20) have been investigated by density functional theory. The calculated lattice constants of Ti2AlC, a and c are just 0.29% and 0.69%, respectively higher than the experimental values[71] and also comparable with the reported values of the parent phase Ti2AlC.[72-75] This demonstrates high level of reproducibility and reliability of our calculations. The nonexistence of energy band gap in the electronic band structure and finite DOS (∼3.0 states per eV) at the Fermi level clearly reveal the metallic nature of the studied solid solutions where the electrical conductivity is expected to be strongly anisotropic. The shifting of dispersion curves in the band structure and peaks in the TDOS towards lower energy and increased accumulation of charge carriers at Ti/Mo sites due to Mo substitution indicate the strengthening of the covalent bonds between unlike atoms (metal-d and C-p states and metal-d and Al-p states). As a consequence various mechanical properties are enhanced. The Fermi level is located in the middle of bonding and antibonding peaks in the TDOS curves (Fig. 2(d)), suggesting a high degree of electronic stability of the studied solid solutions. The computed stiffness constants (C) and polycrystalline elastic constants and ratios (B, G, Y, G/B and υ) are found in good accord with the reported values (for x = 0.0).[56,72,74,81-83] The mechanical stability of the studied solid solutions is validated by Born criteria based on the values of C. The single crystal elastic constants (C11, C33 and C44) and other mechanical parameters of studied solid solutions are found to be increased by a considerable extent owing to Mo substitution for Ti which suggests enhanced mechanical hardness for the same. The values of C11, C33, C44, B, G and Y are increased by 6.5%, 3.3%, 22.4%, 14.0%, 9.7% and 10.5%, respectively for (Ti0.80Mo0.20)2AlC. The Pugh ratio (greater than 0.571), Poisson's ratio (less than 0.26) and Cauchy pressure (−48 to −74 GPa) confirm the brittleness of the studied compositions. The 3D contour plots of elastic moduli (B, G, Y and υ), shear anisotropic factors and universal anisotropy index reveal the anisotropic nature of the considered solid solutions. The values of Debye temperature (quite high compared to many other MAX phases[107-110]) and minimum thermal conductivity decrease (e.g., ΘD: from 723 K for x = 0.00 to 705 K for x = 0.20) while the Grüneisen parameter and melting temperature (γ: from 1.18 to 1.22 and Tm: from 1663 K to 1740 K for x = 0.00 and x = 0.20, respectively) are noted to be increased due to Mo substitution.

To summarize, this study has found that Mo substitution can lead to significant improvement of the mechanical and thermo-physical properties of MAX phase (Ti1−Mo)2AlC solid solutions. Therefore, we hope that the scientific community will be inspired to synthesize and study the Mo substituted solids solutions of 211 MAX phases belonging to different sub-classes in future.

Author contributions

M. A. A. designed the project, performed the analysis and wrote the draft manuscript. S. H. N. finalized the draft paper and supervised the project.

Data availability

The data sets generated and/or analyzed in this study are available from the corresponding authors on reasonable request.

Conflicts of interest

The authors declare no competing interests.