## Giant magnetoelectric coupling in multiferroic PbTi 1 -x V O x 3 from density functional calculations

Lokanath Patra, † ‡ , Ravindran Vidya, ¶ Helmer Fjellv˚g, a § and Ponniah Ravindran* † ‡ ,

† Department of Physics, Central University of Tamil Nadu, Thiruvarur, Tamil Nadu, 610005, India

‡ Simulation Center for Atomic and Nanoscale MATerials, Central University of Tamil Nadu, Thiruvarur, Tamil Nadu, 610005, India

¶ Department of Medical Physics, Anna University, Chennai, 600025, India § Center for Materials Science and Nanotechnology and Department of Chemistry, University of Oslo, Box 1033, Blindern, N -0315, Oslo, Norway

## Abstract

The giant magnetoelectric coupling is a very rare phenomenon which has gained a lot of attention for the past few decades because of fundamental interest as well as practical applications. Here, we have successfully achieved the giant magnetoelectric coupling in PbTi 1 -x V O x 3 ( x =0 -1) with the help of a series of generalized -gradient -corrected (GGA), GGA including on -site coulomb repulsion ( U ) corrected spin polarized calculations based on accurate density functional theory. Our total energy calculations show that PbTi 1 -x V O x 3 stabilizes in C -type antiferromagnetic ground state for x &gt; 0.123. With the substitution of V into PbTiO , the 3 tetragonal distortion is highly enhanced accompanied by a linear increase in polarization. In addition, our band structure analysis shows that for lower x values, the tendency to form 2D magnetism of PbTi 1 -x V O x 3 decreases. The orbital magnetic polarization was calculated with self -consistent field method by including orbital polarization correction in the calculation as well as from the computed X -ray magnetic dichroism spectra. A non -magnetic metallic ground state is observed for the paraelectric phase for V concentration ( x ) = 1 competing with a volume change of 10% showing a large magnetovolume effect. Our orbital projected DOS as well as orbital ordering analysis suggest that the orbital ordering plays a major role in the magnetic to non-magnetic transition when going from ferroelectric to paraelectric phase. The calculated magnetic anisotropic energy shows that the direction [110]is the easy axis of magnetization for x = 1 composition. The partial polarization analysis shows that the Ti/V -O hybridization majorly contribute towards the total electrical polarization. The present study adds a new series of compounds to the magnetoelectric family with rarely existing giant coupling between electric and magnetic order parameters. These results show that such kind of materials can be used for novel practical applications where one can change the magnetic properties drastically (magnetic to non -magnetic as shown here) with external electric field and vice -versa.

## Introduction

The promising coupling between electric and magnetic order parameters and the potential to manipulate one by the application of the other has created much attention in the past few decades. Research on materials with such coupling has grown interest because of their wide range of application in multifunctional devices. 1-5 The ultimate goal of these research is to obtain single phase multiferroics with strong coupling between ferroelectric and magnetic order parameters at room temperature. There has been a considerable recent interest in developing lone pair based magnetoelectrics because of their high value of electrical polarization. Examples for such composites include BiFeO 3 -LaFeO , 3 6 BiFeO 3 -SrTiO , 3 7 BiFeO PbTiO , 3 3 8 BiFeO -BiCoO , 3 3 9 LaFeO 3 -PbTiO , 3 10 BiFeO 3 -BaTiO , 3 11 BiCoO 3 -BaTiO 3 etc. PbTiO 3 (PTO) is a perovskite ferrolectric material with a Curie temperature of 490 ◦ C and a large tetragonal distortion with c/a = 1.06 13 at room temperature. But it is nonmagnetic due to the absence of d electrons. The Ti atom can be substituted with some magnetic ions as a result of which the new compound can be expected to produce both magnetic and ferrolectric behavior. PTO based magnetoelectrics could be much more interesting because of the existence of two types of mechanisms for ferroelectricity i.e. lone pair electrons from Pb 2+ and d 0 -ness from Ti 4+ . So, here we examine the magnetoelectric properties of a PTO based multiferroic series i.e. PbTi 1 -x V O x 3 ( x = 0, 0.25, 0.33, 0.50, 0.67, 0.75, 1) where the 6 s electrons of Pb 2+ and the d 0 -ness of Ti 4+ stabilizes the ferroelectricity and the V induces magnetism.

On the other hand, PbVO 3 (PVO) , which is a member of this series ( x = 1) has stirred a lot of interest in last few years as a strong candidate for multiferroic oxide due to its large electric polarization. Though it is isostructural with PTO, it shows large tetragonal distortion ( c/a = 1.229). Isolated layers of corner shared VO 5 pyramids form a layer type perovskite structure with a space group P mm 4 . 14 PVO has been proposed to have an antiferro-

12

magnetic ordering and a ferroelectric polarization as large as 152 µ C/cm 2 due to the presence of large structural distortions. 14-18 To this point, however, the true magnetic structure and the multifunctional nature of this material is controversial to researchers. Shpanchenko et al. 14 found no long-range magnetic ordering with their neutron powder diffraction measurements down to 1.5 K. Belik et al. 16 proposed that PVO is a two -dimensional spin half square lattice strongly frustrated antiferromagnet due to the antiferromagnetic interactions of the next -nearest -neighbours. With the help of magnetic susceptibility and specific heat measurements as well as band -structure calculations, Tsirlin 19 confirmed that the S = 1/2 square lattice of Vanadium 4+ ions in PVO is strongly frustrated due to the next -nearest -neighbor antiferromagnetic interactions and no long-range magnetic ordering was found down to 1.8 K. Due to the presence of defects or ferromagnetic impurities in the sample, it has been very difficult to clarify the intrinsic magnetic property in PVO. Oka et al. 20 investigated the magnetic properties of PVO by preparing a multidomain single -crystal without any magnetic impurity. The broad maximum centered around 180 K in the temperature dependent magnetization curve indicates the presence of two -dimensional antiferromagnetism. Muon spin rotation ( µ SR) measurement displayed the presence of a long-range order below 43 K. The epitaxial thin films of PVO has been grown 21 using pulsed laser deposition and this brings a step forward to synthesis multiferroic materials outside of high -temperature and high pressure techniques to realize devices with multifunctionalities. In another study on PVO thin films by Kumar et al. , 22 a transition from a ferroelectric only state to a ferroelectric and magnetic state was determined below 100 -130 K using second -harmonic generation and X -ray linear dichroism. Experimental electron energy loss spectroscopy (EELS) investigation on V -L edge show that V in the PVO thin films are in the V 4+ state resulting in a d 1 state. 23

PVO was also investigated computationally by different researchers. First principle cal-

culations were made by Uratani et al. 17 for PVO along with BiCoO . 3 They found the easy axes of spin are different: [110] in PVO and [001] in BiCoO 3 even though both have similar crystal structure. A spin spiral structure was predicted by Solovyev 24 to analyze the absence of long -range magnetic ordering in PVO. Calculations along with experiments were performed by Parveen et al. 25 in 2012 to study the thermal properties of PVO and the results confirmed the observations made by Tsirlin et al. 19 Ming et al. 26 made a comparative study of the structural, electronic, magnetic and phase transition properties in which various exchange -correlation (XC) functionals were used and found that PVO is a 2D C -AF where the d 1 electron of the V 4+ ion occupies the d xy orbital. A ferroelectric to paraelectric phase transition at 1.75 GPa was also noticed. Zhou et al. 27 revisited the structural transitions in PVO with the help of a series of X -ray diffraction (XRD) measurements and first principle calculations. They found that the C -AF insulating and NM metallic states are the ground states for tetragonal and cubic phases, respectively. They have also noticed a noncentrosymmetric tetragonal to centrosymmetric cubic perovskite structural phase transition occurs between a pressure range of 2.7 -6.4 GPa.

Miloˇ sevi´ c et al. 28 performed ab initio calculations to study the electronic structure and optical properties of PVO (and BiCoO ). With the 3 help of first principle DFT calculations Xing et al. 29 showed a first order tetragonal to cubic phase transition with a volume collapse of 11.6% under uniaxial pressure of 1.2 GPa accompanied by a C -AF insulator to a NM metal. Kadiri et al. 30 studied the magnetic properties of PVO with the help of first principle calculations and Monte Carlo simulations and they determined PVO as an S = 1/2 antiferromagnet with a N´ eel temperature T N = 182 K. Recently, Oka et al. 31 examined the cubic phase of PVO under high pressure on experimental and theoretical bases. They determined the transition pressure to be 3 GPa and shown that above which the semiconductor -to -metal transition associated with the structural tran- sition (tetragonal -to -cubic) occurs.

There are some reports about the transition metal substitution at the V site of PVO. Tsuchiya et al. 32 synthesized PbFe 0 5 . V 0 5 . O un3 der high pressure and found the crystal structure as a tetragonal perovskite with c/a = 1.18. The magnetic study revealed that the compound is an antiferromagnet and the electrical polarization is estimated to be as large as 88 µ Ccm -2 . Ar´ evalo et al. 33 studied Ti and Cr substitutions at V site in PVO. A tetragonal to cubic transition was observed as the Cr substitution level reaches 0.4 where as the Ti substitution preserves the P mm 4 symmetry. The interesting result observed by Ar´ evalo et al. 33 was the temperature induced phase transition from tetragonal to cubic in PbTi 0 8 . V 0 2 . O 3 at 730 K. But PVO decomposes to Pb V O 2 2 7 at 570 K before reaching its transition temperature. 15 In another V site Ti substituted Pb M O 3 system (where M stands for V site substituted with Ti), Ricci 34 predicted multiferroicity in PbTi 0 875 . V 0 125 . O 3 with a ferromagnetic -ferroelectric ground state and an electrical polarization of 95 µ C/cm . 2 Recently, Pan et al. 35 experimentally studied the Pb(Ti,V)O 3 system where the V substitution level varies from 0.1 to 0.6. They have found that the whole composition range the P mm 4 symmetry was preserved and the tetragonality was abnormally improved (so the spontaneous polarization) with the replacement of Ti with V. Interestingly, they also observed an intrinsic giant volume contraction of ( ∼ 3.7%) for Pb(Ti 0 7 . V 0 3 . )O 3 during the ferroelectric -to -paraelectric phase transition. These observed interesting findings in the substituted systems motivated us for the present study.

In the present work, we have shown that the V substitution in PTO can induce magnetism which can lead to multiferroicity in the substituted system. Thus our ultimate aim to design magnetoelectric materials from non -magnetic system is satisfied. The coupling between the electric and magnetic order parameters is strong due to the presence of two ferroelectric mechanisms (lone pair of Pb 2+ as well as d 0 -ness of Ti 4+ ). The 2D magnetism arises in

PVO due to the d xy type orbital ordering in the ferroelectric ground state. On the other hand, the disorder induced by the substitution may release the frustration in PVO and that can lead to a long-range magnetic ordering. We hope the substitution could reduce the tendency of the system to form a 2D arrangement of V cations and it may lead to form 3D magnetic ordering.

## Computational details

All the results presented here are obtained from ab-initio density functional calculations using the Vienna Ab initio Simulation Package (VASP) 36 or full potential LAPW method (Wien2k). To obtain the ground state structural parameters for PbTi 1 -x V O , x 3 we have optimized the crystal structure by minimizing the force and stress acting on the system. The generalized gradient approximation (GGA) 37 in the scheme of Perdew -Burke -Ernzerhof (PBE) 38 is employed to treat the exchange -correlation as it gives better equilibrium structural parameters than the local density approximation. The ionic positions and the shape of the crystals were relaxed for different unit cell volumes until the energy as well as force convergence criterion were reached i.e. 10 -6 eV per cell and 1 meV ˚ A -1 /atom, respectively. As the structural parameters are very sensitive to electrical polarization, a very high plane wave cut -off energy of 850 eV 39 was used. A 6x6x6 Monkhorst -Pack k -point 40 mesh was used for ferroelectric PTO and similar density of k -points was used for all other calculations. To find the ground state magnetic ordering, we have considered non -magnetic (NM), ferromagnetic (FM), A -type antiferromagnetic ( A -AF), C -type antiferromagnetic ( C -AF), and G -type antiferromagnetic ( G -AF) configurations. 41 In order to account for the strong electron correlation associated with V 3 d electrons on the electronic structure, GGA+ U method was used with U eff = 3 eV.The Born effective charges were calculated by the Berry phase method with a 8x8x8 k -point mesh per f.u. using the modern theory of polar- ization. 42,43 We have used periodic supercell approach to simulate the substituted systems i.e. for compositions with 33%, and 50% V substitutions 1x1x3, and 1x1x2 supercells were created. Then one Ti atom was replaced with V to mimic the experimental conditions. The supercells are again doubled in each direction to incorporate different magnetic orderings. Similar kind of methodology was used to model other compositions.

The orbital moments at the transition metal sites were calculated using the orbital polarization correction (SO + OP) method implemented in Wien2k. 44 In addition, we have simulated X -ray circular magnetic dichroism (XMCD) 45 spectra using Wien2k and estimated the site selective spin and orbital magnetic moments of V ions using the XMCD sum rules. 46,47 For simulating XMCD, we have used full-potential linearized augmented plane wave (FP-LAPW) method based on density functional theory as implemented in the Wien2k code. The generalized gradient approximation within the PBE scheme was used as the exchange correlation potential. The muffin tin radii for Pb, V and O were chosen as 1.66, 2.49, and 1.50 a.u., respectively. The Brillouin zone integration has been carried out with 5000 kpoints. The energy convergence with a convergence criterion of 10 -5 was obtained. The spin orbit coupling is considered in the calculation. As the ferromagnetic state is a metal, the plasma frequency was first calculated and included in the XMCD calculation.

## Results and Discussions

## Structural phase stability

PVO is a C -type antiferromagnet (See Fig. 1) due to the presence of a d electron with S = 1/2. In oxides, it can either occur in the middle of an octahedron cage (CaVO , SrVO ) 3 3 48 or it can form a strong vanadyl bond with one O with a short bond length ( ∼ 1.55 -1.60 ˚ ). A 49 PVO is reported to exhibit a tetragonal -to -cubic phase transition under pressure but no structural transition was observed from 0 K to its de-

Figure 1: The C -type antiferromagnetic structure of PVO. Pb and V atoms are labelled on the illustration as green (big) and blue (medium) spheres, respectively. Oxygen atoms occur at the corners of the VO 5 square -pyramid and are labelled as red (small) spheres . Two types of oxygen atoms are present in PbVO 3 one at apical (labelled as O1) and the other at the planar (labelled as O2) position. The optimized atomic positions are: Pb [0, 0, 0], V [0.5, 0.5, 0.5708 (0.5677)], O1 [0, 0, 0.2128(0.2087)] and O2 [0, 0, 0.6902(0.6919)]. The values in the bracket are from experimental measurements. 16

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00238v1-figure-8a8bf8c9-77ab-41b7-9942-e4b56b5d1cda.jpg)

composition temperature. 16 Moreover, the minimum (apical) V -O bond length in the VO 5 pyramid in tetragonal PVO is 1.68 ˚ . This relA atively longer bond shows that a comparatively weaker vanadyl bond than usual. The average of the four planar V-O bond lengths is 1.99 ˚ . A In agreement with these findings, the calculated Goldschmidt tolerance factor for PVO also shows a value of 1.036 which usually indicates non-centrosymmetric distortion.

Here, we have made structural optimizations for both the ferrolectric ( P mm 4 ) and paraelectric ( Pm m ¯ 3 ) phases of PbTi 1 -x V O x 3 for x = 0, 0.25, 0.33, 0.50, 0.67, 0.75 and 1 compositions. The substitution of V for Ti in PTO doesn't affect the symmetry of the crystal as it retains the tetragonal symmetry for the whole composition range. 35 With increasing V-content ( x ),

Figure 2: The evolution of lattice parameters in PbTi 1 -x V O x 3 as a function of x . The rapidly increasing c and slowly decreasing a result in unusual enhancement in c/a value. The inset shows the octahedral environment changes into pyramidal one due to increase in tetragonal distortion with increasing x .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00238v1-figure-9383475c-6759-43df-9679-7df09ff06b19.jpg)

the lattice parameter c increases almost linearly and a decreases gradually as shown in Fig. 2. This results in an unusual enhancement in the tetragonality of the PbTi 1 -x V O . x 3 This can be understood on the basis of ionic radii of Ti 4+ and V 4+ ions. The coordination number of Ti 4+ changes from 6 to 5 whilst increasing concentration of vanadium. The ionic radius of V 4+ is 0.53 ˚ ; whereas, the ionic A radii of Ti 4+ are 0.61 ˚ A and 0.51 ˚ for 6 fold A and 5 fold co -ordinations, respectively. 50 So with the decrease in the ionic radii, the lattice parameter a b ( ) decreases in order to accommodate the smaller cation. On the other hand, with the increase of V concentration, the chance of vanadyl bond formation increases which leads to an increase in the length of c -axis. This is also the reason behind the increase in volume with the increase of x . It may be noted that, PTO modified ferroelectrics such as (1x )PTO -x Bi(Ni 1 2 / Ti 1 2 / )O 3 51 and (1x )PTO -x Bi(Mg 1 2 / Ti 1 2 / )O 3 52 usually do not show high value of c/a ratio. But some of the PTO based compounds i.e. (1x )PTO -x BiFeO , 3 53 BiInO 3 -PTO 54 and (1x )PTO -x Bi(Zn 1 2 / Ti 1 2 / )O 3 55 show abnormally high tetragonality. In the present

study, the achieved maximum tetragonality ( c/a ) for a PTO based compound i.e. PbTi 0 25 . V 0 75 . O 3 is 1.24. This value is larger than that of some of the conventional perovskite ferroelectrics such as 0.5PTO -0.5BiFeO 3 ( c/a =1.14), 0.85PTO -0.15BiInO 3 (1.08) and 0.6PTO -0.4Bi(Zn 1 2 / Ti 1 2 / )O 3 ( c/a = 1.11). PTO -based ferroelectric materials with large c/a posses interesting physical properties, such as high value of ferroelectric polarization ( P S ), high Curie temperature ( T C ), and enhanced negative thermal expansion (NTE). In this study, we have shown that PTO -based material with large c/a value possesses strong magnetoelectric coupling as well as giant magnetovolume effect.

The unit cell volume is sensitive to structural transition in PbTi 1 -x V O x 3 i.e. it contracts when a ferroelectric to paraelectric transition takes place. For x = 0.25, the volume contracts by 4.2% during the phase transition. Most interestingly, the volume contraction increases with increase in V content. The volume contraction value increases to 6.3% for x = 0.50 and 8.4% for x = 0.67. When x = 1, the volume contracts by a value as large as 10.4% at the ferroelectric to paraelectric phase transition point. One of the main factors for volume contraction in this series of compounds at the phase transition is the change in bonding between Pb -O and Ti/V -O which is discussed below. PTO based composites with a high value of volume contraction during phase transition can give rise to interesting properties like high NTE. Moreover, a temperature induced phase transition is observed experimentally for x ≤ 0.3. 35 So it can be said that for higher amount of Ti substitution in PbVO , the 3 T C can be measurable without decomposition.

The cubic structure of PVO can be considered as a special case of ideal tetragonal perovskite with c/a =1. The ferroelectric polarization is very sensitive to structural parameters. So, to calculate the c/a ratio for various compositions accurately, we have calculated the total energy as a function of c/a . We have performed this calculation with a 1x1x2 supercell which means we should get a minimum energy for c/a =2 for a cubic system. In consistent with this, our E

∼ c/a graph (See Fig. S1 in supporting information) for PVO ( x = 0) shows a minimum energy for c/a =2 for the paraelectric phase; whereas, for the ferroelectric phase, we get a non -integer c/a . In order to find the possible polarization path, we have calculated the total energy as a function of Pb displacement with respect to V/Ti -O polyhedra. Figure 3 helps us to find the easiest path of the ferroelectric polarization and the role of Pb ions in the tetragonal distortion. We have plotted the displacement of Pb ion with respect to VO 5 pyramid along [111] and [001] directions. The energy vs displacement curves for both the directions have double well shapes. It can be seen that the lowest energy for the off -center displacements is along the [001] direction. It shows that the ferroelectric properties are greatly effected by the polarizability of Pb. The energy difference between the two directions is very small (about 2 meV/f.u.) compared to well -known multiferroic BiFeO . 3 39 So, in principle, the Curie temperature of PVO must be less than that of BiFeO . 3 However, the T C for PVO has not been estimated yet, as it decomposes before reaching to its paraelectric phase. But the T C of PbTi 0 9 . V 0 1 . O 3 was estimated to be 823 K by Zhao Pan et. al. 35 and it showed an increasing trend V concentration (Measured T C for PbTi 0 8 . V 0 2 . O and PbTi 3 0 7 . V 0 3 . O are in be3 tween 823 K to 873 K). Compositions with x &gt; 0.3 decompose before reaching their T C .

We have also performed nudged elastic band (NEB) calculations in order to find out the energetics of ferroelectric to paraelectric phase transition and to identify a minimum energy path for the phase transition. To calculate the minimum energy path, we have considered the tetragonal phase of PbVO 3 as the initial structure whereas the cubic phase as tetragonal structure. Different intermediate images were created and the energies for all the images were calculated considering the volume as sell as shape variation. The energies obtained are plotted versus the images (here mentioned as reaction coordinates).The graph given in Fig. 4 is resulted from NEB calculations and it shows that the ferroelectric-to-paraelectric transition involves an energy barrier of 127 meV which is

Figure 3: (Color online) Total energy as a function of displacement of Pb ion along [001] and [111] directions for the ferroelectric PVO. The structural relaxations were consider through selective dynamics.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00238v1-figure-08fb33d6-653f-4ccd-a330-eaa783841562.jpg)

higher than the difference between the groundstate total energies of these two phases i.e. ∼ 70 meV. It also shows the V-O coordinated polyhedra. In case of ferroelectric phase, the V-O polyhedron is highly distorted and makes an pyramidal arrangement. The intermediate images along the FE-to-PE transition path with higher energies than the FE phase possess comparatively less distorted polyhedra. Finally, the end product i.e. the paraelectric phase exhibits an undistorted octahedron. The high energy barrier for the phase transition is consistent with the high temperature FE-to-PE phase transition.

## Chemical Bonding

The chemical bonding analysis in compounds similar to PVO i.e. BaTiO 3 and PTO shows that the ferroelectric instability arises due to the hybridization interaction between Ti 3 d and O 2 p states. 56 So it is interesting to analyze the chemical bonding between V/Ti and O in PbTi 1 -x V O x 3 in order to reach the depth of the origin of ferroelectricity. Experimental results and theoretical analysis also show that a significant hybridization is present between Pb and O. 57 This also contributes to the Born effective charges (BEC) and consequently to the polarization.

Figure 4: The ferroelectric to paraelectric transition path for PVO obtained from nudged elastic band method. Both volume and shape variation during the transition are considered. The insets show the V and O coordination for the particular images marked by arrows. The atoms colors are same as in Fig. 1

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00238v1-figure-09d79779-7ea3-4b79-865c-8f31a777ebc5.jpg)

The charge density plot for x = 0, 0.50, and 1 compositions are given in Fig. 5. It can be seen from the Fig.5 (a), (b), and (c) that the charge density between the B -site cation and apical oxygen increases as the V concentration increases. In PVO, the presence of strong V -O1 bond along with lone pair electrons at the Pb 2+ site weaken the Pb -O1 bond. Due to these, the Pb atom shifts more towards O2 to maintain the charge balance. Hence, the covalent bonding interaction between Pb and O2 shows increasing trend which is evident from Fig.5 (d), (e) and (f) as the width of the charge density increases with increase in x . These observations indicate that the increase of V concentration in PbTi 1 -x V O x 3 not only increases the covalent bonding between Ti/V -O1 but also that between Pb -O2. As a consequence of these, the tetragonal distortion increases with increase in x .

The charge density plots for PVO in cubic and tetragonal phases are given in Fig. 6(a) upper and lower panel, respectively. It can be seen that, the charge density distribution at the Pb and O sites are essentially isotropic in nature for the paraelectric phase. On the other hand, finite amount of charge can be seen in between Pb and O in the tetragonal phase showing an anisotropic nature of charge density distribu-

Figure 5: Charge density between (a) Ti and O1 in PbTiO , (b) V and O1 in PbTi 3 0 5 . V 0 5 . O , 3 (c) V and O1 bond in PbVO , (d) Pb and O2 3 in PbTiO , 3 (e) in PbTi 0 5 . V 0 5 . O , 3 and (f) in PbVO . 3

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00238v1-figure-a8c66351-cc22-4f6d-bf1b-44913ae109f6.jpg)

tion. So, it can be stated that in cubic phase the Pb-O bond has more ionic character and the tetragonal phase has a mixed iono-covalent bonding character. On the other hand, the bonding between V and O is found to have substantial covalency in both the cubic and tetragonal phases as we can find noticeable charge density between these atomic sites. But, it may be noted that, the covalent interaction between V and O in tetragonal phase is stronger than that in the cubic case. Due to this, the V atom is shifted from the center of the octahedron stabilizing the ferroelectric state..

Charge transfer distribution shows isotropic nature at the Pb and O sites (Fig. 6(b)) in the cubic phase confirming the ionic bonding between Pb -O. The anisotropic nature of charge transfer distribution between V and O indicates the presence of finite covalent bonding between these atoms in the cubic phase. The charge transfer distribution between V and O as well as Pb and O in tetragonal phase shows anisotropic behavior as shown in Fig. 6(b) indicating the presence of finite covalent bonding between them. It may be noted that, though the covalency between Pb -O is weaker than that between V -O in the tetragonal phase, it has equal importance in structural and ferroelectric properties.

Electron localization function (ELF) also provides important information about the bonding interaction between the constituents of a com-

Figure 6: (a) Charge density, (b) Charge transfer, and (c) Electron localization function for the paraelectric (upper panel) and ferroelectric (lower panel) structures.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00238v1-figure-c13664b2-de0e-4dc3-9f7f-2256b0e54e3d.jpg)

pound. It can be seen from Fig. 6(c) that the ELF is maximum at O sites and minimum at Pb and V sites for both the paraelectric and ferroelectric phases. Finite ELF can be seen in between V and O in the tetragonal phase showing the presence of finite covalent bonding between them. The presence of the stereochemically active lone pair electrons at Pb site is also clearly visible in the ELF plot for the tetragonal case. So the bonding interaction between Pb and O as well as V and O in the

Figure 7: (Color online) The valence electron localization function isosurface plotted at a value of 0.7 for PVO in the ferroelectric P mm 4 structure.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00238v1-figure-6eb5d13c-57c6-42ef-82cc-7d26497bda78.jpg)

Figure 8: The site and orbital projected density of states (DOS) for PVO for the respective ground state magnetic structures in the (a) paraelectric and (b) ferroelectric phases. The Fermi level is set to zero.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00238v1-figure-1929eb03-74e0-4d84-b30e-094493c6d34f.jpg)

ferrolectric tetragonal phase can be concluded as having dominant ionic character with finite covalent bonding. The difference in the bonding behaviour of Ti/V -O and Pb -O bonds for paraelectric and ferroelectric phases play a crucial role for the large volume contraction during the ferroelectric-to-paraelectric phase transition. We have also plotted the isosurface of the valence electron localization function with an ELF value of 0.7 and given in Fig. 7. This figure shows the lobe like lone pairs of hybridized 6 s electrons at the Pb site. The charge density and charge transfer plot show that finite amount of electrons from Pb sites covalently interact with O 2 p ensuring that 6s electrons of Pb are not completely chemically inert.

The site and orbital projected DOS for paraelectric and ferroelectric phases for their respective ground state magnetic configurations are given in Fig. 8 (a) and (b), respectively. It can be seen that the V 3 d states and O 2 p states are spread over in the valence band in a range of -7 eV to Fermi level (E F ). The energetically degenerate nature of these two states indicates the presence of strong covalent hybridization between V and O. This hybridization weakens the short range repulsion to lower the energy of the ferroelectric phase. The filled O 2 s states form narrow bands around -18 eV (see Fig. 8(b)). Above that, at around -9 eV, the Pb 6 s lone pair electronic states are present with small contributions from 2 p electrons. This also confirms that the 6 s electrons are stereochemically active and not inert. Though the Pb 6 p states are well separated from its 6 s states and are mostly present in CB, a small amount of Pb 6 p states can be seen which are degenerate with O 2 p states confirming the Pb -O covalency. The hybridization between these states enhances the stability of the ferroelectric phase by lowering the total energy of the system. The partial DOS for PbTi 0 5 . V 0 5 . O 3 calculated with GGA+ U ( U eff = 3 eV) is given as Fig. S2 in the supporting information. Similar kind of DOS distribution can be seen for Pb, V and O atoms in PbTi 0 5 . V 0 5 . O also. The almost empty 3 d states indicates the d 0 -ness of Ti 4+ . Due to the covalent interaction between Ti and O,

the charges redistribute among themselves for which we see very small yet noticeable DOS in Ti valence band.

The nature and strength of the bonding can also be determined by crystal orbital Hamilton population (COHP) analysis where the negative and positive COHP values indicate bonding and anti -bonding interactions, respectively. We have calculated the COHP among the constituents of PVO (Pb -O, Pb -O2, V -O1, V -O2, O1 -O2) and given in Fig. 9. It can be seen that, all the occupied states for Pb -O and V -O interactions have bonding states for both paraelectric and ferroelectric phases. The Pb -O2 interaction for both the paraelectric and ferroelectric phases are almost the same. But the interaction between Pb -O1 in the ferroelectric case has more bonding states which are due to the hybridization between O 2 p and Pb 6 /6 s p states. This indicates the presence of stronger covalency between Pb and O1 in case of the ferroelectric phase than the paraelectric one which is also consistent with our charge density and partial DOS analyses. In the ferroelectric phase the V -O1 bonding interaction is stronger than V -O2 bonding interaction which is due to the off -center displacement of V atom towards the O1 atom.

## Electronic Structure

The total DOS given in Fig. 10 (a) for the paraelectric phase shows a metallic behavior for all magnetic configurations. From the total DOS of the ferroelectric phase (Fig. 10(b)), it can be seen that the nonmagnetic phase shows a metallic state; whereas, the ferromagnetic case is a half metal with a gap of ∼ 0.5 eV in the down spin channel. But, when the antiferromagnetic ordering is introduced, the exchange potential produced by the exchange interaction pushes the V 3 d states towards the lower energy side which results in an insulating state with a narrow gap. So, a combination of spin polarization, magnetic ordering as well as crystal structure play important roles in stabilizing the ferroelectric phase of PVO. The broad band spreading from -7 to -2 eV is mainly due to the strong hybridization between V 3 d O 2 p elec-

Figure 9: Crystal orbital Hamilton population (COHP) for PVO in the paraelectric and ferroelectric phases describing Pb -O, Pb -O2, V -O1, V -O2, O1 -O2. O1 and O2 are the planar and apical oxygen atoms, respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00238v1-figure-ff96e68f-f472-4dd4-a923-ac736ac2401d.jpg)

trons. The localized states at around -9 eV are generated from the hybridized Pb 6 s lone pair electrons. After including the correlation effect through the Hubbard U into the Hamiltonian matrix, the states get more localized to give a larger band gap and this localization of the d electrons produces a larger magnetic moment.

The total DOS for x = 0.33, 0.50 and 0.67 compositions were calculated with the GGA+ U method and given as Fig. S3 in the supporting information. The band gap value decreases with increase in x . As it is known that the bond lengths and lattice parameters have direct impact on band gaps, the decreasing band gap may be attributed to the increase in the average Ti/V -O bond length with x in Ti/V -O 6 polyhedra (The average Ti -O bond length is 2.07 Ain PTO whereas the V ˚ -O bond length is 2.11 ˚ in case of PVO). The lattice parameter A c also increase with x .

The band structure for both the para -and ferroelectric phases of x = 1 composition are given in Fig. 11. The k -path considered in the band structure is given as Fig. S4 in the supporting information. The lowest bands seen around 9 eV below the E F are due to the Pb

Figure 10: Total density of states for PVO in the nonspin polarized (NSP), ferromagnetic, A -, C -, and G -type antiferromagnetic configuration for (a) paraelectric and (b) ferroelectric phases. The Fermi level is set to zero.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00238v1-figure-d9c28bef-5f2b-4891-ae3e-99f3cdc5be1a.jpg)

6 s states. It is to be noted that, these bands are broader in the paraelectric phase; whereas, in the ferroelectric phase, they are more localized in a small energy range. This explains the effect of 6 s lone pair in stabilizing the ferroelectric state. A manifold of occupied O 2 p bands are spread over a range from -7 eV to -2 eV. A pair of bands separated from O 2 p bands and localized just below E F . These are having V 3 d character with no dispersion along Γ -Z direction and a very weak dispersion along X -R as well as M -A directions indicating a 2D characteristic in the ab -plane. 18 Above these band features a narrow gap with the value of 0.2 eV is present in the ferroelectric phase. The conduction bands located in the vicinity of CB-edge are derived from V 3 d states, but the bands at higher energy ( ∼ 3 eV) have mainly Pb 6 p character. We have also calculated the band structure with GGA+ U method. The GGA+ U band structures for x = 0.5 and x = 1 compositions are given in Fig. S5 and S6 in the supporting information. It can be seen that when the V concentration is less, the vanadyl bond weakens and therefore the d xy band gets dispersed along Γ -Z, X -R and M -A directions. This indicates that the tendency of the system to form the 2D magnetism is reduced for lower x values.

## Magnetic Properties

We have calculated total energies for nonmagnetic and all the magnetic orderings mentioned above for PbTi 1 -x V O x 3 in both paraelectric and ferroelectric phases and given in Table S1 in the supporting information. PTO is a well -known non -magnetic compound. Our total energy calculation reveals that the C -AF is the ground state for other compositions. It is to be noted that the C -AF and G -AF states are almost degenerate (See Table S1 in the supporting information). The present observation of C -AF state as a magnetic ground state for PVO is in consistent with experimental observation of a two -dimensional C -AF phase and also with other theoretical studies. 14,17 To identify the exact composition where the non -magnetic to magnetic transition happens, we have plotted ∆E vs x (where ∆ E = E C -AF -E NM ) for the ferroelectric phase and given in the supporting information Fig. S7. The energy difference increases linearly as the concentration of V increases due to the increase in localized d electrons density per cell making the antiferromagnetic state more stable. We have extrapolated the ∆E vs x curve and the

Figure 11: Electronic band dispersion for PVO in the (left) paraelectric and (right) ferroelectric phases. The Fermi level is set to zero.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00238v1-figure-89b7ecf4-d143-4d8b-801f-a9834423f9b5.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00238v1-figure-ba1a0d75-3d07-4efc-ad2a-9201418e6745.jpg)

magnetic transition point is found to be x ∼ 0.123. This is in agreement with previous theoretical finding by Ricci et al. 34 using similar calculations that the magnetic ground state stabilize in PbTi 1 -x V O x 3 for x = 12.5%.

The total energy shows a non -magnetic ground state for the paraelectric phase of PbTi 1 -x V O . So the ferroelectric to paraelecx 3 tric phase transition is associated with a magnetic to non -magnetic transition which shows that a strong coupling between the electric and magnetic order parameters for x &gt; 0.123. Now, the large values of volume contraction during the structural transition can also be attributed to the associated magnetic phase transition. So it can be said that, this series of compounds possess not only giant magnetoelectric effect, but also giant magnetovolume effect. It is to be noted that compounds with strong magnetovolume effect can be used to convert magnetic energy into kinetic energy, or vice-versa, and are used to build actuators and sensors.

empty whereas the t 2 g states are partially filled in the paraelectric phase resulting in a metallic state. The ferroelectric phase of PVO has a square pyramidal coordination (Fig. 1) because of the off -center displacement of V. Due to the presence of VO 5 square pyramidal envi-

The orbital projected DOS for the V 3 d electrons for both paraelectric and ferroelectric phases in PVO are given in Fig. 12 (a) and (b), respectively. In the presence of a cubic octahedral crystal field the d level of transition metal cation splits into t 2 g (with degenerate d xy , d xy and d zx orbitals) and e g (with degenerate d x 2 -y 2 and d z 2 ) orbitals. The e g levels are completely

Figure 12: The orbital projected density of states for the 3 d electrons in the V site of PbVO 3 in the (a) paraelectric and (b) ferroelectric phases for their respective ground state magnetic structures. The Fermi level is set to zero.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00238v1-figure-63f622eb-2f6d-4bb2-9686-feda5501e786.jpg)

Figure 13: Orbital ordering pattern of PVO in (left) the ferroelectric and (right) paraelectric phases as derived from full potential calculations.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00238v1-figure-2ca64432-42fc-4df3-8663-2cdbcb3c7bb7.jpg)

ronment, the t 2 g level of V 3 d electrons splits into b 2 g ( d xy ) and doubly degenerate e g ( d yz , d xz ) levels with an energy difference of ∼ 1 eV. Similarly the degeneracy of e g level breaks to produce separate b 1 g ( d x 2 -y 2 ) and a 1 g ( d z 2 ) levels. 24 The orbital degeneracy is broken because of the second-order-Jahn-Teller ordering. This was also noticed in the YTiO 3 with d 1 configuration of Ti 3+ similar to V 4+ in our study. 58 The orbital projected DOS for the ferroelectric phase given in Fig. 12 (b) shows that the localized V 3 d electron are present in the lowest lying d xy orbital which is separated from other d orbitals by crystal field splitting, making an spin 1/2 antiferromagnetic ordering in the ab plane. Due to the localized nature of this d -electron there is a strong intra -atomic exchange splitting that shifts the unoccupied d xy states to relatively higher energy.

tric and paraelectric phases are given in Fig. 13. It can be seen that the V atom has a d xy pattern OO in the ferroelectric phase. This type of orbital ordering in PVO strengthens the intraplanar V-O-V interactions and also stabilizes the 2D magnetism. Hence, the magnetic and orbital ordering are ingeniously coupled with the lattice distortion and ionic displacements that result in strong magnetoelectric coupling. But, in the cubic paraelectric case, the V - -O V interactions are strong enough to create an itinerant-electron band with all the t 2 g orbitals partially filled. So, the OO pattern for paraelectric phase does not show particular manner of ordering from any of the t 2 g orbitals. This weakens the exchange coupling and stabilizes a non-magnetic ground state.

We know that the electrons present in the vicinity of Fermi level actively participate in electrical conduction and magnetic exchange interactions. So we have used the integrated values of the orbital projected DOS i.e. from E F to -1 eV to study the orbital ordering. We have used same procedure for OO study as given in our previous studies, 59,60 and such analysis not only provides the OO pattern but also a pictorial illustration of the orientation of a particular d orbital. The OO pattern for ferroelec-

The calculated magnetic moment at the V site is 0.93 µ B for the ferroelectric phase. The covalency present between V and O induces a magnetic moment of 0.05 µ B at O site. The orbital magnetic moment at the V site evaluated using SO + OP method is estimated to be -0.045 µ B . This indicates that the orbital magnetic moment is aligned antiparalelly with the spin magnetic moment. The magnitude of the orbital moment of PVO is almost 5 times less than that of BiCoO . 3 61 This is due to the the presence of relatively larger spin moment and stronger spin orbit coupling in BiCoO . 3

We have calculated the magnetic anisotropy energy, 62 63 for x = 0.5 and 1 in their ground state C -AF configurations. For this, we have calculated the variation of total energy by changing the direction of the magnetization with the force theorem for the compositions x = 0.5 and 1. The calculated relative magneto anisotropy energies for x = 0.5 with respect to the easy axis are 0, 0.125, 0.070 and 0.048 meV/f.u. and 0.592, 1.630, 0, and 0.185 meV/f.u. for x = 1 in the [001],[100],[110], and [111] directions, respectively. Therefore, the easy axis for PbTi 0 5 . V 0 5 . O and PVO are [001] 3 and [110], respectively. The change in magnetic easy axis from x = 0.5 to 1 can be attributed to the formation of vanadyl bond with increase in V concentration i.e. the x = 0.5 composition is a BiCoO 3 like C -AF (so it shows similar easy axis) and with increase in V concentration the tendency to form a 2D antiferromagnetism increases (which also reflects in our band structure).

In order to get deeper understanding of the magnetic properties of PVO, we have simulated the XMCD and X-ray absorption spectra (XAS) at the L 2 3 , edges of V in the ferromagnetic configuration of PVO. The XMCD and XAS spectra along with the left and right polarized spectra for V atom in PVO are given in Fig. 14. By applying the sum rules 46,47 given in supporting information, we have evaluated the spin and orbital moments at the V site. The obtained orbital moment of -0.039 µ B for V is consistent with that obtained directly from the self consistent calculations (-0.045 µ B ). The sum rule analysis results a spin moment of at 0.8 µ B /V the V L 2 3 , edges. However, this varies noticeably from the value calculated by SCF method which is ∼ 1 µ B /V. Detailed analysis on different 3 d compounds by Wende et al. 64 shows that the deviation of the sum rule results from the SCF predictions are larger for lighter 3 d elements 65-67 as integral sum rule analysis ignores the spectral shape of the XMCD. So, the detailed fine structure in the XMCD of V should be fitted with the multipole moment analysis as explained in Ref. [67]. 67

Figure 14: Absorption spectra at the V L 2;3 edges for right circular, left circular polarization of the exciting X-rays, X-ray magnetic circular dichroism and X-ray absorption spectra in PVO

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00238v1-figure-269e82da-04ea-469e-92e2-5797bb83931d.jpg)

## Born effective charge and Spontaneous Polarization

We have calculated BECs for all the compositions mentioned above. The diagonal components of the BECs for x =0, 0.5 and 1.0 compositions are given in Table 1. The formal valence for Pb, V, Ti and O in PbTi 1 -x V O x 3 are +2, +4, +4 and -2, respectively. But, the diagonal components of the calculated BEC tensors for the atoms are much larger than their nominal ionic values. This reveals that there is large dynamic contribution superimposed to the static charge. This is due to the strong covalency effect where a large amount of non-rigid delocalized charge flows across this compound during the ionic displacement. This results in additional charges with respect to the nominal ionic values (anomalous contribution) which are obtained at the atomic sites. Alarge BEC value at the oxygen site also suggests that a large force can be generated by applying a small field which can stabilize the polarized state. In an effort to shed light on the polarization of PbTi 1 -x V O x 3 and its origin, we have evaluated the polarization of PbTi 1 -x V O x 3 for x = 0 to 1. PTO and PVO are well studied ferroelectrics with polarization values of 92 µ C/cm 2 and 147 µ C/cm , 2 respectively. The calculated polarization values increase from 92 µ C/cm 2 for pure PTO to 105, 110, 118, 130, 135, and 147 µ C/cm 2 for x = 0.25, 0.33, 0.5, 0.67, 0.75, and 1, respec-

Table 1: Calculated diagonal components of the Born effective charge tensors for x = 0, 0.5 and 1 of PbTi 1 -x V O x 3 in the ferroelectric phase with ground state magnetic ordering.

| x →   | 0      | 0      | 0      | 0.5    | 0.5    | 0.5    | 1      | 1      | 1      |
|-------|--------|--------|--------|--------|--------|--------|--------|--------|--------|
| Atom  | xx     | yy     | zz     | xx     | yy     | zz     | xx     | yy     | zz     |
| Pb    | 3.641  | 3.641  | 3.331  | 3.623  | 3.617  | 3.664  | 3.251  | 3.232  | 3.275  |
| Ti    | 6.166  | 6.166  | 4.961  | 5.480  | 5.406  | 4.778  | -      | -      | -      |
| V     | -      | -      | -      | 4.348  | 4.440  | 4.435  | 4.880  | 4.824  | 3.749  |
| O1    | -1.938 | -1.938 | -4.245 | -1.148 | -1.154 | -3.098 | -1.442 | -1.525 | -3.720 |
| O2    | -2.603 | -5.266 | -2.024 | -2.078 | -3.875 | -1.417 | -2.209 | -4.029 | -1.764 |
| O3    | -5.266 | -2.603 | -2.024 | -3.716 | -1.899 | -1.423 | -3.877 | -2.513 | -1.724 |

tively. This is due to the increase in tetragonal distortion with x as discussed above. The polarization calculated with point charges showed comparatively low values. The reduction in polarization values using point charges indicates that covalency effect plays an important role deciding the polarization of these materials.

The atom -decomposed spontaneous polarization is given for x = 0, 0.5 and 1 compositions in Fig 15. The major contribution towards the polarization is coming from the displacement of the O atom from its centrosymmetric position. It can be seen that the polarization contribution of Ti towards higher x value decreases. In contrast, an opposite trend is seen

Figure 15: Calculated site projected spontaneous electric polarization of PbTi 1 -x V O x 3 for x = 0 (upper panel), 0.5 (middle panel) and 1 (lower panel).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00238v1-figure-df58b0d9-5e0b-4412-967c-ee328dc07759.jpg)

for V atom. This is because of the enhancement of covalency between B site cation and O with V substitution at Ti site.

## Conclusion

PbTi 1 -x V O x 3 possesses strong coupling between magnetic and electric order parameters for V concentration &gt; 12.3%. Lower substituted compositions are non-magnetic ferroelectrics. Also, from the other side, the inclusion of Ti at the V site increases the chance of forming a 3D arrangement of V cations rather than 2D. The strong covalency between Ti/V -O and noticeable covalency between Pb -O bring the non -centrosymmetric distortion and stabilize the ferroelectric ground state. The ferroelectric-to-paraelectric transition is accompanied by a magnetic to non -magnetic transition for x &gt; 0.123. There is a large volume contraction during the above transition, indicating a strong lattice -ferroelectric coupling and strong magnetovolume effect. The calculations show high values of spontaneous electrical polarizations which are mainly due to the displacement of the apical oxygen in the BO 6 octahedra in the paraelectric phase due to strong Ti/V -O covalent interaction. Therefore, the PbTi 1 -x V O x 3 is a multifunctional series of compounds having giant magnetoelectric effect where one can change the magnetic properties drastically (even from magnetic to non -magnetic as shown here) by applying electric field and vice versa, strong magnetovolume effect, induce 2D magnetism as well as orbital ordering by electric field,

strong ferrolectric -lattice coupling etc. Also, recent experimental studies by Zhao Pan et al. 35 show that for lower V concentration, this series shows good NTE response and also a temperature induced tetragonal to cubic transition happens which makes this series more interesting for the researchers to investigate its magneto/electro/multi -caloric properties. With all these properties coexisting in the series, PbTi 1 -x V O x 3 can be useful for many multifunctional device applications.

## Acknowledgements

The authors are grateful to the Research Council of Norway for providing computer time (under the project number NN2875k) at the Norwegian supercomputer consortium (NOTUR). This research was supported by the Indo -Norwegian Cooperative Program (INCP) via Grant No.F. 58-12/2014(IC). L. P. wishes to thank Prof. Anja Olafsen Sj˚stad for her fruita ful discussions. L. P. also thanks Mr Ashwin Kishore MR and Dr. A. Krishnamoorthy for critical reading of the manuscript.

## Associated content

Energy vs. c/a ratio curves, partial dos for x = 0.50 with GGA+ U method, total DOS for x = 0.33, 0.50, and 0.67 with GGA+ U method, the Briilouine zone path used for band calculations, electronic band structure for x = 0.50 with GGA+ U method, electronic band structure for x = 1 with GGA+ U method, variation of energy difference between magnetic and nonmagnetic states with x , and sum rule analysis. This information is available free of charge via the Internet at http://pubs.acs.org.

## Author Information

## Corresponding Author

*Email: raviphy@cutn.ac.in

Tel:

+91 94890 54267

## Notes

The authors declare no competing financial interest.

## References

- (1) Fiebig, M. Revival of the magnetoelectric effect. J. Phys. D 2005 , 38 , R123.
- (2) Spaldin, N. A.; Fiebig, M. The renaissance of magnetoelectric multiferroics. Science 2005 , 309 , 391-392.
- (3) Cheong, S.-W.; Mostovoy, M. Multiferroics: a magnetic twist for ferroelectricity. Nat. Mater. 2007 , 6 , 13-20.
- (4) Ramesh, R.; Spaldin, N. A. Multiferroics: progress and prospects in thin films. Nat. Mater. 2007 , 6 , 21-29.
- (5) Bibes, M.; Barth´ el´my, e A. Multiferroics: Towards a magnetoelectric memory. Nat. Mater. 2008 , 7 , 425-426.
- (6) Gonz´ alez-V´zquez, a O.; Wojdeglyph[suppress], l J. C.; Di´ eguez, O.; ´ ˜iguez, In J. First-principles investigation of the structural phases and enhanced response properties of the BiFeO -LaFeO 3 3 multiferroic solid solution. Phys. Rev. B 2012 , 85 , 064119.
- (7) Ma, Z.; Tian, Z.; Li, J.; Wang, C.; Huo, S.; Duan, H.; Yuan, S. Enhanced polarization and magnetization in multiferroic (1x )BiFeO 3 x SrTiO 3 solid solution. Solid State Sci. 2011 , 13 , 2196-2200.
- (8) Zhu, W.-M.; Guo, H.-Y.; Ye, Z.-G. Structural and magnetic characterization of multiferroic (BiFeO ) 3 1 -x (PbTiO ) 3 x solid solutions. Phys. Rev. B 2008 , 78 , 014401.
- (9) Di´ eguez, O.; ´ ˜iguez, In J. First-Principles Investigation of Morphotropic Transitions and Phase-Change Functional Responses in BiFeO -BiCoO 3 3 Multiferroic Solid Solutions. Phys. Rev. Lett. 2011 , 107 , 057601.

- (10) Singh, A.; Chatterjee, R. Magnetization induced dielectric anomaly in multiferroic LaFeO -PbTiO 3 3 solid solution. Appl. Phys. Lett. 2008 , 93 , 182908.
- (11) Hang, Q.; Xing, Z.; Zhu, X.; Yu, M.; Song, Y.; Zhu, J.; Liu, Z. Dielectric properties and related ferroelectric domain configurations in multiferroic BiFeO 3 BaTiO 3 solid solutions. Ceram. Int. 2012 , 38 , S411-S414.
- (12) Patra, L.; Pan, Z.; Chen, J.; Azuma, M.; Ravindran, P. Metamagnetism stabilized giant magnetoelectric coupling in ferroelectric x BaTiO -(13 x )BiCoO 3 solid solution. Phys. Chem. Chem. Phys. 2018 , 20 , 7021-7032.
- (13) Chen, J.; Hu, L.; Deng, J.; Xing, X. Negative thermal expansion in functional materials: controllable thermal expansion by chemical modifications. Chem. Soc. Rev. 2015 , 44 , 3522-3567.
- (14) Shpanchenko, R. V.; Chernaya, V. V.; Tsirlin, A. A.; Chizhov, P. S.; Sklovsky, D. E.; Antipov, E. V.; Khlybov, E. P.; Pomjakushin, V.; Balagurov, A. M.; Medvedeva, J. E.; Kaul, E. E.; Geibel, C. Synthesis, structure, and properties of new perovskite PbVO . 3 Chem. Mater. 2004 , 16 , 3267-3273.
- (15) Okos, A. Synthesis and characterization of Pb(V 1 -x M )O x 3 compounds. Ph.D. thesis, Universit´ de Grenoble, 2014. e
- (16) Belik, A. A.; Azuma, M.; Saito, T.; Shimakawa, Y.; Takano, M. Crystallographic features and tetragonal phase stability of PbVO , a new member of PbTiO 3 3 family. Chem. Mater. 2005 , 17 , 269-273.
- (17) Uratani, Y.; Shishidou, T.; Oguchi, T. First-principles study on the magnetic anisotropy in multiferroic PbVO 3 and BiCoO . 3 J. Phys. Soc. Jpn. 2009 , 78 , 084709-084709.
- (18) Singh, D. J. Electronic structure and bond competition in the polar magnet PbVO . 3 Phys. Rev. B 2006 , 73 , 094102.
- (19) Tsirlin, A. A.; Belik, A. A.; Shpanchenko, R. V.; Antipov, E. V.; Takayama-Muromachi, E.; Rosner, H. Frustrated spin-1/2 square lattice in the layered perovskite PbVO . 3 Phys. Rev. B 2008 , 77 , 092402.
- (20) Oka, K.; Yamada, I.; Azuma, M.; Takeshita, S.; Satoh, K. H.; Koda, A.; Kadono, R.; Takano, M.; Shimakawa, Y. Magnetic ground-state of perovskite PbVO 3 with large tetragonal distortion. Inorg. Chem. 2008 , 47 , 7355-7359.
- (21) Martin, L. W.; Zhan, Q.; Suzuki, Y.; Ramesh, R.; Chi, M.; Browning, N.; Mizoguchi, T.; Kreisel, J. Growth and structure of PbVO 3 thin films. Appl. Phys. Lett 2007 , 90 , 062903.
- (22) Kumar, A.; Martin, L. W.; Denev, S.; Kortright, J. B.; Suzuki, Y.; Ramesh, R.; Gopalan, V. Polar and magnetic properties of PbVO 3 thin films. Phys. Rev. B 2007 , 75 , 060101.
- (23) Chi, M.; Mizoguchi, T.; Martin, L.; Ramesh, R.; Browning, N. Studying Thinfilm Perovskite Vanadates by EELS. Microsc. Microanal. 2007 , 13 , 1234-1235.
- (24) Solovyev, I. Magnetic structure of the noncentrosymmetric perovskites PbVO 3 and BiCoO : Theoretical analysis. 3 Phys. Rev. B 2012 , 85 , 054420.
- (25) Parveen, A.; Srivastava, A.; Gaur, N. Thermal properties of a novel compound PbVO . 3 J. Phys. Conf. Ser 2012 , 377 , 012059.
- (26) Ming, X.; Yin, J.-W.; Wang, X.-L.; Wang, C.-Z.; Huang, Z.-F.; Chen, G. First-principles comparative study of multiferroic compound PbVO . 3 Solid State Sci. 2010 , 12 , 938-945.

- (27) Zhou, W.; Tan, D.; Xiao, W.; Song, M.; Chen, M.; Xiong, X.; Xu, J. Structural properties of PbVO 3 perovskites under hydrostatic pressure conditions up to 10.6 GPa. J. Phys.: Condens. Matter 2012 , 24 , 435403.
- (28) Miloˇ sevi´, c A. S.; Lali´, c M. V.; Popovi´ c, Z. S.; Vukajlovi´, c F. R. An ab initio study of electronic structure and optical properties of multiferroic perovskites PbVO 3 and BiCoO . 3 Opt. Mater. 2013 , 35 , 1765-1771.
- (29) Ming, X.; Hu, F.; Du, F.; Wei, Y.-J.; Chen, G. First-principles study of multiferroic material PbVO 3 under uniaxial pressure. Eur. Phys. J. B. 2015 , 88 , 212.
- (30) Kadiri, A.; Arejdal, M.; Abbassi, A.; El Amraoui, Y.; Ez-Zahraouy, H.; Benyoussef, A. Calculated Magnetic Properties of the Compound PbVO . 3 J. Supercond. Novel Magn. 2017 , 30 , 2247-2253.
- (31) Oka, K.; Yamauchi, T.; Kanungo, S.; Shimazu, T.; Oh-ishi, K.; Uwatoko, Y.; Azuma, M.; Saha-Dasgupta, T. Experimental and Theoretical Studies of the Metallic Conductivity in Cubic PbVO 3 under High Pressure. J. Phys. Soc. Jpn. 2018 , 87 , 024801.
- (32) Tsuchiya, T.; Katsumata, T.; Ohba, T.; Inaguma, Y. High-pressure synthesis and characterization of a novel perovskite PbFe 1 2 / V 1 2 / O . 3 J. Ceram. Soc. Jpn. 2009 , 117 , 102-105.
- (33) Ar´ evalo-L´pez, o A. M.; AlarioFranco, M. A. Structural Percolation in the PbMO 1 -x MO O x 3 ( M M , = Ti, Cr, and V) Perovskites. Inorg. Chem. 2011 , 50 , 7136-7141.
- (34) Ricci, F.; Fiorentini, V. Multiferroicity in V-doped PbTiO . 3 2013 , 470 , 012013.
- (35) Pan, Z. et al. Colossal Volume Contraction in Strong Polar Perovskites of Pb(Ti,V)O . 3 J. Am. Chem. Soc 2017 , 139 , 14865-14868.
- (36) Kresse, G.; Furthm¨ uller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 1996 , 54 , 11169-11186.
- (37) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996 , 77 , 3865.
- (38) Ernzerhof, M.; Scuseria, G. E. Assessment of the Perdew-Burke-Ernzerhof exchangecorrelation functional. J. Chem. Phys 1999 , 110 , 5029-5036.
- (39) Ravindran, P.; Vidya, R.; Kjekshus, A.; Fjellv˚g, H.; Eriksson, O. Theoretical ina vestigation of magnetoelectric behavior in BiFeO . 3 Phys. Rev. B 2006 , 74 , 224412.
- (40) Monkhorst, H. J.; Pack, J. D. Special points for Brillouin-zone integrations. Phys. Rev. B 1976 , 13 , 5188.
- (41) Patra, L.; Kishore, M. A.; Vidya, R.; Sj˚stad, A. O.; Fjellv˚g, H.; Ravindran, P. a a Electronic and Magnetic Structures of Hole Doped Trilayer La 4 -x Sr Ni O x 3 8 from First-Principles Calculations. Inorg. Chem. 2016 , 55 , 11898-11907.
- (42) Resta, R. Macroscopic polarization in crystalline dielectrics: the geometric phase approach. Rev. Mod. Phys 1994 , 66 , 899.
- (43) King-Smith, R.; Vanderbilt, D. Theory of polarization of crystalline solids. Phys. Rev. B 1993 , 47 , 1651.
- (44) Blaha, P.; Schwarz, K.; Madsen, G. K.; Kvasnicka, D.; Luitz, J. WIEN2k: An Augmented Plane Wave plus Local Orbitals Program for Calculating Crystal Properties, edited by K. Schwarz. (Vienna University of Technology, Austria) 2001 ,
- (45) Erskine, J.; Stern, E. Calculation of the m 23 magneto-optical absorption spectrum of ferromagnetic nickel. Phys. Rev. B 1975 , 12 , 5016.

- (46) Carra, P.; Thole, B.; Altarelli, M.; Wang, X. X-ray circular dichroism and local magnetic fields. Phys. Rev. Lett. 1993 , 70 , 694.
- (47) Thole, B.; Carra, P.; Sette, F.; van der Laan, G. X-ray circular dichroism as a probe of orbital magnetization. Phys. Rev. Lett. 1992 , 68 , 1943.
- (48) Nekrasov, I.; Keller, G.; Kondakov, D.; Kozhevnikov, A.; Pruschke, T.; Held, K.; Vollhardt, D.; Anisimov, V. Comparative study of correlation effects in CaVO 3 and SrVO . 3 Phys. Rev. B 2005 , 72 , 155106.
- (49) Schindler, M.; Hawthorne, F.; Baur, W. Crystal chemical aspects of vanadium: polyhedral geometries, characteristic bond valences, and polymerization of (VO ) n polyhedra. Chem. Mater 2000 , 12 , 12481259.
- (50) Shannon, R. D. Revised effective ionic radii and systematic studies of interatomic distances in halides and chalcogenides. Jpn. J. Appl. Phys 1976 , 32 , 751-767.
- (51) Choi, S. M.; Stringer, C. J.; Shrout, T. R.; Randall, C. A. Structure and property investigation of a Bi-based perovskite solid solution:(1x )Bi(Ni 1 2 / Ti 1 2 / )O 3 x PbTiO . 3 J. Appl. Phys 2005 , 98 , 034108.
- (52) Randall, C.; Eitel, R.; Jones, B.; Shrout, T.; Woodward, D.; Reaney, I. Investigation of a high T c piezoelectric system:(1x )Bi(Mg 1 2 / Ti 1 2 / )O 3 ( x )PbTiO . 3 J. Appl. Phys 2004 , 95 , 3633-3639.
- (53) Chen, J.; Wang, F.; Huang, Q.; Hu, L.; Song, X.; Deng, J.; Yu, R.; Xing, X. Effectively control negative thermal expansion of single-phase ferroelectrics of PbTiO 3 (Bi,La)FeO 3 over a giant range. Sci. Rep 2013 , 3 , 2458.
- (54) Eitel, R. E.; Randall, C. A.; Shrout, T. R.; Rehrig, P. W.; Hackenberger, W.; Park, S.-E. New high temperature morphotropic phase boundary piezoelectrics
- based on Bi(Me)O -PbTiO 3 3 ceramics. Jpn. J. Appl. Phys 2001 , 40 , 5999-6002.
- (55) Suchomel, M. R.; Davies, P. K. Enhanced tetragonality in ( x )PbTiO -(13 x )Bi(Zn 1 2 / Ti 1 2 / )O 3 and related solid solution systems. Appl. Phys. Lett 2005 , 86 , 262905.
- (56) Cohen, R. E. Origin of ferroelectricity in perovskite oxides. Nature 1992 , 358 , 136.
- (57) Kuroiwa, Y.; Aoyagi, S.; Sawada, A.; Harada, J.; Nishibori, E.; Takata, M.; Sakata, M. Evidence for Pb-O covalency in tetragonal PbTiO . 3 Phys. Rev. Lett. 2001 , 87 , 217601.
- (58) Goodenough, J. B.; Zhou, J. Varied roles of Pb in transition-metal PbMO 3 perovskites (M= Ti, V, Cr, Mn, Fe, Ni, Ru). Sci. Technol. Adv. Mater 2015 , 16 , 036003.
- (59) Vidya, R.; Ravindran, P.; Knizek, K.; Kjekshus, A.; Fjellv˚g, H. Density Funca tional Theory Studies of Spin, Charge, and Orbital Ordering in YBaT O 2 2 (T= Mn, Fe, Co). Inorg. Chem. 2008 , 47 , 6608-6620.
- (60) Vidya, R.; Ravindran, P.; Kjekshus, A.; Fjellv˚g, a H. Spin, charge, and orbital ordering in the ferrimagnetic insulator YBaMn O . 2 2 Phys. Rev. B 2002 , 65 , 144422.
- (61) Ravindran, P.; Vidya, R.; Eriksson, O.; Fjellv˚g, H. Magnetic-Instability-Induced a Giant Magnetoelectric Coupling. Adv. Mater 2008 , 20 , 1353-1356.
- (62) Brooks, H. Ferromagnetic anisotropy and the itinerant electron model. Phys. Rev. 1940 , 58 , 909.
- (63) Bozorth, R. Directional ferromagnetic properties of metals. J. Appl. Phys. 1937 , 8 , 575-588.
- (64) Scherz, A.; Wende, H.; Baberschke, K. Fine structure of X-ray magnetic circular

dichroism for early 3 d transition metals. Appl. Phys. A 2004 , 78 , 843-846.

- (65) Wende, H. Recent advances in X-ray absorption spectroscopy. Rep. Prog. Phys 2004 , 67 , 2105.
- (66) Scherz, A.; Wende, H.; Sorg, C.; Baberschke, K.; Minr, J.; Benea, D.; Ebert, H. Limitations of integral XMCD sum-rules for the early 3d elements. Phys. Scr 2005 , 2005 , 586.
- (67) Wende, H.; Scherz, A.; Sorg, C.; Baberschke, K.; Gross, E.; Appel, H.; Burke, K.; Min´ ar, J.; Ebert, H.; Ankudinov, A.; Rehr, J. XMCD analysis beyond standard procedures. AIP Conf. Proc 2007 , 882 , 78-82.
Lokanath Patra Vidya Ravindran Helmer Fjellvag Ponniah Ravindran
Giant magnetoelectric coupling in multiferroic PbTi$_{1-x}$V$_x$O$_{3}$ from density functional calculations
The giant magnetoelectric coupling is a very rare phenomenon which has gained
a lot of attention for the past few decades because of fundamental interest as
well as practical applications. Here, we have successfully achieved the giant
magnetoelectric coupling in PbTi1-xVxO3 (x= 0-1) with the help of a series of
generalized-gradient-corrected (GGA), GGA including on-site coulomb repulsion
(U) corrected spin polarized calculations based on accurate density functional
theory. Our total energy calculations show that PbTi1-xVxO3 stabilizes in
C-type antiferromagnetic ground state for x>0.123. With the substitution of V
into PbTiO3, the tetragonal distortion is highly enhanced accompanied by a
linear increase in polarization. In addition, our band structure analysis shows
that for lower x values, the tendency to form 2Dmagnetism of PbTi1-xVxO3
decreases. A non-magnetic metallic ground state is observed for the
paraelectric phase for V concentration (x) = 1 competing with a volume change
of 10% showing a large magnetovolume effect. Our orbital projected DOS as well
as orbital ordering analysis suggest that the orbital ordering plays a major
role in the magnetic to non-magnetic transition when going from ferroelectric
to paraelectric phase. The calculated magnetic anisotropic energy shows that
the direction [110] is the easy axis of magnetization for x= 1 composition. The
present study adds a new series of compounds to the magnetoelectric family with
rarely existing giant coupling between electric and magnetic order parameters.
These results show that such kind of materials can be used for novel practical
applications where one can change the magnetic properties drastically (magnetic
to non-magnetic as shown here) with external electric fieldand vice-versa.