- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Journal of Materials
Volume 2014 (2014), Article ID 745736, 9 pages
Synthesis of Thin Films by Magnetron Sputtering
1Universidad Nacional Autónoma de México, Centro de Nanociencias y Nanotecnología,
Km 107 Carretera Tijuana-Ensenada, 22800 Ensenada, BCN, Mexico
2CNyN-UNAM, P.O. Box 439036, San Ysidro, CA 92143-9036, USA
Received 9 October 2013; Revised 26 March 2014; Accepted 3 April 2014; Published 18 June 2014
Academic Editor: Roger Narayan
Copyright © 2014 G. Soto et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
In this work films were prepared by reactive magnetron sputtering at room temperature and deposited on a silicon wafer. It was found that the diffractograms of the nitrogen-rich rhenium film are consistent with those produced by high-pressure high-temperature methods, under the assumption that the film is oriented on the substrate. Using density functional calculations it was found that the composition of this compound could be ReN3, instead of ReN2, as stated on previous works. The ReN3 compound fits in the Ama2 (40) orthorhombic space group, and due to the existence of N3 anions between Re layers it should be categorized as an azide. The material is exceptionally brittle and inherently unstable under indentation testing.
In recent years, reports on the synthesis of novel materials by high-pressure and high temperature (HPHT) methods have become common. The use of extreme pressures and temperatures, combined with compositional variables, provides an opportunity to synthesize uncommon materials and/or to tune the physical properties of materials for a wide range of applications . HPHT methods were pioneered by researchers looking to synthesize diamonds, but now are often applied to discover materials that might have mechanical, chemical or thermal resistance comparable to those of diamond. Many transition metal nitrides are poised as promising hard materials principally for machining of ferrous alloys. Using HPHT methods, novel nitrides of heavy metals such as platinum, osmium, and iridium [2, 3] have been synthesized. Due to the observed high bulk modulus of these materials, some HPHT studies have sought to discover novel nitrides of heavy metals. An example of this effort is the work of Kawamura et al., who recently published an article describing the HPHT synthesis of ReN2 .
In direct competition with HPHT methods for the synthesis of artificial diamond are methods using plasma . These have been successful in terms of the quality and purity of the obtained diamond, with the added advantage that these materials can be applied as coatings directly from the synthesis. Similarly, the plasma methods are used to discover alternative materials that can replace diamond. The success of plasmas for producing new materials originates from the high reactivity obtained by breaking molecular bonds after multiple ionization events. Using laser ablation experiments—a very energetic plasma method—it was possible to grow films composed of , , and . The main limitations of plasma synthesis are the low growth rates and the textured crystalline orientations, which make the obtained materials not easy to characterize.
Recently, we published an article describing the synthesis of rhenium nitride using a plasma method . We were successful in describing the experimental conditions and the results when nitrogen is incorporated in low concentrations. Basically, the () consists of a solid solution of nitrogen incorporated into the octahedral interstices of a rhenium lattice. However in this work we failed in interpreting the data when nitrogen is present in high concentrations (). Recently, Kawamura et al. reported the HPHT formation of a rhenium-nitrogen compound, and they conclude to have synthesized ReN2, with structure similar to MoS2 . The report of Kawamura et al. was very useful in resolving the structure of our films in the nitrogen-rich condition. As a result, the objective of this report is to describe the type of compound that is produced when a rhenium target is deposited by DC sputtering in an excess of nitrogen, that is, in target-poisoned state. As a clarification, the target poisoning means that the compound is formed not only on the substrate, as desired, but also on the sputter target, which results in a significantly reduced sputter yield and, therefore, a reduced deposition rate . This regime of deposition is seldom used for technological applications, but it can be used for the preparation of novel compounds.
2.1. Thin Film Preparation and Characterizations
The rhenium nitride films were grown by reactive DC sputtering. The experiment was carried out in a custom-made ultrahigh vacuum (UHV) deposition system. It consists of three stainless steel chambers: sample loading, deposition, and analysis. The analysis chamber is a PHI-548 system equipped with X-ray photoelectron (XPS) and Auger electron (AES) spectroscopies, and it operates at 10−10 Torr, achieved by means of an ion pump. The sample loading and deposition chambers are independently evacuated to pressures 10−8 Torr or lower by using a 130 L s−1 turbomolecular pump. A 2.5 cm diameter Re target (99.9%) is positioned in a magnetron gun, 10 cm above a (111)-silicon substrate. A mechanized rod is used to transfer the sample between chambers, which are isolated by gate valves. With this setup the deposited films can be analyzed without exposure to the outside atmosphere. The rhenium nitride films were deposited at room temperature at different N2 pressures in the following manner. In the deposition chamber, a mass-flow controller fed N2 (ultrahigh purity grade) at increasing flow rates to set pressures ranging from 3 mTorr (1 mL min−1) to 180 mTorr (35 mL min−1). After setting the N2 pressure, an additional 5 mL min−1 flow of Ar (ultrahigh purity grade) was allowed into the chamber. The samples were grown at 120 Watts for 50 minutes, after that they were immediately transferred to the analysis chamber for XPS and AES measurements. A reference sample, metallic Re, was deposited in the absence of N2. The assessment of nitride formation was performed by comparing core-level energy shifts in the XPS lines to the metallic Re sample. The XPS measurements used an Al ( eV) X-ray nonmonochromatic source. The X-ray diffraction (XRD) patterns were collected ex situ in a Philips X’pert diffractometer using the Cu line with a wavelength of nm. The film morphology and roughness were determined by atomic force microscopy using a NanoScope III, Veeco Instruments, in contact mode. The mechanical properties of the samples were measured using a TriboScope nanoindenter, Hysitron.
2.2. First Principle Calculations
Calculations were realized within the framework of density functional theory. We employed the full-potential linearized augmented plane wave method, as implemented in the Wien2k code , using the generalized gradient approximation . We used muffin-tin radii () of 2.0 au for Re atoms and of 1.0 au for N and angular momenta up to k-points were used in the irreducible part of the Brillouin zone. Optimization of volume, (), (), and internal atomic positions were performed recursively. Having arrived at the minimum energy structures, their volumes were varied around the equilibrium value, and their bulk moduli were obtained after fitting the Murnaghan equation of state to the calculated volume-energy data. The compound formation energy, , from molecular nitrogen is computed as the difference between the ab initio total energies of and those of Re metal and molecular N2 . Consider Here and are the number of metal and nitrogen atoms, respectively. From this definition, less positive values mean higher stability of the compound, relative to its spontaneous dissociation to Re and N2 constituents. Similarly, the formation energy from free atomic nitrogen, , is computed as follows: Because the molecular and atomic species cannot be treated directly within Wien2k, the N2 dimer was simulated by constructing a sufficiently large tetragonal unit cell (space group P4/mmm, , and ), locating N2 in the vertices along the z-axis (N-positions (0, 0, ±0.0981)) and performing the calculation for a single k-point at the origin of the first Brillouin zone. The nitrogen atomic species was carried out with analogous calculations for a nitrogen atom placed into the fcc unit cell (Fm3m, ). To confirm that our calculations can be used for nitrogen-rich compounds, we performed equivalent calculations for AgN3. We have minimized the forces acting on all atoms. Our calculations agree with the experimental parameters for AgN3 and with previous calculations for that compound . The results from these calculations are summarized in Table 1.
3.1. Thin Film Deposition
The results for the deposition of films () were discussed in our previous work . In this new work the same experimental setup and method are used for film preparations; however we will focus on the nitrogen-rich conditions, that is, target in poisoned state. The nitrogen-rich rhenium film was prepared at room temperature and a pressure of 180 mTorr N2 ( mL min−1; is the nitrogen flow), with an additional flow of Ar, mL min−1 ( mTorr). Figure 1 presents several XRD patterns of rhenium nitrides grown at different nitrogen flow rates. This figure is important to establish a comparison between compounds with low and high nitrogen contents. The interpretation of these diffractograms follows.
The peaks indicated with an asterisk () represent background from the silicon wafer substrate, which exhibits peaks at 28.40° (111) and 58.85° (222). The sample deposited with and presented in Figure 1(a) shows three main diffraction peaks at 37.51°, 40.4°, and 42.8°. These correspond to the diffraction planes (100), (002), and (101), respectively, of the hexagonal close-packed structure reported for pure Re in the powder diffraction file (PDF) number 89-2935. The (002) peak corresponds to a Re-Re interplanar spacing of 2.23 Å in the -direction. This pure-Re film is to some extent preferentially oriented, with (002)-Re parallel to (111)-Si. For the sample deposited with mL min−1, presented in Figure 1(b), the peaks corresponding to hexagonal Re were not present. Instead, wider peaks at approximately 38.45° (A), 44.69° (B), 65.20° (C), and 78.1° (D) appear in the pattern. This pattern agrees with the face-centered cubic phase reported by Fuchigami et al. , with = . This pattern also agrees with a DFT computer-generated diffractogram (not shown here) for cubic Re4N, with N present in octahedral interstices. For higher values of , the positions of peaks A through D shift to lower angles as the N/Re ratio within the films increases. The shift is more noticeable for the A-peak, which corresponds to the (111)-plane when referenced to the cubic phase. In Figure 1(c), this peak shifts to 38.16° (corresponding to a Re-Re interplanar spacing of 2.35 Å) and then to 36.46° in Figure 1(d) (corresponding a Re-Re interplanar spacing of 2.46 Å). These shifts indicate that the crystalline structure expands with increasing nitrogen content. Taking the most important diffraction peaks from Figure 1, in Figure 2 we show the dependence of interplanar distance on nitrogen flow rate during deposition. Additionally, in Figure 1(d), for the sample prepared with mL min−1, two new broad peaks appear at approximately 51.02° (E) and 65.44° (F). These peaks cannot be indexed to a cubic phase and hereafter are related to a different arrangement of rhenium nitride, given their shape and position. In the same sample, some peaks were present at 16.45° (110) and 25.47° (210); these might be indexed to a hexagonal crystalline phase of ReO3, according to PDF number 40-1155. Keep in mind that the XRD measurements were taken ex situ, and then some amount of oxidation of samples is expected.
At the time that this film was synthesized, we conducted an exhaustive search trying to index this pattern to a known crystalline phase of rhenium compounds. We found that no phase described in the PDF database or in the reviewed literature for rhenium in the metallic, nitride, or oxide forms is in agreement with the diffraction pattern presented in Figure 1(d), except for the peaks at 16.45 and 25.7 related to the formation of ReO3. However, it is clear from Figure 2 that there is a noticeable expansion of the crystal lattice as a function of increasing nitrogen incorporation. The expansion occurs preferentially on the close-packed direction of the rhenium atoms, that is, along the (002)-direction when the hexagonal structure of metallic Re is taken as reference and along the (111)-direction when the cubic structure of is taken as reference.
3.2. Nanoindentation Testing of Nitrogen-Rich Films
Figure 3 shows several load-unload displacement cycles for the sample growth in nitrogen-rich conditions with applied load limits of 50, 500, and 1000 μN. The cycles consist of 5 s of an increasing ramp-up to the maximum load, 5 s of sustaining the maximum load, 5 s of discharge, and 3 s of added acquisition at the contact force. It can be seen that this film shows a very peculiar behavior. For a typical sample during the ramp-up phase, when the indenter contacts the material, it is observed that the sample initially undergoes elastic deformation until the applied load reaches a critical value. When this critical value is exceeded, plastic deformation takes place, which is characterized by a sudden depth incursion or a “pop-in” . The peculiar aspect of this material is the number of pop-ins observed in the load-displacement curve. As we can see from the inset of Figure 3, the elastic deformation is very small before the material experiences plastic deformations. These sudden insertions appear in the plot as vibrations of the Berkovich tip. The deformation is inherently unstable, occurring in bursts of highly localized strains . There are also larger pop-in events that can be observed in the load-displacement curves, . Taken as a whole, the deformation is the sum of the small and large pop-in events that begin at infinitesimal loads. When a material has the tendency to fracture with very little or no detectable deformation beforehand, it is considered brittle. We can consider the current material to have perfectly brittle mechanical behavior because there is no elastic deformation. When the load is held constant, the plastic flow is evident because the penetration depth increases noticeably. Given that the film has a thickness of one micron (estimated by means of X-ray measurements) a creep of 100 nm is significant. Analyzing the initial portion of the unloading response gives an estimate of the elastic modulus of the indented material. As we can see, the initial slope of the unloading portion approaches zero, , which means that this material is very soft. A detailed analysis of the raw data (difficult to detect at the scale pictured) shows extremely small “pop-out” events in the unloading portion of the curve. These “pop-outs” are attributed to the release of internal stresses generated during loading . Some of the deformation is recovered by the pop-outs. The strangest thing is that this material fully recovers from its deformation after removal of the load. To be exact, the hysteresis curve returns to the initial point upon removal of the load, as seen in Figure 3. In experiments using the Berkovich tip as an AFM probe, it can be seen that no vestige of the indentations remains, as shown in Figure 4, which illustrates the film before and after indentation. We have no explanation for this behavior because, as previously stated, this is a brittle material. However, we speculate that this results from chemical changes (oxidation) in the sample rather than a recovery from mechanical deformations of ReN3.
To discuss our results we want to use DFT calculations and the experimental result from the work of Kawamura et al. for a material that was assigned a ReN2 composition . The interpretation of Kawamura was used as the starting point for our DFT calculations. In his article Kawamura explains a metathesis reaction between lithium nitride and rhenium pentachloride in a HPHT environment, which produces materials shaped like platelets. Kawamura provides powder XRD patterns according to Bragg-Brentano and Gandolfi geometries. All diffraction peaks were indexed by Rietveld analysis to a hexagonal P63/mmc (194) space group, which is isostructural with MoS2 and similar to ReB2. The Rietveld analysis gave a lattice constant of , , and Å3 as well as an atomic coordinate of . XPS analysis was used to confirm the Re-N bonding. The authors conclude that their synthesis produces ReN2. However the structural parameters differ from some theoretical calculations for rhenium-nitrogen compounds in 1 : 2 stoichiometric ratios [20–22]. While the reasons for the discrepancy were not clarified, it was assumed that this type of synthesis often results in the development of unusual metastable phases. The unit cell and structural parameters proposed by Kawamura are given in Table 1.
The X-ray pattern for the films growth by sputtering at pressure of 180 mTorr of N2 can be compared with the XRD pattern generated for structure proposed by Kawamura et al. . That is shown in Figure 5, where the experimental diffraction is put side by side with computer-generated patterns from the structure described by Kawamura under particular conditions. Figure 5(a) is the experimental thin film pattern using the Bragg-Brentano geometry for the Cu-Kα radiation. Figure 5(b) presents the generated powder pattern for the Kawamura’s structure using the Cu-Kα wavelength. Figure 5(c) presents only the peaks that can be associated with the experimental pattern. Figure 5(d) presents the diffractogram for the structure proposed by Kawamura with the assumption that this material is highly oriented in the c-direction using the Bragg-Brentano geometry; this pattern agrees with the data reported by Kawamura et al. in the inset of Figure 3 in . It can be observed that the XRD pattern of the sputtered film can be indexed to the material synthesized by Kawamura et al. . However, the film has a texture that is different from that observed by Kawamura. The texture is defined as the distribution of crystallographic orientations of a polycrystalline sample. In the film, the crystallographic orientations are not random but rather assume a preferred orientation of (010)-planes parallel to (111)-planes of the Si substrate. There is some contribution from (014)-planes, indicative of some tilt in the c-direction, and some contribution from (120)- and (124)-planes. Figure 6 presents a sketch illustrating these possibilities. In conclusion, the (010)-planes of the film are parallel to the (111)-Si orientation, with slightly misoriented grains around this preferential orientation. As a consequence of this texture, there are missing reflections from the sample in the Bragg-Brentano geometry and to some extent weak counting rates. In contrast, the material synthesized by Kawamura is highly textured in the c-direction. We conclude that it is very likely that the material that was synthesized by the HPHT method has the same crystalline system as the material that was synthesized by the sputtering method.
To understand more deeply the nature of this compound, we performed DFT calculations using the structure proposed by Kawamura as starting point. The plot of total energy against volume is shown in Figure 7. The derivative at the volume given by the structural parameters reported by Kawamura is not zero. implies that the system has not reached equilibrium. For any metastable point the cancelation of forces must be complete, being that an indispensable condition of metastability. Surprisingly, the structure proposed by Kawamura cannot be a metastable point. This implies that the structure proposed by Kawamura is not fully relaxed. This becomes evident since the energy of ReN2 at the points defined by Kawamura is high in comparison with the energy of the relaxed structure calculated by DFT; see Table 1. In addition, the volume at the points defined by Kawamura is larger than the volume of the relaxed structure. It is unlikely that a structure that was synthesized under pressure would develop a metastable phase that is a “tensioned variation” of the relaxed structure. Very recently it was proposed, using DFT structural optimizations, that the MoS2-like structure of ReN2 will reform due to an rearrangement of the N2 dimers, in such a way that the hexagonal symmetry is lost . Evidently, the diffraction pattern of the DFT-improved structure will not match the X-ray diffraction reported by Kawamura and, therefore, will not match either of the diffractograms obtained in this work. Additionally, in , the MoS2-like structure of ReN2 was found to be mechanically unstable according to the Born-Huang criterion for a hexagonal structure because a negative value for the C44 elastic constant was found. A negative value for the C44 constant implies that the structure is tensioned and that the related strain tensor needs to be relaxed in the opposite direction. This is identical to our findings and is illustrated in Figure 7, where the relaxation of to lower ratios is a crucial step toward finding the lowest-energy point. In brief, the structure as it was proposed by Kawamura lacks certain features in the -direction. The DFT calculations reveal that the empty space in -direction cannot be justified by chemical bonds or by Van der Waals forces. With the goal of elucidating this irregularity, it was conjectured that there was sufficient space between the Re layers to insert (N3)-azide anions instead of (N2)2-diazenide anions. The azide functional group is larger than the diazenide group; the larger anion could therefore match the interplanar spacing between the Re layers. Table 1 lists the structural parameters resulting from the DFT optimization of ReN3 constrained within the P63/mmc (194) space group. It can be seen that the relaxed structure calculated by DFT for the ReN3 composition is a better match for the cell parameters given by Kawamura. This suggests that the compound of interest might be ReN3 rather than ReN2. However, the atomic positions of nitrogen and rhenium are highly constrained when the P63/mmc (194) space group is used. To avoid this, the structural degrees of freedom were increased by reducing the symmetry of the space group. Instead of the initial choice of P63/mmc (194), the Am2a (40) space group was used. Even using the lower-symmetry Am2a space group, the diffraction pattern generated is close to the pattern observed by Kawamura et al., as shown in Figure 8. More DFT refinements are recommended to improve the quality of the match.
From the experimental conditions, where a rhenium target was sputtered in nitrogen environments at different pressures, it can be concluded that incorporation of nitrogen causes important expansion of the rhenium lattice, mainly in c-direction. The material synthesized by magnetron sputtering is similar to those prepared by the HPHT method, but with a different crystallographic texture. A diazenide group, as it was proposed by Kawamura, inserted within rhenium interstices could not explain the diffraction patterns of the recently prepared rhenium-nitrogen compounds. We propose that the azide group between the rhenium layers provides a better match for these diffraction patterns. The diffractograms given by DFT with ReN3 agree with the experimental evidence. Given that, the material should be ReN3 instead of ReN2. The ReN3 is a brittle material and also it is mechanically unstable under indentation. The possibility of forming of metal azides by sputtering deposition is a new field that must be further studied.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
The authors would like to give thanks to D. Dominguez for technical assistance. Financial support from DGAPA-UNAM Grants IN114209, IT00314, and CONACYT 83275 is acknowledged. A. Reyes appreciates the sabbatical scholarship from PASPA-DGAPA-UNAM at Department of Physics, The Pennsylvania State University.
- E. Horvath-Bordon, R. Riedel, A. Zerr et al., “High-pressure chemistry of nitride-based materials,” Chemical Society Reviews, vol. 35, no. 10, pp. 987–1014, 2006.
- J. C. Crowhurst, A. F. Goncharov, B. Sadigh et al., “Synthesis and characterization of the nitrides of platinum and iridium,” Science, vol. 311, no. 5765, pp. 1275–1278, 2006.
- A. F. Young, C. Sanloup, E. Gregoryanz, S. Scandolo, R. J. Hemley, and H.-K. Mao, “Synthesis of novel transition metal nitrides IrN2 and OsN2,” Physical Review Letters, vol. 96, no. 15, Article ID 155501, 2006.
- F. Kawamura, H. Yusa, and T. Taniguchi, “Synthesis of rhenium nitride crystals with MoS2 structure,” Applied Physics Letters, vol. 100, no. 25, Article ID 251910, 2012.
- S. Vepřek, “The search for novel, superhard materials,” Journal of Vacuum Science and Technology A, vol. 17, no. 5, pp. 2401–2420, 1999.
- A. P. Caricato, M. Fernàndez, G. Leggieri et al., “Reactive pulsed laser deposition of gold nitride thin films,” Applied Surface Science, vol. 253, no. 19, pp. 8037–8040, 2007.
- G. Soto, “Synthesis of PtNx films by reactive laser ablation,” Materials Letters, vol. 58, no. 16, pp. 2178–2180, 2004.
- M. G. Moreno-Armenta, J. Diaz, A. Martinez-Ruiz, and G. Soto, “Synthesis of cubic ruthenium nitride by reactive pulsed laser ablation,” Journal of Physics and Chemistry of Solids, vol. 68, no. 10, pp. 1989–1994, 2007.
- G. Soto, A. Rosas, M. H. Farias, W. De la Cruz, and J. A. Diaz, “Characterization of rhenium nitride films produced by reactive pulsed laser deposition,” Materials Characterization, vol. 58, no. 6, pp. 519–526, 2007.
- D. Güttler, B. Abendroth, R. Grötzschel, W. Möller, and D. Depla, “Mechanisms of target poisoning during magnetron sputtering as investigated by real-time in situ analysis and collisional computer simulation,” Applied Physics Letters, vol. 85, no. 25, pp. 6134–6136, 2004.
- J. L. P. Blaha, K. Schwarz, G. Madsen, and D. Kvasnicka, Wien2k, an Augmented Plane Wave + Local Orbitals Program for Calculating Crystal Properties, Technische Universität Wien, Wien, Austria, 2001.
- J. P. Perdew, K. Burke, and M. Ernzerhof, “Generalized gradient approximation made simple,” Physical Review Letters, vol. 77, no. 18, pp. 3865–3868, 1996.
- C. Stampfl, W. Mannstadt, R. Asahi, and A. J. Freeman, “Electronic structure and physical properties of early transition metal mononitrides: density-functional theory LDA, GGA, and screened-exchange LDA FLAPW calculations,” Physical Review B, vol. 63, no. 15, Article ID 155106, 11 pages, 2001.
- A. B. Gordienko, Y. N. Zhuravlev, and A. S. Poplavnoi, “Energy band structure of silver azide (AgN3),” Soviet Physics Journal, vol. 35, no. 2, pp. 130–132, 1992.
- G. Soto, H. Tiznado, A. Reyes, and W. De La Cruz, “First principles calculations of interstitial and lamellar rhenium nitrides,” Journal of Alloys and Compounds, vol. 514, pp. 127–134, 2012.
- M. Fuchigami, K. Inumaru, and S. Yamanaka, “Interstitial binary nitride ReNx phases prepared by pulsed laser deposition: structure and superconductivity dependence on nitrogen stoichiometry,” Journal of Alloys and Compounds, vol. 486, no. 1-2, pp. 621–627, 2009.
- A. C. Fischer-Cripps, Nanoindentation, Springer, New York, NY, USA, 2002.
- B. Moser, J. Kuebler, H. Meinhard, W. Muster, and J. Michler, “Observation of instabilities during plastic deformation by in-situ SEM indentation experiments,” Advanced Engineering Materials, vol. 7, no. 5, pp. 388–392, 2005.
- C. A. Schuh, “Nanoindentation studies of materials,” Materials Today, vol. 9, no. 5, pp. 32–40, 2006.
- G. Soto, “Computational study of Hf, Ta, W, Re, Ir, Os and Pt pernitrides,” Computational Materials Science, vol. 61, pp. 1–5, 2012.
- X. P. Du, Y. X. Wang, and V. C. Lo, “Investigation of tetragonal ReN2 and WN2 with high shear moduli from first-principles calculations,” Physics Letters A, vol. 374, no. 25, pp. 2569–2574, 2010.
- E. Zhao and Z. Wu, “Structural, electronic and mechanical properties of ReN2 from first principles,” Computational Materials Science, vol. 44, no. 2, pp. 531–535, 2008.
- Y. Wang, T. Yao, J.-L. Yao, J. Zhang, and H. Gou, “Does the real ReN2 have the MoS2 structure?” Physical Chemistry Chemical Physics, vol. 15, no. 1, pp. 183–187, 2013.