Arxiv:1212.1837V2 [Cond-Mat.Mtrl-Sci] 14 Feb 2013 0.35), That Remains Controversial Because Muon Spin Ro- of the Magnetic Sites in the Ab Plane
Total Page:16
File Type:pdf, Size:1020Kb
Square-lattice magnetism of diaboleite Pb2Cu(OH)4Cl2 Alexander A. Tsirlin,1, 2, ∗ Oleg Janson,2 Stefan Lebernegg,2 and Helge Rosner2, y 1National Institute of Chemical Physics and Biophysics, 12618 Tallinn, Estonia 2Max Planck Institute for Chemical Physics of Solids, N¨othnitzerStr. 40, 01187 Dresden, Germany We report on the quasi-two-dimensional magnetism of the natural mineral diaboleite Pb2Cu(OH)4Cl2 with a tetragonal crystal structure, which is closely related to that of the frustrated 1 spin- 2 magnet PbVO3. Magnetic susceptibility of diaboleite is well described by a Heisenberg spin model on a diluted square lattice with the nearest-neighbor exchange of J ' 35 K and about 5 % of non-magnetic impurities. The dilution of the spin lattice reflects the formation of Cu vacancies that are tolerated by the crystal structure of diaboleite. The weak coupling between the magnetic planes triggers the long-range antiferromagnetic order below TN ' 11 K. No evidence of magnetic frustration is found. We also analyze the signatures of the long-range order in heat-capacity data, and discuss the capability of identifying magnetic transitions with heat-capacity measurements. PACS numbers: 75.30.Et,75.50.Ee,75.10.Jm,91.60.Pn I. INTRODUCTION The crystal structure of PbVO3 is a tetragonal deriva- tive of the perovskite type. The tetragonal distor- tion, which is essential for the quasi-two-dimensional Quantum spin systems show intricate low-temperature (2D) square-lattice magnetism, is only rarely observed phenomena of fundamental1{4 and even applied5,6 inter- in perovskites. Similar structures are found in BiCoO3 est. Despite strong quantum fluctuations that tend to (Ref. 21) and in the \supertetragonal" polymorph of impair and eventually destroy ordered spin states, most 22 BiFeO3, which is stabilized in thin films. However, none quantum magnets develop long-range magnetic order at of these compounds is suitable as a model system for the sufficiently low temperatures. In two and three dimen- 1 +3 quantum spin- 2 square-lattice model because both Co sions, the lack of the ordered state at zero temperature and Fe+3 bear much higher and largely classical spins. is only possible in frustrated magnets, where the com- Surprisingly, an appropriate structural analog of PbVO3 petition between magnetic couplings amplifies quantum can be found in a completely different family of materi- fluctuations. For example, systems of interest are those als. The Cu-based mineral diaboleite, Pb Cu(OH) Cl , based on spin- 1 in the square-lattice or kagom´e-lattice 2 4 2 2 has tetragonal crystal structure23,24 that can be derived 3 geometries. The simple square lattice with nearest- from the perovskite structure type. neighbor Heisenberg interactions (J1) is non-frustrated. The structure of diaboleite features isolated CuO4 pla- However, a second-neighbor coupling J2 frustrates the quettes that form magnetic layers in the ab plane. These system and leads to a rich phase diagram of the J1 − J2 layers are interleaved by Pb Cl slabs (Fig. 1). The re- frustrated square lattice model that was extensively stud- 2 2 7{11 lation to the parent ABX3 perovskite structure can be ied in the past. understood by writing the chemical composition of dia- Experimental search for the J1 − J2 systems on a boleite as Pb2(Cu)X6, where X = OH, Cl. This way, square lattice remains a challenging problem. While half of the B-type positions are vacant, with the result- − materials with ferromagnetic (FM) J1 can be prepared ing voids filled by the bulky OH anions. While this in a rather systematic fashion using the building blocks non-trivial transformation leaves little chemical similar- of V+4 phosphates,12,13 systems with antiferromagnetic ity to perovskites, the crystal structure of diaboleite still (AFM) J1 are less studied. In the limit of J2=J1 −! 0, features the space group P 4mm, same as in PbVO3 and +2 14{17 Cu -based coordination compounds and VOMoO4 PbTiO3. This tetragonal crystal structure is perfectly (Ref. 18) are excellent material prototypes of square- suited for a quasi-2D magnetism, because the magnetic lattice systems with weak magnetic frustration. The layers of Cu+2 atoms are well separated, with the short- stronger frustration with J2=J1 approaching 0.5 has been est interlayer distance of 5.5 A.˚ The robust tetragonal so far observed only in one material, PbVO3 (J2=J1 ' symmetry entails the perfect square-lattice arrangement arXiv:1212.1837v2 [cond-mat.mtrl-sci] 14 Feb 2013 0:35), that remains controversial because muon spin ro- of the magnetic sites in the ab plane. tation experiments detect the magnetic ordering below In the following, we present an experimental and com- 43 K,19 but neither thermodynamic measurements nor putational study of diaboleite. We show that this com- 1 neutron powder diffraction are capable of observing this pound is indeed a good material prototype of the spin- 2 20 magnetic transition. Latest theoretical results suggest square lattice. However, its J2=J1 ratio is very small, that above J2=J1 = 0:35 − 0:4 the J1 − J2 square-lattice hence the magnetic frustration is weak. We discuss the system enters the spin-liquid regime, and the magnetic origin of this low J2, and pinpoint another peculiarity, 7,9{11 order vanishes. While PbVO3 may lie on (or be close the dilution of the spin lattice through the formation of to) the boundary of this spin-liquid region, further sys- Cu vacancies that may have implications for experimen- 1 tems with J2=J1 ' 0:5 are highly desirable. tal studies of the diluted square lattices with spin- 2 . 2 a a b H O c Pb b Cl J ^ J' ^ Cu J2 2mm J1 FIG. 1. (Color online) An overall view of the diaboleite crystal structure (left panel) and its projection on the ab plane (right panel). Large spheres depict the Pb (violet) and Cl (brown) atoms located between the [Cu(OH)4] magnetic layers. Green spheres inside the plaquettes are Cu atoms. Small spheres show the O (orange) and H (gray) atoms. The middle panel shows the natural sample of diaboleite used in the present study. Note the translucent blue crystals of diaboleite and opaque white crystals of laurionite and phosgenite (see Sec. II for details). The outline of this paper is as follows. After summariz- tures were visualized using the VESTA software.27 ing experimental and computational methods in Sec. II, Magnetization measurements were performed with we report the detailed characterization of the natural Quantum Design MPMS SQUID magnetometer in the mineral sample (Sec. III A) and elaborate on possible de- temperature range 2 − 380 K in fields up to 5 T. Heat viations from the ideal crystal structure (Sec. III B). We capacity was measured with Quantum Design PPMS in further develop microscopic magnetic models for both the the temperature range 1:8 − 100 K in zero field and in ideal and distorted crystal structures (Sec. III C), and the fields of 3 T, 6 T, and 9 T. All measurements were apply them to thermodynamic properties of diaboleite performed on polycrystalline samples. (Sec. III D). Our results are discussed in Sec. IV, which Magnetic couplings in diaboleite, as well as differ- is followed by a summary and an outlook. ent aspects of its structural distortions, were analyzed by density-functional theory (DFT) band-structure cal- culations performed in the FPLO28 and VASP29 codes II. METHODS that implement basis sets of local orbitals and pro- jected augmented waves,30 respectively. Local density A natural sample of diaboleite (Mammoth-Saint approximation (LDA)31 and generalized gradient approx- Antony Mine, Pinal Co., Arizona, USA) was provided imation (GGA)32 for the exchange-correlation potential by the mineralogical collection of the Department of Ma- were used. Band-structure results are obtained for well- terials Research and Physics at Salzburg University. A converged k meshes with 336 points in the symmetry- visual inspection of this sample (Fig. 1, middle panel) irreducible part of the first Brillouin zone for the crys- identified translucent blue crystals of diaboleite mixed tallographic unit cell, and about 100 points for the su- with opaque white crystals of other minerals. Both blue percells. Residual forces in optimized crystal structures and white crystals were detached from the support and were below 0.01 eV/A.˚ Effects of strong electronic cor- separated manually. This way, a 100 mg batch of the relations in the Cu 3d shell were treated by introducing diaboleite crystals was obtained. the LDA band structure into an effective Hubbard model Phase composition of the samples was determined by (model approach), or within the mean-field DFT+U laboratory x-ray diffraction (Huber G670 Guinier cam- procedure, with the on-site Coulomb repulsion param- ◦ era, CuKα1 radiation, ImagePlate detector, 2θ = 3−100 eter Ud = 6:5 eV and the Hund's exchange parame- angle range). The batch of the blue crystals mostly ter Jd = 1 eV applying the around-mean-field double- contained diaboleite, whereas white crystals were iden- counting correction scheme.33 Details of the computa- tified as a mixture of phosgenite (Pb2Cl2CO3) and lau- tional procedure are described in Sec. III C. rionite (PbOHCl). Considering possible contamination Thermodynamic properties and N´eeltemperatures of of the natural samples with variable { both crystalline the spin models relevant to diaboleite were obtained from and amorphous { impurities, sample quality was further quantum Monte-Carlo (QMC) simulations based on the checked by bulk chemical analysis and high-resolution loop34 and dirloop_sse35 algorithms of the ALPS36 sim- x-ray diffraction (XRD). The chemical analysis was per- ulation package. Simulations were performed on 2D formed with the ICP-OES method.25 The high-resolution (L × L) and 3D (L × L × L=5) finite lattices with L ≤ 80 XRD data were collected at the ID31 beamline of the Eu- and L ≤ 50, respectively.