Physics Research International

Physics Research International / 2011 / Article

Research Article | Open Access

Volume 2011 |Article ID 536183 | 7 pages | https://doi.org/10.1155/2011/536183

First Principles Computations of Second-Order Elastic Constants (SOEC) and Equations of State of Rutile TiO2

Academic Editor: Ravindra R. Pandey
Received09 Apr 2011
Revised20 Jun 2011
Accepted21 Jun 2011
Published06 Sep 2011

Abstract

First principles computations of second-order elastic constants (SOECs) and bulk moduli (B) are carried out by ELASTCON and equation of state (EOS) programs. Computational results of lattice parameters, elastic constants, and bulk moduli are obtained with a wide variety of potentials and a limited combination of basis sets and are compared with computational and experimental results by other researchers in the field. DFT hybrid potentials provided the best comparison with the experiment.

1. Introduction

Titanium dioxide (TiO2) is an important transition metal oxide, which exists in anatase, rutile, and brookite phases [1]. There is a considerable interest in the fabrication [2] of this material. It is employed extensively in semiconductors, optical devices, photovoltaic cells, gas sensors, and electrochemical storage devices [3]. Due to its technological importance, the crystalline structure and properties of rutile TiO2 phases have remained a major focus of researchers [1, 4, 5].

Ab-initio Hartree Fock (HF) and density functional theory (DFT) techniques have been employed by various research groups to compute the optimized electronic structure, band gap, and charge density of rutile TiO2 [1, 6โ€“8].

The ab-initio computation of the elastic constants and bulk moduli of rutile TiO2 is presented in this work using the CRYSTAL09 code [9]. The availability of precise and accurate values of experimental lattice parameters, elastic properties, [10] and chemical bonding [11] has provided a considerable challenge to ab-initio computational codes. Computational codes based upon the linear combination of atomic orbitals (LCAO) and plane waves (PW) were employed to compute the optimized crystalline structure of rutile TiO2 [1, 6โ€“8][5]. However, the inadequacy of the existing potentials [1, 6, 7] to accurately predict the physics of TiO2 has motivated a detailed analysis of a wide range of existing and new potentials [6, 12]. Research on the structural and electronic properties has confirmed the relative accuracy of the hybrid potentials [1, 6, 7, 12].

There has been no systematic effort to explore the elastic properties of rutile TiO2 using the LCAO code CRYSTAL09 [9]. CRYSTAL09 possesses a combination of geometry optimization techniques, basis sets, potentials, and algorithms such as ELASTCON [9] (for second-order elastic constants) and EOS (equation of state) [9]. The term ELASTCON represents a computational algorithm that calculates second-order elastic constants and bulk moduli in a highly automated manner. The algorithm possesses various parameters which can be explored to enhance the computational accuracy for a wide variety of crystalline systems. The experimental values of lattice parameters, elastic constants and bulk moduli [10] of the TiO2 rutile phase provided an opportunity to confirm the relative merit of the self-consistent field (SCF) process, the ELASTCON, and EOS programs [13, 14].

2. Potentials and Basis Sets

The HF, local and non-local DFT potentials have consistently resulted in inadequate results of the lattice parameters. First principles computations using PW codes [15โ€“20] can not employ the hybrid mixing of HF exchange and DFT correlations.

Therefore, we employed hybrid potentials in our computations of elastic properties as suggested in [1, 6, 7]. The hybridization between the Ti๐‘‘-orbitals and O๐‘-orbitals requires that an adequate percentage of exchange and correlation is introduced in the potential while computing the electronic structure of a material [21]. The DFT-PWGGA and DFT-PBE potentials lack the exchange part and HF lacks the correlation part, barring them from being as accurate as hybrid potentials.

In addition to using hybrid potentials, the proper choice of basis sets, SCF tolerances, ELASTCON, and EOS parameters can achieve the optimum efficiency as well as accuracy (see Sections 5 and 7). Moreover, it is necessary to employ adequate SCF tolerances consistently.

We selected the O-8411d1 [22] and O-6311d1 [23] basis sets for O atoms and defined these basis sets as basis set 1 and 2, respectively. Both basis sets used a pseudopotential basis set for the Ti atom [9]. The aim of employing these combinations was to separate the basis set dependency from other factors such as potentials, SCF parameters, ELASTCON and EOS tolerances.

The sensitivity of the elastic properties with respect to deviations in the lattice parameters [1] suggests that a well-chosen HF, DFT, and hybrid potentials must be employed to compute the lattice parameters and understand their influences on the elastic constant values. The impact of a particular potential and basis set on the computational values of lattice parameters, elastic constants, and bulk moduli values can be understood by careful comparisons of results achieved for a variety of these basis sets and potentials.

Further, we have also compared our computational results for each potential and basis set combination with [1, 7] to confirm their findings.

3. Computational Procedure

We employed two different basis sets and a variety of potentials to compute the optimized lattice constants and elastic properties (see Section 2 for the discussion of potentials and basis sets). The ELASTCON [9] and EOS [9] algorithms are employed to compute the elastic constants and bulk moduli of rutile TiO2 in an automated manner.

The two basis sets employed are named as basis set 1 and 2 (see Section 2 for basis set definitions). Optimized lattice parameters are computed with HF, DFT, and hybrid potentials using the two unique basis sets. The computations of elastic constants and bulk moduli are performed by employing the ELASTCON algorithm [9]. The computations of the equation of state and bulk moduli are carried out by the EOS algorithm. Finally, the bulk moduli obtained with ELASTCON and EOS algorithms are compared.

Optimized lattice parameters, elastic constants, and bulk moduli, obtained with basis sets 1 and 2, are also compared with the experimental values where possible. A possible contribution of this research work is to assist a general reader in understanding the complex dependance of the elastic properties on the quality of basis sets, potentials, SCF process, ELASTCON, and EOS parameters. Experimental values of lattice parameters, elastic constants, and bulk moduli of rutile TiO2 provide additional assistance for implementation of ELASTCON and EOS programs.

A significant number of computations and experiments have been conducted on rutile TiO2. It is also noticed that a high precision has been achieved in the experimental lattice parameters of rutile TiO2 [24]. The availability of the experimental values of lattice constants, elastic constants, and bulk moduli of rutile TiO2 provides a highly valuable resource to conduct new research. In contrast, computations of elastic constants and bulk moduli have been obtained in an isolated and nonsystematic manner.

A vast majority of ab-initio computations have already tried DFT exchange and correlation potentials [1, 7]. The need for hybrid potentials arose because HF underestimated and DFT potentials overestimated the lattice parameters of rutile TiO2. It has been found that the higher exchange and correlation associated with the transition metals requires new hybrid potentials with variable exchange and correlation [1, 7]. The HF, DFT, and hybrid potentials are employed in this work [9] to fill the gap with regard to the efficient, accurate, and systematic computation of elastic properties of rutile TiO2.

4. Computational Algorithms

4.1. ELASTCON

The computation of elastic constants and bulk moduli is an automated procedure within the ELASTCON algorithm. The computational process begins with determining the crystalline symmetry of rutile TiO2. The analytic first derivative and the numerical second derivative of the total energy are carried out for each specific strain. Levenberg Marquardt (LM) curve fitting [25] is used to compute the required elastic constants for the rutile TiO2 symmetry.

The elastic constants can be extracted from the second derivative of the total energy as๐ถ๐›ผ๐›ฝ=1๐‘‰๐œ•2๐ธ๐œ•๐œ–๐›ผ๐œ•๐œ–๐›ผ||||0.(1) In (1), the terms ๐ถ๐›ผ๐›ฝ, ๐ธ, and ๐‘‰ express the elastic constant tensor, energy, and volume of the crystalline structure, respectively.

The appropriate number of strains is applied in a systematic manner, the elastic constants are calculated, and the compliance coefficients are computed from (3). The compliance coefficients are then utilized for the computation of the bulk modulus as shown in the following:1๐ต=๐‘†11+๐‘†22+๐‘†33๎€ท๐‘†+212+๐‘†13+๐‘†23๎€ธ,[๐‘†]=[๐ถ](2)โˆ’1.(3)

The terms ๐‘†๐‘–๐‘— and ๐ต in (2) express the compliance tensor elements and bulk modulus, respectively. Whereas the term ๐ถ represents the second-order elastic constant.

4.2. EOS

The EOS algorithm [13] employs changes in the optimized volume of rutile TiO2 in a systematic manner. Additionally, the EOS algorithm allows the user to select a range of volumes and number of volumes within that range. For each of the volumes in the range, the CRYSTAL09 optimizes the internal coordinates and lattice parameters while keeping the volume constant. The energy versus volume results are then curve-fitted to an EOS such as the Murnaghan EOS [26] shown in (4) (see [13] for the CVOLOPT option employed during the constant volume geometry optimization for EOS computations).

The EOS algorithm in CRYSTAL09 contains a diversity of equations of state such as Birch Murnaghan, third-order Birch Murnaghan, logarithmic, Vinet, and polynomial. The third-order Birch Murnaghan equation of state algorithm utilized for computing the bulk moduli from the energy versus volume computations for rutile TiO2 is๐ธ(๐‘‰)=๐ต0๐‘‰0๎ƒฌ1๐ต๎…ž๎€ท๐ต๎…ž๎€ธ๎‚ต๐‘‰โˆ’10๐‘‰๎‚ถ๐ตโ€ฒโˆ’1+๐‘‰๐ต๎…ž๐‘‰0โˆ’1๐ต๎…ž๎‚นโˆ’1+๐ธ0.(4) In (4), ๐‘‰0 represents the volume at the minimum energy, ๐ต0 is the bulk modulus at pressure ๐‘ƒ=0, ๐ต๎…ž is the derivative of the bulk modulus at ๐‘ƒ=0, and ๐ธ0 is the minimum energy.

The bulk moduli results are obtained with Levenberg-Marquardt curve fitting of the E versus V computations. The detailed discussion about the ELASTCON [9] and EOS [9] algorithms can be seen in [13, 14].

5. Computational versus Experimental Attempts

Due to its technological importance, a significant number of experiments have been conducted on rutile TiO2 [10, 20, 27]. The experimental results of the electronic structure, band gap, and optic and elastic properties are available for rutile TiO2. The lattice parameters of the rutile TiO2 have been determined [24] precisely up to the fifth significant figure. In addition, the pressure and temperature dependence of elastic constants and bulk moduli have also been explored through experimental means.

The experimental and computational values of elastic constants and bulk moduli [10] of rutile TiO2 have special significance as they serve the purpose of linking these vital branches of research. The importance of achieving a consistency in the lattice parameters of rutile TiO2 among diverse experimental techniques and computational codes is a remarkable achievement which has been repeatedly observed [1, 6, 7].

However, there are deviations in the experimental values of elastic constants and bulk moduli due to their dependence on the experimental details, pressure, and temperature conditions [28]. Therefore, the sources that cause the variations in the experimental values [10, 28โ€“30] must also be considered. However, the dependence of the experimental results on the nature of the experimental setup and ambient conditions can be explored and understood by robust and reliable ab-initio computational methods.

An increase in pressure has shown an increase in the elastic constants and bulk moduli values [31] for rutile TiO2. Moreover, the values of ๐ถ11, ๐ถ33, ๐ถ66, ๐ถ12, and ๐ถ13 increase with increasing pressure, whereas the pressure dependence of ๐ถ44 is not clear. Moreover, the elastic constants are nonlinearly dependent on the temperature.

Another important area where experiments were performed is the volume charge density and chemical bonding of rutile TiO2 [11]. The quantitative convergent beam electron diffraction (QCBED) technique was employed to determine the experimental volume charge density and chemical bonding. The experimental data was utilized to confirm the contribution of ionic and covalent bonding in rutile TiO2. The experimental charge density maps predicted the p-d hybridization between Ti 3-๐‘‘ electrons and O ligands.

The charge density influenced by highly localized ๐‘‘-orbitals of Ti atoms impacts the selection of proper basis sets and potentials in ab-initio computations (see Sections 2 and 7).

Unlike the precision in the experimental lattice parameters, the computational values of lattice parameters vary in the second significant figures. The variation in the lattice parameter values is partially due to the complex nature of the chemical bonding of rutile TiO2. Moreover, the sensitivity of the computational values of lattice parameters is attributed to the choice of potentials and basis sets. However, the hybrid potentials can map the chemical bonding and charge density of rutile TiO2 with a considerable accuracy.

In general, the ab-initio computations are lacking the level of accuracy of the experiments. However, the extensive computational effort by [1, 6, 7, 34] has provided guidance for the present study.

6. Computational Parameters

For the computation of lattice parameters and elastic properties, the SCF tolerances and other computational parameters were carefully chosen. The ELASTCON, EOS, and SCF tolerances were adjusted due to the highly localized nature of transition metal Ti๐‘‘-orbitals. The ELASTCON and other parameters were chosen as STEPSIZE = 0.01, NUMDERIV=7, LGRID = (75, 434), and SHRINK = 9ร—9. The SCF tolerances were fixed as TOLINTEG = 9 9 9 9 18 and TOLDEE = 9 [20].

7. Discussion of Results

Table 1 shows the optimized lattice parameters computed with basis sets 1 and 2, respectively. The unit cell of rutile TiO2, computed with CRYSTAL09 code, is shown in Figure 1. The percentage deviation between the computational and experimental values of lattice parameters and volume is also shown for basis sets 1 and 2. The experimental values from [2, 4, 10] are provided at the bottom of Table 1. Agreement between the computational and experimental values of lattice parameters can be observed for basis sets 1 as well as 2.


๐‘Ž (โ„ซ) ๐‘ (โ„ซ) V o l . ( ร… 3 )

HF4.568 (โˆ’.54)2.980 (.74)62.11 (0)
4.561 (โˆ’.70)2.991 (1.10)62.24 (0)
LDA 4.559 (.74)2.932 (.87) 60.98 ( 2 )
4.539 (โˆ’1.18)2.904 (โˆ’1.85)59.84 ( 4 )
PWGGA4.640 ( 1 ) 2.976 (.33) 64.08 ( 2 )
4.619 (.56)2.946 (โˆ’.40)62.82 ( 1 )
PBE 4.647 (1.17)2.978 (.67)64.32 ( 3 )
4.625 (.69)2.949 (โˆ’.30)63.05 (1.03)
BLYP4.66 (1.45)3.01 (1.75)65.67 (5)
4.657 (1.37)2.971 (.43)64.41 (3.12)
B3LYP4.629 (.78)2.976 (.60)63.78 ( 2 )
4.607 (.21)2.957 (โˆ’.03)62.75 (.55)
B3PW4.599 (.13)2.961 (.10) 62.63 (0)
4.583 (โˆ’.21)2.942 (-.54)61.81 (.95)
PBE04.627 (.73)2.973 (.50)63.69 ( 2 )
4.571 (โˆ’.48)2.940 (โˆ’.61)61.46 (1.52)
Exps. [2, 4, 10, 32]4.5932.95862.40

Elastic constants and bulk moduli computed with ELASTCON are presented in Tables 2 and 3. Significant agreement between the computational and experimental values of elastic constants and bulk moduli is found for basis sets 1 and 2 which can be observed in Figures 2 and 3. The percentage deviation between the computational and experimental values of elastic constants and bulk moduli is shown in Tables 2 and 3.


๐ถ 1 1 ๐ถ 1 2 ๐ถ 1 3 ๐ถ 3 3 ๐ถ 4 4 ๐ถ 6 6 ๐ต

HF364. (36) 216. (23)184. (25) 626. (29)164. (32) 277. (456)270. (27,17)
LDA 311. (16) 211. (21)175. ( 4 ) 504. (28) 159. (33)253. (15)243. (15,6)
PWGGA266. (โˆ’1)177. ( 1 ) 149. ( 1 ) 462. (โˆ’5)138. (11) 223. (17)208. (โˆ’2,โˆ’10)
PBE262. (โˆ’2)176. ( 1 ) 146. (โˆ’1) 458. (โˆ’6) 135. (9)221. (16)205. (โˆ’3,โˆ’11)
BLYP256. (โˆ’5)133. (โˆ’24)143. ( 3 ) 478. (โˆ’1)154. (โˆ’19)206. (8)188. (โˆ’11,โˆ’18)
B3LYP 281. (5) 186. (6)157. (7)506. (5)140. (13)236. (24)220. (4,โˆ’4)
B3PW294. (10)194. (11)164. (12)517. (7)147. (19)246. (29)229. (8,โˆ’1)
PBE0277. ( 3 ) 185. (6)155. (5)494. ( 4 ) 139. (12)234. (23)217. (2,โˆ’6)
Exps. [10, 29, 30, 33]268.00175.00147.00484.00124.00190.00212.00, 230.00


๐ถ 1 1 ๐ถ 1 2 ๐ถ 1 3 ๐ถ 3 3 ๐ถ 4 4 ๐ถ 6 6 ๐ต

HF 393. (47)237. (35)210. (43)662. (37) 171. (38) 303. (60)295. (39,28)
LDA 171. (โˆ’36) 385. (120) 205. (39) 577. (19)130. (5)279. (47) 266. (23,16)
PWGGA 223. (โˆ’17)273. (56) 182. (24) 522. (8) 123. (โˆ’1) 246. (29)237. (12,3)
PBE220. (โˆ’18)270. (54)181. (23)517. (7)122. (โˆ’2)244. (28) 235. (11,2)
BLYP260. (โˆ’3)226. (29)175. (19)509. (5)123. (โˆ’1)236. (24) 231. (9,1)
B3LYP 295. (10)240. (37)193. (31) 562. (16)137. (10) 261. (37)255. (20,11)
B3PW 269. (0)269. (54)197. (34) 569. (18)137. (10) 267. (41)257. (21,12)
PBE0 267. (โˆ’1)282. (61)203. (38)579. (20)139. (12)273. (44) 263. (24,14)
Exps. [10, 29, 30, 33]268.00175.00147.00484.00124.00190.00212.00, 230.00

Table 4 and Figure 3 show the comparisons between the bulk moduli values computed with ELASTCON and EOS algorithms. Excellent agreement is observed between the computational values of bulk moduli achieved with ELASTCON and EOS programs. The agreement between the bulk moduli values points at the computational accuracy of the ELASTCON and EOS programs. Moreover, it is important that the computational results should also be crosschecked with the experimental values of the bulk moduli given at the bottom of Table 4.


๐ต E O S (GPa) ๐‘‰ ๐‘œ ( ร… 3 ) ๐ต E L (GPa)

HF 267. 62.29270.
293. 62.31295.
LDA 241.60.98243.
279.59.84283.
PWGGA 206. 64.09208.
239. 62.87237.
PBE 202. 64.33205.
236.63.10235.
BLYP 191. 65.93188.
229. 64.45 231.
B3LYP 217.63.81220.
252. 62.81255.
B3PW 226. 62.66229.
256. 61.83257.
PBE0232.3662.24217.12
263.61.49263.
Exps.[10, 29, 30, 33] โ€”โ€”212.00, 230.00

However, there is a considerable disagreement between the computational and experimental values of the elastic constants and ๐ต for basis set 1 and 2 with HF and DFT-LDA potentials as shown in Tables 2 and 3. On the other hand, the computational values of the elastic constants and ๐ต provide a significant agreement with the experimental values computed with DFT-PWGGA, DFT-PBE, and hybrid potentials for basis set 1 as shown in Table 2.

The hybrid potentials have shown better agreement due to the adequate percentage of exchange and correlation contributions to total energy of the crystal specifically important for the highly correlated physics of the Ti transition metal. The localized nature of the Ti atom ๐‘‘-orbitals contributes to the higher exchange and correlation effects.

The HF, local DFT, and nonlocal DFT potentials cannot predict results as effectively as hybrid potentials. It can be easily seen in Tables 1โ€“4. In fact, the lack of correlation in HF and the localized nature of DFT-LDA potentials make these potentials less effective for transition metal oxides which possess the covalent as well as ionic nature of chemical bonding.

It must be mentioned that DFT-PWGGA and DFT-PBE exhibit better agreement with experimental results due to the nonlocalized nature of the rutile TiO2 volume charge density. However, the agreement between the DFT-PWGGA, DFT-PBE, and experimental results of elastic constants may not be adequate which can be confirmed by values in Tables 1, 3, and 4.

8. Conclusions

The technological applications of titanium dioxide (TiO2) have generated a significant research activity in experimental and computational sciences. Computational research has resulted in testing the merit of basis sets, potentials, and new programs. The employment of basis sets 1 and 2 has resulted in a considerable agreement between the computational and experimental values of lattice parameters, elastic constants, and bulk moduli. The elastic properties of rutile TiO2 are computed and compared with experimental values. The dependance of experimental values on experimental setup, temperature, and pressure conditions cannot be ignored. The computations of elastic constants and bulk moduli by a wide variety of ab-initio techniques provide results to generate new experiments on this material.

We have separated the factors that determine the quality of computational results of the lattice parameters and elastic properties. The nonlocal DFT and hybrid potentials present better agreement with the experimental values of lattice parameters, elastic constants, and bulk moduli. However, the disagreement between the computational and experimental values of the elastic constants and bulk moduli for HF and DFT-LDA potentials [15, 16, 18] is significant.

The results presented for rutile TiO2 with different potentials are motivated by a variety of challenges related to the existing potentials. For a crystalline system such as rutile TiO2, which has considerable visibility in the experimental research arena, an integrated set of computational results have significant utility. These computations can guide researchers to appreciate the subtle influences of the charge density on lattice parameters and elastic properties of rutile TiO2.

Acknowledgment

One of the authors (W. F. Perger) gratefully acknowledges the support of the Office of Naval Research Grant N00014-01-1-0802 through the MURI program.

References

  1. F. Labat, P. Baranek, C. Domain, C. Minot, and C. Adamo, โ€œDensity functional theory analysis of the structural and electronic properties of TiO2 rutile and anatase polytypes: performances of different exchange-correlation functionals,โ€ The Journal of Chemical Physics, vol. 126, no. 15, pp. 1โ€“12, 2007. View at: Publisher Site | Google Scholar
  2. S. Tanemuraa, L. Miaoa, W. Wunderlicha et al., โ€œFabrication and characterization of anatase/rutileโ€”TiO2 thin films by magnetron sputtering: a review,โ€ Science and Technology of Advanced Materials, vol. 6, no. 1, pp. 11โ€“17, 2005. View at: Publisher Site | Google Scholar
  3. S. D. Mo and W. Y. Ching, โ€œElectronic and optical properties of three phases of titanium dioxide: rutile, anatase, and brookite,โ€ Physical Review B, vol. 51, no. 19, pp. 13023โ€“13032, 1995. View at: Publisher Site | Google Scholar
  4. H. Yao, L. Ouyang, and W. Y. Ching, โ€œAb initio calculation of elastic constants of ceramic crystals,โ€ Journal of the American Ceramic Society, vol. 90, no. 10, pp. 3194โ€“3204, 2007. View at: Publisher Site | Google Scholar
  5. D. I. Bilc, R. Orlando, R. Shaltaf, G. M. Rignanese, J. รรฑiguez, and P. Ghosez, โ€œHybrid exchange-correlation functional for accurate prediction of the electronic and structural properties of ferroelectric oxides,โ€ Physical Review B, vol. 77, article 165107, no. 16, pp. 1โ€“13, 2008. View at: Publisher Site | Google Scholar
  6. F. Corรก, โ€œThe performance of hybrid density functionals in solid state chemistry: the case of BaTiO3,โ€ Molecular Physics, vol. 103, no. 18, pp. 2483โ€“2496, 2005. View at: Publisher Site | Google Scholar
  7. C. Adamo and V. Barone, โ€œToward reliable density functional methods without adjustable parameters: the PBE0 model,โ€ Journal of Chemical Physics, vol. 110, no. 13, pp. 6158โ€“6170, 1999. View at: Google Scholar
  8. M. E. Stizaumainis, T. Ejima, and W. J. James, โ€œThe TiO2 phase explored by the lattice constant and density method,โ€ Acta Crystallographica, vol. 8, pp. 493โ€“497, 1961. View at: Google Scholar
  9. R. Dovesi, V. R. Saunders, C. Roetti et al., CRYSTAL09 User's Manual, University of Torino, Torino, Italy, 2009.
  10. M. Iuga, G. Steinle-Neumann, and J. Meinhardt, โ€œAb-initio simulation of elastic constants for some ceramic materials,โ€ European Physical Journal B, vol. 58, no. 2, pp. 127โ€“133, 2007. View at: Publisher Site | Google Scholar
  11. B. Jiang, J. M. Zuo, N. Jiang, M. O'Keeffe, and J. C. H. Spence, โ€œCharge density and chemical bonding in rutile, TiO2,โ€ Acta Crystallographica, vol. 59, no. 4, pp. 341โ€“350, 2003. View at: Publisher Site | Google Scholar
  12. F. Corรก, M. Alfredsson, G. Mallia et al., Principles and Applications of Density Functional Theory in Inorganic Chemistry II, 2004.
  13. W. F. Perger, J. Criswell, B. Civalleri, and R. Dovesi, โ€œAb-initio calculation of elastic constants of crystalline systems with the CRYSTAL code,โ€ Computer Physics Communications, vol. 180, no. 10, pp. 1753โ€“1759, 2009. View at: Publisher Site | Google Scholar
  14. W. F. Perger, โ€œFirst-principles calculation of second-order elastic constants and equations of state for lithium azide, LiN3, and lead azide, Pb(N3)2,โ€ International Journal of Quantum Chemistry, vol. 110, no. 10, pp. 1916โ€“1922, 2010. View at: Publisher Site | Google Scholar
  15. M. Mikami, S. Nakamura, O. Kitao, H. Arakawa, and X. Gonze, โ€œFirst-principles study of titanium dioxide: rutile and anatase,โ€ Japanese Journal of Applied Physics, vol. 39, no. 8 B, pp. L847โ€“L850, 2000. View at: Google Scholar
  16. X. Ma, P. Liang, L. Miao et al., โ€œPressure-induced phase transition and elastic properties of TiO2 polymorphs,โ€ Physica Status Solidi (B), vol. 246, no. 9, pp. 2132โ€“2139, 2009. View at: Publisher Site | Google Scholar
  17. Y. Liu, L. Ni, Z. Ren, G. Xu, C. Song, and G. Han, โ€œNegative pressure induced ferroelectric phase transition in rutile TiO 2,โ€ Journal of Physics Condensed Matter, vol. 21, article 275901, no. 27, 2009. View at: Publisher Site | Google Scholar
  18. J. Muscat, V. Swamy, and N. M. Harrison, โ€œFirst-principles calculations of the phase stability of TiO2,โ€ Physical Review B, vol. 65, no. 22, Article ID 224112, pp. 2241121โ€“22411215, 2002. View at: Google Scholar
  19. E. Shojaee and M. R. Mohammadizadeh, โ€œFirst-principles elastic and thermal properties of TiO2: a phonon approach,โ€ Journal of Physics Condensed Matter, vol. 22, no. 1, pp. 1โ€“8, 2010. View at: Publisher Site | Google Scholar
  20. J.-Z. Zhao, G.-T. Wang, and Y.-C. Liang, โ€œMechanical properties and electronic structures of cotunnite TiO2,โ€ Chinese Physics Letters, vol. 25, no. 12, pp. 4356โ€“4359, 2008. View at: Publisher Site | Google Scholar
  21. G. Francois and K. Winfried, โ€œLattice dynamics of oxides with rutile structure and instabilities at the metal-semiconductor phase transitions of NbO2 and VO2,โ€ Physical Review B, vol. 31, no. 8, pp. 4809โ€“4814, 1985. View at: Publisher Site | Google Scholar
  22. T. Bredow, K. Jug, and R. A. Evarestov, โ€œElectronic and magnetic structure of ScMnO3,โ€ Physica Status Solidi, vol. 243, no. 2, pp. R10โ€“R12, 2006. View at: Publisher Site | Google Scholar
  23. C. Gatti, V. R. Saunders, and C. Roetti, โ€œCrystal field effects on the topological properties of the electron density in molecular crystals: the case of urea,โ€ The Journal of Chemical Physics, vol. 101, no. 12, pp. 10686โ€“10696, 1994. View at: Google Scholar
  24. B. H. O'Connor and S. Pratapa, โ€œImproving the accuracy of Rietveld-derived parameters by an order of magnitude,โ€ Advances in X-ray Analysis, vol. 45, pp. 158โ€“165, 2002. View at: Google Scholar
  25. W. Marquardt, โ€œAn algorithm for least-squares estimation of nonlinear parameters,โ€ The SIAM Journal on Applied Mathematics, vol. 11, pp. 431โ€“444, 1963. View at: Google Scholar
  26. F. D. Murnaghan, โ€œThe compressibility of media under extreme pressures,โ€ Proceedings of the National Academy of Sciences of the United States of America, vol. 30, no. 9, pp. 244โ€“247, 1944. View at: Google Scholar
  27. I. J. Fritz, โ€œPressure and temperature dependences of the elastic properties of rutile (TiO2),โ€ Journal of Physics and Chemistry of Solids, vol. 35, no. 7, pp. 817โ€“826, 1974. View at: Google Scholar
  28. R. F. S. Hearmon, โ€œThe elastic constants of anisotropic materials,โ€ Reviews of Modern Physics, vol. 18, no. 3, pp. 409โ€“440, 1946. View at: Publisher Site | Google Scholar
  29. L. Gerward and J. S. Olsen, โ€œPost-rutile high-pressure phases in TiO2,โ€ Journal of Applied Crystallography, vol. 30, no. 3, pp. 259โ€“264, 1997. View at: Google Scholar
  30. L. Kรณci, D. Y. Kim, J. S. D. Almeida, M. Mattesini, M. M. E. Isaev, and R. Ahuja, โ€œMechanical stability of TiO2 polymorphs under pressure: ab initio calculations,โ€ Journal of Physics Condensed Matter, vol. 20, no. 34, Article ID 345218, 2008. View at: Publisher Site | Google Scholar
  31. Z. Jun, Y. J. Xin, W. Y. Ju, C. X. Rong, and J. F. Qian, โ€œFirst-principles calculations for elastic properties of rutile TiO2 under pressure,โ€ Chinese Physics B, vol. 17, no. 6, pp. 2216โ€“2221, 2008. View at: Google Scholar
  32. J. K. Burdett, T. Hughbanks, G. J. Miller, J. W. Richardson, and J. V. Smith, โ€œStructural-electronic relationships in inorganic solids: powder neutron diffraction studies of the rutile and anatase polymorphs of titanium dioxide at 15 and 295 K,โ€ Journal of the American Chemical Society, vol. 109, no. 12, pp. 3639โ€“3646, 1987. View at: Google Scholar
  33. S. Ogata, H. Iyetomi, K. Tsuruta et al., โ€œVariable-charge interatomic potentials for molecular-dynamics simulations of TiO2,โ€ Journal of Applied Physics, vol. 86, no. 6, pp. 3036โ€“3041, 1999. View at: Google Scholar
  34. K. Held, V. Anisimov, V. Eyert et al., โ€œLDA+DMFT investigations of transition metal oxides and f-electron materials,โ€ in Advances in Solid State Physics, vol. 43, pp. 267โ€“286, Springer, Berlin, Germany, 2003. View at: Google Scholar

Copyright ยฉ 2011 Ghous Narejo and Warren F. Perger. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

1175ย Views | 732ย Downloads | 2ย Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at help@hindawi.com to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.