Helicoidal magnetic structure and ferroelectric polarization in. Cu 3 Nb 2 O 8
|
|
- John Jackson
- 5 years ago
- Views:
Transcription
1 Helicoidal magnetic structure and ferroelectric polarization in Cu 3 Nb 2 O 8 Zheng-Lu Li, 1 M.-H. Whangbo, 2 X. G. Gong, 1, * and H. J. Xiang 1, * 1 Key Laboratory of Computational Physical Sciences (Ministry of Education), State Key Laboratory of Surface Physics, and Department of Physics, Fudan University, Shanghai , People s Republic of China 2 Department of Chemistry, North Carolina State University, Raleigh, North Carolina , USA ABSTRACT We investigate the origin of the coplanar helicoidal magnetic structure and the ferroelectric polarization in Cu 3 Nb 2 O 8 by combining first-principles calculations and our spin-induced ferroelectric polarization model. The coplanar helicoidal spin state comes from the competition between the isotropic exchange interactions, and the ferroelectric polarization from the symmetric exchange striction with slight spin canting. However, the direction of the polarization is not determined by the orientation of the spin rotation plane. PACS numbers: t, b, Et, e I. INTRODUCTION In recent years, multiferroics 1 with ferroelectric polarization P induced by a magnetic order have attracted much research interests. 1-3 In these materials the ferroelectric and magnetic orders can coexist to present a strong magnetoelectric effect, leading to their potential applications to spintronics. 4 The ferroelectric polarization P driven by spin spiral order are often explained by the spin current model 5 or equivalently the inverse Dzyaloshinskii-Moriya (DM) interaction model, 6 but these models fail in some cases, 7-10 particularly, in predicting the polarization induced by helical magnetic structure. 7,8 Recently, we developed a general model for describing P induced by a helical magnetic structure, and have successfully explained the origin of P in those multiferroics that are not accounted for by the abovementioned two models. 7,8 Johnson et al. have found some puzzling multiferroics where P is perpendicular to the spin rotation plane, which is not explained by the spin current model or the inverse DM interaction model. 9,10 To describe the microscopic origin of P, they proposed a phenomenological ferroaxial coupling mechanism, namely, P = γσa, where γ is a coupling constant, σ represents the macroscopic chirality and A is the macroscopic axial vector that depends on the crystal structure 9,10 This mechanism indicates that P is related to both the crystal and the magnetic structures. Recently Johnson et al. reported that the multiferroic 1
2 CaMn 7 O 12 has giant polarization and used the ferroaxial mechanism to explain the origin of P. 9 However, by using our general model, we showed that P originates from the symmetric exchange striction rather than spin-orbit coupling (SOC) in CaMn 7 O Johnson et al. also reported another multiferroic Cu 3 Nb 2 O 8 and showed that Cu 3 Nb 2 O 8 has a non-collinear coplanar helicoidal magnetic structure with the propagation vector k = (0.4876, , ) below a transition temperature T 24 K, where the polarization with a magnitude of μc/cm 2 exists, being almost perpendicular to the spin rotation plane, with an angle ~14 to the vector normal to the spin rotation plane. 10 However, the polarization P was not perpendicular to k, contradicting the prediction of the spin current or inverse DM interaction model, and thus is suggested to originate from the ferroaxial coupling mechanism of Johnson et al. In the present work, by combining first-principles calculations with our spin-induced ferroelectric polarization model, 7,8 we extract the exchange parameters and the polarization coefficients using the mapping analysis. 7,8,11,12 The magnetic structure of Cu 3 Nb 2 O 8, including the spin rotation plane and the propagation vector k, is determined by the competition between isotropic exchange interactions. Instead of SOC, the P is determined by the symmetric exchange striction with the slight spin canting due to the anisotropy of the system. The direction of P is not determined by the orientation of the spin rotation plane, in contradiction to the ferroaxial explanation. 9,10 Our study has successfully explained the origin of the magnetic structure and the ferroelectric polarization in Cu 3 Nb 2 O 8. This work is organized as follows. In Sec. II, we describe our first-principles computational methods. In Sec. III, we present our results on the magnetic structure and the ferroelectric polarization. We demonstrate that the competition between the isotropic exchange interactions is responsible for the magnetic structure, and the symmetric exchange striction with spin canting is the origin of the ferroelectric polarization. The main conclusions are summarized in Sec. IV. II. CALCULATION DETAILS AND GEOMETRICAL STRUCTURES We performed first-principles calculations based on density functional theory (DFT) by using the Vienna ab initio simulation package 13 (VASP) with the projector-augmented-wave method, 14,15 the generalized gradient approximation (GGA) by Perdew, Burke, and Ernzerhof 16 for exchange-correlation functional. We employed the GGA plus on-site repulsion (U) method 17 (GGA+U) to describe the strong electron correlation of Cu 3d orbitals, with U = 6 ev and J = 1 ev. 18,19 The plane-wave cut-off energy set 500 ev and the total energy was converged to 10-6 ev. The Berry phase method 20 was used to calculate the ferroelectric polarization. We also include the SOC effect in some calculations (GGA+U+SOC), to study the origin of the polarization, and in these cases, the spin direction are fixed but their magnitude are allowed to relax. The crystal structure of Cu 3 Nb 2 O 8 has the space group of P1, the lattice parameters and the atomic coordinates used for our calculations were taken from the experiment. 10 There are two non-equivalent Cu atoms, Cu1 and Cu2. In the unit cell, there is one Cu1 atom at the inversion center and the other two Cu2 atoms (marked as Cu2_1 and Cu2_2) are equivalent. The structure of Cu 3 Nb 2 O 8 is generally made up of Cu 3 O 8 units, and in one Cu 3 O 8 unit, the eleven atoms are almost in one plane [Fig. 1(a)]. The Cu 3 O 8 units connect to form the staggered Cu-O chains along the a axis, and further the chains extend to become the Cu-O layers, separated by the non-magnetic Nb atoms in between these layers [Fig. 1(b)]. Our mapping analysis 7,8,11,12 requires 2
3 the use of large supercells to avoid the interactions between spin dimers in the periodic cells, and thus we use 2 2 2, 2 2 3, and supercells to extract the parameters of exchange interactions and the ferroelectric polarization. The total polarization simulating the experimental magnetic structures from first-principles calculations was derived with supercell. In all abovementioned supercells, we use the k-points mesh. Convergence tests have been performed to assure the high accuracy of the energy term we are concerned with, which is the energy difference between different spin states. The spin exchange parameter, for example, J 1 equals mev (see Table I) using the k-point mesh, and mev using the k-point mesh. Thus, for all systems we use the k-point mesh or even a larger k-point set to keep our results reliable. III. RESULTS AND DISCUSSIONS A. Magnetic Structure We first consider the formation of non-collinear magnetic structure. We have considered totally ten exchange paths of the Cu pairs [Fig. 2(a)], i.e., all those Cu pair distances less than 6 Å in Cu-O chains and those less than 5 Å between the Cu-O chains and layers. The energy of exchange interactions can be written as follows: E ex = J ij S i S j <i,j> (1) where we take S i = 1, and J ij is the effective exchange constant between spins at the i and j sites (i.e., J ij S i S j ). We employ the energy mapping analysis 11,12 to determine the value of the exchange interactions J values from first-principles GGA+U calculations by using large enough supercells. Table I lists the Cu Cu distances (from short to long) of the 10 spin exchange paths and their J values extracted from the mapping analysis. J 1 and J 10 are the interactions just in one Cu 3 O 8 unit; J 3, J 4, J 8 and J 9 are the interactions between different Cu 3 O 8 units in one Cu-O chain; J 2 and J 7 are the interactions between different Cu-O chains in one Cu-O layer; J 5 and J 6 are the interactions between different Cu-O layers. Considering the symmetry, J 2, J 3, J 6, J 8 and J 10 have the inversion symmetry with the inversion center locates at the mid-point of the pairs, and the other pairs have no inversion symmetry [Fig. 2(a)]. Our results show that the intrachain interactions are dominant and those interchain and interlayer ones play a minor role. The magnetic structure of Cu 3 Nb 2 O 8 10 shows that the three Cu atoms in one Cu 3 O 8 unit have almost ferromagnetic (FM) spin configuration (hereafter the Cu 3 O 8 unit in FM state is referred to as FM-u), and have almost an antiferromagnetic (AFM) spin arrangement with those in adjacent Cu 3 O 8 units, leading to the AFM Cu-O chain (hereafter AFM-c) [Fig. 1(a)]. Let us first look at the spin structure of the Cu-O chain qualitatively. If we consider the interactions only within the Cu 3 O 8 unit, i.e., the nearest neighbor (NN) FM J 1 = mev and the next nearest neighbor (NNN) AFM J 10 = 2.78 mev, the three spins in the Cu 3 O 8 unit would show a completely non-collinear spin state, rather than a FM arrangement found in experiment. 10 Thus, we have to take other intra-chain interactions into consideration. We consider J 1, J 4, J 9 and J 10 in determining the spin structure of the Cu 3 O 8 unit and the Cu-O chain, and these four interactions are the strongest among the ten (see Table I). The arrangement of these four interactions is shown in Fig. 3
4 2(b) together with the Cu1, Cu2_1, Cu2_2 and Cu2_1 atoms. In the triangle made up of J 1, J 4 and J 9, the Cu1 and Cu2_1 spins would be coupled ferromagnetically not only through J 1, but also through J 4 and J 9. Because both J 4 and J 9 are AFM, this provide an enhanced effective FM coupling between Cu1 and Cu2_1 with J FM = J 1 - J 4 - J 9 = mev (see APPENDIX for details). If J FM > 2 J 10 = 5.56 mev, the spins of the three Cu atoms in one Cu 3 O 8 unit would have a FM arrangement to achieve the lowest energy, while the spin of Cu2_1 would be antiparallel to that of Cu1 and Cu2_1 without any constraints due to strong AFM J 4 and J 9 leading to an AFM arrangement between neighboring Cu 3 O 8 units (see APPENDIX for details). Consequently, the magnetic structure acquires the FM-u and the AFM-c arrangements, in agreement with experiment. In fact, the smaller J 2 and J 7 would even strengthen this result in a similar manner, leading to the effective J FM = J 1 - J 4 - J 9 - J 2 - J 7 = mev. To quantitatively investigate the magnetic structure of Cu 3 Nb 2 O 8, we performed Monte Carlo (MC) simulations using the J values derived from the GGA+U calculations (Table I). Similar spin structures (i.e., FM-u and AFM-c) are obtained when we consider all the ten J s, when the interlayer exchanges are deleted, and when all the interlayer and interchain spin exchanges are deleted. Thus, our qualitatively analysis based only on the strongest intrachain interactions is reasonable, and the spin structure of the Cu-O chain is basically determined by the intrachain interactions. Furthermore, by considering all the ten J s, with the supercell, we derive perfect the spin rotation plane from MC simulations. Given the a axis parallel to the x axis with the b axis in the xy plane, the normal vector of the spin rotation plane is specified as n expr. = (θ, φ). Experimentally, it was found that θ = 75.5 and φ = Our MC simulations led to two normal vectors n 1 = (85.7, -7.2 ) and n 2 = (19.7, -1.3 ), which are different from each other and from the experimental values. This is expected because the J s, being isotropic interactions, cannot induce any preferred spin direction in the space. After determining the spin rotation plane, we focus on the propagation vector k. The experiment shows k = (k a*, k b*, k c* ) = (0.4876, , ), which shows an important relation that k a* k b* + k c* 0.5. We reproduce this relation in the following on the basis of the FM-u and AFM-c magnetic structures. We employ the notation (n a, n b, n c ) to mark the unit cells in the direct lattice, e.g., (0, 0, 0) for the unit cell at the origin, and (1, 0, 0) for the next unit cell along the a axis, etc. For the Cu1 atom, which is at the inversion center in the (0, 0, 0) cell, the corresponding Cu 3 O 8 unit has the Cu2_1 in the (-1, 0, 0) cell and Cu2_2 in the (0, -1, -1) cell [see Fig. 2(a)]. The spins in the spiral state of Cu 3 Nb 2 O 8 are described by the expression 10 S i = R cos(k R L + ξ i ) + I sin(k R L + ξ i ) (2) where i denotes the type of the Cu atom, R and I are two orthogonal vectors determining the spin rotation plane, R L = (n a, n b, n c ) represents the lattice vector in the real space, and ξ i is a relative phase (ξ Cu1 = 0, ξ Cu2_1 = 1.03π and ξ Cu2_2 = 1.05π from experiment). 10 Now, we introduce the in-plane spin rotation angle φ i = k R L + ξ i = 2πn a k a* +2πn b k b* + 2πn c k c* +ξ i. In this Cu 3 O 8 unit, we have φ Cu1 = ξ Cu1, φ Cu2_1 = 2πk a* + ξ Cu2_1, φ Cu2_2 = 2πk b* 2πk c* + ξ Cu2_2. Because ξ i just represents a relative phase, we choose ξ Cu1 = 0 and thus φ Cu1 = 0. The FM-u magnetic structure has inversion symmetry, so we have ξ Cu2_1 = ξ Cu2_2, and the FM nature requires φ Cu1 = φ Cu2_1 = φ Cu2_2 = 0, so that k a* = ξ Cu2_1 /2π, and k b* + k c* = ξ Cu2_2 /2π. Because of the AFM-c spin arrangement, k a* = 0.5, as a result, ξ Cu2_1 = ξ Cu2_2 = π, and k b* + k c* = 0.5, consistent with the experimental values. So far, the above results and analyses show that the coplanar magnetic 4
5 structure (including the spin rotation plane, the FM-u and the AFM-c) is a result of the competition between the isotropic exchange interactions J s, and that the propagation vector k consistent with the experiment is derived by analyzing the crystal and magnetic structures of Cu 3 Nb 2 O 8. B. Ferroelectric Polarization Now we turn to the ferroelectric polarization P. To reduce the computational task, we employ the propagation vector k = (1/2, 1/3, 1/3) with a supercell, and specify the spins by using Eq. (2) and the experimental ξ i 10 to simulate the experimental helicoidal magnetic structure from first-principles calculations. Our GGA+U calculation gives P = (138, 294, -52) 10-4 μc/cm 2 in the xyz coordinates, and the GGA+U+SOC calculation gives almost an identical result with P = (138, 302, -48) 10-4 μc/cm 2. Both are substantially greater than the experimental value of μc/cm 2 in magnitude. The experiment 10 shows the angle between P and n expr. (<P, n expr. >) to be about 14, and our calculation shows that <P, n expr. > = 22.8 from the GGA+U calculations. The above results indicate that the polarization of Cu 3 Nb 2 O 8 basically originates from the exchange striction between Cu pairs rather than SOC because SOC does not play an important role on the polarization. This is not in support of the ferroaxial mechanism 10 because the direction of the polarization does not depend on the spin rotation plane of the magnetic structure. To see which Cu pairs are mainly responsible for the ferroelectric polarization, w e apply our general spin-induced ferroelectric polarization model to this Cu 3 Nb 2 O 8 case. Our general model is illustrated in Table II. Since the P comes from the exchange striction, we have (see Table II): P = P es ij (S i S j ) /V <i,j> (3) where Pes ij represents the exchange striction polarization coefficients of spins at i and j sites, and V is the volume of the cell containing the pairs in the summation. We again employ the mapping analysis 7,8 to derive the exchange striction polarization coefficients from first-principles GGA+U calculations, as listed in Table I. We note that those pairs with inversion symmetry have P es = 0, and those without inversion symmetry have comparatively large P es (Table I). We use our model [Eq. (3)] to sum up all the pair contributions to the total polarization P with the polarization coefficients P es derived from the GGA+U calculations. We adopt the experimental values k = (0.4876, , ), ξ Cu1 = 0, ξ Cu2_1 = 1.03π and ξ Cu2_2 = 1.05π to obtain the summed up result P = (30, 32, 9) 10-4 μc/cm 2 (see Table III), which has the magnitude of μc/cm 2 and <P, n expr. > = 8.6. The obtained P is closer (see below for further discussions on the magnitude of P) to the experimental value, assuring that the total P comes from the exchange striction. We now examine why the P calculated directly from first-principles calculations using k = (1/2, 1/3, 1/3) is so large. With the experimental ξ i values and k = (1/2, 1/3, 1/3), we find P = (119, 206, -22) 10-4 μc/cm 2 from our model (see Table III), which is consistent with our first-principles result. This indicates that the discrepancy in the polarization between the direct first-principles results and the experimental value is due to the improper spin arrangement and the improper k value adopted in the GGA+U calculation. It should be point out that the total polarization P is very sensitive to the propagation vector k and the relative phase ξ i. We define the commensurate magnetic structure, i.e., purely FM-u and 5
6 AFM-c through all crystal, with strictly k a* = k b* + k c* = 0.5, ξ Cu1 = 0, and ξ Cu2_1 = ξ Cu2_2 = π. However, the real magnetic structure from experiment 10 is slightly incommensurate, with k a* = , k b* + k c* = , ξ Cu1 = 0, ξ Cu2_1 = 1.03π and ξ Cu2_2 = 1.05π. Because the crystal structure of Cu 3 Nb 2 O 8 has inversion symmetry, when k is commensurate and ξ Cu2_1 = ξ Cu2_2 = π, the magnetic structure cannot break the inversion symmetry of the crystal structure, and thus would lead no ferroelectric polarization, i.e., P = 0 (see last row of Table III). Only when the k and ξ i are incommensurate (even if slightly), there would be the canting of spins, and the pair contribution to the P would not cancel out, resulting in non-zero macroscopic observable polarization. Specifically, the pairs with large P es, i.e., P 1 es, P 4 es, P 5 es and P 7 es, would dominate the P due to the spin canting. However, P 9 es is an exception, because the contributions of the two J 9 pairs in one unit cell strictly cancel out due to symmetry. The incommensurate behaviors of k and ξ i may come from the perturbation of the interchain and interlayer J s, the antisymmetric DM interactions and the single-ion anisotropy. 8,11,12 These effects and hence the associated spin canting are small, but are critical in determining the magnitude of P. The above explanation needs further validation. Results from our spin-induced ferroelectric polarization model, presented in Fig. 3 (a) and (b), show how the propagation vector k and the relative phases ξ i affect the total ferroelectric polarization P. The experimental values are k = (k a*, k b*, k c* ) = (0.4876, , ), ξ Cu1 = 0, ξ Cu2_1 = 1.03π and ξ Cu2_2 = 1.05π [10]. We first consider the effects of k. We fix two of the k i (i = a*, b*, c*) as the experimental values, and vary the remaining component k i from 0 to 1 to obtain the three curves in Fig. 3(a). Explicitly, P could take a value in a vast range, i.e., from almost 0 to more than μc/cm 2, owing to the largest polarization coefficients in Table I, as mentioned before. We note the three curves take the form of sine function, and this is consistent with Eq. (2). In Fig. 3(b), we keep k and ξ Cu1 as the experimental values and see the impact to P brought about by ξ i (i = Cu2_1, Cu2_2). Again we see the curves are in the form of sine function except for the vicinity of π, where there exist fluctuations that may come from the incommensurate property of the experimental k and the other ξ i. The experimentally measured polarization is μc/cm 2. If we take all experimental parameters for the spin structures, the predicted polarization is μc/cm 2 (see Table III), which is larger than the experimental value by a factor of 2.5. Several reasons may account for this discrepancy. First, the ferroelectric polarization depends sensitively on the spin structure. We note that ξ i (i = Cu2_1, Cu2_2) deviate from π very slightly, and the small deviation takes a relatively large error; as shown in Ref. 10, the ξ Cu2_1 and ξ Cu2_2 are 1.03(7) π and 1.05(7) π, respectively. Note that, as shown in the inset of Fig. 3(b), one can either reproduce the experimental P or even underestimate it within the error range. Second, it is very common (see for example Ref. 8) in multiferroics that the theoretical polarization is larger than the experimental value due to the small magnitude and the possible electric leakage. Third, it was found that ionic displacement may have some effects on the polarization, 21,22 which is however beyond the scope of the current work. The above discussions are based on our polarization model. Now we show how well the spin-induced polarization model agrees with the first-principles results. As summarized in Table IV, for our first-principles calculations for the total ferroelectric polarization P, we have selected several k values. Results indicate that the model give reliable results as do first principles calculations. As for the ξ i (i = Cu2_1, Cu2_2), Fig. 3(c) shows that our model generally also provides similar results as do GGA+U calculations, but the results of our model are smaller than those of first-principles calculations. This may due to some longer-range interactions that we do 6
7 not consider. Furthermore, our GGA+U and GGA+U+SOC results are close, strongly indicating the P originates from the exchange striction. Taken together, our polarization model can predict similar and reliable results as do first principles calculations, and allows one to study much more configurations of k and ξ i when direct first-principles calculations cannot be afforded. All the results shown in Fig. 3 indicate that the polarization P, which originates from the exchange striction, is very sensitive to both k and ξ i and can take an extremely large value, but that the experimental P is relatively rather small because the magnetic structure of Cu 3 Nb 2 O 8 is almost collinear and slightly incommensurate. C. Effects of Spin-Orbit Coupling Finally, we investigate in more detail the effects of SOC on the magnetic and ferroelectric properties of Cu 3 Nb 2 O 8 by performing GGA+U+SOC calculations. Let us first check the isotropy of J, and we take J 1 and J 4 as examples. It is recalled that J 1 = mev and J 4 = 2.13 mev from GGA+U calculations. When SOC is considered, 12 we have Jxx 1 = mev, J yy 1 = mev, Jzz 1 = mev and Jxx 4 = 2.15 mev, J yy 4 = 2.14 mev, Jzz 4 = 2.14 mev, showing good isotropy and consistency with our GGA+U result. Given SOC as perturbation, the antisymmetric DM interaction is a second-order perturbation term and exists only when inversion symmetry absent. 23 The DM energy terms is E DM = D ij (S i S j ) <i,j> (4) and D ij is a vector. We again perform energy mapping analysis, 11,12 and find D 1 = (-0.15, 0.19, 0.15) mev and D 4 = (0.09, -0.08, 0.13) mev in the xyz coordinates, with the ratio D/J about 0.10 and 0.08, respectively. The comparatively small DM interaction 8,12,24 would not change the basic magnetic structure, but may determine the direction of the spin-rotation plane as a consequence of gaining more energy from it. As to the effect of SOC on P, we follow the third column of Table II and calculate all P αβ 12 terms of P 12 int for P 1 (corresponding to J 1 ) as an example, using mapping analysis 7,8 with GGA+U+SOC calculations. We find P ( 196, 334, 150) ( 23, 40, 12) ( 16, 17, 7) 1 int = [ (23,41,11 ) ( 197, 336,150) ( 8, 17, 3) ] 10 5 ea (20,20, 7) (8, 17, 4) ( 195, 336,150) where the tensor corresponds to the spin represented in the xyz coordinates, and the diagonal terms are quite consistent with our DFT+U calculated P 1 es = (-205, -337, 126) 10 5 eå (see Table I), and the contribution of SOC (mainly off-diagonal terms) are indeed small comparing with the effect of exchange striction. These results ensure that SOC has no obvious impact on both the magnetic structure and the ferroelectric polarization, and confirms the validity of our general spin-induced ferroelectric polarization model. 7,8 (5) IV. SUMMARY 7
8 We have performed mapping analysis based on first-principles calculations to extract the exchange parameters and the polarization coefficients of Cu 3 Nb 2 O 8. The magnetic structure of Cu 3 Nb 2 O 8 originates from the competition isotropic exchange interactions and is basically described by the FM-u and AFM-c arrangements leading to a coplanar helicoidal spin spiral order. Our MC simulations lead to the magnetic structure similar to the one experimentally observed. The total polarization P of Cu 3 Nb 2 O 8 is induced by the exchange striction rather than by SOC. The magnitude of P is governed by a delicate spin canting arising from the anisotropy of the incommensurate magnetic structure, and the direction of P is not determined by the orientation of the spin rotation plane. ACKNOWLEDGEMENTS This work is supported by NSF of China, the Special Funds for Major State Basic Research and the Research Program of Shanghai Municipality (Pujiang, Eastern Scholar), Foundation for the Author of National Excellent Doctoral Dissertation (FANEDD) of P. R. China. APPENDIX: THE MAGNETIC GROUND STATE OF THE Cu-O CHAIN In our main text, we have qualitatively analyzed the origin of the FM-u and AFM-c magnetic structure, and now we discuss it in more detail. We have defined the effective J FM competing with the antiferromagnetic J 10, denoting J 10 = J AFM in this APPENDIX. Precisely, J FM = J 1 + J 4 S Cu1 S Cu2_1 + J 9 S Cu2_1 S Cu2_1, i.e., the value of J FM varies with the spin directions [see Fig. 2(b)]. When the magnetic structure consists of the FM-u and AFM-c arrangements, the J FM would become the strongest ferromagnetic coupling with J FM = J 1 - J 4 - J 9. Because the exchange interactions are isotropic, we can fix the spin of Cu1 along the z axis, and that of Cu2_1 is in the xz plane, and thus the spin direction of the Cu2_1 is decided by a variable θ 1. The spin direction of Cu2_2 is decided by θ 2 and φ 2. For any θ 1 and θ 2, to lower the energy of spin exchange between Cu2_1 and Cu2_2, we have φ 2 = 0 or π, i.e., all three spins are in the xz plane. Considering the periodicity, the energy of the Cu-O chain can be written as E = J FM S Cu1 S Cu2_1 + J FM S Cu1 S Cu2_2 + J AFM S Cu2_1 S Cu2_2 = J FM cosθ 1 + J FM cosθ 2 + J AFM cosθ 12 (A1) where θ 12 = min (θ 1 + θ 2, 2π θ 1 θ 2 ) [min(a, B) equals to the smaller one between the two numbers A and B]. We rewrite Eq. (A1) E = J FM cosθ 1 + J FM cosθ 2 + J AFM cos(θ 1 + θ 2 ) We make a = cosθ 1, b = cosθ 2, and thus = J FM cosθ 1 + J FM cosθ 2 + J AFM (cosθ 1 cosθ 2 sinθ 1 sinθ 2 ) E = J FM (a + b) + J AFM [ab (1 a 2 )(1 b 2 )] Note that J AFM > 0 and the basic inequality that xy x2 +y 2, where x and y are any real numbers. Thus, we have 8 2 (A2) (A3)
9 E J FM (a + b) + J AFM [ab 2 a2 b 2 (a + b)2 ] = J 2 FM (a + b) + J AFM [ 1] 2 where the equality is valid when a = b. If we take t = (a + b) as a variable, the right hand side of Eq. (A4) is a quadratic function J AFM t 2 + J 2 FM t J AFM. Note that 2 t 2, J AFM > 0, 2 J FM J AFM > 0, and the minimum occurs at t = min ( J FM J AFM, 2). Thus, when J FM 2J AFM < 1, the three spins form a non-collinear spiral state. When J FM 2J AFM 1, t = 2 and cos θ 1 = cos θ 2 = 1 (i.e., θ 1 = θ 2 = 0), the lowest energy results with the three spins in FM arrangement. When J FM = J 1 - J 4 - J 9 (i.e., the strongest ferromagnetic coupling), we have J FM = J 1 J 4 J 9 = 5.72 mev > 2 J AFM = 5.56, and at this time, we have the AFM-c arrangement. Note that any other spin configuration cannot lead to a lower energy, which confirms that the FM-u and AFM-c arrangements form the magnetic ground state of the Cu-O chain. (A4) * xggong@fudan.edu.cn * hxiang@fudan.edu.cn REFERENCES 1 S.-W. Cheong and M. Mostovoy, Nature Mater. 6, 13 (2007); R. Ramesh and N. Spaldin, Nature Mater. 6, 21 (2007); S. Picozzi and C. Ederer, J. Phys. Condens. Matter 21, (2009). 2 T. Kimura, Annu. Rev. Mater. Res. 37, 387 (2007). 3 Y. Tokura and S. Seki, Adv. Mater. 22, 1554 (2010). 4 M. Bibes, J. E. Villegas, and A. Barthélémy, Adv. Phys. 60, 5 (2011). 5 H. Katsura, N. Nagaosa, A. V. Balatsky, Phys. Rev. Lett. 95, (2005). 6 I. A. Sergienko and E. Dagotto, Phys. Rev. B 73, (2006). 7 H. J. Xiang, E. J. Kan, Y. Zhang, M.-H. Whangbo, and X. G. Gong, Phys. Rev. Lett. 107, (2011) 8 X. Z. Lu, M.-H. Whangbo, Shuai Dong, X. G. Gong, and H. J. Xiang, Phys. Rev. Lett. 108, (2012). 9 R. D. Johnson, L. C. Chapon, D. D. Khalyavin, P. Manuel, P. G. Radaelli, and C. Martin, Phys. Rev. Lett. 108, (2012). 10 R. D. Johnson, Sunil Nair, L. C. Chapon, A. Bombardi, C. Vecchini, D. Prabhakaran, A. T. Boothroyd, and P. G. Radaelli, Phys. Rev. Lett. 107, (2011). 11 H. J. Xiang, E. J. Kan, Su-Huai Wei, M.-H. Whangbo, and X. G. Gong, Phys. Rev. B 84, (2011). 12 Z. L. Li, J. H. Yang, G. H. Chen, M.-H. Whangbo, H. J. Xiang, and X. G. Gong, Phys. Rev. B 85, (2012). 9
10 13 G. Kresse and J. Furthmüller, Phys. Rev. B 54, (1996). 14 P. E. Blöchl, Phys. Rev. B 50, (1994). 15 G. Kresse and D. Joubert, Phys. Rev. B 59, 1758 (1999). 16 J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996). 17 A. I. Liechtenstein, V. I. Anisimov, and J. Zaanen, Phys. Rev. B 52, R5467 (1995). 18 H. J. Xiang and M.-H. Whangbo, Phys. Rev. Lett. 99, (2007). 19 H. J. Xiang, C. Lee, and M.-H. Whangbo, Phys. Rev. B 76, (R) (2007). 20 R. D. King-Smith and D. Vanderbilt, Phys. Rev. B 47, 1651 (1993); R. Resta, Rev. Mod. Phys. 66, 899 (1994). 21 H. J. Xiang, Su-Huai Wei, M.-H. Whangbo, and Juarez L. F. Da Silva, Phys. Rev. Lett. 101, (2008). 22 S. Picozzi, K. Yamauchi, B. Sanyal, I. A. Sergienko, and E. Dagotto, Phys. Rev. Lett. 99, (2007). 23 D. Dai, H. J. Xiang, M.-H. Whangbo, J. Comput. Chem. 29, 13 (2008). 24 H. J. Xiang, E. J. Kan, M.-H. Whangbo, C. Lee, Su-Huai Wei, and X. G. Gong, Phys. Rev. B 83, (2011). 10
11 Table I. The exchange interactions and the exchange striction polarization coefficients derived from first-principles calculations with the mapping methods. Path Pair Type a Pair distance ( Å) J (mev) P es [10-5 (e Å)] 1 Cu1-Cu2 SE (-205, -337, 126) 2 Cu2-Cu2 SE-I (0, 0, 0) 3 Cu2-Cu2 SE-I (0, 0, 0) 4 Cu1-Cu2 SE (-174, -243, -96) 5 Cu1-Cu2 SSE (9, 4, -10) 6 Cu2-Cu2 SSE-I (0, 0, 0) 7 Cu1-Cu2 SSE (28, 22, 23) 8 Cu1-Cu1 SSE-I (1, 1, 0) 9 Cu2-Cu2 SSE (-69, -92, -26) 10 Cu2-Cu2 SSE-I (0, 0, 0) a SE represents superexchange, SSE super-superexchange. The exchange pair with inversion symmetry is indicated by adding I. 11
12 Table II. Explaining the ferroelectric polarization induced by a helical magnetic structure using our spin-induced ferroelectric polarization model. In general, the polarization P has two origins, i.e., the exchange striction and SOC. Under different situations, our model is reduced to simpler forms. P 12 (S 1,S 2 ) = P αβ 12 S 1α S 2β = (S 1x,S 1y,S 1z ) ( αβ Exchange striction 8 Any rotation is allowed, so the Pxx xy 12 P 12 yx P 12 yy P 12 zy Pzx 12 P 12 With inversion symmetry 7 Pxz 12 yz P 12 Pzz 12 Spin-orbit coupling (SOC) S 2x )( S 2y ) = S 1 P 12 int S 2 S 2z With no inversion symmetry diagonal terms of P 12 int are equal and the off-diagonal ones are zero: P 12 = P es 12 (S 1 S 2 ) P es 12 = 0 if the spin dimer has inversion symmetry. P 12 int reduces to an antisymmetric tensor. Generalized spin current model: P 12 = M(S 1 S 2 ) where M is a matrix. All P αβ 12 terms of P 12 int are needed from first-principles. CaMn 7 O 12, Cu 3 Nb 2 O 8 MnI 2 12
13 Table III. The total polarization derived from our spin-induced ferroelectric polarization model for the case of the exchange striction, with different propagation vector k and relative phases. ξ Cu1 = 0 is fixed in this table. Spin orientation S i = R cos(k R L + ξ i ) + I sin(k R L + ξ i ) Polarization P = P es ij (S i S j ) /V <i,j> k = (k a*, k b*, k c* ) ξ Cu2_1 / π ξ Cu2_2 / π P ( 10-4 μc/cm 2 ) (0.4876, , ) (30, 32, 9) (0.4876, , ) (10, 10, 4) (1/2, 1/3, 1/3) (119, 206, -22) (1/2, 1/3, 1/3) (165, 283, -30) (1/2, 1/4, 1/4) (13, 12, 3) (1/2, 1/4, 1/4) (0, 0, 0) 13
14 Table IV. Values of P obtained from the first principles and model calculations for various k with ξ Cu1 = 0, ξ Cu2_1 = 1.03π and ξ Cu2_2 = 1.05π. P varies k drastically, showing that the model can predict results consistent with first-principles calculations. Propagation vector k Polarization P ( 10-4 μc/cm 2 ) GGA+U GGA+U+SOC Model (1, 1/2, 1/2) (5, 8, -5) (5, 7, -3) (-1, -1, 2) (1/2, 1, 1/2) (2, 6, -2) (3, 4, -4) (3, 5, -1) (1/2, 1/2, 1) (3, 5, -1) (5, 6, -4) (3, 5, 0) (1/2, 1/2, 1/2) (790, 1141, -125) (806, 1151, -120) (732, 1179, -91) (1/2, 1/2, 1/3) (567, 898, -146) (563, 904, -106) (498, 804, -63) (1/2, 1/3, 1/2) (523, 953, -173) (532, 989, -134) (499, 805, -62) (1/3, 1/2, 1/2) (538, 727, -39) (535, 754, -22) (511, 823, -69) (1/2, 1/3, 1/3) (138, 294, -52) (138, 302, -48) (119, 206, -22) 14
15 (a) (b) FIG. 1 (Color online). Crystal structure of Cu 3 Nb 2 O 8. (a) The Cu 3 O 8 units connect to form the Cu-O chain. The Cu1 atoms that are at the inversion center are in brown, and the Cu2 atoms are in blue. The magnetic structure is approximately given by the FM-u and AFM-c arrangements. (b) Cu-O chains stacked to form Cu-O layers, separated by non-magnetic Nb atoms. The Nb-O bonds are not shown. 15
16 (a) (b) FIG. 2 (Color online). (a) The ten spin exchange paths. The dashed line represents one Cu 3 O 8 unit. The empty arrow points to the ith atom and the red solid arrow the jth atom. The inversion centers of those pairs with inversion symmetry are marked as green cross. (b) The main intrachain interactions leading to the FM-u and AFM-c magnetic structure, which is determined not only by the intra-trimer FM J 1 and AFM J 10 but also by the inter-trimer AFM J 4 and J 9. 16
17 (a) (b) (c) 17
18 FIG. 3 (Color online). (a) The effect of the propagation vector k on the magnitude of P predicted from our model using the ferroelectric polarization coefficients derived from the first-principles calculations. We take the experimental ξ i and fix two components of k, to see how the third component affects P. Results indicate that P varies in a wide range as a function of k. The experimental magnetic structure is almost collinear so that k satisfies k a* k b* + k c* 0.5, leading to a small P. (b) The effect of the relative phase ξ i to P, indicated by the polarization model. Results show that P is very sensitive to a small change in ξ i, given the experimental value of μc/cm 2 (see inset). (c) The effect of ξ Cu2_2 on P obtained from the first-principles and the model calculations for a simple system with k = (1, 1, 1), for which both the GGA+U and the GGA+U+SOC calculations can be performed without much computation task. The DFT and the model calculations produce similar results, indicating the correctness of the model. See Table IV for the effect of k on P. 18
Realizing Magnetoelectric Coupling with Hydroxide as a Knob J.Y. Ni 1,2, P.S. Wang 1,2, J. L. Lu 1,2, and H. J. Xiang 1,2*
Realizing Magnetoelectric Coupling with Hydroxide as a Knob J.Y. Ni 1,2, P.S. Wang 1,2, J. L. Lu 1,2, and H. J. Xiang 1,2* 1 Key Laboratory of Computational Physical Sciences (Ministry of Education), State
More informationDefects in TiO 2 Crystals
, March 13-15, 2013, Hong Kong Defects in TiO 2 Crystals Richard Rivera, Arvids Stashans 1 Abstract-TiO 2 crystals, anatase and rutile, have been studied using Density Functional Theory (DFT) and the Generalized
More informationSupporting information. Realizing Two-Dimensional Magnetic Semiconductors with. Enhanced Curie Temperature by Antiaromatic Ring Based
Supporting information Realizing Two-Dimensional Magnetic Semiconductors with Enhanced Curie Temperature by Antiaromatic Ring Based Organometallic Frameworks Xingxing Li and Jinlong Yang* Department of
More informationSnO 2 Physical and Chemical Properties due to the Impurity Doping
, March 13-15, 2013, Hong Kong SnO 2 Physical and Chemical Properties due to the Impurity Doping Richard Rivera, Freddy Marcillo, Washington Chamba, Patricio Puchaicela, Arvids Stashans Abstract First-principles
More informationHigh-pressure cupric oxide: a room-temperature
1 High-pressure cupric oxide: a room-temperature multiferroic Xavier Rocquefelte 1*, Karlheinz Schwarz 2, Peter Blaha 2 Sanjeev Kumar 3, Jeroen van den Brink 4 1. Institut des Matériaux Jean Rouxel, UMR
More informationSpin density wave as a superposition of two magnetic states of opposite chirality and its implications
Spin density wave as a superposition of two magnetic states of opposite chirality and its implications Elijah E. Gordon a, Shahab Derakhshan b, Corey M. Thompson c, and Myung-Hwan Whangbo a,d,e, * a Department
More informationBerry s phase under the Dzyaloshinskii-Moriya interaction
PHYSICAL REVIEW A 77, 8 erry s phase under the Dzyaloshinskii-Moriya interaction M. K. Kwan, Zeynep ilhan Gurkan, and L. C. Kwek, ational Institute of Education, anyang Technological University, anyang
More informationChromium Cluster on Defected Graphene
Chromium Cluster on Defected Graphene Yuhang Liu June 29, 2017 Abstract In this work, diffusion process of Cr atoms on two types of defected graphene and structure and magnetic properties of Cr cluster
More informationDirect visualization of the Jahn Teller effect coupled to Na ordering in Na 5/8 MnO 2
Direct visualization of the Jahn Teller effect coupled to Na ordering in Na 5/8 MnO 2 Xin Li 1, Xiaohua Ma 1, Dong Su 2, Lei Liu 1, Robin Chisnell 3, Shyue Ping Ong 1, Hailong Chen 1, Alexandra Toumar
More informationExplaining the apparent arbitrariness of the LDA-1/2 self-energy. correction method applied to purely covalent systems
Explaining the apparent arbitrariness of the LDA-1/2 self-energy correction method applied to purely covalent systems Kan-Hao Xue, 1,2 Leonardo R. C. Fonseca, 3 and Xiang-Shui Miao 1,2 1 School of Optical
More informationSupplementary Information. Structural Transition and Unusually Strong Antiferromagnetic Superexchange Coupling in Perovskite KAgF3
Supplementary Information Structural Transition and Unusually Strong Antiferromagnetic Superexchange Coupling in Perovskite KAgF3 Dominik Kurzydłowski, *a Zoran Mazej, b Zvonko Jagličić, c Yaroslav Filinchuk
More informationQuantum Chemical Study of Defective Chromium Oxide
, March 13-15, 2013, Hong Kong Quantum Chemical Study of Defective Chromium Oxide Richard Rivera, Soraya Jácome, Frank Maldonado, Arvids Stashans 1 Abstract Through the use of first-principles calculations
More information100 Tesla multishot. 60 Tesla long pulse. Los Alamos branch of the Magnet Lab Pulsed magnetic fields
Los Alamos branch of the Magnet Lab Pulsed magnetic fields 100 Tesla multishot 100 80 60 40 20 Magnetic field (T) 0 0 0.5 1 1.5 2 2.5 3 time (s) 60 Tesla long pulse 60 40 20 0 0 1 2 3 time (s) Magnetization,
More informationAntiferromagnetic Textures
Antiferromagnetic Textures This image cannot currently be displayed. ULRICH K. RÖSSLE IFW DRESDEN SPICE Workshop Antiferromagnetic Spintronics 26.-30/09/2016 Schloss Waldthausen u.roessler@ifw-dresden.de
More informationarxiv: v2 [cond-mat.mtrl-sci] 7 Oct 2010
Orbital degrees of freedom as origin of magnetoelectric coupling in magnetite Kunihiko Yamauchi and Silvia Picozzi Consiglio Nazionale delle Ricerche (CNR-SPIN), 671 L Aquila, Italy (Dated: July 14, 218)
More informationSupporting Information
Supporting Information The Origin of Active Oxygen in a Ternary CuO x /Co 3 O 4 -CeO Catalyst for CO Oxidation Zhigang Liu, *, Zili Wu, *, Xihong Peng, ++ Andrew Binder, Songhai Chai, Sheng Dai *,, School
More informationCrystalline and Magnetic Anisotropy of the 3d Transition-Metal Monoxides
Crystalline and of the 3d Transition-Metal Monoxides Institut für Festkörpertheorie und -optik Friedrich-Schiller-Universität Max-Wien-Platz 1 07743 Jena 2012-03-23 Introduction Crystalline Anisotropy
More informationAnalysis of Spin Exchange Interactions in (C 2 N 2 H 10 )[Fe(HPO 3 )F 3 ] on the Basis of Electronic Structure Calculations
Analysis of Spin Exchange Interactions in (C N H 10 )[Fe(HPO 3 )F 3 ] Bull. Korean Chem. Soc. 011, Vol. 3, No. 467 DOI 10.501/bkcs.011.3..467 Analysis of Spin Exchange Interactions in (C N H 10 )[Fe(HPO
More informationStructural and Electronic Effects on the Properties of Fe 2 (dobdc) upon Oxidation with N 2 O
Supporting information for paper in Inorganic Chemistry, April 11, 016, page S-1 Structural and Electronic Effects on the Properties of Fe (dobdc) upon Oxidation with N O oshua Borycz, 1, oachim Paier,
More informationMgO-decorated carbon nanotubes for CO 2 adsorption: first principles calculations
MgO-decorated carbon nanotubes for CO 2 adsorption: first principles calculations Zhu Feng( ), Dong Shan( ), and Cheng Gang( ) State Key Laboratory for Superlattices and Microstructures, Institute of Semiconductors,
More informationarxiv: v2 [cond-mat.mtrl-sci] 8 Jun 2018
Relationship between crystal structure and multiferroic orders in orthorhombic perovskite manganites arxiv:1805.02172v2 [cond-mat.mtrl-sci] 8 Jun 2018 Natalya S. Fedorova, 1, Yoav William Windsor, 2 Christoph
More informationMagnetoelectricity and multiferroics. Charles Simon Laboratoire CRISMAT, CNRS and ENSICAEN, F14050 Caen.
Magnetoelectricity and multiferroics Charles Simon Laboratoire CRISMAT, CNRS and ENSICAEN, F14050 Caen. Introduction : The possibility for a material to be both ferromagnetic and ferroelectric was predicted
More informationPolarization Enhancement in Perovskite Superlattices by Oxygen. Octahedral Tilts
Polarization Enhancement in Perovskite Superlattices by Oxygen Octahedral Tilts X. Z. Lu a), X. G. Gong, and H. J. Xiang b) Key Laboratory of Computational Physical Sciences (Ministry of Education), State
More informationSupplementary information
Supplementary Figures Supplementary information Supplementary Figure 1: Microscopic magnetic model and the effective model. a, microscopic magnetic model is based on five leading couplings: JW FM, J S,
More informationSupporting Information
Supporting Information Controlled Growth of Ceria Nanoarrays on Anatase Titania Powder: A Bottom-up Physical Picture Hyun You Kim 1, Mark S. Hybertsen 2*, and Ping Liu 2* 1 Department of Materials Science
More informationarxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 13 Nov 2003
1. 14 August 1996 (final accepted version arxiv:cond-mat/0311297v1 [cond-mat.mtrl-sci] 13 Nov 2003 2. Non-collinear magnetism in distorted perovskite compounds 3. I.V.Solovyev a,, N.Hamada b, K.Terakura
More informationExperiment Section Fig. S1 Fig. S2
Electronic Supplementary Material (ESI) for ChemComm. This journal is The Royal Society of Chemistry 2018 Supplementary Materials Experiment Section The STM experiments were carried out in an ultrahigh
More informationMultiferroic BiFeO 3. Sang-Wook Cheong. Seongsu Lee. Partially supported by NSF-MRSEC
Multiferroic BiFeO 3 Seongsu Lee Chenglin Zhang YoungJai Choi Seongsu Lee Soonyong Park Y. Horibe Valery Kiryukhin Y. J. Cho L. Balicas Florida SU. X. S. Xu J. Musfeldt U. of Tennessee J. G. Park SKKU
More informationFacet engineered Ag 3 PO 4 for efficient water photooxidation
Supporting Information Facet engineered Ag 3 PO 4 for efficient water photooxidation David James Martin, Naoto Umezawa, Xiaowei Chen, Jinhua Ye and Junwang Tang* This file includes the following experimental/theoretical
More informationExchange-induced negative-u charge order in N-doped WO 3 : A spin-peierls-like system
Exchange-induced negative-u charge order in N-doped WO 3 : A spin-peierls-like system Muhammad N. Huda,*, Yanfa Yan, Su-Huai Wei, and Mowafak M. Al-Jassim National Renewable Energy Laboratory, Golden,
More informationTunable Band Gap of Silicene on Monolayer Gallium Phosphide Substrate
2017 International Conference on Energy Development and Environmental Protection (EDEP 2017) ISBN: 978-1-60595-482-0 Tunable Band Gap of Silicene on Monolayer Gallium Phosphide Substrate Miao-Juan REN
More informationarxiv: v1 [cond-mat.mtrl-sci] 29 Sep 2010
Spin-spiral inhomogeneity as origin of ferroelectricity in orthorhombic manganites I. V. Solovyev National Institute for Materials Science, -- Sengen, Tsukuba, Ibaraki 0-0047, Japan (Dated: December, 08)
More informationElectron Spin Resonance and Quantum Dynamics. Masaki Oshikawa (ISSP, University of Tokyo)
Electron Spin Resonance and Quantum Dynamics Masaki Oshikawa (ISSP, University of Tokyo) Electron Spin Resonance (ESR) E-M wave electron spins H measure the absorption intensity Characteristic of ESR single
More informationSupporting Information Tuning Local Electronic Structure of Single Layer MoS2 through Defect Engineering
Supporting Information Tuning Local Electronic Structure of Single Layer MoS2 through Defect Engineering Yan Chen, 1,2,,$, * Shengxi Huang, 3,6, Xiang Ji, 2 Kiran Adepalli, 2 Kedi Yin, 8 Xi Ling, 3,9 Xinwei
More informationLaboratory of Surface Physics, Collaborative Innovation Center of Advanced Microstructures,
Predicting a Ferrimagnetic Phase of Zn2FeOsO6 with Strong Magnetoelectric Coupling P. S. Wang, 1 W. Ren, 2 L. Bellaiche, 3 and H. J. Xiang* 1 1 Key Laboratory of Computational Physical Sciences (Ministry
More informationSupporting Information for 2D Intrinsic. Ferromagnets from van der Waals. Antiferromagnets
Supporting Information for 2D Intrinsic Ferromagnets from van der Waals Antiferromagnets Naihua Miao,, Bin Xu, Linggang Zhu,, Jian Zhou, and Zhimei Sun,, School of Materials Science and Engineering, Beihang
More informationCompeting Ferroic Orders The magnetoelectric effect
Competing Ferroic Orders The magnetoelectric effect Cornell University I would found an institution where any person can find instruction in any study. Ezra Cornell, 1868 Craig J. Fennie School of Applied
More informationSupporting Information for. Electric-magneto-optical Kerr effect in a. hybrid organic-inorganic perovskite
Supporting Information for Electric-magneto-optical Kerr effect in a hybrid organic-inorganic perovskite Feng-Ren Fan, Hua Wu,,, Dmitrii Nabok, Shunbo Hu, Wei Ren, Claudia Draxl, and Alessandro Stroppa,,
More informationGiant biquadratic interaction-induced magnetic anisotropy in the iron-based superconductor A x Fe 2 y Se 2
PHYSICAL REVIEW B 93, 024511 (2016) Giant biquadratic interaction-induced magnetic anisotropy in the iron-based superconductor A x Fe 2 y Se 2 Hai-Feng Zhu, 1,2 Hai-Yuan Cao, 1,2 Yun Xie, 1,2 Yu-Sheng
More informationSupplementary Figures
Supplementary Figures 8 6 Energy (ev 4 2 2 4 Γ M K Γ Supplementary Figure : Energy bands of antimonene along a high-symmetry path in the Brillouin zone, including spin-orbit coupling effects. Empty circles
More informationSkyrmions and Anomalous Hall Effect in a Dzyaloshinskii-Moriya Magnet
Skyrmions and Anomalous Hall Effect in a Dzyaloshinskii-Moriya Magnet Jung Hoon Han (SungKyunKwanU, Suwon) Su Do Yi SKKU Shigeki Onoda RIKEN Naoto Nagaosa U of Tokyo arxiv:0903.3272v1 Nearly ferromagnetic
More informationFirst Principles Calculation of Defect and Magnetic Structures in FeCo
Materials Transactions, Vol. 47, No. 11 (26) pp. 2646 to 26 Special Issue on Advances in Computational Materials Science and Engineering IV #26 The Japan Institute of Metals First Principles Calculation
More informationA New Electronic Orbital Order Identified in Parent Compound of Fe-Based High-Temperature Superconductors
A New Electronic Orbital Order Identified in Parent Compound of Fe-Based High-Temperature Superconductors Cooperative Research Team on Predictive Capability for Strongly Correlated Systems Summary: The
More informationCurvature-enhanced Spin-orbit Coupling and Spinterface Effect in Fullerene-based Spin Valves
Supplementary Information Curvature-enhanced Spin-orbit Coupling and Spinterface Effect in Fullerene-based Spin Valves Shiheng Liang 1, Rugang Geng 1, Baishun Yang 2, Wenbo Zhao 3, Ram Chandra Subedi 1,
More informationLuigi Paolasini
Luigi Paolasini paolasini@esrf.fr LECTURE 4: MAGNETIC INTERACTIONS - Dipole vs exchange magnetic interactions. - Direct and indirect exchange interactions. - Anisotropic exchange interactions. - Interplay
More informationSUPPLEMENTARY INFORMATION
Supplementary Information Anisotropic conductance at improper ferroelectric domain walls D. Meier 1,, *, J. Seidel 1,3, *, A. Cano 4, K. Delaney 5, Y. Kumagai 6, M. Mostovoy 7, N. A. Spaldin 6, R. Ramesh
More informationCover Page. The handle holds various files of this Leiden University dissertation.
Cover Page The handle http://hdl.handle.net/1887/49403 holds various files of this Leiden University dissertation. Author: Keesman, R. Title: Topological phases and phase transitions in magnets and ice
More informationNovel Magnetic Properties of Carbon Nanotubes. Abstract
Novel Magnetic Properties of Carbon Nanotubes Jian Ping Lu Department of Physics and Astronomy, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599 jpl@physics.unc.edu arxiv:cond-mat/94779v1
More informationshows the difference between observed (black) and calculated patterns (red). Vertical ticks indicate
Intensity (arb. unit) a 5 K No disorder Mn-Pt disorder 5 K Mn-Ga disorder 5 K b 5 K Observed Calculated Difference Bragg positions 24 28 32 2 4 6 8 2 4 2θ (degree) 2θ (degree) Supplementary Figure. Powder
More informationSTRUCTURAL AND MECHANICAL PROPERTIES OF AMORPHOUS SILICON: AB-INITIO AND CLASSICAL MOLECULAR DYNAMICS STUDY
STRUCTURAL AND MECHANICAL PROPERTIES OF AMORPHOUS SILICON: AB-INITIO AND CLASSICAL MOLECULAR DYNAMICS STUDY S. Hara, T. Kumagai, S. Izumi and S. Sakai Department of mechanical engineering, University of
More informationarxiv: v1 [cond-mat.mtrl-sci] 21 May 2012
Magnetic structure of hexagonal YMnO 3 and LuMnO 3 from a microscopic point of view I. V. Solovyev arxiv:1205.4478v1 [cond-mat.mtrl-sci] 21 May 2012 Computational Materials Science Unit, National Institute
More informationSpin Superfluidity and Graphene in a Strong Magnetic Field
Spin Superfluidity and Graphene in a Strong Magnetic Field by B. I. Halperin Nano-QT 2016 Kyiv October 11, 2016 Based on work with So Takei (CUNY), Yaroslav Tserkovnyak (UCLA), and Amir Yacoby (Harvard)
More informationarxiv: v3 [cond-mat.mtrl-sci] 4 Mar 2017
Transition between strong and weak topological insulator in ZrTe 5 and HfTe 5 Zongjian Fan 1, Qi-Feng Liang 2, Y. B. Chen 3, Shu-Hua Yao 1, and Jian Zhou 1,4,* arxiv:1611.04263v3 [cond-mat.mtrl-sci] 4
More informationAnisotropic phase diagram of the frustrated spin dimer compound Ba 3 Mn 2 O 8
Anisotropic phase diagram of the frustrated spin dimer compound Ba 3 Mn 2 O 8 E. C. Samulon, 1 K. A. Al-Hassanieh, 2 Y.-J. Jo, 3,4 M. C. Shapiro, 1 L. Balicas, 3 C. D. Batista, 2 and I. R. Fisher 1 1 Geballe
More informationSupporting Online Material for
www.sciencemag.org/cgi/content/full/325/5948/1670/dc1 Supporting Online Material for Coordinatively Unsaturated Al 3+ Centers as Binding Sites for Active Catalyst Phases of Platinum on γ-al 2 O 3 Ja Hun
More informationPhysics. Physics Research Publications. Purdue University Year 2009
Physics Physics Research Publications Purdue University Year 2009 First-Principles Calculations of the Electronic Structure of Tetragonal alpha-fete and alpha-fese Crystals: Evidence for a Bicollinear
More informationarxiv: v2 [cond-mat.str-el] 17 Jun 2009
Structural distortion and frustrated magnetic interactions in the layered copper oxychloride (CuCl)LaNb 2 O 7 arxiv:0901.0154v2 [cond-mat.str-el] 17 Jun 2009 Alexander A. Tsirlin 1, 2, and Helge Rosner
More informationAntiferromagnetic Spin Coupling between Rare Earth Adatoms and Iron Islands Probed by Spin-Polarized Tunneling
Supplementary information Antiferromagnetic Spin Coupling between Rare Earth Adatoms and Iron Islands Probed by Spin-Polarized Tunneling David Coffey, *,+,,# José Luis Diez-Ferrer, + David Serrate, +,#
More informationOutline. Introduction: graphene. Adsorption on graphene: - Chemisorption - Physisorption. Summary
Outline Introduction: graphene Adsorption on graphene: - Chemisorption - Physisorption Summary 1 Electronic band structure: Electronic properties K Γ M v F = 10 6 ms -1 = c/300 massless Dirac particles!
More informationEFFECTIVE MAGNETIC HAMILTONIANS: ab initio determination
ICSM212, Istanbul, May 3, 212, Theoretical Magnetism I, 17:2 p. 1 EFFECTIVE MAGNETIC HAMILTONIANS: ab initio determination Václav Drchal Institute of Physics ASCR, Praha, Czech Republic in collaboration
More informationMultiferroicity due to Charge Ordering
Multiferroicity due to Charge Ordering Jeroen van den Brink ESRF 8/2/2011 Multiferroicity due to Charge Ordering Jeroen van den Brink Gianluca Giovannetti & Sanjeev Kumar Sylvia Picozzi & Daniel Khomskii
More informationarxiv: v1 [cond-mat.str-el] 17 Jun 2013
MnSb 2 O 6 : A polar magnet with a chiral crystal structure arxiv:1306.3854v1 [cond-mat.str-el] 17 Jun 2013 R. D. Johnson, 1,2, K. Cao, 3 L. C. Chapon, 4 F. Fabrizi, 1 N. Perks, 1 P. Manuel, 2 J. J. Yang,
More informationT C SUH. Tuning multiferroics under extreme conditions: Effects of high pressure, magnetic fields, and substitutions. Bernd Lorenz
Tuning multiferroics under extreme conditions: Effects of high pressure, magnetic fields, and substitutions Bernd Lorenz Texas Center for Superconductivity and Department of Physics, University of Houston
More informationTheoretical Concepts of Spin-Orbit Splitting
Chapter 9 Theoretical Concepts of Spin-Orbit Splitting 9.1 Free-electron model In order to understand the basic origin of spin-orbit coupling at the surface of a crystal, it is a natural starting point
More informationThe Low Temperature Conversion of Methane to Methanol on CeO x /Cu 2 O catalysts: Water Controlled Activation of the C H Bond
The Low Temperature Conversion of Methane to Methanol on CeO x /Cu 2 O catalysts: Water Controlled Activation of the C H Bond Zhijun Zuo, a Pedro J. Ramírez, b Sanjaya Senanayake, a Ping Liu c,* and José
More informationAnisotropic Magnetic Structures in Iron-Based Superconductors
Anisotropic Magnetic Structures in Iron-Based Superconductors Chi-Cheng Lee, Weiguo Yin & Wei Ku CM-Theory, CMPMSD, Brookhaven National Lab Department of Physics, SUNY Stony Brook Another example of SC
More informationUniversity of Chinese Academy of Sciences, Beijing , People s Republic of China,
SiC 2 Siligraphene and Nanotubes: Novel Donor Materials in Excitonic Solar Cell Liu-Jiang Zhou,, Yong-Fan Zhang, Li-Ming Wu *, State Key Laboratory of Structural Chemistry, Fujian Institute of Research
More informationH alf metallic system is of particular interest in spintroincs because one spin channel is completely suppressed
OPEN SUBJECT AREAS: MAGNETIC PROPERTIES AND MATERIALS TWO-DIMENSIONAL MATERIALS Received 29 April 2014 Accepted 22 July 2014 Published 14 August 2014 Correspondence and requests for materials should be
More informationTowards Direct-Gap Silicon Phases by the Inverse Band Structure Design Approach. P. R. China. Jiangsu , People s Republic of China
Towards Direct-Gap Silicon Phases by the Inverse Band Structure Design Approach H. J. Xiang 1,2*, Bing Huang 2, Erjun Kan 3, Su-Huai Wei 2, X. G. Gong 1 1 Key Laboratory of Computational Physical Sciences
More informationSpin-wave dispersion in half-doped La3/2Sr1/2NiO4
Physics Physics Research Publications Purdue University Year 2007 Spin-wave dispersion in half-doped La3/2Sr1/2NiO4 D. X. Yao E. W. Carlson This paper is posted at Purdue e-pubs. http://docs.lib.purdue.edu/physics
More informationSelectivity in the initial C-H bond cleavage of n-butane on PdO(101)
Supporting Information for Selectivity in the initial C-H bond cleavage of n-butane on PdO(101) Can Hakanoglu (a), Feng Zhang (a), Abbin Antony (a), Aravind Asthagiri (b) and Jason F. Weaver (a) * (a)
More informationThe Nature of the Interlayer Interaction in Bulk. and Few-Layer Phosphorus
Supporting Information for: The Nature of the Interlayer Interaction in Bulk and Few-Layer Phosphorus L. Shulenburger, A.D. Baczewski, Z. Zhu, J. Guan, and D. Tománek, Sandia National Laboratories, Albuquerque,
More information[100] directed Cu-doped h-coo Nanorods: Elucidation of. Growth Mechanism and Application to Lithium-Ion Batteries
Supplementary Information [100] directed Cu-doped h-coo Nanorods: Elucidation of Growth Mechanism and Application to Lithium-Ion Batteries Ki Min Nam, Young Cheol Choi, Sung Chul Jung, Yong-Il Kim, Mi
More informationUnderstanding the Magnetic Ground States for Improper Multiferroic Materials
Wright State University CORE Scholar Physics Seminars Physics 4-11-2013 Understanding the Magnetic Ground States for Improper Multiferroic Materials Jason T. Haraldsen Follow this and additional works
More informationSupplementary Materials for Oxygen-induced self-assembly of quaterphenyl molecule on metal surfaces
Electronic Supplementary Material (ESI) for ChemComm. This journal is The Royal Society of Chemistry 2014 Supplementary Materials for Oxygen-induced self-assembly of quaterphenyl molecule on metal surfaces
More informationSpinons and triplons in spatially anisotropic triangular antiferromagnet
Spinons and triplons in spatially anisotropic triangular antiferromagnet Oleg Starykh, University of Utah Leon Balents, UC Santa Barbara Masanori Kohno, NIMS, Tsukuba PRL 98, 077205 (2007); Nature Physics
More informationStability of skyrmion lattices and symmetries of Dzyaloshinskii-Moriya magnets. Alexey A. Kovalev Utkan Güngördü Rabindra Nepal
Stability of skyrmion lattices and symmetries of Dzyaloshinskii-Moriya magnets Alexey A. Kovalev Utkan Güngördü Rabindra Nepal Outline Discuss possible 2D Dzyaloshinskii-Moriya magnets Study phase diagram
More informationOur first-principles calculations were performed using the Vienna Ab Initio Simulation
Supplementary Note 1: Computational details First-principles calculations Our first-principles calculations were performed using the Vienna Ab Initio Simulation Package (VASP) 1, which is based on density
More informationEnergy Stabilities, Magnetic Properties, and Electronic Structures of Diluted Magnetic Semiconductor Zn 1 x Mn x S(001) Thin Films
CHINESE JOURNAL OF CHEMICAL PHYSICS VOLUME 24, NUMBER 1 FEBRUARY 27, 2011 ARTICLE Energy Stabilities, Magnetic Properties, and Electronic Structures of Diluted Magnetic Semiconductor Zn 1 x Mn x S(001)
More informationMassively parallel Monte Carlo simulation of a possible topological phase transition in two-dimensional frustrated spin systems
Massively parallel Monte Carlo simulation of a possible topological phase transition in two-dimensional frustrated spin systems Tsuyoshi OKUBO Institute for Solid State Physics, University of Tokyo Kashiwa-no-ha,
More informationSupporting Information. Don-Hyung Ha, Liane M. Moreau, Clive R. Bealing, Haitao Zhang, Richard G. Hennig, and. Richard D.
Supporting Information The structural evolution and diffusion during the chemical transformation from cobalt to cobalt phosphide nanoparticles Don-Hyung Ha, Liane M. Moreau, Clive R. Bealing, Haitao Zhang,
More informationSupplementary Figure 1 Representative sample of DW spin textures in a
Supplementary Figure 1 Representative sample of DW spin textures in a Fe/Ni/W(110) film. (a) to (d) Compound SPLEEM images of the Fe/Ni/W(110) sample. As in Fig. 2 in the main text, Fe thickness is 1.5
More informationObservation of ferroelectric domains in bismuth ferrite using coherent diffraction techniques
Observation of ferroelectric domains in bismuth ferrite using coherent diffraction techniques James Vale October 25, 2011 Abstract Multiferroic materials have significant potential for both the scientific
More informationMultiwalled nanotube faceting unravelled
Multiwalled nanotube faceting unravelled Itai Leven, Roberto Guerra, Andrea Vanossi, Erio Tosatti, and Oded Hod NATURE NANOTECHNOLOGY www.nature.com/naturenanotechnology 1 1. Achiral DWCNTs optimized using
More informationToward Electric Control of Magnetism
Toward Electric Control of Magnetism Collaorators Yoshi Tokura Dept of Appl hys, Univ Tokyo Correlated Electron Research Center (CERC), AIST and ERATO Multiferroics roject,jst Univ. Tokyo Y. Yamasaki,
More informationMagnetism in transition metal oxides by post-dft methods
Magnetism in transition metal oxides by post-dft methods Cesare Franchini Faculty of Physics & Center for Computational Materials Science University of Vienna, Austria Workshop on Magnetism in Complex
More informationHydrogenation of Penta-Graphene Leads to Unexpected Large. Improvement in Thermal Conductivity
Supplementary information for Hydrogenation of Penta-Graphene Leads to Unexpected Large Improvement in Thermal Conductivity Xufei Wu, a Vikas Varshney, b,c Jonghoon Lee, b,c Teng Zhang, a Jennifer L. Wohlwend,
More informationFirst-Principles Calculation of Exchange Interactions
Chapter 2 First-Principles Calculation of Exchange Interactions Before introducing the first-principles methods for the calculation of exchange interactions in magnetic systems we will briefly review two
More informationPeriodic DFT Study of Molecular Crystals
, March 13-15, 2013, Hong Kong Periodic DFT Study of Molecular Crystals Richard Rivera, Soraya Jácome, Darwin Castillo, Arvids Stashans 1 Abstract Two molecular crystals have been studied using the first-principles
More informationStructural Effect on the Oxygen Evolution Reaction in the Electrochemical Catalyst FePt
New Physics: Sae Mulli, Vol. 65, No. 9, September 2015, pp. 878 882 DOI: 10.3938/NPSM.65.878 Structural Effect on the Oxygen Evolution Reaction in the Electrochemical Catalyst FePt Wonseok Jeong Gijae
More informationUnconventional magnetic order in 3D Kitaev materials revealed by resonant x-ray diffraction Radu Coldea
Unconventional magnetic order in 3D Kitaev materials revealed by resonant x-ray diffraction Radu Coldea Oxford Collaborators Alun Biffin (Oxford->PSI) Roger D. Johnson S. Choi P. Manuel A. Bombardi Sample
More informationMonte Carlo Study of Planar Rotator Model with Weak Dzyaloshinsky Moriya Interaction
Commun. Theor. Phys. (Beijing, China) 46 (2006) pp. 663 667 c International Academic Publishers Vol. 46, No. 4, October 15, 2006 Monte Carlo Study of Planar Rotator Model with Weak Dzyaloshinsky Moriya
More informationHIGH-PRESSURE CHARACTERISTICS OF α Fe 2 O 3 USING DFT+U
HIGH-PRESSURE CHARACTERISTICS OF α Fe 2 O 3 USING DFT+U G. ROLLMANN and P. ENTEL Institute of Physics, University of Duisburg-Essen, 47048 Duisburg, Germany A. ROHRBACH and J. HAFNER Center for Computational
More informationSupplementary Information: Construction of Hypothetical MOFs using a Graph Theoretical Approach. Peter G. Boyd and Tom K. Woo*
Electronic Supplementary Material ESI) for CrystEngComm. This journal is The Royal Society of Chemistry 2016 Supplementary Information: Construction of Hypothetical MOFs using a Graph Theoretical Approach
More informationSupporting information. The Unusual and the Expected in the Si/C Phase Diagram. Guoying Gao, N. W. Ashcroft and Roald Hoffmann.
Supporting information The Unusual and the Expected in the Si/C Phase Diagram Guoying Gao, N. W. Ashcroft and Roald Hoffmann Table of Contents Computational Methods...S1 Hypothetical Structures for Si
More informationMagnetoelectric effect
Department of Physics Seminar Magnetoelectric effect The challenge of coupling magnetism and ferroelectricity Luka Vidovič Mentor: prof. dr. Denis Arčon Ljubljana, december 2009 Abstract Magnetism and
More informationCritical entanglement and geometric phase of a two-qubit model with Dzyaloshinski Moriya anisotropic interaction
Chin. Phys. B Vol. 19, No. 1 010) 010305 Critical entanglement and geometric phase of a two-qubit model with Dzyaloshinski Moriya anisotropic interaction Li Zhi-Jian 李志坚 ), Cheng Lu 程璐 ), and Wen Jiao-Jin
More informationDepartment of Chemistry, North Carolina State University, Raleigh, North Carolina 27695,
Most spin-1/2 transition-metal ions do have single ion anisotropy Jia Liu 1, Hyun-Joo Koo 2, Hongjun Xiang 3, *, Reinhard K. Kremer 4 and Myung-Hwan Whangbo 1, * 1 Department of Chemistry, North Carolina
More informationQuasi-1d Frustrated Antiferromagnets. Leon Balents, UCSB Masanori Kohno, NIMS, Tsukuba Oleg Starykh, U. Utah
Quasi-1d Frustrated Antiferromagnets Leon Balents, UCSB Masanori Kohno, NIMS, Tsukuba Oleg Starykh, U. Utah Outline Frustration in quasi-1d systems Excitations: magnons versus spinons Neutron scattering
More informationSelf-compensating incorporation of Mn in Ga 1 x Mn x As
Self-compensating incorporation of Mn in Ga 1 x Mn x As arxiv:cond-mat/0201131v1 [cond-mat.mtrl-sci] 9 Jan 2002 J. Mašek and F. Máca Institute of Physics, Academy of Sciences of the CR CZ-182 21 Praha
More information