Neutron diffraction has been a very prominent tool to investigate high-temperature superconductors, in particular through the discovery of an incommensurate magnetic signal known as stripes. We here report the findings of a neutron diffraction experiment on the superconductor (La,Sr)2CuO4, where a spurious signal appeared to be magnetic stripes. The signal strength was found to be strongly dependent on the neutron energy, peaking at . We therefore attribute the origin of this signal to be a combination of multiple scattering and crystal twinning. A forward calculation of the scattering intensity including these two effects almost completely recovers our experimental observations. We emphasize the need for employing such analysis when searching for ways to avoid spurious scattering signals.
Since their discovery more than 30 years ago [2], considerable scientific effort has been put into the quest of understanding the emergence of superconductivity in ceramic cuprates. During this formidable amount of research, much evidence points to the interplay between magnetism and superconductivity as playing an important role for the mechanisms that form the superconducting Cooper pairs [8].
Lanthanum-based cuprate superconductors possess a very complex electronic phase diagram with intertwined orders ranging from an antiferromagnetic Mott insulator to a metal. The long-range antiferromagnetic order in these compounds is destroyed by as little as strontium doping which introduces one electron hole per substituted atom [7]. Upon hole doping, these materials are known to form so called magnetic stripes, also known as incommensurate antiferromagnetic (IC AFM) order, depicted as insulating domains of antiferromagnetically arranged spins separated by one dimensional “rivers of charge” within the CuO2 layers. Spin stripes were first proposed by Tranquada et al. in [20] and later confirmed in (LSCO), as summarized by Yamada et al. [24]. In neutron scattering experiments, spin stripes are observed as a quartet of incommensurate peaks around the AFM reflection with a wavevector transfer , or crystallographical equivalent positions, such as .
Since their initial discovery, numerous neutron scattering experiments have focused on characterizing the static and dynamic spin correlations under different conditions in the attempt of determining their relationship with superconductivity. The evidence so far points, in broad terms, to a competition between superconductivity and the static magnetic order, as the second was seen to vanish at optimal (and higher) doping [10] and is enhanced by an applied magnetic field, which is known to suppress superconductivity [11,14].
From the early magnetic scattering studies performed on cuprates, multiple scattering was acknowledged as a possible contaminant of the stripe signal [20]. Here we present the results of a systematic neutron study aimed at determining the magnitude and origin of the double scattering events in a highly underdoped superconducting La1.93Sr0.07CuO4 single crystal. The findings presented here can be taken into account as guidelines in the planning phase of future experiments (i.e. by fine tuning the neutron beam energy) in order to obtain optimal measurements of magnetic stripes.
Experimental method
The sample used throughout this study is a highly underdoped crystal with nominal doping . The large single crystal has been grown using the Travelling Solvent Floating Zone method and, later on, cut into two pieces, one annealed under oxygen gas flow as described in Ref. [9] (3.4 g) which is the main sample of this paper and a second which was kept as-grown. The true (average) doping has been determined, before annealing, in a neutron diffraction experiment on the Morpheus diffractometer at the Paul Scherrer Institute (PSI), by following the temperature dependence of the strongly hole-doping dependent structural phase transition [22], from the high temperature tetragonal phase (I4/mmm space group number 139) to the low temperature othorombic one (Bmab space group number 64). The measurements revealed a transition temperature which corresponds to a doping of . The superconducting critical temperature of the as-grown and annealed crystals was established, through magnetic susceptibility measurements, to be for both samples which is in agreement with measurements from literature on samples of similar doping [13].
Elastic neutron scattering investigations were carried out on the cold-neutron triple axis spectrometer RITA II at PSI [15]. The set-up contained a vertically focusing graphite monochromator, 80′ horizontal collimation before the sample, a cooled beryllium filter after the sample in order to minimize the effects of higher order scattering contamination, followed by a graphite analyzer consisting of 9 blades, each 25 mm wide, and a coarse cross-talk-avoiding radial collimator before a position-sensitive detector. The sample was aligned with the CuO2 layers in the scattering plane, making it possible to access wavevectors in the -plane.
Throughout this paper the orthorhombic notation is used, meaning that a commensurate antiferromagnetic reflection would be observed at , and crystallographically equivalent positions, and the stripe signal, due to its incommensurate magnetic modulation, should be visible at , and equivalent positions, where in this low-doping regime [24]. In this notation, the 10 K lattice parameters of a doped crystal (a value very similar to the doping level of our crystal) were earlier reported as , , [17].
To increase the efficiency of grid scanning around the antiferromagnetic reflection, we made use of the monochromatic imaging mode of the multi-blade analyzer. In this configuration, all 9 blades elastically scatter neutrons of the same energy, but due to the difference in scattering angle, each blade maps a slightly different portion of reciprocal space [1]. In this configuration, the effective collimation is determined by the width of the sample, the (narrow) analyzer, and the distance between them. This leads to a collimation sequence of open–––open. However, the results presented throughout the paper are not particular to this set-up. The same observations can be reproduced on a standard triple axis spectrometer with a single-blade analyzer.
Results
Initial neutron diffraction measurements were performed with incident and outgoing energies of 4.6 meV, corresponding to , where and are the wave vectors of the initial and final neutron beams, respectively. This is a standard configuration of cold triple-axis instruments that allows for an optimal filtering of higher order scattering by cooled polycrystalline beryllium. The measurements revealed the presence of a peak at the antiferromagnetic position, which has previously not, to our knowledge, been reported in homogeneously doped systems of . In addition to this central peak, a number of satellite peaks at incommensurate positions were also measured, as depicted in Fig. 1(a), with incommensurability , which is much smaller than the value expected from the dopant concentration: .
(a) Elastic scattering grid scan around the reflection taken at 2.5 K on LSCO with (annealed crystal). Temperature dependence of the (b) peak intensity obtained from Gaussian fits, (c) incommensurability for both annealed and as-grown crystals of the same doping and (d) peaks position. The incommensurability is defined as half the distance between the two side peaks marked as ‘left’ and ‘right’ in sub-figure (a).
The temperature dependence of the signal, shown in Fig. 1, was however, significantly different from that of static magnetic stripes. First of all, the intensity persisted up to room temperature with only a slight decrease as a function of increased temperature in contrast with genuine spin stripes signal that usually correlates with [6]. Secondly, a decrease in incommensurability with increasing temperature was measured (Figs 1(c, d)), whereas the spin modulation in stripe signals was earlier found to be invariable to temperature variations. These two temperature effects were invariant to sample annealing (Fig. 1(c)).
Because the temperature dependence of the peaks resembles that of the tetragonal-to-orthorhombic phase transition around 360 K, multiple scattering involving reflections allowed only in the orthorhombic phase was soon proposed as the origin of the signal. Following the textbook by Shirane, Shapiro, and Tranquada [19], the strict condition for multiple scattering to occur is that more than one reciprocal lattice point lies on the surface of the Ewald sphere. Its presence can be tested by either rotating the crystal around the contaminated scattering vector or by changing the neutron energy of the diffraction experiment (i.e. changing the radius of the Ewald sphere). It should be mentioned that in a real experiment, where the incident beam is not perfectly monochromatized, the Ewald sphere should be seen rather as a spherical shell, the thickness of which corresponds to the spread in wave vectors of the neutron beam.
The method of changing wavelength of the diffraction experiment was employed in our experiments at RITA-II and the data is shown in Fig. 2. Our findings clearly support the hypothesis of multiple scattering, since the incommensurate side peaks were only found at certain energies (4.5 meV to 4.7 meV) and the intensity of the central peak greatly decreased when the initial and final energies were modified. The more than one order of magnitude difference in intensity can be visualised in Fig. 3 where 1 dimensional cuts thought the grids are plotted for two different energies.
Figure 2 also shows the result of a virtual experiment by the McStas ray-tracing simulation package [16,23] of a hypothetical crystal showing single scattering at the reflection, in effect providing the resolution function of the instrument. Through the same McStas simulation, it was found that the width of incoming energy, , on the sample was 0.13 meV (FWHM). This value is the relevant number to determine the effective thickness of the Ewald sphere and is not equal to the energy resolution of the spectrometer (which is 0.20 meV FWHM), because the energy resolution of the secondary spectrometer is here not taken into account.
Diffraction measurements performed around the reflection at 125 K. The caption in each figure indicates the value of the initial and final energy of the neutron beam. Sub-figure (a) shows in red the direction of the two diagonal scans plotted in 2D in Fig. 3. Panel (h) shows the simulated resolution of the instrument as presented in the text.
Elastic diagonal scans taken at different values of the initial and final energy, noted in the caption, in the direction exemplified in Fig. 2(a). The solid line represents a Gaussian fit to the data.
Discussion
Double scattering
Due to the fact that the signal appears in a very narrow energy range and because we know the space groups of the crystallographic phases allowed in this sample, we were able to determine the reflections involved in a potential double scattering event. Following the treatment by Shirane, Shapiro, and Tranquada [19], we construct the Ewald sphere corresponding to our experiment as shown in Fig. 4(a). In our particular case, the condition for double scattering is fulfilled if another allowed reflection of the low temperature orthorombic (LTO) phase lies on the surface of the sphere. This poses two requirements: (1) the existence of a secondary reflection wavevector of the same length as ; (2) the two scattering wavevectors and should correspond to allowed reflections in the LTO phase (space group no. 64 in the unconventional Bmab setting [5]). is geometrically defined as:
Using the scattering angle to calculate the angle and knowing the lattice parameters corresponding to the LTO phase, we are able to rewrite the components of and as:
where , and are the Miller indices of the allowed reflections in the LTO phase. We have investigated indices with values between and 3. This decomposition allows us to calculate the corresponding in the energy range 4.5–4.6 meV for all the allowed reflections of the LTO phase:
Checking the above mentioned condition (1), meaning comparing the values of and , two pairs of reflections appear to lie on the surface of the Ewald sphere (equivalent to ) and (equivalent to ). However, condition (2) is only fulfilled for the first pair since both and require the disallowed reflections ( and respectively) in order to produce intensity at . In contrast, and reflections pair with and respectively, which are only allowed in the orthorombic phase and not in the tetragonal one. This explanation is thus in direct agreement with the specific temperature dependence of the signal as presented in Fig. 1.
Figure 4(b) shows the excellent correlation between the fulfilment of condition (1) for the reflection and the peak of the measured integrated intensity. The red symbols close to correspond to the fulfilment of the first condition for multiple scattering: for the which takes place in the same energy range as the observed increase in scattered intensity. The reflection behaves in an equivalent manner.
(a) Graphical representation of an instance of double scattering within the Ewald sphere. and define the incoming and outgoing neutron beam, and denote the tested reflections allowed in the LTO phase and is the secondary reflection wavevector. (b) Blue symbols show the integrated intensity of the measured signal (see Fig. 2) and follow the left side y-axis. The red symbols follow the right side y-axis and show the calculated difference between the magnitude of the two scattering wavevectors and for the case and . The x-axis shows the incoming and outgoing energy of the neutron beam ().
The effect of twinning
We attribute the complex pattern of multiple peaks around the reflection, as it can be observed in Fig. 1(a), to a combination of the double scattering mechanism, presented above, and twinning of the reflections. It is well known that LSCO crystals in the low temperature orthorombic phase are composed of two sets of two twin orientations sharing the or planes [3]. This co-existence of a 4-phase structure with equal lattice parameters but different orientations of the planes is observed in neutron scattering experiments as a peak splitting of the Bragg reflections along the h and k directions (as exemplified in Fig. 1 from Ref. [3]).
Thus, if we consider all possible combinations of , where and are the wavevectors of the 4 satellite peaks around and respectively (i.e. and ), we obtain a total of 16 peaks appearing around the reflection (as depicted by the red circles in Fig. 5). Note that the incommensurability (t) of the satellite peaks is related to the twin splitting ( [3]) as , in the small angle approximation, and it does not affect the l value since there is no splitting along the c-axis. In Fig. 5(b), we plot the calculated pattern obtained as Gaussian distributions centred at the calculated peak positions and with amplitudes tuned to match the measured data shown in Fig. 5(a). The twin splitting used in the calculations corresponds to a lattice parameters ratio , a value slightly higher than the one from the literature [17]. A similar figure could likely be produced from a true simulation of the scattering intensities by modifying the volume fractions of the 4 twin phases. However, we have chosen not to address this involved optimisation scheme, since we do not think it would add to the general understanding. The good agreement between the measured data and the calculated reconstruction implies that the observed multiple peaks pattern can be regarded as a superposition of the twinning patterns of all the reflections that simultaneously intersect the surface of the Ewald sphere.
(a) Elastic scattering grid scan around the reflection taken at 2.5 K on LSCO with annealed. (b) Calculated intensity patterns of the 16 peaks originating from the twinning of the sample. The red circles indicate the position of all reflections. The intensity of each Gaussian distribution around the calculated position is tuned to match the measured data.
Earlier studies of stripes in cuprates
Although numerous studies from the literature follow the behaviour of static spin stripes in LSCO superconductors as a function of doping, temperature and applied magnetic field, this spurious process has not been addressed in any of the publications to our knowledge. There are two main ways in which the additional double scattering signal went unreported: (a) at higher doping () the true magnetic signal appears at higher values compared to the spurious one, thus, depending on the scanning direction, one is able to measure the stripe signal without detecting the spurion (for example by performing scans of the sample orientation, A3). Secondly, as we have argued throughout this paper, by tuning the incoming and outgoing neutron beam energy the contamination can easily be avoided. Indeed, many of the experiments presented in the literature have been carried out at energies outside the 4.5–4.7 meV regime [4,12–14,18].
That multiple scattering can make forbidden peaks appear was already mentioned by the first discovery of stripes [20], where Tranquada et al. reported a neutron energy (13.9 meV) almost free of spurious scattering at . In this respect, not all findings in the present article are new as such. However, the effect of multiple scattering, and in particular the effect that twinning causes IC peaks to appear, is in our work thoroughly documented, as a warning and explanation to future experimentalists.
Conclusion
Our neutron diffraction measurements and calculations support the hypothesis of a double scattering event, originating from the - and -type reflections, corroborated with twinning of the crystal, as the cause for the observed varying intensity of the reflection and the appearance of incommensurate peaks around it.
Because the intensity of the static magnetic stripe signal is very low in cuprates in general and in LSCO in particular [21], the contribution of multiple scattering should always be carefully considered and the neutron energy should be tuned as to minimize this intrinsic effect, which could easily lead to detrimental pollution of the experimental data. In particular for cold-neutron experiment on the LSCO-family of compounds, energies around 4.6 meV should be avoided. As alternatives, energies of 3.7 meV (the BeO filter edge) or 5.0 meV (the Be-filter edge) can be strongly recommended. Our findings explain why, historically, best diffraction results on stripes in cuprate superconductors were found for exactly these energies. Although the multiple scattering explanation of spurious signals was noted early by Tranquada et al., we have here presented a comprehensive investigation of this effect that we now believe to be fully understood.
We emphasize the need for applying these considerations for the planning of neutron diffraction experiments on similar systems, in order to minimize contamination by spurious signals.
Footnotes
Acknowledgements
We thank Jean-Claude Grivel and Maria Retuerto for valuable assistance with the preparation of the single crystal samples.
We are grateful to Philippe Bourges for assisting us with a related test experiment at the LLB and for bringing the idea of multiple scattering to our attention.
We thank the students from the University of Copenhagen Neutron Scattering Course as well as Linda Udby, Ursula B. Hansen, and Sonja Holm-Dahlin for participating in the initial measurements.
We thank Anne Bartholdy and Henrik Jacobsen for assisting in the neutron scattering characterization of the structural phase transition of the crystals.
This work was supported by the Danish Agency of Science and Technology through DANSCATT. Ana Elena Ţuţueanu was supported by a Ph.D. grant from the Institute Laue-Langevin.
The neutron scattering data were obtained at the SINQ neutron source, Paul Scherrer Institute, Switzerland.
References
1.
C.Bahl, K.Lefmann, A.Abrahamsen, H.Ronnow, F.Saxild, T.Jensen, L.Udby, N.Andersen, N.Christensen, H.Jakobsen, T.Larsen, P.Hafliger, S.Streule and C.Niedemayer, Inelastic neutron scattering experiments with the monochromatic imaging mode of the RITA-II spectrometer, Nuclear Instruments and Methods B246(2) (2006), 452–462. doi:10.1016/j.nimb.2006.01.023.
2.
J.G.Bednorz and K.A.Müller, Possible high Tc superconductivity in the Ba–La–Cu–O system, in: Ten Years of Superconductivity: 1980–1990, Springer, 1986, pp. 267–271. doi:10.1007/978-94-011-1622-0_32.
3.
M.Braden, G.Heger, P.Schweiss, Z.Fisk, K.Gamayunov, I.Tanaka and H.Kojima, Characterization and structural analysis of twinned crystals by neutron diffraction, Physica C: Superconductivity191(3–4) (1992), 455–468. doi:10.1016/0921-4534(92)90944-8.
4.
J.Chang, C.Niedermayer, R.Gilardi, N.Christensen, H.Rønnow, D.McMorrow, M.Ay, J.Stahn, O.Sobolev, A.Hiesset al., Tuning competing orders in cuprate superconductors by the application of an external magnetic field, Physical Review B78(10) (2008), 104525. doi:10.1103/PhysRevB.78.104525.
5.
J.K.Cockcroft, Crystallographic space group diagrams and tables.
6.
T.Croft, C.Lester, M.Senn, A.Bombardi and S.Hayden, Charge density wave fluctuations in and their competition with superconductivity, Physical Review B89(22) (2014), 224513. doi:10.1103/PhysRevB.89.224513.
7.
E.Dagotto, Correlated electrons in high-temperature superconductors, Reviews of Modern Physics66(3) (1994), 763–840. doi:10.1103/RevModPhys.66.763.
8.
E.Fradkin, S.A.Kivelson and J.M.Tranquada, Colloquium: Theory of intertwined orders in high temperature superconductors, Reviews of Modern Physics87 (2015), 457–482. doi:10.1103/RevModPhys.87.457.
9.
M.Fujita, K.Yamada, H.Hiraka, P.Gehring, S.Lee, S.Wakimoto and G.Shirane, Static magnetic correlations near the insulating-superconducting phase boundary in , Physical Review B65(6) (2002), 064505. doi:10.1103/PhysRevB.65.064505.
10.
M.-H.Julien, Magnetic order and superconductivity in : A review, Physica B: Condensed Matter329 (2003), 693–696. doi:10.1016/S0921-4526(02)01997-X.
11.
S.Katano, M.Sato, K.Yamada, T.Suzuki and T.Fukase, Enhancement of static antiferromagnetic correlations by magnetic field in a superconductor with , Physical Review B62(22) (2000), R14677. doi:10.1103/PhysRevB.62.R14677.
12.
B.Khaykovich, S.Wakimoto, R.Birgeneau, M.Kastner, Y.Lee, P.Smeibidl, P.Vorderwisch and K.Yamada, Field-induced transition between magnetically disordered and ordered phases in underdoped , Physical Review B71(22) (2005), 220508. doi:10.1103/PhysRevB.71.220508.
13.
M.Kofu, S.-H.Lee, M.Fujita, H.-J.Kang, H.Eisaki and K.Yamada, Hidden quantum spin-gap state in the static stripe phase of high-temperature superconductors, Physical Review Letters102(4) (2009), 047001. doi:10.1103/PhysRevLett.102.047001.
14.
B.Lake, H.Rønnow, N.Christensen, G.Aeppli, K.Lefmann, D.McMorrow, P.Vorderwisch, P.Smeibidl, N.Mangkorntong, T.Sasagawaet al., Antiferromagnetic order induced by an applied magnetic field in a high-temperature superconductor, Nature415(6869) (2002), 299–302. doi:10.1038/415299a.
15.
K.Lefmann, C.Niedemayer, A.Abrahamsen, C.Bahl, N.Christensen, H.Jacobsen, T.Larsen, P.Hafliger, U.Filges and H.Ronnow, Realizing the full potential of a RITA spectrometer, Physica B385–386 (2006), 1083–1085. doi:10.1016/j.physb.2006.05.372.
16.
K.Lefmann, P.K.Willendrup, L.Udby, B.Lebech, K.Mortensen, J.O.Birk, K.Klenø, E.Knudsen, P.Christiansen, J.Saroun, J.Kulda, U.Filges, M.Konnecke, P.Tregenna-Piggott, J.Peters, K.Lieutenant, G.Zsigmond, P.Bentley and E.Farhi, Virtual experiments: The ultimate aim of neutron ray-tracing simulations, J. Neutron Research16 (2008), 97–111. doi:10.1080/10238160902819684.
17.
P.Radaelli, D.Hinks, A.Mitchell, B.Hunter, J.Wagner, B.Dabrowski, K.Vandervoort, H.Viswanathan and J.Jorgensen, Structural and superconducting properties of as a function of Sr content, Physical Review B49(6) (1994), 4163–4175. doi:10.1103/PhysRevB.49.4163.
18.
A.T.Rømer, J.Chang, N.B.Christensen, B.Andersen, K.Lefmann, L.Mähler, J.Gavilano, R.Gilardi, C.Niedermayer, H.M.Rønnowet al., Glassy low-energy spin fluctuations and anisotropy gap in La1.88Sr0.12CuO4, Physical Review B87(14) (2013), 144513. doi:10.1103/PhysRevB.87.144513.
19.
G.Shirane, S.M.Shapiro and J.M.Tranquada, Neutron Scattering with a Triple-Axis Spectrometer: Basic Techniques, Cambridge University Press, 2002. doi:10.1017/CBO9780511534881.
20.
J.Tranquada, B.Sternlieb, J.Axe, Y.Nakamura and S.Uchida, Evidence for stripe correlations of spins and holes in copper oxide superconductors, Nature375(6532) (1995), 561–563. doi:10.1038/375561a0.
21.
A.E.Ţuţueanuet al., Magnetic stripes and fluctuations in strongly underdoped (La,Sr)2CuO4, 2020 (in preparation).
22.
S.Wakimoto, S.Lee, P.M.Gehring, R.J.Birgeneau and G.Shirane, Neutron scattering study of soft phonons and diffuse scattering in insulating La1.95Sr0.05CuO4, Journal of the Physical Society of Japan73(12) (2004), 3413–3417. doi:10.1143/JPSJ.73.3413.
23.
P.K.Willendrup and K.Lefmann, McStas (i): Introduction, use, and basic principles for ray-tracing simulations, Journal of Neutron ResearchPre-press (2019), 1–16. doi:10.3233/JNR-190108.
24.
K.Yamada, C.Lee, K.Kurahashi, J.Wada, S.Wakimoto, S.Ueki, H.Kimura, Y.Endoh, S.Hosoya, G.Shiraneet al., Doping dependence of the spatially modulated dynamical spin correlations and the superconducting-transition temperature in , Physical Review B57(10) (1998), 6165–6172. doi:10.1103/PhysRevB.57.6165.