# Phonons and Thermal Expansion Behavior of NiSi and NiGe

^{1}Solid State Physics Division, Bhabha Atomic Research Centre, Mumbai, India^{2}Homi Bhabha National Institute, Mumbai, India

We have carried out first principles calculations of the vibrational and thermodynamic behavior in NiSi and isostructural compound NiGe. Phonon density of states has also been measured in NiSi using inelastic neutron scattering techniques. We find that the vibrational spectra of the two compounds are very different, due to the difference in the size and mass of Si and Ge. Interesting anomalous thermal behavior of NiSi due to anharmonic phonons is brought out well in our calculations, particularly the negative thermal expansion (NTE) along the b-axis of the orthorhombic unit cell. Large difference in thermal expansion behavior of NiSi and NiGe is very well reproduced by the calculations. Additionally, calculations enable to identify the phonon modes which lend major contribution to the negative thermal expansion behavior in NiSi, and reasons for negligible NTE in NiGe. Such typical representative modes at the zone-boundary along b-axis involve transverse vibrations of Si/Ge along c-axis.

**PACS numbers:** 78.70.Nx, 63.20.-e, 65.40.–b.

## Introduction

Transition metal silicides are considered important in microelectronics industry. These materials are used to reduce contact resistance at source, drain and gate terminals in standard complementary metal oxide semiconductor (CMOS) devices and very-large-scale integration (VLSI) technology. There are significant experimental and theoretical studies (Wu et al., 2005; Kim, 2015; Lin et al., 2015; Parshin et al., 2016; Stern and Madsen, 2016; Tang et al., 2017; Sidorova et al., 2018; Zhang et al., 2018a,b) on silicides to investigate their electronic and transport properties. Si-Ge alloys are the backbone of the metal-oxide semiconductor industry (Von Känel et al., 1992; Dahal et al., 2016; Lee et al., 2017; Simon et al., 2017; Kim et al., 2018). Lower resistivity, lower silicon consumption and higher stability make them better alternative to conventional silicides, like TiSi_{2} and CoSi_{2}, for CMOS devices. Apart from its electronic industry applications, the compounds are reported to be good candidates for infrared sensor and ultraviolet photo detector applications. Detailed studies on the thermal behavior has been carried out on NiSiGe (Perrin et al., 2007) and FeNiSi (Simon et al., 2017) non-stoichiometric compositions, for fruitful integration of these derived compounds in the field of microelectronics.

NiSi (Franciosi et al., 1982; Lee et al., 2000; Detavernier et al., 2003; Jin et al., 2013), FeSi (Damascelli et al., 1997a,b), CoSi (Acun and Soyalp, 2012), and similar compounds (Kim, 2015) have been extensively studied with above applications in mind, and also for their various unusual physical and chemical properties (Weber et al., 2006; Lin et al., 2015; Geenen et al., 2018). NiSi (Connétable and Thomas, 2009) is structurally similar to NiGe. FeSi, a fascinating material, has been studied for its unusual magnetic and thermodynamic properties like thermal expansion, elastic behavior, heat capacity, Seebeck coefficients, magnetic susceptibility, insulator-metal transition with increasing temperature (Mandrus et al., 1995; Caracas and Wentzcovitch, 2004; Delaire et al., 2011; Parshin et al., 2016; Stern and Madsen, 2016). It is also described as a strongly correlated system, a possible Kondo insulator. Several electronic structure calculations have been reported on FeSi to understand its anomalous thermal expansion and to elucidate a link between its electronic density of states with the thermal disorder (Vočadlo et al., 2002). Phonon density of states have been studied in FeSi (Damascelli et al., 1997a,b; Racu et al., 2008), both experimentally and using *ab initio* molecular dynamics, to investigate softening of phonons (Delaire et al., 2013) and coupling between phonons and electronic structure in order to explain anomalous thermal expansion. Equation of state (Sarrao et al., 1994; Dobson et al., 2002) and elasticity of FeSi (Vočadlo et al., 2012) have been studied extensively to gain better understanding of the earth's core, which is mainly comprised of an alloy of iron and nickel, with a few light alloying elements like Si, C, S, H, and O.

Although NiSi and FeSi are not structurally similar at ambient conditions, their thermal expansion behavior is comparable. A cumulative study of the phonon density of states in NiSi using neutrons and first principles lattice dynamics calculations of thermal expansion would fill the gaps in complete understanding of these monosilicides and monogermanides. NiSi exhibits unusual thermal anisotropy (Donthu et al., 2004) with an expansion along the **a-** and **c-** direction while along **b-** direction there is significant contraction with increasing temperature. Although the overall thermal expansion coefficient has a positive value, yet the behavior along **b** direction makes it an interesting study. This behavior is consistent in both bulk and in thin films. In case of NiGe, there is a contraction along b- axis but it is not as large as in NiSi and FeSi. It would be very interesting to understand the differences in the behavior of NiSi, FeSi, and NiGe. The calculation of the elastic properties of NiSi and NiGe in comparison to the reported studies on FeSi would be useful to know the strength of these materials. High pressure theoretical studies on NiSi up to 500 GPa have already been reported (Lord et al., 2012; Vočadlo et al., 2012) for unraveling the mysteries surrounding the composition and behavior of earth's core. NiSi_{x}Ge_{(1−x)} compounds are found to be promising in future microelectronics industry.

The crystal structure of NiSi and NiGe at ambient temperature is found to be orthorhombic (Pauling and Soldate, 1948), with space group Pnma. There are four Ni and four Si/Ge atoms in a unit cell. Although Ni is ferromagnetic, NiSi shows no magnetic ordering even up to lowest temperature around 0 K (Dahal et al., 2016). First principle studies to understand its high pressure behavior have been extensively reported (Jin et al., 2013). Anisotropy of thermal expansion in NiSi thin films have been studied (Perrin et al., 2007). Molar specific heat up to 400 K has been measured (Acker et al., 1999). Transport properties and neutron scattering to probe its magnetic ordering at 0.48 K have been reported. FeSi has a cubic P2_{1}3 (B20) structure at ambient conditions. It transforms to the CsCl structure (Pm3m) at around 24 GPa and beyond 1950 K. NiSi adopts a B20 structure above 10 GPa and at around 900 K. It finally adopts a CsCl structure above 46 GPa and 1900 K.

This paper is a combined study of the vibrational properties of ambient phase of NiSi using neutron inelastic scattering experiments and theoretical calculations to understand the anomalous behavior of the thermal contraction along b-axis. First principles calculations have been carried out to derive the phonon density of states and thermal expansion behavior is obtained using *ab initio* lattice dynamics simulations. We also performed the calculation in NiGe. These studies will throw light on the specific phonon modes which are responsible for the anisotropic thermal behavior in NiSi, and the changes which come into play as we move from NiSi to NiGe. This is particularly interesting as the anisotropic thermal behavior is much less pronounced in NiGe.

## Experimental and Computational Details

Inelastic neutron scattering to study the phonon density of states in polycrystalline NiSi sample has been carried out using the Triple Axis Spectrometer at Dhruva reactor, Trombay. A sample of about 20 g of NiSi polycrystalline powder was used. All the measurements are carried out in the energy loss mode with constant momentum transfer (** Q**) geometry over a range of

*Q*-values. The energy of the analyzer in various scans was kept fixed at 25 and 34.4 meV. In the incoherent one-phonon approximation (Carpenter and Price, 1985; Price and Skold, 1986), the measured scattering function

*S*(

*Q*,

*E*), as observed in the neutron experiments, is related to the phonon density of states g

^{(n)(}E) as follows:

Where the + or – signs correspond to energy loss or gain of the neutrons respectively and $n(E,T)={\left[exp(E/{k}_{B}T)-1\right]}^{-1}$. *A* and *B* are normalization constants.

The observed neutron-weighted phonon density of states is a sum of the partial components of the density of states due to the various atoms, weighted by their scattering length squares.

Here *b*_{k}, *m*_{k}, and *g*_{k}(*E*) are, respectively, the neutron scattering length, mass, and partial density of states of the *k*^{th} atom in the unit cell. The quantity between < > represents suitable average over all *Q* values at a given energy. 2*W*(*Q*) is the Debye-Waller factor averaged over all the atoms. The weighting factors $\frac{4\pi {b}_{k}^{2}}{{m}_{k}}$ for the various atoms in the units of barns/amu are: Ni: 0.331, and Si: 0.077.

The density functional theory calculations using VASP (Kresse and Furthmüller, 1996; Kresse and Joubert, 1999) along with PHONON software (Parlinksi, PHONON Software, 2003) is used to calculate the phonon frequencies of NiSi and NiGe in entire Brillouin zone. The software requires the forces on various atoms after displacement along various directions. These are obtained by computing the Hellman-Feynman forces on various atom in various configuration in a supercell with various (±x, ±y, ±z) atomic displacement patterns. The Hellman-Feynman forces on the atoms in the super cell have been calculated using density functional theory as implemented in the VASP software (Kresse and Furthmüller, 1996; Kresse and Joubert, 1999). NiSi and NiGe in orthorhombic phase (space group 62- Pnma) have 2 symmetrically inequivalent atoms in the unit cell, and total number of atoms in the primitive unit cell is 8. The force calculations is performed by displacing the symmetrically inequivalent atoms along (±x, ±y, ±z) directions. The energy cutoff is 1,000 eV and a 16 × 24 × 16 k point mess has been used to obtain energy convergence in total energy of the order of meV, which is sufficient to obtain the required accuracy in phonon energies. The Monkhorst Pack method is used for k point sampling (Monkhorst and Pack, 1976). The exchange-correlation contributions were included within PBE generalized gradient approximation (GGA) (Perdew et al., 1996). The structure is optimized with atomic pseudopotential with d^{8}s^{2}p^{2}, s^{2}p^{2} and s^{2}p^{2} valence electrons in Ni, Si, and Ge atom respectively. The convergence criteria for the total energy and ionic forces were set to 10^{−8} eV and 10^{−3} eV Å^{−1}, respectively.

## Results and Discussion

### Phonon Spectra

The calculated total density of states in NiSi and NiGe along with respective partial contributions of various atoms is given in Figure 1. Replacing lighter Si in NiSi with Ge brings about a change in the lattice parameters and in the bond lengths. The lattice parameter as well as volume of NiGe is larger in comparison to that of NiSi. This is reflected in the phonon density of states in the two compounds. In case of NiSi, the phonon spectra extend to 50 meV, while for NiGe it extends only up to 34 meV. The spectral range of the partial contribution of Ni is different in both the compounds. Ge (72.63 amu) is heavier in comparison to Si (28.095 amu). The low energy peak in the partial density of states of Ge in NiGe is at 8 meV in comparison to that at 16 meV in Si spectrum of NiSi. The shift in the peak position may be mainly due to the mass difference of both the atoms in respective compounds. Further, Ni-Ge bond lengths are longer in comparison to that of Ni-Si and hence the phonon frequencies are lower. As a result, the span of the vibrational spectrum of NiGe is up to 34 meV is lower in comparison to that of 50 meV in NiSi.

**Figure 1**. (Color online) The calculated total and partial density of states in NiSi (black line) and NiGe (red line).

In order to explain the difference in the range of phonon spectrum in NiSi and NiGe, we have plotted the force constant matrix components (Figure 2) as a function of the interatomic distance. Typically, the force constant Φ(*k, k*′) between a pair of atoms (k, k′) is determined from the force on the atom k' due to a small displacement of the atom k. On the other hand, the self force constant Φ(*k, k*) occurs due the force on the atom k when the same atom is displaced. The values at the zero distance (Figure 2) refer to the self force constants. The self force constants for NiGe are about 70% of the values for NiSi. The first nearest neighbor distances occur in the range of 2.2–2.6 Å in NiSi and NiGe, which correspond to Ni-Si and Ni-Ge bond lengths. We can see that the magnitude of force constant matrix elements in NiSi is larger than that in NiGe, which reflects in the phonon spectrum of these compounds. The phonon spectrum of NiSi extends up to 50 meV, while in NiGe the range is only up to 34 meV.

**Figure 2**. (Color online) The calculated force constant matrix components as a function of the interatomic distance in NiSi(open circle) and NiGe (solid circle).

The measured neutron weighted phonon density of states using inelastic neutron scattering in NiSi has been shown in Figure 3. The experimental spectra compare very well with the calculations. The phonon spectra extend up to about 50 meV. The calculated neutron cross-section weighted partial phonon density of states of various atoms (Figure 3) clearly shows that the first peak in the phonon spectrum at 20 meV has mainly contribution from the Ni atom, while above 30 meV both the Ni and Si atoms contribute equally. The general characteristics of the experimental features are very well reproduced by the calculations. We may note that the peak heights in the calculated spectrum do not match very well with the measured inelastic spectrum. This may be due to the known inherent limitations both in experiment (e.g., incoherent approximation) and theory.

**Figure 3**. (Color online) Experimental and computed neutron weighted phonon density of states of NiSi. The calculated neutron cross-section weighted partial phonon density of states of various atoms is also shown. The calculated spectrum has been broadened by the energy resolution in the experiment.

Phonon dispersion relation in NiSi and NiGe has been calculated (Figure 4). In accordance to the calculated phonon density of states, the spectral range of phonons in NiSi extends up to 50 meV, while the highest mode in NiGe is at around 34 meV. Comparison of the dispersion curve shows that the acoustic modes in NiSi extend up to 20 meV while the energy range in NiGe is up to 12 meV. These will give rise to difference in the elastic properties which arises from the slope of acoustic phonon branches along various high symmetry directions in the Brillion zone of both the compounds. The optic modes in NiSi are well separated and give rise to a three peak structure in the phonon spectrum (Figure 1), while in NiGe the modes are not separated and they give rise to a nearly broad peak (Figure 1).

**Figure 4**. The calculated phonon dispersion curves in NiSi and NiGe along the high symmetry points in the Brillioun zone.

The group theoretical decomposition of phonons at the zone center is given as:

One each of the B _{1u}, B_{2u}, and B _{3u} modes are zero energy acoustic modes. The calculated Raman and infrared modes in NiSi and NiGe is given in Table 1. The comparison of zone center modes also shows that energies of various modes in NiSi extend up to about 50 meV while the maximum range in NiGe is up to 31 meV.

The calculated elastic constants for both the compounds are given in Table 2. It can be seen that all the elastic constants in NiSi have larger values in comparison to that in NiGe. This compares very well with the observed slopes of the dispersion curves in both the compounds. Similarly the bulk modulus value of NiSi (162.6 GPa) is found to be larger in comparison to that of NiGe (119.5 GPa). The bulk modulus of NiSi has been measured using synchrotron-based X-ray powder diffraction method (Lord et al., 2012) and found to be 165 ± 3 GPa, which is in very good agreement with our calculation.

### Thermal Expansion Behavior

Both NiSi and NiGe are reported to show negative thermal expansion along the **b**-direction, while expansion coefficient is positive along **a-** and **c**-axis. The overall volume thermal expansion coefficient is positive. The anisotropic thermal expansion behavior can be computed under the quasiharmonic approximation. Each phonon mode of energy *E*_{qj} (*j*^{th} phonon mode at point q in the Brillouin zone) contributes to the thermal expansion coefficient, which is given by following relations for an orthorhombic (Mittal et al., 2018) system:

Where V_{0} is the unit cell volume, and *s*_{ij} are elements of elastic compliances matrix *s* = *C*^{−1}. *C*_{v} *(q, j, T)* is the specific-heat contribution from the phonon of energy *E*_{q, j}. Γ_{a}, Γ_{b} and Γ_{c} are the anisotropic mode Grüneisen parameters of phonon of energy *E*_{q, j} as given by (Grüneisen and and Goens, 1924),

The volume thermal expansion coefficient for a orthorhombic system is given by:

The calculation of anisotropic thermal expansion requires calculation of anisotropic Grüneisen parameters. The calculations of Γ_{l}(*E*) are performed by applying an anisotropic stress by changing the lattice constant “a” and keeping the “b” and “c” parameters constant; and vice versa following eq.(6). Figure 5 gives the calculated variation of anisotropic Grüneisen parameters with energy averaged over all the phonon modes in the Brillouin zone in NiSi and NiGe along a, b, and c directions.

**Figure 5**. (Color online) Calculated variation of anisotropic Grüneisen parameters with energy averaged over all the phonon modes in the Brillouin zone in NiSi and NiGe along a, b, and c directions.

As shown in Figure 5, for NiSi, the acoustic modes up to about 2.5 meV have negative values of Γ_{l}(*E*) (*l* = *a, b.c*) of about −12, while in NiGe the Γ_{l}(*E*) have relatively smaller negative values of −6 along b- and c- axis and positive values of +8 along a- axis. Further for NiSi, the modes in the energy range of 12–16 meV show small negative values of −0.2 along b-axis. The Γ_{l}(*E*) (*l* = *a, b.c*) in the remaining part of the phonon spectrum have positive values. In case of NiGe, the Γ_{b} shows a dip in the energy range of 12–16 meV with positive values of +0.5. Above 2.5 meV, Γ_{l}(*E*) in NiGe along the three axes have positive values.

The calculated thermal expansion coefficients as obtained from equations (3)–(5) are shown in Figure 6. We find that for NiSi, the thermal expansion coefficients are calculated to be positive along a and c-axes, while it is negative along b-axis. Similarly for NiGe, the α_{b} has small negative values up to about 100 K and it is positive above 100 K. The α_{l} (*l* = a, c) along a and c-axis are found to be positive in the whole temperature range. The experimental thermal expansion coefficient α_{b} has large negative value of −16.1 × 10^{−6} K^{−1} for NiSi, while in NiGe α_{b} has small negative values of ~ −0.4 × 10^{−6} K^{−1}. Our calculations are not able to quantitatively reproduce the negative values along b-axis in both the compounds, while qualitatively they agree very well with the experimental observations of a large difference in α_{b} in both the compounds.

**Figure 6**. (Color online) Calculated anisotropic behavior of the linear thermal expansion coefficients along a, b, and c directions as a function of temperature in NiSi and NiGe.

The comparison between the experimental (Perrin et al., 2007) and calculated change in the lattice parameters with increasing temperature is given in Figure 7, along with change in unit cell volume of the two compounds. The change in a- and c- lattice parameters as computed in this work is in good agreement with reported experimental data (Perrin et al., 2007). However, along the b- axis the contraction with increasing temperature is found to be underestimated in our calculations than from the experiments (Perrin et al., 2007). The change in b- lattice parameter is found to be more abrupt in the experimental observations for NiSi than that in our calculated values. The calculations are found to underestimate the contraction along the b-axis in NiSi, while for NiGe the comparison between the experiments and calculations seems to show good agreement with the reported data. Further, we find that the change in the volume of NiSi (Figure 7) from our calculations is found to be underestimated in comparison to the experimental data. The calculated volume is found to be slightly overestimated in NiGe (Figure 7).

**Figure 7**. (Color online) Change in lattice parameters and volume of NiSi and NiGe with increasing temperature. Our computed values (solid lines) in comparison with experimental data (Perrin et al., 2007) (solid symbols).

As described above, the thermal expansion coefficient at a given temperature is an outcome of contribution of phonon modes of different energies in entire Brillouin zone. To get a better insight, we have calculated (Figure 8) the contribution of phonons of energy E, averaged over Brillouin zone, to the linear thermal expansion coefficient at 300 K. In case of NiSi, phonons from entire Brillouin zone contribute to positive expansion along both a- and c- axis. Interestingly, phonons below 30 meV contribute to negative expansion coefficient along b- axis. We find that phonon modes below 30 meV contribute in different ways along the a-, c-, and b- axis. Further, in case of NiGe, we find similar behavior of phonon contribution along a- and c- axis i.e., positive contribution to thermal expansion coefficient. However, along the b- axis, phonons up to 15 meV contribute to negative expansion coefficient. We find that the range of phonon energies contributing to NTE along b- axis is significantly different for NiSi and NiGe; and this leads to significant difference in the magnitude of negative thermal expansion along b- axis in these compounds. The important inference of the above analysis is that the modes below 15 meV strongly contribute to NTE along b- axis in both the compounds.

**Figure 8**. (Color online) The contribution to anisotropic linear and volume thermal expansion coefficients at 300 K, from phonon mode of energy E averaged over the Brillouin zone.

As mentioned in the computation details above, the thermal expansion coefficient along different axes is given by the sum of the products of anisotropic Γ and the elastic compliance components relevant to the given direction. These values have been calculated and given in Table 3. For NiSi and NiGe, thermal expansion α_{b}(T) is given by including the values of the individual compliance values from Table 3:

In case of NiSi, the values of Γ_{a} and Γ_{c} are comparable for both the compounds, and they cancel out: so only 0.129Γ_{b} term contributes. As Γ_{b} values are negative, the net value of α_{b} becomes negative. In case of NiGe, sum of the first and the third term is less than the second term, and since Γ_{b} values are positive, the net value of α_{b} is positive. The Γ_{b} values for different modes in the range of 12–18 meV (NiSi) are almost zero or negligibly small negative, while in case of NiGe they are small positive. As a result, the value of the thermal expansion coefficient in case of NiSi becomes more negative as compared to NiGe. Hence negative thermal expansion is more pronounced in NiSi.

In case of NiSi, the specific modes exhibiting negative values in α_{b} and the vibrational amplitude of the atoms in these low energy modes would throw light on the physical mechanism involved in NTE. This in turn would eventually lead us to better understanding of the observed negative thermal expansion along b-direction in NiSi. The substitution of Ge in place of Si, brings about changes in the lattice which cancels out the possibility of similar behavior in NiGe. We have calculated the mean square amplitudes of various atom in NiSi and NiGe along a-, b-, and c- axis as a function of phonon energy to understand the mechanism of NTE along b- axis in both the compound (Figure 9). In Figure 8, we have observed that the NTE along b- axis in NiSi (NiGe) is mainly contributed by phonon modes of energy 8–20 meV. However, along a- and c- axis all the phonon modes contribute to positive expansion behavior (Figure 8). The density of states around 8–20 meV energy is very large, hence these modes play the dominating role in thermal expansion behavior of NiSi(NiGe). In Figure 9, the mean square displacements of Ni and Si below 20 meV show similar amplitude along a- and c- axis. However, along b- axis the magnitude of Ni is about 4 times than that of Si. In NiGe compound, the mean square displacements of Ge along b-axis are large in comparison to that of Si in NiSi. The ratio of u^{2} values, along b-axis, for Ni to Ge is only about 1.5. These results of mean-square amplitudes clearly bring out the similarity between the dynamics along the **a**- and **c**- axes and a different behavior along the b-axis, particularly for the modes in the 8–20 meV range. This dynamical behavior is also reflected in the nature of the anisotropic thermal expansion behavior along the b-axis.

**Figure 9**. (Color online) The calculated anisotropic mean square displacement of various atom in NiSi and NiGe as a function of phonon energy E averaged over Brillouin zone at 300 K.

The NiSi structure consists of hexagonal sheets formed by Ni and Si atoms stacked along the **a**- axis. In Figure 10, we have given two of such representative modes at Y(0 1/2 0) point of 16.5 and 20.9 meV, which lead to negative expansion along **b**- axis. From Eq. (4), the 16.5/20.9 meV mode gives significant contraction along b- axis (α_{b} = −6.0 × 10^{−6} K^{−1}/−3.69 × 10^{−6} K^{−1}, Γ_{a} = 1.19/1.73, Γ_{b} = −1.31/0.00, Γ_{c} = 3.53/2.90) and expansion in the a–c plane. The same representative modes in NiGe appear at an energy of 11.3/13.8 meV and contribute to NTE along b- axis (α_{b} = −5.5/−3.02 × 10^{−6} K^{−1}, Γ_{a} = 1.68/1.03, Γ_{b} = −0.59/−0.24, Γ_{c} = 5.12/2.98). The displacement pattern of 16.5 meV mode shows that, two of the three neighboring Ni atoms in a given hexagon move anti parallel to each other along b- axis, while the third Ni is almost at rest; and similarly two of the Si atoms vibrate along the c-axis. In case of 20.9 meV mode, all the three Ni atoms and the Si atoms in the hexagon have significant displacements. Both of these modes propagate along the b-axis, and involve transverse Si vibrations largely along the c-axis. It may be noted that transverse vibrations may contribute to NTE since it causes a negative stress in the longitudinal direction.

**Figure 10**. (Color online) The displacement pattern of two representative modes at the Brillouin zone boundary (0 1/2 0) responsible for negative thermal expansion in NiSi and NiGe. The numbers below the figures give the phonon energy, α_{b} at 300 K, Γ_{a}, Γ_{b}, and Γ_{c}, respectively. Key: Ni–Red and Si/Ge–Blue.

The displacement patterns of these modes at 16.5/20.9 meV in NiSi are very similar to the modes at 11.3/13.8 meV for NiGe. However, in NiGe, other higher energy modes which give contribution to positive expansion are well below 30 meV and significantly contribute at ambient temperature, thus nullifying the NTE along b- axis from the above phonon mode, while in NiSi the high-energy modes are distributed up to 50 meV, hence positive contributions from these modes are less significant.

## Conclusions

Inelastic neutron scattering study of the phonon density of states in NiSi is reported. The measured data are in very good agreement with our calculations. We also report calculations on isostructural compound NiGe. Due to the heavier and bigger Ge atom, the phonon density of states of the two compounds differs in their total span. The contribution of Ni to the phonon density of states in both the compounds differs. NiSi shows large NTE along b-axis, while NiGe compound has very small NTE behavior. The difference in the negative thermal expansion along b-axis in both compounds is brought out very well in our calculations. The different thermal expansion along the b-axis for the two compounds could be understood in terms of the respective values of the Grüneisen parameters as well as the elastic compliances. We have identified representative phonon modes involving transverse vibration of Si/Ge atoms, which contribute significantly to the phenomenon of NTE along b-axis.

## Author Contributions

PG contributed in inelastic neutron scattering experiment, analysis and interpretation of experiment and ab-initio calculation, and writing of the manuscript. MG and BS contributed in inelastic neutron scattering experiment, ab-initio calculations, interpretation of experiment and ab-initio calculation, and writing of the manuscript. SM contributed in inelastic neutron scattering experiment, analysis and interpretation of experiment and ab-initio calculation and writing of the manuscript. RM formulated the problem, contributed in inelastic neutron scattering experiment, interpretation of experiment and ab-initio calculation, and writing of the manuscript. SC formulated the problem, contributed in interpretation of experiment and ab-initio calculation, and writing of the manuscript.

## Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The reviewer LW and handling Editor declared their shared affiliation.

## Acknowledgments

SC would like to thank the Department of Atomic Energy, India for the award of Raja Ramanna Fellowship. The use of ANUPAM super-computing facility at BARC is acknowledged.

## References

Acker, J., Van Den Berg, G., Bohmhammel, K., Kloc, C., and Van Miltenburg, J. (1999). Thermodynamic properties of the nickel silicide NiSi between 8 and 400 K. *Thermochim. Acta* 339, 29–33. doi: 10.1016/S0040-6031(99)00293-2

Acun, A. D., and Soyalp, F. (2012). Elastic and phonon properties of FeSi and CoSi in the B2 structure. *Philos. Mag.* 92, 635–646. doi: 10.1080/14786435.2011.617323

Caracas, R., and Wentzcovitch, R. (2004). Equation of state and elasticity of FeSi. *Geophys. Res. Lett.* 31:L20603. doi: 10.1029/2004GL020601

Carpenter, J. M., and Price, D. L. (1985). Correlated motions in glasses studied by coherent inelastic neutron scattering. *Phys. Rev. Lett.* 54, 441–443. doi: 10.1103/PhysRevLett.54.441

Connétable, D., and Thomas, O. (2009). First-principles study of the structural, electronic, vibrational, and elastic properties of orthorhombic NiSi. *Phys. Rev. B* 79:094101. doi: 10.1103/PhysRevB.79.094101

Dahal, A., Gunasekera, J., Harringer, L., Singh, D. K., and Singh, D. J. (2016). Metallic nickel silicides: experiments and theory for NiSi and first principles calculations for other phases. *J. Alloys Compd.* 672, 110–116. doi: 10.1016/j.jallcom.2016.02.133

Damascelli, A., Schulte, K., Van Der Marel, D., Fäth, M., and Menovsky, A. (1997a). Optical phonons in the reflectivity spectrum of FeSi. *Physica B: Condensed Matter* 230, 787–789.

Damascelli, A., Schulte, K., Van der Marel, D., and Menovsky, A. (1997b). Infrared spectroscopic study of phonons coupled to charge excitations in FeSi. *Phys. Rev. B* 55:R4863.

Delaire, O., Al-Qasir, I. I., Ma, J., dos Santos, A. M., Sales, B. C., Mauger, L., et al. (2013). Effects of temperature and pressure on phonons in FeSi${}_{1{-}x}$Al${}_{x}$. *Phys. Rev. B* 87:184304. doi: 10.1103/PhysRevB.87.184304

Delaire, O., Marty, K., Stone, M. B., Kent, P. R., Lucas, M. S., Abernathy, D. L., et al. (2011). Phonon softening and metallization of a narrow-gap semiconductor by thermal disorder. *Proc. Natl. Acad. Sci. U.S.A.* 108, 4725–4730. doi: 10.1073/pnas.1014869108

Detavernier, C., Lavoie, C., and d'Heurle, F. (2003). Thermal expansion of the isostructural PtSi and NiSi: negative expansion coefficient in NiSi and stress effects in thin films. *J. Appl. Phys.* 93, 2510–2515. doi: 10.1063/1.1545156

Dobson, D. P., Vocadlo, L., and Wood, I. G. (2002). A new high-pressure phase of FeSi. *Am. Mineral.* 87, 784–787. doi: 10.2138/am-2002-5-623

Donthu, S., Tripathy, S., Chi, D., and Chua, S. (2004). Raman scattering probe of anharmonic effects in NiSi. *J. Raman Spectr.* 35, 536–540. doi: 10.1002/jrs.1164

Franciosi, A., Weaver, J., and Schmidt, F. (1982). Electronic structure of nickel silicides Ni 2 Si, NiSi, and Ni Si 2. *Phys. Rev. B* 26:546. doi: 10.1103/PhysRevB.26.546

Geenen, F., Solano, E., Jordan-Sweet, J., Lavoie, C., Mocuta, C., and Detavernier, C. (2018). The influence of alloying on the phase formation sequence of ultra-thin nickel silicide films and on the inheritance of texture. *J. Appl. Phys.* 123:185302. doi: 10.1063/1.5022070

Jin, Y.-Y., Kuang, X.-Y., Wang, Z.-H., and Huang, X.-F. (2013). Effect of pressure on structural phase transition and elastic properties in NiSi. *High Press. Res.* 33, 15–26. doi: 10.1080/08957959.2012.761213

Kim, A., Choi, E., Kwon, S. H., Zheng, L., Kang, K., Nam, M., et al. (2018). Nickel silicide phase distribution of electroless deposited contacts on silicon P+ junction substrate using nano-indented atomic force microscopy. *Sci. Adv. Materials* 10, 542–546. doi: 10.1166/sam.2018.3059

Kim, H. (2015). First-principles calculations of the lattice instability and the symmetry-lowering modulation of PtSi. *J. Korean Phys. Soc.* 66, 612–616. doi: 10.3938/jkps.66.612

Kresse, G., and Furthmüller, J. (1996). Efficient iterative schemes for *{*ab initio} total-energy calculations using a plane-wave basis set. *Phys. Rev. B* 54, 11169–11186. doi: 10.1103/PhysRevB.54.11169

Kresse, G., and Joubert, D. (1999). From ultrasoft pseudopotentials to the projector augmented-wave method. *Phys. Rev. B* 59, 1758–1775. doi: 10.1103/PhysRevB.59.1758

Lee, P. S., Mangelinck, D., Pey, K. L., Shen, Z., Ding, J., Osipowicz, T., et al. (2000). Micro-Raman spectroscopy investigation of nickel silicides and nickel (platinum) silicides. *Electrochem. Solid-State Lett.* 3, 153–155. doi: 10.1149/1.1390986

Lee, S. K., Seok, K. H., Chae, H. J., Lee, Y. H., Han, J. S., Jo, H. A., et al. (2017). Bottom-gate poly-Si thin-film transistors by nickel silicide seed-induced lateral crystallization with self-aligned lightly doped layer. *Solid State Electron.* 129, 6–9. doi: 10.1016/j.sse.2016.11.019

Lin, J.-Y., Hsu, H.-M., and Lu, K.-C. (2015). Growth of single-crystalline nickel silicide nanowires with excellent physical properties. *CrystEngComm* 17, 1911–1916. doi: 10.1039/C4CE02513J

Lord, O. T., Vočadlo, L., Wood, I. G., Dobson, D. P., Clark, S. M., and Walter, M. J. (2012). High-pressure phase transitions and equations of state in NiSi. II. Experimental results. *J. Appl. Crystallogr.* 45, 726–737. doi: 10.1107/S0021889812016809

Mandrus, D., Sarrao, J., Migliori, A., Thompson, J., and Fisk, Z. (1995). Thermodynamics of FeSi. *Phys. Rev. B* 51, 4763–4767. doi: 10.1103/PhysRevB.51.4763

Mittal, R., Gupta, M. K., and Chaplot, S. L. (2018). Phonons and anomalous thermal expansion behaviour in crystalline solids. *Prog. Mater. Sci.* 92, 360–445. doi: 10.1016/j.pmatsci.2017.10.002

Monkhorst, H. J., and Pack, J. D. (1976). Special points for Brillouin-zone integrations. *Phys. Rev. B* 13, 5188–5192. doi: 10.1103/PhysRevB.13.5188

Parshin, P., Chumakov, A., Alekseev, P., Nemkovski, K., Perßon, J., Dubrovinsky, L., et al. (2016). Experimental observation of phonons as spectators in FeSi electronic gap formation. *Phys. Rev. B* 93:081102. doi: 10.1103/PhysRevB.93.081102

Pauling, L., and Soldate, A. (1948). The nature of the bonds in the iron silicide, FeSi, and related crystals. *Acta Crystallogr.* 1, 212–216. doi: 10.1107/S0365110X48000570

Perdew, J. P., Burke, K., and Ernzerhof, M. (1996). Generalized gradient approximation made simple. *Phys. Rev. Lett.* 77, 3865–3868. doi: 10.1103/PhysRevLett.77.3865

Perrin, C., Nemouchi, F., Clugnet, G., and Mangelinck, D. (2007). Anisotropy of the thermal expansion of the Ni(Si1–xGex) phases investigated by high-temperature x-ray diffraction. *J. Appl. Phys.* 101:073512. doi: 10.1063/1.2713992

Racu, A.-M., Menzel, D., Schoenes, J., Marutzky, M., Johnsen, S., and Iversen, B. B. (2008). Phonon properties in narrow gap FeSi and FeSb2 single crystals. *J. Appl. Phys.* 103:07C912. doi: 10.1063/1.2885067

Sarrao, J., Mandrus, D., Migliori, A., Fisk, Z., and Bucher, E. (1994). Elastic properties of FeSi. *Physica B: Condensed Matter* 199, 478–479. doi: 10.1016/0921-4526(94)91875-9

Sidorova, M. V., Kozorezov, A. G., Semenov, A. V., Korneeva, Y. P., Mikhailov, M. Y., Devizenko, A. Y., et al. (2018). Nonbolometric bottleneck in electron-phonon relaxation in ultrathin WSi films. *Phys. Rev. B* 97:184512. doi: 10.1103/PhysRevB.97.184512

Simon, A. H., Elizabeth, T. H. W., David, P. D., Lindunka, V., and Ian, G. W. (2017). The thermal expansion of (Fe 1– y Ni y)Si. *J. Phys. Condensed Matter* 29:335701. doi: 10.1088/1361-648X/aa797d

Stern, R., and Madsen, G. K. (2016). *Ab initio* investigation of the anomalous phonon softening in FeSi. *Phys. Rev. B* 94:144304. doi: 10.1103/PhysRevB.94.144304

Tang, P., Zhou, Q., and Zhang, S.-C. (2017). Multiple types of topological fermions in transition metal silicides. *Phys. Rev. Lett.* 119:206402. doi: 10.1103/PhysRevLett.119.206402

Vočadlo, L., Knight, K., Price, G., and Wood, I. (2002). Thermal expansion and crystal structure of FeSi between 4 and 1173 K determined by time-of-flight neutron powder diffraction. *Phys. Chem. Minerals* 29, 132–139. doi: 10.1007/s002690100202

Vočadlo, L., Wood, I. G., and Dobson, D. P. (2012). High-pressure phase transitions and equations of state in NiSi. I. *Ab initio* simulations*. J. Appl. Crystallogr.* 45, 186–196. doi: 10.1107/S0021889812000337

Von Känel, H., Mäder, K., Müller, E., Onda, N., and Sirringhaus, H. (1992). Structural and electronic properties of metastable epitaxial FeSi 1+ x films on Si (111). *Phys. Rev. B* 45:13807.

Weber, W. M., Geelhaar, L., Graham, A. P., Unger, E., Duesberg, G. S., Liebau, M., et al. (2006). Silicon-nanowire transistors with intruded nickel-silicide contacts. *Nano Lett.* 6, 2660–2666. doi: 10.1021/nl0613858

Wu, H., Kratzer, P., and Scheffler, M. (2005). First-principles study of thin magnetic transition-metal silicide films on Si(001). *Phys. Rev. B* 72:144425. doi: 10.1103/PhysRevB.72.144425

Zhang, T., Song, Z., Alexandradinata, A., Weng, H., Fang, C., Lu, L., et al. (2018a). Double-Weyl phonons in transition-metal monosilicides. *Phys. Rev. Lett.* 120:016401. doi: 10.1103/PhysRevLett.120.016401

Keywords: inelastic, neutron, phonons, ab-initio, negative thermal expansion, grunesien parameter

Citation: Goel P, Gupta MK, Mishra SK, Singh B, Mittal R and Chaplot SL (2018) Phonons and Thermal Expansion Behavior of NiSi and NiGe. *Front. Chem*. 6:331. doi: 10.3389/fchem.2018.00331

Received: 05 June 2018; Accepted: 16 July 2018;

Published: 14 August 2018.

Edited by:

Jun Chen, University of Science and Technology Beijing, ChinaReviewed by:

Qiang Sun, Zhengzhou University, ChinaLei Wang, University of Science and Technology Beijing, China

Copyright © 2018 Goel, Gupta, Mishra, Singh, Mittal and Chaplot. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Ranjan Mittal, rmittal@barc.gov.in