[Paper Review] First-principles calculations atomic structure and elastic properties of Ti-Nb alloys

A. N. Timoshevskiĭ, S. O. Yablonovskii+1 more|arXiv (Cornell University)|Nov 30, 2011

Titanium Alloys Microstructure and Properties3 citations

TL;DR

This study uses first-principles density functional theory (DFT) calculations to investigate the atomic structure and elastic properties of Ti-Nb and Ti-Zr-Nb β-phase alloys. It reveals non-monotonic concentration-dependent Young’s modulus in binary Ti₁₋ₓNbₓ alloys due to Nb atom distribution effects, and predicts that partial Zr substitution in Ti-Nb alloys significantly reduces Young’s modulus while stabilizing the single-phase β-phase, identifying optimal compositions for low-modulus biomedical implants.

ABSTRACT

Elastic properties of Ti based β-alloy were studied by the method of the model structure first principle calculations. Concentrational dependence of Young modulus for the binary β-alloy Ti-Nb was discovered. It is shown that peculiarities visible at 15-18% concentrations can be related to the different Nb atoms distribution. Detailed comparison of the calculation results with the measurement results was done. Young modulus for the set of the ordered structures with different Nb atoms location, which simulate triple β-alloys Ti-29.7%Zr-18.5%Nb and Ti-51.8%Zr-18.5%Nb have been calculated. The results of these calculations allowed us to suggest the concentration region for single-phase ternary β-phase alloys possessing low values of Young's modulus.

Motivation & Objective

• To determine the concentration dependence of Young’s modulus in binary Ti₁₋ₓNbₓ β-phase alloys using first-principles calculations.
• To investigate how Nb atom distribution (ordered vs. disordered) affects elastic properties in these alloys.
• To evaluate the effect of Zr substitution on the electronic structure and elastic moduli of ternary Ti-Zr-Nb alloys.
• To identify composition regions in Ti-Zr-Nb alloys that exhibit both single-phase β-phase stability and low Young’s modulus.
• To provide theoretical guidance for designing low-modulus, biocompatible β-titanium alloys for medical implants.

Proposed method

• First-principles calculations were performed using the full-potential linearized augmented plane wave (FP-LAPW) method in the Wien2k package.
• The generalized gradient approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) was used for the exchange-correlation functional.
• Structural optimization was performed under both homogeneous and inhomogeneous deformations, including full relaxation of lattice parameters and atomic positions.
• Elastic constants (C₁₁, C₁₂, C′) and Young’s modulus were calculated using the Voigt-Reuss-Hill (VRH) approximation.
• Model supercells of 54 atoms (3×3×3 BCC unit cells) were used to simulate ordered solid solutions with controlled Nb and Zr concentrations.
• The calculations were conducted at 0 K, with high k-point sampling (2000 k-points) and high angular momentum basis sets (lmax=12, Lmax=6) for accuracy.

Experimental results

Research questions

• RQ1 How does the concentration of Nb in Ti₁₋ₓNbₓ alloys affect the Young’s modulus, and is the dependence monotonic?
• RQ2 What is the role of Nb atom distribution (e.g., clustering vs. random) in influencing the elastic properties of β-Ti-Nb alloys?
• RQ3 How does partial substitution of Ti with Zr affect the elastic modulus and phase stability in Ti-Zr-Nb ternary alloys?
• RQ4 Can first-principles calculations predict compositions of single-phase β-Ti-Zr-Nb alloys with low Young’s modulus suitable for biomedical implants?
• RQ5 What is the theoretical Young’s modulus of ordered Ti-Zr-Nb structures with Nb content fixed at 18.5 at.% and varying Zr content?

Key findings

• The concentration dependence of Young’s modulus in binary Ti₁₋ₓNbₓ alloys shows non-monotonic behavior with anomalies at 15–18 at.% Nb, attributed to specific Nb atom distribution configurations.
• For the ordered Ti₄₄₋ₘZrₘNb₁₀ structure with m=16 (Ti₀.₅₁₉Zr₀.₂₉₆Nb₀.₁₈₅), the calculated Young’s modulus was 43.04 GPa, significantly lower than in the binary Ti-Nb system.
• Further increasing Zr content to m=28 (Ti₀.₂₉₆Zr₀.₅₁₉Nb₀.₁₈₅) resulted in a Young’s modulus of 49.75 GPa, indicating a moderate increase due to higher Zr content despite phase stabilization.
• The A-type structures (with Nb atoms isolated from each other) showed lower Young’s modulus than B-type structures (with Nb-Nb neighbors), confirming the importance of Nb distribution.
• The study predicts that Zr substitution stabilizes the single-phase β-phase at lower Nb contents, with the hatched region in the phase diagram (Fig. 4) indicating optimal compositions for low modulus.
• Theoretical predictions were validated by comparison with experimental data, showing good agreement for elastic constants and Young’s modulus in similar systems.