Abstract

Inelastic neutron scattering experiments were performed on polycrystalline samples of NdFeAsO0.85F0.15 over a wide temperature range (3 K–250 K). Based on the analysis of the experimental data, a Nd3+ CF energy level scheme is proposed to give a consistent explanation about the observed CF transitions. The observed extra ground-state CF transitions could not be simply explained by the transitions between five Kramers doublets split from the Nd3+   ground state in the point symmetry. A reliable explanation would be a superposition of crystal fields due to different local symmetries around the Nd3+ ions induced by the fluorine doping.

1. Introduction

For materials with rare-earth (R) elements, the interactions between f-electrons and the crystal field potentials produced by local lattice environment play important roles in determining the electronic ground state of R ions, as well as in understanding the magnetic and thermodynamic properties of these compounds. Inelastic neutron scattering has been widely accepted as the most reliable technique in determining the crystal field potentials in rare-earth based materials. Particularly, for superconducting materials, in either high- cuprates or other R based heavy fermion superconductors, important research progress has been made in investigating the neutron spectroscopy of the CF excitations [1, 2].

The discovery of Fe-based superconductor has caused a tremendous upsurge of scientific interest since 2008 [3]. Shortly after the report of 26 K superconductivity in , through the substitution of La by other R elements with f-electrons (such as Ce, Pr, Nd, and Sm), superconducting critical temperature above 50 K could be observed in [46]. So far, these so-called FeAs-1111 compounds mark the highest in all Fe-based superconductors, which also provided a unique opportunity in studying the correlation between CF excitations and physical properties in these exotic superconducting families. Up to now, there have been several publications about crystal field studies in FeAs-1111 compounds through inelastic neutron scattering technique. For example, in , the Ce CF levels were determined and the intrinsic linewidth of 18.7 meV was found to have a clear anomaly at [7]; inelastic neutron scattering measurements of CF transitions in (RE = Pr, Nd) provided evidences that there were two distinct charge environments around the RE ions arising from a random distribution of fluorine ions [8]. The CF level structures for both paramagnetic phase and antiferromagnetic phase of NdFeAsO were constructed by Xiao et al. [9]. However, for , there are still insufficient neutron scattering data on the CF transitions and the CF energy level scheme of Nd ion is still not clear.

In this paper, we report inelastic neutron scattering studies on superconducting samples. Based on the analysis of inelastic neutron scattering spectra at different temperatures, we propose a crystal field energy level scheme for Nd ion, which could explain the CF transitions observed in the experimental data. The additional CF transitions beyond the requirements of point symmetry indicate the superposition of crystal field and inhomogeneous ligand environments around the rare-earth ions in .

2. Experimental Results and Discussions

About 20 grams of polycrystalline samples were synthesized using the solid state reaction method. Resistivity measurements show that the onset superconducting transition temperature of our samples is about 47 K while the state of zero resistivity is reached at about 43 K, as shown in the inset of Figure 1(a). Inelastic neutron scattering measurements on were carried out on thermal triple-axis spectrometer TAIPAN at the Australian research reactor OPAL [10]. The final neutron energy was set to be 14.7 meV, with pyrolytic graphite (PG) as the monochromator and analyzer. The collimation was 48′-40′-40′-120′.

Figure 1 shows the raw data of the inelastic neutron spectra for at different temperatures. From 3 K to 250 K, constant- scans were performed every 50 Kelvin. The momentum transfer was fixed to be 3.1 Å in the whole measurements. First of all, the strong peaks at about 45 meV and 40 meV (marked with red stars) could be easily distinguished as the spurious peaks, which originated from the elastic neutron scattering. (The best energy windows for PG filter are 14.7 meV and 13.5 meV. In our experimental setup, was fixed to be 14.7 meV; when  meV and 54 meV, the energy transfers of neutron (14.7 meV and 13.5 meV) were corresponding to the strong elastic scattering signals. That is the origin of the spurious peaks at the energy transfers equaling 44.1 meV and 40.5 meV.) At below 30 meV plenty of CF excitation peaks could be clearly observed in the spectra at 3 K (marked with green arrows). With increasing temperature, the intensities of these peaks gradually decreased, which indicates that the corresponding CF excitations are ground-state excitations. In high temperature regions, some new CF peaks gradually emerge at around 11 meV and their absolute intensities got higher with increasing temperature. So these peaks should correspond to the excitations between excited CF states.

In order to fit and carefully analyze the CF peaks, after subtracting the spurious peaks, the next thing needed is to clarify the phonon contributions from the inelastic neutron spectrum. From both experimental data and theoretical calculations, it is known that the phonon spectra are almost the same for both parent compound and F-doped superconductors in FeAs-1111 system [1113]. So with the generalized phonon density of states (GDOS) of NdFeAsO in reference paper [9], after transformation with Bose factor and the neutron cross section equation [14, 15]we can finally calculate the phonon background for each temperature in our inelastic neutron spectra. The results calculated for  K are shown in Figure 1(f) (red hollow scatters).

After subtracting phonons, we noticed that the spectrum at 3 K almost overlaps with the spectrum at high temperatures in energy transfer above 36 meV except for notable additional intensity contributions between 40 meV and 48 meV. These additional neutron intensities gradually disappeared with increasing temperatures, which indicates the existence of ground-state CF excitations in this energy range. So it is reasonable to choose the spectrum above 36 meV at high temperature as instrument background because these intensities could not come from CF excitations and most likely result from the neutron contributions. After subtracting both the phonon and instrument background, the spectrum obtained at 3 K could be well fitted by seven CF excitation peaks with the method of Gaussian fitting, as shown in Figure 2(a) (labelled as A–G, resp.). The similar crystal field excitations at energy transfer below 30 meV were also observed by Goremychkin et al. [8]. Because our data included higher energy transfer region, we could observe the CF excitations at 41.5 meV and 44.5 meV. From both Figures 1 and 2(b), it can be seen that the intensities of all these seven CF peaks decrease with increasing temperatures, which is clearly the characteristic of ground-state excitations.

Therefore based on our inelastic neutron data, we proposed a CF energy level scheme for , which is illustrated in Figure 3. The crystal field energy levels of NdFeAsO were obtained from another publication [9]. As shown in Figure 3, under crystal field with the point symmetry, the ground state splits into five doublets in NdFeAsO at temperatures above the antiferromagnetic and structural phase transition [9]. For , the ground-state CF transitions marked as “A” to “G” in Figure 3 were observed in the inelastic neutron spectra as shown in both Figures 1 and 2(a). From the comparison between the 150 K and 3 K data in Figure 2(b), we could also notice that, in high temperature neutron spectra, there were some new CF peaks emerging at around 11 meV and 26 meV. These new CF peaks were corresponding to the excitations between excited states, as illustrated in Figure 3 (i.e., transitions “H” and “I”). The consistence between the observed CF transitions and the CF energy level scheme proved the validity of our CF model proposed in this paper.

It is clear that the CF level scheme of is significantly different from that of tetragonal NdFeAsO. In order to understand this difference, first of all we need to check whether there were CF states splitting caused by magnetic phase transitions. In the inset of Figure 1(a), the typical temperature dependence of resistances was measured on our samples. There were no traces of phase transitions because there was no anomaly in the normal state resistance which should be the typical features in the samples with structural or magnetic phase transitions. The resistance measurements were also performed in different batches of our samples which show similar features. Also based on our previous neutron measurements on samples with similar doping [16], we could safely rule out the possibility of structural or magnetic phases transitions induced additional CF states.

All the CF peaks should correspond to certain transitions of Nd f-electrons between different CF states. In the weak crystal field approximation, which applies well to the R elements, the multiplets of different ’s (total angular moment) are well separated in energy compared to the intramultiplet splittings; therefore, the interactions between multiplets can be ignored. In such a case, the crystal field states can be considered based solely on the ground-state splitting [17].

According to Hund’s rules, the ground state of ion, which contains three f-electrons, is . The splitting of CF levels is directly related to the local lattice symmetry. samples have a tetragonal crystal structure with space group. The Nd atoms are located at the crystallographic site therefore having point symmetry. This gives the crystal field Hamiltonian as follows:where are the CF parameters and are Steven’s operator equivalents. This CF symmetry splits the ground-state multiplet of the ion into five Kramers doublets.

However, suppose the CF peaks at  K in our samples all came from the transitions between the five doublets as in the high temperature tetragonal phase of NdFeAsO [9]; then the number of ground-state excitation peaks should be at most four. According to our data, the observed CF transitions for outnumbered the transitions allowed for with point symmetry. Therefore the Hamiltonian in (2) is insufficient to describe the CF states in F-doped samples.

A reliable interpretation of the experimental data would be the occurrence of a superposition of crystal fields due to different local symmetries around the ions induced by the fluorine doping, similar to the case of oxygen defects in the R compounds [18] and [19]. There are four nearest-neighboring oxygen ions around the ion; the ways of substitution of fluorine for oxygen could not all be the same. So different local charge environments and CF potentials around ions are expected. As in the analysis of [19], in this situation because of the lowered point symmetry, additional terms are needed to enter the Hamiltonian in (2) which would produce more CF states as in our experimental data. Although it is difficult to give a parameterization in terms of the current complex CF Hamiltonian based on our data, the clear observation of seven ground-state CF transitions would provide strong evidence for the above analysis.

Besides , the additional CF states induced by F-doping were also observed in [8]. However it is interesting to mention that, for another typical FeAs-1111 compound CeFeAs(O,F) [7], there is not much difference between the CF states in tetragonal CeFeAsO phase and phase. The reason could be that the line widths of major CF transitions in the Ce compounds are of the order of 10 meV, which is significantly larger than that in Nd compounds (mainly below 3 meV). So the presumably small effect of fluorine doping could not be detected.

3. Concluding Remarks

Inelastic neutron scattering measurements were carried out on polycrystalline samples of at different temperatures. Through the analysis of inelastic neutron spectra, up to seven ground-state CF excitations were identified. A crystal field energy level scheme for was proposed which shows quite different features comparing with that in tetragonal NdFeAsO phase. This indicates the occurrence of superposition of crystal fields due to the inhomogeneous ligand charge environments and lowered local point symmetries around the upon fluorine doping.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors appreciate the useful discussions with R. Osborn, Songxue Chi, Ninghua Tong, Huining Dong, Tianlong Xia, and Yinguo Xiao. This work is supported by NSFC (no. 11204373), the Fundamental Research Funds for the Central Universities, and the Research Funds of Renmin University of China.