Jiajun Li(李佳俊), Giao Ngoc Phan, Xingyu Wang(王兴玉), Fazhi Yang(杨发枝), Quanxin Hu(胡全欣),Ke Jia(贾可),4, Jin Zhao(赵金), Wenyao Liu(刘文尧), Renjie Zhang(张任杰), Youguo Shi(石友国),4,5,Shiliang Li(李世亮),5, Tian Qian(钱天), and Hong Ding(丁洪)
1Beijing National Laboratory for Condensed Matter Physics and Institute of Physics,Chinese Academy of Sciences,Beijing 100190,China
2School of Physics,University of Chinese Academy of Sciences,Beijing 100190,China
3CAS Center for Excellence in Topological Quantum Computation,University of Chinese Academy of Sciences,Beijing 100190,China
4Center of Materials Science and Optoelectronics Engineering,University of Chinese Academy of Sciences,Beijing 100049,China
5Songshan Lake Materials Laboratory,Dongguan 523808,China
6Tsung-Dao Lee Institute,Shanghai Jiao Tong University,Shanghai 201210,China
Keywords: iron-based superconductor,angle-resolved photoelectron spectroscopy,uniaxial strain,symmetry breaking
BaFe2As2,well known for its rich phase diagram and interplay among various orders, provides a good material system for studying many novel properties, such as the superconducting mechanism of iron-based superconductors.[1–7]In BaFe2As2, as the temperature decreases, the sample undergoes an antiferromagnetic phase transition as well as a structural phase transition.[8,9]With the structural phase transition,the sample changes from an original single-domain state to two orthogonal twinned crystal phase compositions and exhibits a new electronic nematic order near the phase transition temperature,which has been observed in previous studies on iron-based superconductors by a scanning tunneling microscope(STM)and angle-resolved photoemission spectroscopy(ARPES) at low temperatures.[10–13]Considering the simultaneous appearance of nematic order, antiferromagnetic order and superconductivity, it is challenging to independently study these strongly correlated intertwined properties.Therefore,finding a method to independently study the nature of the nematic order without introducing the influence of other orders is highly desirable.In previous studies, a single domain state has been obtained in detwinned BaFe2As2by applying a small uniaxial strain wherein an electronic nematic order was observed.[14–16]By further applying a considerable uniaxial strain to BaFe2As2, the in-plane anisotropy of the resistance,which is closely related to the nematic order, would be further enhanced.[17,18]On the other hand,while a study on FeSe has revealed that the band positions corresponding to different electric orbitals shift differently under a larger strain,[19]it is still worth studying how band positions shift when a larger uniaxial strain is applied to BaFe2As2.
In this paper, we present the evolution of the electronic properties of BaFe2As2under uniaxial strains.The observed transport results demonstrate that uniaxial strains can enhance the nematic order.On the other hand,at theMYpoint,the upper electron band is drastically shifted toward the lower one.Consequently,the energy splitting between these two electron bands diminishes under a large uniaxial strain,even though the nematic order is enhanced.These surprising behaviors suggest that the band splitting at theMYpoint is not induced by the nematic order in BaFe2As2.
Single crystals of BaFe2As2were grown by a self-flux method.The samples with better single-crystal quality were selected,which could withstand larger strains and show clearer band dispersions.Considering that a strain would be applied along the[1 1 0]direction in theabplane,crystal structure of BaFe2As2in theabplane is shown in Fig.1(a).The upper inset in Fig.1(a)shows the strain device, which can continuously apply uniaxial strains to the mounted sample by turning the screw.With this design, anin situstrain along one desirable direction can be applied to the sample.Here,we use the Brillouin zone(BZ)notations corresponding to the true crystallographic unit cell of two Fe atoms,in whichΓ–MYis along the[1 1 0]direction andΓ–MXis along the[1−1 0]direction,as shown in Fig.1(b).
High-resolution ARPES measurements were carried out at the Institute of Physics,Chinese Academy of Sciences,with an R4000 analyzer with a He discharge lamp(hν=21.2 eV).The angular and momentum resolutions were set to 0.2°and 4 meV, respectively.To clearly study the electron band characteristics at the boundary of the BaFe2As2Brillouin zone,the direction of applied strain was set to be perpendicular to the electron-receiving slit of the ARPES analyzer,as shown in Fig.1(c).In order to keep the same cleave surfaces, all samples in this paper were cleavedin situat room temperature in a chamber with a vacuum of 1×10−10Torr, and As dimer surfaces were obtained.[20–22]For all of the ARPES measurements described below, the temperature is 14 K and the vacuum is better than 3×10−11Torr.
The in-plane uniaxial strains were applied to a single crystal of BaFe2As2by turning the screw of a strain device on which the sample is mounted, and thus the screw movements could be deemed to be the deformationsδLof the sample in the direction of applied strains.[23]The uniaxial strain strengthδL/Lis defined as the ratio of the screw movement to the original lengthLof the sample, which was measured directly before mounting it.
The changes in lattice constants indicated by XRD results prove that it is reliable to regard the screw movements as the deformations of the sample.The specific calculation method is as follows.Considering that the lattice constant along the[2 2 0] direction of the sample cannot be directly measured by XRD due to the sample holder,we indirectly calculated the deformation and lattice constant along the[2 2 0]direction by measuring the lattice constants along the[0 0 12]and[2 2 12]directions.Figures 1(d)and 1(e)show the corresponding XRD results.
It is known that the interplanar spacing for the strain-free sample, denoted byd, is 1.08618 ˚A along the[0 0 12]direction and 0.8573 ˚A along the[2 2 12]direction.Strains S1,S2,S3 and S4,which represent four gradually increasing uniaxial strains,were applied to the samples in turn.Taking strain S1 as an example,by substituting diffraction anglesθbefore and after applying uniaxial strain into the Bragg diffraction formula,2d×sin(θ)=nλ,the interplanar spacing for the strained sample,denoted bydS1,is 1.08497 ˚A along the[0 0 12]direction and 0.8551 ˚A along the[2 2 12]direction.
Fig.1.(a)Crystal structure of BaFe2As2 in the ab plane with a uniaxial strain applied along the[1 1 0]direction.Upper inset: picture of the in situ strain device(with sample).A red arrow indicates the direction in which strain is applied.Strain can be sequentially enhanced by turning the screw.(b)Schematic of the corresponding in-plane Brillouin zone for two Fe atoms per unit cell.(c)Schematic of the ARPES experiment.The direction of applied strain is perpendicular to the electron-receiving slit of the ARPES analyzer.[(d),(e)]XRD results at room temperature for(d)[0 0 12]and(e)[2 2 12]Bragg peaks under different applied strains of 0,S1 and S2 at room temperature.Shifts of the Bragg peak indicate the changes in lattice constants.(f)XRD results at 90 K for the sample under zero strain and under the minimum uniaxial strain S1.
According to the formula
whereh,kandldenote Miller indices,cS1for the[0 0 1]direction is derived as 13.0196 ˚A whendS1[0012]is 1.08497 ˚A,and thenaS1for the[1 1 0]direction is derived as 5.5574 ˚A.Thus,the change in lattice constant(aS1−a)/ais−0.49%.On the other hand,for strain S1,the screw movement is 0.0175 mm,the original lengthLof the sample is 3.5 mm and the ratio of the two is 0.5%, which is consistent with (aS1−a)/a0of−0.49%calculated from the XRD results.
Similarly, for strain S2, (aS2−a)/acalculated from the XRD results is−0.84%.The screw movement for strain S2 is 0.025 mm, the original lengthLof the sample is 2.8 mm and the ratio of the two is 0.89%, which is consistent with(aS2−a)/aof−0.84%.
Therefore, it is reliable to define uniaxial strain strengthδL/Las the ratio of screw movement to the original length of the sample at room temperature.For cases at low temperature we take the thermal expansion coefficients into account,which are 17.0×10−6for BeCu (C17200) in the strain device and 18.9×10−6for BaFe2As2; we found the difference caused by thermal expansion coefficients is 0.995 at 14 K.Thus,δL/Lat 14 K is−0.498%under strain S1 and−0.889%under strain S2.By the same method, strain strengthsδL/Lcorresponding to strains S3 and S4 can be obtained as 1.357%and 2.001%,respectively.
XRD at a low temperature of 90 K,for both the strain-free sample and the S1 strained sample,was measured as shown in Fig.1(f).It can be clearly observed that under zero strain,the[2 2 12] Bragg peak splits due to the twinning effect at low temperatures,showing two structural domains in a strain-free sample.For the S1 strained sample,there is only one[2 2 12]Bragg peak at low temperature,indicating that the sample remains in a single-domain state under the minimum uniaxial strain S1.It is worth mentioning that the strains in this report are much stronger than those in previous reports for detwinned samples that exhibit nematic order.[24]
With increasing strains,the in-plane transport anisotropy of BaFe2As2becomes stronger,as shown in Fig.2,which implies a stronger nematic order in strained BaFe2As2, in good agreement with previous reports.[17,18]In Fig.2,RbandRacorrespond to the resistances measured along the [1 1 0] direction and the [1−1 0] direction respectively on the same samples and normalized by the resistances at 300 K.
Fig.2.The in-plane resistance anisotropy indicated by(a)Rb −Ra and(b)Rb/Ra for five typical strained samples.The more heavily strained samples S4 and S3 exhibit stronger anisotropy than the moderately strained samples S2 and S1.
Fig.3.ARPES intensity plots measured around the Γ point for(a)strain-free and(b)–(e)increasingly strained BaFe2As2 samples.S1,S2,S3 and S4 represent four gradually increasing uniaxial strains.The three hole bands shift downward under the increasing uniaxial strains.
Fig.4.ARPES intensity plots measured around the MY point for(a)strain-free and(b)–(e)increasingly strained BaFe2As2 samples.The two electron bands shift downward under the increasing uniaxial strains.The upper electron band(α band)is more sensitive to the strains,while the lower one(β band)shifts slightly under the applied strains.The black dashed lines are guides for the eyes.
Figure 3 shows the evolution of electronic structure with increasing uniaxial strains around theΓpoint.A strain-free sample exhibits three hole bands near the Fermi level similar to previous reports.[15,25]With increasing uniaxial strains,the three hole bands shift downward by less than 13 meV,similar to the previous report on the variation of hole bands near theΓpoint under uniaxial strains in iron-based superconductors.[19]It is well known that the splitting of hole bands at theΓpoint is independent of the nematic order.[13]Therefore, the downward shifts of the three hole bands here can be considered to be a result of lattice distortions caused by the uniaxial strains.
Figure 4 shows the evolution of electronic structure with increasing uniaxial strains around theMYpoint.A strain-free sample exhibits two electron bands near the Fermi level,similar to previous reports.[15,25]With increasing uniaxial strains,the electron bands shift downward at different rates.The upper electron band(αband)shifts downward significantly while the lower one(βband)shifts slightly.As a result,the upper band shifts toward the lower one, and eventually these two bands merge into each other at theMYpoint.
In previous studies of the detwinned sample under a smaller uniaxial strain, it was reported that the positions of the bottoms of the two electron bands at theMYpoint do not shift significantly before or after the detwinning.[14,15]However, in this paper, when a larger strain is applied to the sample,it is directly observed from the spectra that the bottoms of two electron bands gradually merge under the increasing uniaxial strains.Moreover, the transport results in Fig.2 imply that a larger uniaxial strain leads to a stronger nematic order.That is,with the increasing strains,the splitting of two electron bands decreases although the nematic order becomes stronger.Therefore, it can be concluded that the splitting of these two electron bands is not induced by the nematic order.
In order to quantitatively analyze the effect of strains at theMYpoint,we extract the energy distribution curves(EDCs)and energy positions of the bottoms of the two electron bands under different strains,as shown in Figs.5(a)and 5(b),respectively.As the uniaxial strain increases from S1 to S4,the bottom of theβband shifts from−60 meV to−72 meV,which is consistent with the trend at theΓpoint as shown in Figs.3(a)–3(e).However,the energy position of the bottom of theαband shows a different shift when the uniaxial strain is applied to the sample.It gradually shifts from−10 meV to−33 meV — a change of 23 meV—when the uniaxial strain increases from S1 to S3.Finally, under the uniaxial strain of S4, the bottom of theαband merges with the bottom of theβband.It can be seen thatαband exhibits significantly more sensitivity to strain.
Fig.5.(a)Comparison of EDCs extracted from Figs.4(a)–4(e).Different colored lines indicate different applied strains of 0,S1,S2,S3 and S4,as shown in the diagram.(b)The energy positions of the bottoms of the α band and the β band and the energy distance δE between two bottoms of the electron band are plotted as functions of strain strength defined as the ratio between the deformation and the sample’s size δL/L in BaFe2As2.(c) Schematic of the band structure at the MY point in BaFe2As2.Also depicted is a diagram of the evolution in the two electron bands as the strain increases.
The difference between the bottoms of the two electron bands and their energy distanceδEare also extracted in Fig.5(b).When no strain was applied to the sample, the energy distanceδEbetween the two bottoms of the two electron bands was about 50 meV.When the strain is large enough,δEgradually decreases with increasing strains.This also proves that the two electron bands do not change uniformly under increasing strains.
In previous studies of uniaxial strains in the BaFe2As2system, it was found that the energy positions of different bands along two in-plane directions before and after detwinning changed only modestly.In this paper, varying uniaxial strains were applied to BaFe2As2samples.Under increasing strains, three hole bands at theΓpoint show similar shifts,which are affected by the lattice distortions.However, two electron bands at theMYpoint exhibit different properties: the shift of theβband is similar to that of the three hole bands at theΓpoint.In contrast,a more significant shift of theαband is observed.Under a large strain,the two electron bands tend to close together even though the nematic order is enhanced.This clearly suggests that the splitting of these two electron bands is not induced by the nematic order.
Acknowledgments
Project supported by the National Natural Science Foundation of China (Grant Nos.11888101 and U1832202), the Chinese Academy of Sciences (Grant Nos.QYZDB-SSWSLH043, XDB28000000, and XDB33000000), the K.C.Wong Education Foundation(Grant No.GJTD-2018-01),and the Informatization Plan of Chinese Academy of Sciences(Grant No.CAS-WX2021SF-0102).