ACS Publications. Most Trusted. Most Cited. Most Read
Investigation into Structural Phase Transitions in Layered Titanium-Oxypnictides by a Computational Phonon Analysis
My Activity

Figure 1Loading Img
  • Open Access
Article

Investigation into Structural Phase Transitions in Layered Titanium-Oxypnictides by a Computational Phonon Analysis
Click to copy article linkArticle link copied!

View Author Information
Research Center for Advanced Computing Infrastructure, JAIST, Asahidai 1-1, Nomi, Ishikawa 923-1292, Japan
PRESTO, Japan Science and Technology Agency (JST), Kawaguchi, Saitama 332-0012, Japan
§ Center for Materials Research by Information Integration, Research and Services Division of Materials Data and Integrated System, National Institute for Materials Science, Tsukuba 305-0047, Japan
School of Information Science, JAIST, Asahidai 1-1, Nomi, Ishikawa 923-1292, Japan
Open PDF

Inorganic Chemistry

Cite this: Inorg. Chem. 2017, 56, 22, 13732–13740
Click to copy citationCitation copied!
https://doi.org/10.1021/acs.inorgchem.7b01709
Published November 2, 2017

Copyright © 2017 American Chemical Society. This publication is licensed under these Terms of Use.

Abstract

Click to copy section linkSection link copied!

We applied ab initio phonon analysis to layered titanium-oxypnictides, Na2Ti2Pn2O (Pn = As and Sb), and found a clear contrast between the cases with lighter/heavier pnictogen in comparison with experiments. The result completely explains the experimental structure at low temperature, C2/m for Pn = As, within the conventional charge density wave, while there arise discrepancies when the pnictogen gets heavier. Our phonon calculation using the GGA-PBE functional predicts that a Cmce polymorph is more stable than the experimentally observed one (Cmcm) for Pn = Sb. On the basis of further quantitative analysis, we suggest the possibility that the GGA-PBE functional does not properly reproduce the electron correlation effects for Pn = Sb, and this could be the reason for the present discrepancy.

Copyright © 2017 American Chemical Society

Synopsis

We applied ab initio phonon analysis to layered titanium-oxypnictides, Na2Ti2Pn2O (Pn = As and Sb), and found a clear contrast between the cases with lighter/heavier pnictogen in comparisons with experiments. Our phonon calculations completely explain the experimental structure, C2/m for Pn = As, within the conventional charge density wave, while discrepancies arise when the pnictogen gets heavier, with the same trend for BaTi2Pn2O (Pn = As, Sb, and Bi).

Introduction

Click to copy section linkSection link copied!

Recently discovered layered titanium-oxypnictides, ATi2Pn2O [A = Na2, Ba, (SrF)2, and (SmO)2; Pn = As, Sb, and Bi], (1-13) comprise a new superconducting family, and this superconducting mechanism has attracted intensive attentions because the crystal and electronic structures are similar to those of exotic superconductors such as cuprates (14) or iron arsenides. (15) The superconducting mechanism has been expected to be electron–phonon driven, (16-19) and the accompanying singularities in resistivity and susceptibility at low temperature are regarded as being due to conventional charge density wave (CDW). (17, 19-21) Though our latest phonon evaluation using GGA-PBE without assuming any magnetic orderings revealed that the singularities could be explained within the conventional CDW for Pn = As, (22) it is still unknown if this is the case for Pn = Sb and Bi. (16, 22) Interestingly, recent works are pointing out the possibility of exotic mechanism of the structural transition for these compounds: For BaTi2Pn2O (Pn = As and Sb), Frandsen et al. reported the breaking of the 4-fold symmetry at low temperature without any superstructure peaks by neutron and electron diffractions (23) and suggested the possibility of intra-unit-cell nematic CDW, which is observed in cuprates and iron arsenides superconductors, (24-26) to account for the breaking. A theoretical suggestion was made that such an exotic CDW could be realized by the orbital ordering mediated by spin fluctuations. (27) Thus, the mechanism of structural transition in layered titanium-oxypnictides is still under debate.
Davies et al. (28) have just synthesized large single crystals of Na2Ti2Pn2O (Pn = As and Sb; Figure 1) in order to achieve reliable diffaction measurements and observed clear superstructure peaks of 2 × 2 × 2 and 2 × 2 × 1 for Pn = As and Sb, respectively. They did not detect any indication of spin orderings, so they concluded that the structural transition could be attributed to the conventional phonon-assisted CDW. From the theoretical side, examining whether or not the observed superstructures are explained by ab initio phonon calculations can provide a critical clue to the mechanism of structural transition. In the present study, we performed phonon calculations for Na2Ti2Pn2O (Pn = As and Sb) using the GGA-PBE functional and compared the superstructures predicted by the calculations with the experimentally observed ones. The result shows that the phonon calculation can explain the observed C2/m (monoclinic No.12) structure for Pn = As, indicating that the present phonon calculation works well for this system. For Pn = Sb, however, we obtained a Cmce (orthorhombic No. 64) polymorph, being inconsistent with the observed one [Cmcm (orthorhombic No. 63)]. On the basis of these results combined with preceding studies (16, 22, 29-32) and further quantitative analysis, we propose a new possibility that a heavier Pn causes stronger electron correlations and induces the discrepancy.

Figure 1

Figure 1. Crystal structure of Na2Ti2Pn2O (Pn = As and Sb). The space group is I4/mmm.

Results

Click to copy section linkSection link copied!

Figure 2 shows band structures of undistorted Na2Ti2Pn2O (Pn = As and Sb). The results are consistent with the previous calculations, (32, 33) meaning that there is no artifact due to the choice of cutoff, k-meshes, smearing parameters, and pseudo potentials. Figure 3a shows the phonon dispersions for Pn = As, giving imaginary frequencies appearing around N = (2π/a)·(1/2, 0, a/2c) and Σ′ = (2π/a)·(1/2, 0, 0), where a and c are conventional lattice constants for the undistorted structure. The compatible modes with the experimental observation (28) where the twice larger periodicity along c-axis is realized are N = (2π/a)·(1/2, 0, a/2c) and N′ = (2π/a)·(0, 1/2, a/2c). We note that the previous phonon calculation (34) could not reproduce this doubled periodicity. We therefore took these for further lattice relaxations from the original I4/mmm structure along the mode displacements to get the C2/m (monoclinic, No. 12) superstructure (Figures 4a and 5a), being consistent with the experiments. (28) The results of geometry optimizations for the undistorted (I4/mmm) and distorted (C2/m) structures are shown in Tables 1 and 2. The resultant optimized geometry of the C2/m superstructure gives fairly good agreement with the experiments (28) within deviations at most 2.5%. (35) Although we could not confirm whether or not the negative modes disappear in the C2/m superstructure even at high-symmetry points because of intractable computational costs for the enlarged unit cell, it should be confirmed in the future.

Figure 2

Figure 2. Electronic band structures for undistorted Na2Ti2Pn2O structure ((a) Pn = As, (b) Pn = Sb). The Fermi energy is set to 0 eV.

Figure 3

Figure 3. Phonon dispersions and phonon DOS for undistorted (I4/mmm) Na2Ti2Pn2O (Pn = (a) As, (b) Sb). The red cross marks in (b) represent the frequencies with spin–orbit coupling.

Figure 4

Figure 4. In-plane superstructures obtained by our calculations for Pn = (a) As and (b) Sb, and (c) the observed in-plane superstructure for Pn = Sb. (28) Pn is located above and below the Ti2O plane. Solid lines represent the original unit cells for I4/mmm. Thin dash lines represent the unit cells for 2 × 2 superstructures. Solid dash lines represent redefined unit cells for distorted superstructures. The conventional lattice vectors of the superstructure for (a) Pn = As (C2/m) are redefined as , , and , those for (b) Pn = Sb (Cmce) are redefined as , , and , and those for (c) Pn = Sb (Cmcm) are redefined as , , and , where , , and are conventional lattice vectors of the undistorted structures.

Figure 5

Figure 5. 3D superstructures obtained from our calculations for Pn = (a) As (C2/m) and (b) Sb (Cmce).

Table 1. Results of Geometry Optimization for Undistorted Na2Ti2Pn2O Structure (Pn = As and Sb)a
labelxyzwyckoff
Pn = As (I4/mmm)
Na0.500000.500000.181464e
Ti0.500000.000000.000004c
As0.000000.000000.117524e
O0.500000.500000.000002b
Pn = Sb (I4/mmm)
Na0.500000.500000.183694e
Ti0.500000.000000.000004c
Sb0.000000.000000.120154e
O0.500000.500000.000002b
a

The relaxed conventional lattice parameters are a = 4.07 Å, c = 15.44 Å (a = 4.14 Å, c = 16.98 Å) for Pn = As (Sb), giving good agreement with the experiments. (2, 54) The height of Sb toward Ti2O plane obtained from the present calculation is 2.04 Å, which is consistent with the experimental value at room temperature (2.01 Å). (54) For Pn = As, the height of As toward Ti2O plane has not been determined in experiments, but the obtained value (1.81 Å) is quite consistent with the experimental results of BaTi2As2O (1.77 Å) (4) and (SrF)2Ti2As2O (1.81 Å). (3)

Table 2. Results of Geometry Optimization for Na2Ti2Pn2O Superstructures (Pn = As and Sb) Obtained by Phonon Calculationsa
labelxyzwyckoff
Pn = As (C2/m)
Na0.909070.000000.636994i
Na0.409430.000000.636944i
Na0.159320.249810.637258j
Ti0.130310.119790.000008j
Ti0.880190.369690.000008j
As0.808770.000000.236254i
As0.309330.000000.236254i
As0.059050.249690.236238j
O0.000000.000000.000002a
O0.250000.250000.000004e
O0.000000.500000.000002b
Pn = Sb (Cmce)
Na0.184110.000000.000008d
Ti0.000000.232660.268128f
Sb0.378830.000000.000008d
O0.000000.000000.000004a
a

The conventional lattice parameters after relaxation are a′=11.51 Å, b′ = 11.51 Å, c′ = 8.23 Å, and β = 110.46° (a′ = 17.04 Å, b′ = 5.84 Å, and c′ = 5.83 Å) for Pn = As (Sb).

For Pn = Sb, we obtained phonon dispersions as shown in Figure 3b. Imaginary frequencies appear widely around X = (2π/a)·(1/2, 1/2, 0) and P = (2π/a)·(1/2, 1/2, a/2c). The experiment (28) reported that the twice larger super periodicity is realized only in ab-plane, not along c-axis. Therefore, the X mode is likely to realize the instability toward the superstructure. By optimizing the geometry along the mode displacement, we obtained a Cmce (orthorhombic, No. 64, Figures 4b and 5b), being inconsistent with the observed Cmcm (orthorhombic, No. 63, Figure 4c) in the experiment. (28) To confirm the present conclusion further, we examined the stability by seeing if the imaginary frequencies disappear at the relaxed structure, as shown in Figure 6. In addition to the stability check, we compared enthalpies between our Cmce and the observed Cmcm. The optimized geometry with Cmce is found to be more stable by 1.1 (0.7) mRy per formula unit than that with Cmcm when using Marzari–Vanderbilt (Gaussian) smearing. These results indicate that the experimentally observed Cmcm is unfavorable at 0 K at the GGA level of DFT. It is, however, expected that GGA+U could reverse the order as discussed later. A comparison of Helmholtz free energies (ETS) under quasiharmonic approximation (36) is necessary to discuss the stability more precisely, where the entropic term is calculated from phonon DOS. (36) It was, however, intractable because the large unit cell and the low symmetry of the Cmcm polymorph requires much more computational time to calculate its entropic term (TS). The results of geometry optimizations for the undistorted (I4/mmm) and distorted (Cmce) structures, and those for the experimentally observed superstructure (Cmcm) are shown in Tables 13.

Figure 6

Figure 6. Phonon dispersions for Na2Ti2Sb2O superstructure (Cmce).

Table 3. Results of Geometry Optimization for Experimentally Observed Na2Ti2Sb2O Superstructurea
labelxyzwyckoff
Pn = Sb (Cmcm)
Na0.000000.374980.932938f
Na0.000000.125030.432948f
Na0.749990.125000.9327516h
Ti0.374170.000830.250008g
Ti0.625840.249180.250008g
Ti0.87581–0.000850.250008g
Ti0.124180.250820.250008g
Sb0.000000.375330.129098f
Sb0.000000.125320.629088f
Sb0.750330.124990.1290716h
O0.000000.625000.250004c
O0.000000.125000.250004c
O0.750000.375000.250008g
a

The conventional lattice parameters after relaxation are a′ = 11.74 Å, b′ = 11.74 Å, and c′ = 16.83 Å.

Discussion

Click to copy section linkSection link copied!

The present results confront the sharp contrast between Pn = As and Sb in terms of whether or not the present phonon calculations can reproduce the experimental structural phase transitions. Since it is clearly reported (28) that the measured XRD results cannot be identified by the × × 1 superstructure (Figure 4b), the discrepancy would be intrinsic so that it should be accounted for by the effects that are not taken into account in the present calculations. A spin–orbit effect is likely to be the origin of the discrepancy because Sb is heavier than As. To confirm, we again performed phonon calculations with fully relativistic pseudopotentials at high symmetry positions (Figure 3b). Although the frequencies slightly change because of the spin–orbit effect, the X point that induces the Cmce superstructure still shows the largest imaginary frequency. This result indicates that the spin–orbit effect is not the origin of the discrepancy.
We suggest the possibility that the GGA-PBE functional does not properly reproduce the electron correlation effects for Pn = Sb, and this could be the reason for the present discrepancy. The conclusion was obtained by the following discussion: In iron arsenide superconductors, the strength of electron correlation is well-captured by a trend of h, a vertical distance between Fe layer and Pn/Ch (Pn = P and As; Ch = Se and Te). (30) Miyake et al. applied constrained random-phase approximation (cRPA) to iron arsenide superconductors and revealed that the covalency of Fe-Pn/Ch bonding estimated by DFT-LDA calculation and the resultant Wannier orbitals can be used for a measure of the magnitude of electron correlation in the family. (30) When h gets larger, the covalency gets weaker to make Fe 3d more localized, and the spin/orbital polarizations get enhanced as one of the electron correlation effects. This measure probably works well for layered titanium-oxypnictides because they have crystal and electronic structures similar to those of iron arsenides. (8, 15, 23, 30, 31) In order to apply the measure to the layered titanium-oxypnictides, we quantitatively analyzed the covalencies of Ti–Pn bonding in Na2Ti2Pn2O (Pn = As and Sb). The ratio of partial density of states (pDOS) of each atom to total density of states (tDOS) can be used as a measure of covalency, (37, 38) which is one of the common methods to evaluate covalency within DFT framework. (37-41) In this method, it is necessary to calculate a center of mass (CM) of atomic orbital for each atom. The CM is calculated by the following equation: (38)(1)where A is an element, ε is an energy, gA(ε) is pDOS of A at ε, and CMA is a center of mass of the element A. Note that the energy window [ε2, ε1] in the equation has to be chosen in such a way that it encompasses all relevant bands contributing to the bond. (38) The difference in CM can be used as a measure of covalency of the bond. (38) If the difference in the CM between X and Y (|CMX – CMY|) is small, then the covalency of the bond is strong, and vice versa. Figure 7 shows the obtained pDOS and center of masses of Ti 3d, As 4p, and Sb 5p orbitals for Na2Ti2Pn2O (Pn = As and Sb). The pDOS indicates the existence of two types of bonds: One is Ti–O bonding contributing to [− 8.0 to −4.9 (4.5) eV], another is Ti–Pn bonding contributing to [− 4.9 (4.5) to 0.0 eV (Fermi energy)] for Pn = As (Sb). Note that the result that Ti–Pn bonding is related to Fermi surface is consistent with the previous study by Singh. (31) The lack of overlap between Pn and O states indicates that no Pn–O bonding is formed in the compound (Figure 7). Figure 7 shows that Ti–O bonding does not change between As and Sb. On the other hand, there is a significant difference in pDOS between As and Sb, indicating that the nature of Ti–Pn bonding changes between them. To quantify the covalencies of Ti–Pn bonding, we calculated the center of masses according to eq 1. The results are shown in Figure 7 as solid lines. We adopted [ε2, ε1] = [− 4.9 eV, 0.0 eV] and [− 4.5 eV, 0.0 eV] for As and Sb, respectively. As a result, we obtained |CMTi3d – CMSb5p| = 0.77 eV and |CMTi3d – CMAs4p| = 0.68 eV. The smaller value of |CMTi3d – CMAs4p| indicates that Ti–As bonding is more covalent than Ti–Sb bonding. According to electronegativity, which is one of the simplest measures of covalency, Ti–Sb bonding seems more covalent than Ti–As one. Our results, however, show that the simple measure is invalid for the layered titanium-oxypnicitdes. In the compounds, the distance of Ti–Pn (the height of Pn toward Ti2O plane) governs the covalency as well as electronegativity. This is the same trend as in iron arsenide superconductors. To confirm our conclusion, we also calculated electron density of Na2Ti2Pn2O and visualize (0 2 0) plane including Ti–Pn bondings. Figure 8 shows that the electron density between Ti–Sb (Ti–As) bonding is 0.03–0.04 electron/bohr3 (0.04–0.05 electron/bohr3). The result indicates that Ti–As bonding is more covalent than Ti–Sb one again, which is consistent with the above pDOS analysis. Our analysis implies the correlation effect gets more enhanced for Sb than As in layered titanium-oxypnictides. The trend is again consistent with our previous study for BaTi2Pn2O (Pn = As, Sb, and Bi), where we showed that the phonon calculation using GGA-PBE could explain experiments only for Pn = As, but not for the other heavier Pn. (22) It might be a general tendency also applicable to the layered titanium-oxypnictides that the larger h enhances the electron correlation as the origin of the present discrepancy. Calculating the absolute values of electron correlation (e.g., effective Coulombic interaction, U) of the layered titanium-oxypnictides is an important future challenge.

Figure 7

Figure 7. Partial density of states and center of masses of Na2Ti2Pn2O (a) Pn = As, (b) Pn = Sb.

Figure 8

Figure 8. Electron densities of Na2Ti2Pn2O (a) Pn = As, (b) Pn = Sb along Ti–Pn bonds. The electron densities are plotted in the same range (0.03–0.07 electron/bohr3), and the aspects are consistent with each crystal structure.

Putting the discrepancy for Pn = Sb aside for a while, we can provide a plausible explanation for the origin of the lattice instability (Figure 3) from the above analysis: The lattice instability mainly derives from in-plane Ti displacements (Figure 4). To account for the interatomic forces to restore the lattice displacements, we can ignore (at least for the discussion of the instability) the contribution from Na because there are no bonds between Na and Ti (Figure 1). Starting with X mode for Pn = As and Sb, the imaginary frequency appears only for Pn = Sb. Figure 9a shows the change in potential energy accompanying the displacement of Ti atoms according to the phonon showing the imaginary (lowest) frequency at X point in Pn = Sb (As). The result shows that the original position is the most stable for Pn = As, but it is not for Pn = Sb. A double-well potential is realized for Pn = Sb. These phonons correspond to the Ti displacement toward Pn on the Ti2O plane as shown in Figure 9b. One of possible restoring forces would come from the Coulombic attraction from the valence electrons between Ti and O atoms, but the force is almost orthogonal to the displacement and hence we can ignore it. (42) In contrast, Pn is located in the same direction of Ti displacement, being responsible for the restoring forces. The displacement heads for the intermediate point of Pn located above and below the Ti2O plane (Figure 9c). Ti atom is surrounded by these Pn (Figure 9c), so forces in the z direction were canceled. However, the forces in the y direction could be responsible for the restoring forces (Figure 9c). To discuss the difference between As and Sb, we focus on the high electron density region surrounded by Ti and Pn atoms (Figure 8). When we focus on the Ti atom shown in Figure 9c, it approaches (recedes from) the left (right) high electron density region. The Ti atom feels the Coulombic repulsion from the left region, because another Ti atom on the opposite side simultaneously approaches the region (Figure 9b). They compress the high electron density region, resulting in the increase of the Coulomb energy. The magnitude of the high density region is 0.04 electron/bohr3 (0.03 electron/bohr3) for Pn = As (Sb). This implies that the Coulombic repulsion effect in Pn = As is larger than that in Pn = Sb. In contrast, the Ti atom feels Coulombic attraction from the right region. The effect of Coulombic attraction is again larger for Pn = As because of the higher electron density of the region. These are the reasons for the difference that the original position is the most stable for Pn = As, but it is not for Pn = Sb. For N mode, the imaginary frequency appears commonly for Pn = As and Sb. This is also accountable along the above Discussion section. In this mode, the Ti–Pn interaction has little effect on the restoring forces because Ti and Pn move in the same direction (Figure 4a). The magnitude of the Coulombic repulsion and attraction from the high electron density region are expected to be much smaller because Pn moves in the same direction of Ti. In this case, the possible restoring force would originate from the nearest oxygen; again, this is orthogonal to the displacement, giving too weak a contribution and leading to the instability. This is the reason why the imaginary frequencies appear at N mode for both Pn.

Figure 9

Figure 9. (a) Change in potential energy accompanying the displacement of Ti atoms according to the phonon showing the imaginary (lowest) frequency at X point in Pn = Sb (As). The horizontal axis corresponds to the magnitude of the displacement from the original positon. (b) Two-dimensional image of the displacement parttern of Ti atoms. (c) Three-dimensional image of the displacement parttern of the Ti atom enclosed in a red broken line in (b).

We return to the discussion of the discrepancy for Pn = Sb. If the discrepancy is attributed to the electron correlation as we suggest, then calculations using exchange-correlation functionals beyond GGA such as GGA+U and hybrid functionals are necessary to reproduce the experimentally observed superstructure of Pn = Sb within DFT framework. Though we could not perform further phonon calculations with GGA+U or hybrid functional due to the limitation of the implementation, we note that an interesting case has been reported for BaTi2Sb2O that only structural optimization with GGA+U could reproduce the observed transition from P4/mmm to Pmmm, otherwise without +U it cannot be explained. (43) GGA+U has also been applied to Na2Ti2Sb2O (32) to examine possible magnetic orderings, identifying a bicollinear antiferromagnetic ordering is the most stable. The ordering is unfortunately not able to solely account for the observed superstructure, (28) but when combined with a phonon evaluation under the supposed ordering, it might give further instability toward the enlarged superstructure, which is possibly consistent with the observed symmetry. A model Hamiltonian can be used to explicitly take the electron correlation effect into consideration. In fact, plausible explanations for the observed intra-unit-cell nematic CDW transition, (23) which could not be explained by the present phonon calculations, have been proposed by using a model Hamiltonian. (27, 43) The ideas are based on spin fluctuations to form magnetic orderings at low temperature, triggering another structural instability over the conventional prediction. Such a mechanism is not able to be described within the present calculation, (16, 22) and could be a possible origin for the discrepancy. For the key issue (i.e. spin fluctuations), there are still a small number of experiments (NMR and μSR) giving the controversy about the existence of the fluctuations. (17-19, 28, 44) Inelastic magnetic neutron scattering can be expected to provide conclusive information about this.

Conclusions

Click to copy section linkSection link copied!

We applied ab initio phonon analysis to layered titanium-oxypnictides, Na2Ti2Pn2O (Pn = As and Sb), and found a clear contrast between the cases with lighter/heavier pnictogen in comparison with experiments. The result completely explains the experimental structure at low temperature, C2/m for Pn = As, within the conventional charge density wave, while the calculation predicts that a Cmce polymorph is more stable than the experimentally observed superstructure (Cmcm) for Pn = Sb. To reveal the origin of discrepancy, we focus on the height of Pn toward Ti2O plane in layered titanium-oxypnictides. This is because the strength of electron correlation is well-captured by a trend of h, a vertical distance between Fe layer and Pn or Ch in iron arsenide superconductors that have similar electronic structures to our layered titanium-oxypnictides. On the basis of quantitative analysis, we concluded that Ti–As bonding is more covalent than Ti–Sb one and the distance of Ti–Pn (the height of Pn toward Ti2O plane) governs the covalency. The result indicates that the covalency of Ti–Pn bonding cannot be explained only by electronegativity in layered titanium-oxpnictides, which is the same trend as in iron arsenide superconductors. The result supports the similarity between iron arsenide and layered titanium-oxypnictide superconductors that has been proposed so far and probably stimulates further experimental and theoretical studies of layered titanium-oxypnictides. Our analysis could support the correlation effect gets more enhanced for Sb than for As in layered titanium-oxypnictides. Thus, we suggest the possibility that the GGA-PBE functional does not properly reproduce the electron correlation effects for Pn = Sb. This could cause the present discrepancy. If the discrepancy is attributed to the electron correlation as we suggest, then calculations using exchange-correlation functionals beyond GGA such as GGA+U and hybrid functionals are necessary to reproduce the experimentally observed superstructure of Pn = Sb. As shown in the present study, examining whether or not the observed superstructures are explained by ab initio phonon calculations can provide a clue to the mechanism of structural transition. The combination of diffraction studies and phonon calculations will be more important in the field of inorganic chemistry.

Computational Details

Click to copy section linkSection link copied!

All the calculations were done within DFT using the GGA-PBE exchange-correlation functional, (45) implemented in Quantum Espresso package. (46) We chose the PBE functional because it has been used for the layered titanium-oxypnictides so for and then no critical discrepancy has been reported, for example, compared to the ARPES measurement. (21) We adopted PAW (47) pseudopotentials in most calculations. The present PAW implementation takes into account the scalar relativistic effects. (48, 49) Only for spin–orbit case, we adopted Ultrasoft (US) pseudo potentials that take into account the fully relativistic effects. (48, 49) We restricted ourselves to spin unpolarized calculations. The valence electron configurations used for the US and PAW pseudopotentials are 2s22p63s1 (Na), 3s23p63d24s2 (Ti), 4s24p3 (As), 5s25p3 (Sb), and 2s22p4 (O). Lattice instabilities were detected by the negative (imaginary) phonon dispersions evaluated for undistorted and distorted structures. Taking each of the negative phonon modes, the structural relaxations along the mode were evaluated by the BFGS optimization scheme with the structural symmetries fixed to C2/m for Na2Ti2As2O and Cmce for Na2Ti2Sb2O. For phonon calculations, we used the linear response theory implemented in Quantum Espresso package. (50) Crystal structures and Fermi surfaces were depicted by using VESTA (51) and XCrySDen, (52) respectively.
To deal with all the compounds systematically, we checked the convergence of plane-wave cutoff energies (Ecut), k-meshes, q-meshes, and smearing parameters. The most strict condition among the compounds was taken to achieve the convergence within ±1.0 mRy per formula unit in the ground state energy, resulting in Ecut(WF) = 90 Ry for wave function and Ecut(ρ) = 800 Ry for charge density. For undistorted Na2Ti2As2O and Na2Ti2Sb2O, (6 × 6 × 6) k-meshes were used for the Brillouin zone integration. Phonon dispersions were calculated on (6 × 6 × 6) q-meshes. For Na2Ti2As2O (C2/m) and Na2Ti2Sb2O (Cmce), (4 × 4 × 4) k-meshes were used. For Na2Ti2Sb2O, (4 × 4 × 4) q-meshes were used to calculate phonon dispersions. We also calculated a ground state energy of the experimentally observed superstructure for Na2Ti2Sb2O (Cmcm) (28) with (6 × 6 × 3) k-meshes. The Marzari–Vanderbilt cold smearing scheme (53) with a broadening width of 0.01 Ry was applied to all the compounds. All the calculations were performed with primitive Brillouin zones.
The primitive Brillouin zone for undistorted Na2Ti2Pn2O (I4/mmm) is shown in Figure 10a. The special k and q points are Γ = (2π/a)·(0, 0, 0), N = (2π/a)·(1/2, 0, a/2c), P = (2π/a)·(1/2, 1/2, a/2c), X = (2π/a)·(1/2, 1/2, 0), Σ′ = (2π/a)·(1/2, 0, 0), and Z = (2π/a)·(0, 0, a/c) in Cartesian axis, where a and c are conventional lattice constants for the undistorted structures. The primitive Brillouin zone for Na2Ti2Sb2O superstructure (Cmce) is also shown in Figure 10b. The special k and q points are Γ = (2π/a′)·(0, 0, 0), X = (2π/a′)·(1, 0, 0), S = (2π/a′)·(1/2, a′/2b′, 0), Z = (2π/a′)·(0, 0, a′/2c′), A = (2π/a′)·(1, 0, a′/2c′), R = (2π/a′)·(1/2, a′/2b′, a′/2c′) in Cartesian axis, where a′, b′, and c′ are conventional lattice constants for the superstructure.

Figure 10

Figure 10. Primitive Brillouin zones (a) for undistorted Na2Ti2Pn2O structure (I4/mmm) and (b) for Na2Ti2Sb2O superstructure (Cmce).

Author Information

Click to copy section linkSection link copied!

  • Corresponding Authors
  • Author
    • Kenta Hongo - Research Center for Advanced Computing Infrastructure, JAIST, Asahidai 1-1, Nomi, Ishikawa 923-1292, JapanPRESTO, Japan Science and Technology Agency (JST), Kawaguchi, Saitama 332-0012, JapanCenter for Materials Research by Information Integration, Research and Services Division of Materials Data and Integrated System, National Institute for Materials Science, Tsukuba 305-0047, JapanOrcidhttp://orcid.org/0000-0002-2580-0907
  • Notes
    The authors declare no competing financial interest.

Acknowledgment

Click to copy section linkSection link copied!

The computation in this work has been performed using the facilities of the Research Center for Advanced Computing Infrastructure (RCACI) at JAIST. R.M. and K.H. are also grateful to MEXT-FLAGSHIP2020 (hp170269, hp170220) for their computational resources. K.H. is grateful for financial support from a KAKENHI grant (JP17K17762), a Grant-in-Aid for Scientific Research on Innovative Areas “Mixed Anion” project (JP16H06439) from MEXT, PRESTO (JPMJPR16NA) and the Materials research by Information Integration Initiative (MI2I) project of the Support Program for Starting Up Innovation Hub from Japan Science and Technology Agency (JST). R.M. is grateful to MEXT-KAKENHI (26287063, 17H05478), a Grant-in-Aid for Scientific Research on Innovative Areas "Mixed Anion" project (JP16H06440), a grant from the Asahi glass Foundation, and US-AFOSR-AOARD for their financial supports.

References

Click to copy section linkSection link copied!

This article references 54 other publications.

  1. 1
    Axtell, E. A.; Ozawa, T.; Kauzlarich, S. M.; Singh, R. R. Phase Transition and Spin-gap Behavior in a Layered Tetragonal Pnictide Oxide J. Solid State Chem. 1997, 134, 423 DOI: 10.1006/jssc.1997.7715
  2. 2
    Ozawa, T. C.; Kauzlarich, S. M.; Bieringer, M.; Greedan, J. E. Possible Charge-Density-Wave/Spin-Density-Wave in the Layered PnictideOxides:Na2Ti2Pn2O (Pn = As, Sb) Chem. Mater. 2001, 13, 1804 DOI: 10.1021/cm010009f
  3. 3
    Liu, R.; Song, Y.; Li, Q.; Ying, J.; Yan, Y.; He, Y.; Chen, X. Structure and Physical Properties of the Layered Pnictide-Oxides:(SrF)2Ti2Pn2O (Pn = As, Sb) and (SmO)2Ti2Sb2O Chem. Mater. 2010, 22, 1503 DOI: 10.1021/cm9027258
  4. 4
    Wang, X. F.; Yan, Y. J.; Ying, J. J.; Li, Q. J.; Zhang, M.; Xu, N.; Chen, X. H. Structure and physical properties for a new layered pnictide-oxide: BaTi2As2O J. Phys.: Condens. Matter 2010, 22, 075702 DOI: 10.1088/0953-8984/22/7/075702
  5. 5
    Doan, P.; Gooch, M.; Tang, Z.; Lorenz, B.; Moeller, A.; Tapp, J.; Chu, P. C. W.; Guloy, A. M. Ba1–xNaxTi2Sb2O (0.0 ≤ x ≤ 0.33) A Layered Titanium-Based Pnictide Oxide Superconductor J. Am. Chem. Soc. 2012, 134, 16520 DOI: 10.1021/ja3078889
  6. 6
    Yajima, T.; Nakano, K.; Takeiri, F.; Ono, T.; Hosokoshi, Y.; Matsushita, Y.; Hester, J.; Kobayashi, Y.; Kageyama, H. Superconductivity in BaTi2Sb2O with a d1 Square Lattice J. Phys. Soc. Jpn. 2012, 81, 103706 DOI: 10.1143/JPSJ.81.103706
  7. 7
    Yajima, T.; Nakano, K.; Takeiri, F.; Hester, J.; Yamamoto, T.; Kobayashi, Y.; Tsuji, N.; Kim, J.; Fujiwara, A.; Kageyama, H. Synthesis and Physical Properties of the New Oxybismuthides BaTi2Bi2O and (SrF)2Ti2Bi2O with a d1 Square Net J. Phys. Soc. Jpn. 2013, 82, 013703 DOI: 10.7566/JPSJ.82.013703
  8. 8
    Yajima, T.; Nakano, K.; Takeiri, F.; Nozaki, Y.; Kobayashi, Y.; Kageyama, H. Two Superconducting Phases in the Isovalent Solid Solutions BaTi2Pn2O (Pn = As, Sb, and Bi) J. Phys. Soc. Jpn. 2013, 82, 033705 DOI: 10.7566/JPSJ.82.033705
  9. 9
    Zhai, H.-F.; Jiao, W.-H.; Sun, Y.-L.; Bao, J.-K.; Jiang, H.; Yang, X.-J.; Tang, Z.-T.; Tao, Q.; Xu, X.-F.; Li, Y.-K.; Cao, C.; Dai, J.-H.; Xu, Z.-A.; Cao, G.-H. Superconductivity, charge- or spin-density wave, and metal-nonmetal transition in BaTi2(Sb1-xBix)2O Phys. Rev. B: Condens. Matter Mater. Phys. 2013, 87, 100502 DOI: 10.1103/PhysRevB.87.100502
  10. 10
    Nakano, K.; Yajima, T.; Takeiri, F.; Green, M. A.; Hester, J.; Kobayashi, Y.; Kageyama, H. Tc Enhancement by Aliovalent Anionic Substitution in Superconducting BaTi2(Sb1–xSnx)2O J. Phys. Soc. Jpn. 2013, 82, 074707 DOI: 10.7566/JPSJ.82.074707
  11. 11
    Pachmayr, U.; Johrendt, D. Superconductivity in Ba1-xKxTi2Sb2O (0 ≤ x ≤ 1) controlled by the layer charge Solid State Sci. 2014, 28, 31 DOI: 10.1016/j.solidstatesciences.2013.12.005
  12. 12
    von Rohr, F.; Nesper, R.; Schilling, A. Superconductivity in rubidium-substituted Ba1-xRbxTi2Sb2O Phys. Rev. B: Condens. Matter Mater. Phys. 2014, 89, 094505 DOI: 10.1103/PhysRevB.89.094505
  13. 13
    Ji, Q.; Ma, Y.; Hu, K.; Gao, B.; Mu, G.; Li, W.; Hu, T.; Zhang, G.; Zhao, Q.; Zhang, H.; Huang, F.; Xie, X. Synthesis, Structural, and Transport Properties of Cr-Doped BaTi2As2O Inorg. Chem. 2014, 53, 13089 DOI: 10.1021/ic502192h
  14. 14
    Bednorz, J.; Müller, K. Possible high Tc superconductivity in the Ba-La-Cu-O system Z. Phys. B: Condens. Matter 1986, 64, 189 DOI: 10.1007/BF01303701
  15. 15
    Kamihara, Y.; Watanabe, T.; Hirano, M.; Hosono, H. Iron-Based Layered Superconductor La[O1-xFx]FeAs (x = 0.05–0.12) with Tc = 26 K J. Am. Chem. Soc. 2008, 130, 3296 DOI: 10.1021/ja800073m
  16. 16
    Subedi, A. Electron-phonon superconductivity and charge density wave instability in the layered titanium-based pnictide BaTi2Sb2O Phys. Rev. B: Condens. Matter Mater. Phys. 2013, 87, 054506 DOI: 10.1103/PhysRevB.87.054506
  17. 17
    von Rohr, F.; Schilling, A.; Nesper, R.; Baines, C.; Bendele, M. Conventional superconductivity and charge-density-wave ordering in Ba1-xNaxTi2Sb2O Phys. Rev. B: Condens. Matter Mater. Phys. 2013, 88, 140501 DOI: 10.1103/PhysRevB.88.140501
  18. 18
    Kitagawa, S.; Ishida, K.; Nakano, K.; Yajima, T.; Kageyama, H. s-wave superconductivity in superconducting BaTi2Sb2O revealed by 121/123Sb-NMR/nuclear quadrupole resonance measurements Phys. Rev. B: Condens. Matter Mater. Phys. 2013, 87, 060510 DOI: 10.1103/PhysRevB.87.060510
  19. 19
    Nozaki, Y.; Nakano, K.; Yajima, T.; Kageyama, H.; Frandsen, B.; Liu, L.; Cheung, S.; Goko, T.; Uemura, Y. J.; Munsie, T. S. J.; Medina, T.; Luke, G. M.; Munevar, J.; Nishio-Hamane, D.; Brown, C. M. Muon spin relaxation and electron/neutron diffraction studies of BaTi2(As1-xSbx)2O: Absence of static magnetism and superlattice reflections Phys. Rev. B: Condens. Matter Mater. Phys. 2013, 88, 214506 DOI: 10.1103/PhysRevB.88.214506
  20. 20
    Gooch, M.; Doan, P.; Tang, Z.; Lorenz, B.; Guloy, A. M.; Chu, P. C. W. Weak coupling BCS-like superconductivity in the pnictide oxide Ba1-xNaxTi2Sb2O (x = 0 and 0.15) Phys. Rev. B: Condens. Matter Mater. Phys. 2013, 88, 064510 DOI: 10.1103/PhysRevB.88.064510
  21. 21
    Tan, S. Y.; Jiang, J.; Ye, Z. R.; Niu, X. H.; Song, Y.; Zhang, C. L.; Dai, P. C.; Xie, B. P.; Lai, X. C.; Feng, D. L. Photoemission study of the electronic structure and charge density waves of Na2Ti2Sb2O Sci. Rep. 2015, 5, 9515 DOI: 10.1038/srep09515
  22. 22
    Nakano, K.; Hongo, K.; Maezono, R. Phonon dispersions and Fermi surfaces nesting explaining the variety of charge ordering in titanium-oxypnictides superconductors Sci. Rep. 2016, 6, 29661 DOI: 10.1038/srep29661
  23. 23
    Frandsen, B. A.; Bozin, E. S.; Hu, H.; Zhu, Y.; Nozaki, Y.; Kageyama, H.; Uemura, Y. J.; Yin, W.-G.; Billinge, S. J. L. Intra-unit-cell nematic charge order in the titanium-oxypnictide family of superconductors Nat. Commun. 2014, 5, 5761 DOI: 10.1038/ncomms6761
  24. 24
    Lawler, M. J.; Fujita, K.; Lee, J.; Schmidt, A. R.; Kohsaka, Y.; Kim, C. K.; Eisaki, H.; Uchida, S.; Davis, J. C.; Sethna, J. P.; Kim, E.-A. Intra-unit-cell electronic nematicity of the high-Tc copper-oxide pseudogap states Nature 2010, 466, 347 DOI: 10.1038/nature09169
  25. 25
    Kasahara, S.; Shi, H.; Hashimoto, K.; Tonegawa, S.; Mizukami, Y.; Shibauchi, T.; Sugimoto, K.; Fukuda, T.; Terashima, T.; Nevidomskyy, A. H.; Matsuda, Y. Electronic nematicity above the structural and superconducting transition in BaFe2(As1-xPx)2 Nature 2012, 486, 382 DOI: 10.1038/nature11178
  26. 26
    Fujita, K.; Hamidian, M. H.; Edkins, S. D.; Kim, C. K.; Kohsaka, Y.; Azuma, M.; Takano, M.; Takagi, H.; Eisaki, H.; Uchida, S.-i.; Allais, A.; Lawler, M. J.; Kim, E.-A.; Sachdev, S.; Davis, J. C. S. Direct phase-sensitive identification of a d-form factor density wave in underdoped cuprates Proc. Natl. Acad. Sci. U. S. A. 2014, 111, E3026 DOI: 10.1073/pnas.1406297111
  27. 27
    Nakaoka, H.; Yamakawa, Y.; Kontani, H. Theoretical prediction of nematic orbital-ordered state in the Ti oxypnictide superconductor BaTi2(As, Sb)2O Phys. Rev. B: Condens. Matter Mater. Phys. 2016, 93, 245122 DOI: 10.1103/PhysRevB.93.245122
  28. 28
    Davies, N. R.; Johnson, R. D.; Princep, A. J.; Gannon, L. A.; Ma, J.-Z.; Qian, T.; Richard, P.; Li, H.; Shi, M.; Nowell, H. Coupled commensurate charge density wave and lattice distortion in Na2Ti2Pn2O (Pn = As, Sb) determined by X-ray diffraction and angle-resolved photoemission spectroscopy Phys. Rev. B: Condens. Matter Mater. Phys. 2016, 94, 104515 DOI: 10.1103/PhysRevB.94.104515
  29. 29
    Kuroki, K.; Usui, H.; Onari, S.; Arita, R.; Aoki, H. Pnictogen height as a possible switch between high-Tc nodeless and low-Tc nodal pairings in the iron-based superconductors Phys. Rev. B: Condens. Matter Mater. Phys. 2009, 79, 224511 DOI: 10.1103/PhysRevB.79.224511
  30. 30
    Miyake, T.; Nakamura, K.; Arita, R.; Imada, M. Comparison of ab initio low-energy models for LaFePO, LaFeAsO, BaFe2As2, LiFeAs, FeSe, and FeTe: electron correlation and covalency J. Phys. Soc. Jpn. 2010, 79, 044705 DOI: 10.1143/JPSJ.79.044705
  31. 31
    Singh, D. J. Electronic structure, disconnected Fermi surfaces and antiferromagnetism in the layered pnictide superconductor NaxBa1-xTi2Sb2O New J. Phys. 2012, 14, 123003 DOI: 10.1088/1367-2630/14/12/123003
  32. 32
    Yan, X.-W.; Lu, Z.-Y. Layered pnictide-oxide Na2Ti2Pn2O (Pn = As, Sb): a candidate for spin density waves J. Phys.: Condens. Matter 2013, 25, 365501 DOI: 10.1088/0953-8984/25/36/365501
  33. 33
    Suetin, D.; Ivanovskii, A. Structural, electronic properties, and chemical bonding in quaternary layered titanium pnictide-oxides Na2Ti2Pn2O and BaTi2Pn2O (Pn = As, Sb) from FLAPW-GGA calculations J. Alloys Compd. 2013, 564, 117 DOI: 10.1016/j.jallcom.2013.02.155
  34. 34
    Chen, D.; Zhang, T.-T.; Song, Z.-D.; Li, H.; Zhang, W.-L.; Qian, T.; Luo, J.-L.; Shi, Y.-G.; Fang, Z.; Richard, P.; Ding, H. New phase transition in Na2Ti2As2O revealed by Raman scattering Phys. Rev. B: Condens. Matter Mater. Phys. 2016, 93, 140501 DOI: 10.1103/PhysRevB.93.140501
  35. 35

    Note that we need to take special care when comparing our calculated geometry with an experimental one due to Davis et al. (28) We found that the atomic positions listed in their Table 1 are incompatible with their superstructure vectors, , , and , described in their main text: (28) (1) If the description in the main text is assumed to be correct, then Na and Sb positions in the superstructure are respectively located unlikely far from their undistorted positions. (2) O is not located at (0.0, 0.0) in their Figure 6a, being inconsistent with their Table 1. Accordingly, we speculate that they again define the superstructure vectors compatible with the atomic positions in their Table 1, which differs from that in their main text. We also found y = 0.0 for As3 in their Table 1 must be a typo and correctly y = 0.25 because As3 should occupy 8j Wyckoff site and its undistorted position is y = 0.5. Assuming our speculation is correct, our optimized geometry parameters agree well with their experimental values.

  36. 36
    Van De Walle, A.; Ceder, G. The effect of lattice vibrations on substitutional alloy thermodynamics Rev. Mod. Phys. 2002, 74, 11 DOI: 10.1103/RevModPhys.74.11
  37. 37
    Whangbo, M.-H.; Koo, H.-J.; Villesuzanne, A.; Pouchard, M. Effect of Metal-Oxygen Covalent Bonding on the Competition between Jahn-Teller Distortion and Charge Disproportionation in the Perovskites of High-Spin d4 Metal Ions LaMnO3 and CaFeO3 Inorg. Chem. 2002, 41, 1920 DOI: 10.1021/ic0110427
  38. 38
    Cammarata, A.; Rondinelli, J. M. Covalent dependence of octahedral rotations in orthorhombic perovskite oxides J. Chem. Phys. 2014, 141, 114704 DOI: 10.1063/1.4895967
  39. 39
    Hughbanks, T.; Hoffmann, R. Chains of trans-edge-sharing molybdenum octahedra: metal-metal bonding in extended systems J. Am. Chem. Soc. 1983, 105, 3528 DOI: 10.1021/ja00349a027
  40. 40
    Dronskowski, R.; Bloechl, P. E. Crystal orbital Hamilton populations (COHP): energy-resolved visualization of chemical bonding in solids based on density-functional calculations J. Phys. Chem. 1993, 97, 8617 DOI: 10.1021/j100135a014
  41. 41
    Grechnev, A.; Ahuja, R.; Eriksson, O. Balanced crystal orbital overlap population—a tool for analysing chemical bonds in solids J. Phys.: Condens. Matter 2003, 15, 7751 DOI: 10.1088/0953-8984/15/45/014
  42. 42

    The angle of O–Ti–Pn becomes 87° even at the most stable point for Pn = Sb. The corresponding displacement is 0.125 Å (Figure 9a). In this case, the restoring force between Ti and the nearest O contributing to the direction of the displacement is negligible because cos(87°) is only 0.05.

  43. 43
    Zhang, G.; Glasbrenner, J. K.; Flint, R.; Mazin, I. I.; Fernandes, R. M. Double-stage nematic bond ordering above double stripe magnetism: Application to BaTi2Sb2O Phys. Rev. B: Condens. Matter Mater. Phys. 2017, 95, 174402 DOI: 10.1103/PhysRevB.95.174402
  44. 44
    Fan, G.; Zhang, X.; Shi, Y.; Luo, J. Charge density wave transition in Na2Ti2Sb2O probed by 23Na NMR Sci. China: Phys., Mech. Astron. 2013, 56, 2399 DOI: 10.1007/s11433-013-5347-3
  45. 45
    Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple Phys. Rev. Lett. 1996, 77, 3865 DOI: 10.1103/PhysRevLett.77.3865
  46. 46
    Giannozzi, P. QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials J. Phys.: Condens. Matter 2009, 21, 395502 DOI: 10.1088/0953-8984/21/39/395502
  47. 47
    Blöchl, P. E. Projector augmented-wave method Phys. Rev. B: Condens. Matter Mater. Phys. 1994, 50, 17953 DOI: 10.1103/PhysRevB.50.17953
  48. 48
    Jollet, F.; Torrent, M.; Holzwarth, N. Generation of Projector Augmented-Wave atomic data: A 71 element validated table in the XML format Comput. Phys. Commun. 2014, 185, 1246 DOI: 10.1016/j.cpc.2013.12.023
  49. 49
    Kucukbenli, E.; Monni, M.; Adetunji, B.; Ge, X.; Adebayo, G.; Marzari, N.; de Gironcoli, S.; Corso, A. D. Projector augmented-wave and all-electron calculations across the periodic table: a comparison of structural and energetic properties. 2014, arXiv:cond-mat.mtrl-sci/1404.3015. arXiv.org e-Print archive. https://arxiv.org/abs/1404.3015 (accessed July 6, 2017).
  50. 50
    Baroni, S.; de Gironcoli, S.; Dal Corso, A.; Giannozzi, P. Phonons and related crystal properties from density-functional perturbation theory Rev. Mod. Phys. 2001, 73, 515 DOI: 10.1103/RevModPhys.73.515
  51. 51
    Momma, K.; Izumi, F. VESTA3 for three-dimensional visualization of crystal, volumetric and morphology data J. Appl. Crystallogr. 2011, 44, 1272 DOI: 10.1107/S0021889811038970
  52. 52
    Kokalj, A. XCrySDen a new program for displaying crystalline structures and electron densities J. Mol. Graphics Modell. 1999, 17, 176 DOI: 10.1016/S1093-3263(99)00028-5
  53. 53
    Marzari, N.; Vanderbilt, D.; De Vita, A.; Payne, M. C. Thermal Contraction and Disordering of the Al(110) Surface Phys. Rev. Lett. 1999, 82, 3296 DOI: 10.1103/PhysRevLett.82.3296
  54. 54
    Ozawa, T. C.; Pantoja, R.; Axtell, E. A.; Kauzlarich, S. M.; Greedan, J. E.; Bieringer, M.; Richardson, J. W. Powder Neutron Diffraction Studies of Na2Ti2Sb2O and Its Structure Property Relationships J. Solid State Chem. 2000, 153, 275 281 DOI: 10.1006/jssc.2000.8758

Cited By

Click to copy section linkSection link copied!
Citation Statements
Explore this article's citation statements on scite.ai

This article is cited by 5 publications.

  1. Genki Imam Prayogo, Andrea Tirelli, Keishu Utimula, Kenta Hongo, Ryo Maezono, Kousuke Nakano. Shry: Application of Canonical Augmentation to the Atomic Substitution Problem. Journal of Chemical Information and Modeling 2022, 62 (12) , 2909-2915. https://doi.org/10.1021/acs.jcim.2c00389
  2. Peng Song, Zhufeng Hou, Pedro Baptista de Castro, Kousuke Nakano, Kenta Hongo, Yoshihiko Takano, Ryo Maezono. High-Tc Superconducting Hydrides Formed by LaH24 and YH24 Cage Structures as Basic Blocks. Chemistry of Materials 2021, 33 (24) , 9501-9507. https://doi.org/10.1021/acs.chemmater.1c02371
  3. Keishu Utimula, Tom Ichibha, Ryo Maezono, Kenta Hongo. Ab Initio Search of Polymer Crystals with High Thermal Conductivity. Chemistry of Materials 2019, 31 (13) , 4649-4656. https://doi.org/10.1021/acs.chemmater.9b00020
  4. Shibghatullah Muhammady, Adelya S. Erlyanti, Rena Widita, Yudi Darma. Investigation of structural and electronic properties by pnictogen substitution in the layered oxypnictides (LaO)Zn Pn ( Pn = P, As, Sb). International Journal of Quantum Chemistry 2020, 120 (3) https://doi.org/10.1002/qua.26090
  5. Takuya Yasunaga, Makoto Kobayashi, Kenta Hongo, Kotaro Fujii, Shunsuke Yamamoto, Ryo Maezono, Masatomo Yashima, Masaya Mitsuishi, Hideki Kato, Masato Kakihana. Synthesis of Ba1−Sr YSi2O5N and discussion based on structure analysis and DFT calculation. Journal of Solid State Chemistry 2019, 276 , 266-271. https://doi.org/10.1016/j.jssc.2019.05.012

Inorganic Chemistry

Cite this: Inorg. Chem. 2017, 56, 22, 13732–13740
Click to copy citationCitation copied!
https://doi.org/10.1021/acs.inorgchem.7b01709
Published November 2, 2017

Copyright © 2017 American Chemical Society. This publication is licensed under these Terms of Use.

Article Views

1073

Altmetric

-

Citations

Learn about these metrics

Article Views are the COUNTER-compliant sum of full text article downloads since November 2008 (both PDF and HTML) across all institutions and individuals. These metrics are regularly updated to reflect usage leading up to the last few days.

Citations are the number of other articles citing this article, calculated by Crossref and updated daily. Find more information about Crossref citation counts.

The Altmetric Attention Score is a quantitative measure of the attention that a research article has received online. Clicking on the donut icon will load a page at altmetric.com with additional details about the score and the social media presence for the given article. Find more information on the Altmetric Attention Score and how the score is calculated.

  • Abstract

    Figure 1

    Figure 1. Crystal structure of Na2Ti2Pn2O (Pn = As and Sb). The space group is I4/mmm.

    Figure 2

    Figure 2. Electronic band structures for undistorted Na2Ti2Pn2O structure ((a) Pn = As, (b) Pn = Sb). The Fermi energy is set to 0 eV.

    Figure 3

    Figure 3. Phonon dispersions and phonon DOS for undistorted (I4/mmm) Na2Ti2Pn2O (Pn = (a) As, (b) Sb). The red cross marks in (b) represent the frequencies with spin–orbit coupling.

    Figure 4

    Figure 4. In-plane superstructures obtained by our calculations for Pn = (a) As and (b) Sb, and (c) the observed in-plane superstructure for Pn = Sb. (28) Pn is located above and below the Ti2O plane. Solid lines represent the original unit cells for I4/mmm. Thin dash lines represent the unit cells for 2 × 2 superstructures. Solid dash lines represent redefined unit cells for distorted superstructures. The conventional lattice vectors of the superstructure for (a) Pn = As (C2/m) are redefined as , , and , those for (b) Pn = Sb (Cmce) are redefined as , , and , and those for (c) Pn = Sb (Cmcm) are redefined as , , and , where , , and are conventional lattice vectors of the undistorted structures.

    Figure 5

    Figure 5. 3D superstructures obtained from our calculations for Pn = (a) As (C2/m) and (b) Sb (Cmce).

    Figure 6

    Figure 6. Phonon dispersions for Na2Ti2Sb2O superstructure (Cmce).

    Figure 7

    Figure 7. Partial density of states and center of masses of Na2Ti2Pn2O (a) Pn = As, (b) Pn = Sb.

    Figure 8

    Figure 8. Electron densities of Na2Ti2Pn2O (a) Pn = As, (b) Pn = Sb along Ti–Pn bonds. The electron densities are plotted in the same range (0.03–0.07 electron/bohr3), and the aspects are consistent with each crystal structure.

    Figure 9

    Figure 9. (a) Change in potential energy accompanying the displacement of Ti atoms according to the phonon showing the imaginary (lowest) frequency at X point in Pn = Sb (As). The horizontal axis corresponds to the magnitude of the displacement from the original positon. (b) Two-dimensional image of the displacement parttern of Ti atoms. (c) Three-dimensional image of the displacement parttern of the Ti atom enclosed in a red broken line in (b).

    Figure 10

    Figure 10. Primitive Brillouin zones (a) for undistorted Na2Ti2Pn2O structure (I4/mmm) and (b) for Na2Ti2Sb2O superstructure (Cmce).

  • References


    This article references 54 other publications.

    1. 1
      Axtell, E. A.; Ozawa, T.; Kauzlarich, S. M.; Singh, R. R. Phase Transition and Spin-gap Behavior in a Layered Tetragonal Pnictide Oxide J. Solid State Chem. 1997, 134, 423 DOI: 10.1006/jssc.1997.7715
    2. 2
      Ozawa, T. C.; Kauzlarich, S. M.; Bieringer, M.; Greedan, J. E. Possible Charge-Density-Wave/Spin-Density-Wave in the Layered PnictideOxides:Na2Ti2Pn2O (Pn = As, Sb) Chem. Mater. 2001, 13, 1804 DOI: 10.1021/cm010009f
    3. 3
      Liu, R.; Song, Y.; Li, Q.; Ying, J.; Yan, Y.; He, Y.; Chen, X. Structure and Physical Properties of the Layered Pnictide-Oxides:(SrF)2Ti2Pn2O (Pn = As, Sb) and (SmO)2Ti2Sb2O Chem. Mater. 2010, 22, 1503 DOI: 10.1021/cm9027258
    4. 4
      Wang, X. F.; Yan, Y. J.; Ying, J. J.; Li, Q. J.; Zhang, M.; Xu, N.; Chen, X. H. Structure and physical properties for a new layered pnictide-oxide: BaTi2As2O J. Phys.: Condens. Matter 2010, 22, 075702 DOI: 10.1088/0953-8984/22/7/075702
    5. 5
      Doan, P.; Gooch, M.; Tang, Z.; Lorenz, B.; Moeller, A.; Tapp, J.; Chu, P. C. W.; Guloy, A. M. Ba1–xNaxTi2Sb2O (0.0 ≤ x ≤ 0.33) A Layered Titanium-Based Pnictide Oxide Superconductor J. Am. Chem. Soc. 2012, 134, 16520 DOI: 10.1021/ja3078889
    6. 6
      Yajima, T.; Nakano, K.; Takeiri, F.; Ono, T.; Hosokoshi, Y.; Matsushita, Y.; Hester, J.; Kobayashi, Y.; Kageyama, H. Superconductivity in BaTi2Sb2O with a d1 Square Lattice J. Phys. Soc. Jpn. 2012, 81, 103706 DOI: 10.1143/JPSJ.81.103706
    7. 7
      Yajima, T.; Nakano, K.; Takeiri, F.; Hester, J.; Yamamoto, T.; Kobayashi, Y.; Tsuji, N.; Kim, J.; Fujiwara, A.; Kageyama, H. Synthesis and Physical Properties of the New Oxybismuthides BaTi2Bi2O and (SrF)2Ti2Bi2O with a d1 Square Net J. Phys. Soc. Jpn. 2013, 82, 013703 DOI: 10.7566/JPSJ.82.013703
    8. 8
      Yajima, T.; Nakano, K.; Takeiri, F.; Nozaki, Y.; Kobayashi, Y.; Kageyama, H. Two Superconducting Phases in the Isovalent Solid Solutions BaTi2Pn2O (Pn = As, Sb, and Bi) J. Phys. Soc. Jpn. 2013, 82, 033705 DOI: 10.7566/JPSJ.82.033705
    9. 9
      Zhai, H.-F.; Jiao, W.-H.; Sun, Y.-L.; Bao, J.-K.; Jiang, H.; Yang, X.-J.; Tang, Z.-T.; Tao, Q.; Xu, X.-F.; Li, Y.-K.; Cao, C.; Dai, J.-H.; Xu, Z.-A.; Cao, G.-H. Superconductivity, charge- or spin-density wave, and metal-nonmetal transition in BaTi2(Sb1-xBix)2O Phys. Rev. B: Condens. Matter Mater. Phys. 2013, 87, 100502 DOI: 10.1103/PhysRevB.87.100502
    10. 10
      Nakano, K.; Yajima, T.; Takeiri, F.; Green, M. A.; Hester, J.; Kobayashi, Y.; Kageyama, H. Tc Enhancement by Aliovalent Anionic Substitution in Superconducting BaTi2(Sb1–xSnx)2O J. Phys. Soc. Jpn. 2013, 82, 074707 DOI: 10.7566/JPSJ.82.074707
    11. 11
      Pachmayr, U.; Johrendt, D. Superconductivity in Ba1-xKxTi2Sb2O (0 ≤ x ≤ 1) controlled by the layer charge Solid State Sci. 2014, 28, 31 DOI: 10.1016/j.solidstatesciences.2013.12.005
    12. 12
      von Rohr, F.; Nesper, R.; Schilling, A. Superconductivity in rubidium-substituted Ba1-xRbxTi2Sb2O Phys. Rev. B: Condens. Matter Mater. Phys. 2014, 89, 094505 DOI: 10.1103/PhysRevB.89.094505
    13. 13
      Ji, Q.; Ma, Y.; Hu, K.; Gao, B.; Mu, G.; Li, W.; Hu, T.; Zhang, G.; Zhao, Q.; Zhang, H.; Huang, F.; Xie, X. Synthesis, Structural, and Transport Properties of Cr-Doped BaTi2As2O Inorg. Chem. 2014, 53, 13089 DOI: 10.1021/ic502192h
    14. 14
      Bednorz, J.; Müller, K. Possible high Tc superconductivity in the Ba-La-Cu-O system Z. Phys. B: Condens. Matter 1986, 64, 189 DOI: 10.1007/BF01303701
    15. 15
      Kamihara, Y.; Watanabe, T.; Hirano, M.; Hosono, H. Iron-Based Layered Superconductor La[O1-xFx]FeAs (x = 0.05–0.12) with Tc = 26 K J. Am. Chem. Soc. 2008, 130, 3296 DOI: 10.1021/ja800073m
    16. 16
      Subedi, A. Electron-phonon superconductivity and charge density wave instability in the layered titanium-based pnictide BaTi2Sb2O Phys. Rev. B: Condens. Matter Mater. Phys. 2013, 87, 054506 DOI: 10.1103/PhysRevB.87.054506
    17. 17
      von Rohr, F.; Schilling, A.; Nesper, R.; Baines, C.; Bendele, M. Conventional superconductivity and charge-density-wave ordering in Ba1-xNaxTi2Sb2O Phys. Rev. B: Condens. Matter Mater. Phys. 2013, 88, 140501 DOI: 10.1103/PhysRevB.88.140501
    18. 18
      Kitagawa, S.; Ishida, K.; Nakano, K.; Yajima, T.; Kageyama, H. s-wave superconductivity in superconducting BaTi2Sb2O revealed by 121/123Sb-NMR/nuclear quadrupole resonance measurements Phys. Rev. B: Condens. Matter Mater. Phys. 2013, 87, 060510 DOI: 10.1103/PhysRevB.87.060510
    19. 19
      Nozaki, Y.; Nakano, K.; Yajima, T.; Kageyama, H.; Frandsen, B.; Liu, L.; Cheung, S.; Goko, T.; Uemura, Y. J.; Munsie, T. S. J.; Medina, T.; Luke, G. M.; Munevar, J.; Nishio-Hamane, D.; Brown, C. M. Muon spin relaxation and electron/neutron diffraction studies of BaTi2(As1-xSbx)2O: Absence of static magnetism and superlattice reflections Phys. Rev. B: Condens. Matter Mater. Phys. 2013, 88, 214506 DOI: 10.1103/PhysRevB.88.214506
    20. 20
      Gooch, M.; Doan, P.; Tang, Z.; Lorenz, B.; Guloy, A. M.; Chu, P. C. W. Weak coupling BCS-like superconductivity in the pnictide oxide Ba1-xNaxTi2Sb2O (x = 0 and 0.15) Phys. Rev. B: Condens. Matter Mater. Phys. 2013, 88, 064510 DOI: 10.1103/PhysRevB.88.064510
    21. 21
      Tan, S. Y.; Jiang, J.; Ye, Z. R.; Niu, X. H.; Song, Y.; Zhang, C. L.; Dai, P. C.; Xie, B. P.; Lai, X. C.; Feng, D. L. Photoemission study of the electronic structure and charge density waves of Na2Ti2Sb2O Sci. Rep. 2015, 5, 9515 DOI: 10.1038/srep09515
    22. 22
      Nakano, K.; Hongo, K.; Maezono, R. Phonon dispersions and Fermi surfaces nesting explaining the variety of charge ordering in titanium-oxypnictides superconductors Sci. Rep. 2016, 6, 29661 DOI: 10.1038/srep29661
    23. 23
      Frandsen, B. A.; Bozin, E. S.; Hu, H.; Zhu, Y.; Nozaki, Y.; Kageyama, H.; Uemura, Y. J.; Yin, W.-G.; Billinge, S. J. L. Intra-unit-cell nematic charge order in the titanium-oxypnictide family of superconductors Nat. Commun. 2014, 5, 5761 DOI: 10.1038/ncomms6761
    24. 24
      Lawler, M. J.; Fujita, K.; Lee, J.; Schmidt, A. R.; Kohsaka, Y.; Kim, C. K.; Eisaki, H.; Uchida, S.; Davis, J. C.; Sethna, J. P.; Kim, E.-A. Intra-unit-cell electronic nematicity of the high-Tc copper-oxide pseudogap states Nature 2010, 466, 347 DOI: 10.1038/nature09169
    25. 25
      Kasahara, S.; Shi, H.; Hashimoto, K.; Tonegawa, S.; Mizukami, Y.; Shibauchi, T.; Sugimoto, K.; Fukuda, T.; Terashima, T.; Nevidomskyy, A. H.; Matsuda, Y. Electronic nematicity above the structural and superconducting transition in BaFe2(As1-xPx)2 Nature 2012, 486, 382 DOI: 10.1038/nature11178
    26. 26
      Fujita, K.; Hamidian, M. H.; Edkins, S. D.; Kim, C. K.; Kohsaka, Y.; Azuma, M.; Takano, M.; Takagi, H.; Eisaki, H.; Uchida, S.-i.; Allais, A.; Lawler, M. J.; Kim, E.-A.; Sachdev, S.; Davis, J. C. S. Direct phase-sensitive identification of a d-form factor density wave in underdoped cuprates Proc. Natl. Acad. Sci. U. S. A. 2014, 111, E3026 DOI: 10.1073/pnas.1406297111
    27. 27
      Nakaoka, H.; Yamakawa, Y.; Kontani, H. Theoretical prediction of nematic orbital-ordered state in the Ti oxypnictide superconductor BaTi2(As, Sb)2O Phys. Rev. B: Condens. Matter Mater. Phys. 2016, 93, 245122 DOI: 10.1103/PhysRevB.93.245122
    28. 28
      Davies, N. R.; Johnson, R. D.; Princep, A. J.; Gannon, L. A.; Ma, J.-Z.; Qian, T.; Richard, P.; Li, H.; Shi, M.; Nowell, H. Coupled commensurate charge density wave and lattice distortion in Na2Ti2Pn2O (Pn = As, Sb) determined by X-ray diffraction and angle-resolved photoemission spectroscopy Phys. Rev. B: Condens. Matter Mater. Phys. 2016, 94, 104515 DOI: 10.1103/PhysRevB.94.104515
    29. 29
      Kuroki, K.; Usui, H.; Onari, S.; Arita, R.; Aoki, H. Pnictogen height as a possible switch between high-Tc nodeless and low-Tc nodal pairings in the iron-based superconductors Phys. Rev. B: Condens. Matter Mater. Phys. 2009, 79, 224511 DOI: 10.1103/PhysRevB.79.224511
    30. 30
      Miyake, T.; Nakamura, K.; Arita, R.; Imada, M. Comparison of ab initio low-energy models for LaFePO, LaFeAsO, BaFe2As2, LiFeAs, FeSe, and FeTe: electron correlation and covalency J. Phys. Soc. Jpn. 2010, 79, 044705 DOI: 10.1143/JPSJ.79.044705
    31. 31
      Singh, D. J. Electronic structure, disconnected Fermi surfaces and antiferromagnetism in the layered pnictide superconductor NaxBa1-xTi2Sb2O New J. Phys. 2012, 14, 123003 DOI: 10.1088/1367-2630/14/12/123003
    32. 32
      Yan, X.-W.; Lu, Z.-Y. Layered pnictide-oxide Na2Ti2Pn2O (Pn = As, Sb): a candidate for spin density waves J. Phys.: Condens. Matter 2013, 25, 365501 DOI: 10.1088/0953-8984/25/36/365501
    33. 33
      Suetin, D.; Ivanovskii, A. Structural, electronic properties, and chemical bonding in quaternary layered titanium pnictide-oxides Na2Ti2Pn2O and BaTi2Pn2O (Pn = As, Sb) from FLAPW-GGA calculations J. Alloys Compd. 2013, 564, 117 DOI: 10.1016/j.jallcom.2013.02.155
    34. 34
      Chen, D.; Zhang, T.-T.; Song, Z.-D.; Li, H.; Zhang, W.-L.; Qian, T.; Luo, J.-L.; Shi, Y.-G.; Fang, Z.; Richard, P.; Ding, H. New phase transition in Na2Ti2As2O revealed by Raman scattering Phys. Rev. B: Condens. Matter Mater. Phys. 2016, 93, 140501 DOI: 10.1103/PhysRevB.93.140501
    35. 35

      Note that we need to take special care when comparing our calculated geometry with an experimental one due to Davis et al. (28) We found that the atomic positions listed in their Table 1 are incompatible with their superstructure vectors, , , and , described in their main text: (28) (1) If the description in the main text is assumed to be correct, then Na and Sb positions in the superstructure are respectively located unlikely far from their undistorted positions. (2) O is not located at (0.0, 0.0) in their Figure 6a, being inconsistent with their Table 1. Accordingly, we speculate that they again define the superstructure vectors compatible with the atomic positions in their Table 1, which differs from that in their main text. We also found y = 0.0 for As3 in their Table 1 must be a typo and correctly y = 0.25 because As3 should occupy 8j Wyckoff site and its undistorted position is y = 0.5. Assuming our speculation is correct, our optimized geometry parameters agree well with their experimental values.

    36. 36
      Van De Walle, A.; Ceder, G. The effect of lattice vibrations on substitutional alloy thermodynamics Rev. Mod. Phys. 2002, 74, 11 DOI: 10.1103/RevModPhys.74.11
    37. 37
      Whangbo, M.-H.; Koo, H.-J.; Villesuzanne, A.; Pouchard, M. Effect of Metal-Oxygen Covalent Bonding on the Competition between Jahn-Teller Distortion and Charge Disproportionation in the Perovskites of High-Spin d4 Metal Ions LaMnO3 and CaFeO3 Inorg. Chem. 2002, 41, 1920 DOI: 10.1021/ic0110427
    38. 38
      Cammarata, A.; Rondinelli, J. M. Covalent dependence of octahedral rotations in orthorhombic perovskite oxides J. Chem. Phys. 2014, 141, 114704 DOI: 10.1063/1.4895967
    39. 39
      Hughbanks, T.; Hoffmann, R. Chains of trans-edge-sharing molybdenum octahedra: metal-metal bonding in extended systems J. Am. Chem. Soc. 1983, 105, 3528 DOI: 10.1021/ja00349a027
    40. 40
      Dronskowski, R.; Bloechl, P. E. Crystal orbital Hamilton populations (COHP): energy-resolved visualization of chemical bonding in solids based on density-functional calculations J. Phys. Chem. 1993, 97, 8617 DOI: 10.1021/j100135a014
    41. 41
      Grechnev, A.; Ahuja, R.; Eriksson, O. Balanced crystal orbital overlap population—a tool for analysing chemical bonds in solids J. Phys.: Condens. Matter 2003, 15, 7751 DOI: 10.1088/0953-8984/15/45/014
    42. 42

      The angle of O–Ti–Pn becomes 87° even at the most stable point for Pn = Sb. The corresponding displacement is 0.125 Å (Figure 9a). In this case, the restoring force between Ti and the nearest O contributing to the direction of the displacement is negligible because cos(87°) is only 0.05.

    43. 43
      Zhang, G.; Glasbrenner, J. K.; Flint, R.; Mazin, I. I.; Fernandes, R. M. Double-stage nematic bond ordering above double stripe magnetism: Application to BaTi2Sb2O Phys. Rev. B: Condens. Matter Mater. Phys. 2017, 95, 174402 DOI: 10.1103/PhysRevB.95.174402
    44. 44
      Fan, G.; Zhang, X.; Shi, Y.; Luo, J. Charge density wave transition in Na2Ti2Sb2O probed by 23Na NMR Sci. China: Phys., Mech. Astron. 2013, 56, 2399 DOI: 10.1007/s11433-013-5347-3
    45. 45
      Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple Phys. Rev. Lett. 1996, 77, 3865 DOI: 10.1103/PhysRevLett.77.3865
    46. 46
      Giannozzi, P. QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials J. Phys.: Condens. Matter 2009, 21, 395502 DOI: 10.1088/0953-8984/21/39/395502
    47. 47
      Blöchl, P. E. Projector augmented-wave method Phys. Rev. B: Condens. Matter Mater. Phys. 1994, 50, 17953 DOI: 10.1103/PhysRevB.50.17953
    48. 48
      Jollet, F.; Torrent, M.; Holzwarth, N. Generation of Projector Augmented-Wave atomic data: A 71 element validated table in the XML format Comput. Phys. Commun. 2014, 185, 1246 DOI: 10.1016/j.cpc.2013.12.023
    49. 49
      Kucukbenli, E.; Monni, M.; Adetunji, B.; Ge, X.; Adebayo, G.; Marzari, N.; de Gironcoli, S.; Corso, A. D. Projector augmented-wave and all-electron calculations across the periodic table: a comparison of structural and energetic properties. 2014, arXiv:cond-mat.mtrl-sci/1404.3015. arXiv.org e-Print archive. https://arxiv.org/abs/1404.3015 (accessed July 6, 2017).
    50. 50
      Baroni, S.; de Gironcoli, S.; Dal Corso, A.; Giannozzi, P. Phonons and related crystal properties from density-functional perturbation theory Rev. Mod. Phys. 2001, 73, 515 DOI: 10.1103/RevModPhys.73.515
    51. 51
      Momma, K.; Izumi, F. VESTA3 for three-dimensional visualization of crystal, volumetric and morphology data J. Appl. Crystallogr. 2011, 44, 1272 DOI: 10.1107/S0021889811038970
    52. 52
      Kokalj, A. XCrySDen a new program for displaying crystalline structures and electron densities J. Mol. Graphics Modell. 1999, 17, 176 DOI: 10.1016/S1093-3263(99)00028-5
    53. 53
      Marzari, N.; Vanderbilt, D.; De Vita, A.; Payne, M. C. Thermal Contraction and Disordering of the Al(110) Surface Phys. Rev. Lett. 1999, 82, 3296 DOI: 10.1103/PhysRevLett.82.3296
    54. 54
      Ozawa, T. C.; Pantoja, R.; Axtell, E. A.; Kauzlarich, S. M.; Greedan, J. E.; Bieringer, M.; Richardson, J. W. Powder Neutron Diffraction Studies of Na2Ti2Sb2O and Its Structure Property Relationships J. Solid State Chem. 2000, 153, 275 281 DOI: 10.1006/jssc.2000.8758