# Theory of PbTiO3, BaTiO3, and SrTiO3 Surfaces

## I Introduction

The surfaces of insulating cubic perovskite materials such as PbTiO<sub>3</sub>, BaTiO<sub>3</sub>, and SrTiO<sub>3</sub> are of interest from several points of view. First, some of these materials (notably SrTiO<sub>3</sub>) are very widely used as substrates for growth of other oxide materials (e.g., layered high-$`T_c`$ superconductors and “colossal magnetoresistance” materials). Second, this class of materials is of enormous importance for actual and potential applications that make use of their unusual piezoelectric, ferroelectric, and dielectric properties (e.g., for piezoelectric transducers, non-volatile memories, and wireless communications applications, respectively). Many of these applications are increasingly oriented towards thin-film geometries, where surface properties are of growing importance. Third, the bulk materials display a variety of structural phase transitions; the ferroelectric (FE) structural phases are of special interest, but antiferroelectric (AFE) or antiferrodistortive (AFD) transitions can also take place. It is then of considerable fundamental interest to consider how these structural distortions couple to the surface, e.g., whether the presence of the surface acts to enhance or suppress the structural distortion. The ferroelectric properties are well known to degrade in thin-film and particulate geometries, and it is very important to understand whether such behavior is intrinsic to the presence of a surface, or whether it arises from extrinsic factors such as compositional non-uniformities or structural defects in the surface region. Finally, the cubic perovskites can serve as model systems for the study of transition-metal oxide surfaces more generally.

In the last decade, there has been a surge of activity in the application of first-principles computational methods based on density-functional theory (DFT) to the study of the bulk properties, and especially the ferroelectric transitions, in bulk perovskite oxides. (For a recent review, see Ref.  or .) The importance of these methods was recently underlined by the award of the Nobel Prize in Chemistry to Walter Kohn, the primary originator of DFT. In the materials theory community, these methods have been widely used for two decades to predict properties of semiconductors and simple metals. However, recent advances in computational algorithms and computer power now allow these methods to be applied to more complex materials (e.g., perovskites) and more complex geometries (e.g., defects and surfaces). In particular, pioneering studies of BaTiO<sub>3</sub>  and SrTiO<sub>3</sub>  surfaces have recently appeared.

Experimental investigations of the surface structure of cubic perovskites have not been very extensive. Such studies are hindered by the difficulties of preparing clean and defect-free surfaces, and of overcoming charging effects associated with many experimental probes. Even for SrTiO<sub>3</sub>, the best-studied of these surfaces, there is a disappointing level of agreement among experimental results  and between experiment and theory. We are not aware of comparable studies of BaTiO<sub>3</sub> and PbTiO<sub>3</sub> surfaces.

The purpose of the present contribution is to present new theoretical work on the structural properties of the PbTiO<sub>3</sub> (001) surface, and to compare and contrast these results with the previous work of our group on BaTiO<sub>3</sub> and SrTiO<sub>3</sub> surfaces. As regards bulk properties, lead-based compounds such as PbTiO<sub>3</sub> and PbZrO<sub>3</sub> are known to behave quite differently from alkaline-earth based perovskites such as BaTiO<sub>3</sub> and SrTiO<sub>3</sub>. Previous theoretical work has shown that the FE distortion is typically larger and that Pb atoms participate much more strongly in (and sometimes even dominate) the FE distortion, compared with non-Pb perovskites. Moreover, the Pb-based compounds are generally more susceptible to more complex AFD and AFE instabilities involving tilting of the oxygen octahedra, and the ground-state structures often involve the formation of some quite short Pb–O bonds. All of these effects point to a strong and active involvement of the Pb atoms in the bonding, most naturally interpreted in terms of the formation of partially covalent Pb–O bonds with the closest oxygen neighbors. Finally, a focus on Pb-based materials is motivated by the fact that these are the leading candidates for many practical piezoelectric and switching applications, especially in the form of solid solutions such as PZT (PbZr<sub>x</sub>Ti<sub>1-x</sub>O<sub>3</sub>), PMN (PbMg<sub>1/3</sub>Nb<sub>2/3</sub>O<sub>3</sub>), and PZN (PbZn<sub>1/3</sub>Nb<sub>2/3</sub>O<sub>3</sub>).

The manuscript is organized as follows. Section II contains a brief account of the technical details of the work, including the theoretical methods used, the slab geometries studied, and the formulation of the surface energy. In Sec. III we present the computed structural relaxations of the PbTiO<sub>3</sub> surfaces, and compare these to the previous results on BaTiO<sub>3</sub> and SrTiO<sub>3</sub> surfaces. Additionally, we discuss the surface energetics (surface energies and surface relaxation energies), and point out some characteristic differences in the surface electronic structure of the three compounds. Finally, the paper ends with a summary in Sec. IV.

## II Preliminaries

### A Theoretical Methods

We carried out self-consistent plane-wave pseudopotential calculations within Kohn-Sham density-functional theory using a conjugate-gradient technique. Exchange and correlation were treated using the Ceperley-Alder form. Vanderbilt ultrasoft pseudopotentials were employed, with semicore Pb $`5d`$, Ba $`5s`$ and $`5p`$, Sr $`4s`$ and $`4p`$, and Ti $`3s`$ and $`3p`$ orbitals included as valence states. A plane-wave cutoff of 25 Ry has been used throughout. Relaxations of atomic coordinates are iterated until the forces are less than 0.01 eV/Å. Justification of the convergence and accuracy of this approach can be found in the previously published work.

### B Surface and Slab Geometries

In this work we consider only II-IV cubic perovskites, i.e., ABO<sub>3</sub> perovskites in which atoms A and B are divalent and tetravalent, respectively. In this case, two non-polar (001) surface terminations are possible: the AO–terminated surface, and the BO<sub>2</sub>–terminated surface.

We have studied both types of surface termination for all three materials (PbTiO<sub>3</sub>, BaTiO<sub>3</sub>, and SrTiO<sub>3</sub>) using a repeated slab geometry. The slabs are symmetrically terminated and typically contain seven layers (17 or 18 atoms), as illustrated in Fig. 1. The vacuum region was chosen to be two lattice constants thick. The calculations were done using a (4,4,2) Monkhorst-Pack mesh, corresponding to three or four k-points in the irreducible Brillouin zone for cubic and tetragonal surfaces respectively. The convergence of the calculations has been very carefully checked for PbTiO<sub>3</sub> by repeating some of the calculations with asymmetrically terminated eight-layer slabs and symmetrically terminated nine-layer slabs. Additionally, we have enlarged the vacuum region to a thickness of three lattice constants, and we have checked the convergence of the Brillouin zone integration by going to a (6,6,2) k-point mesh. In all cases, the results for the structural properties of the surfaces given in the Tables I to V change by less than 0.2%.

For all three materials, we first computed the relaxations for the “cubic” surface, i.e., for the case where there is no symmetry lowering relative to a slab of ideal cubic material. In this case we preserved $`M_x`$, $`M_y`$, and $`M_z`$ mirror symmetries relative to the center of the slab, and set the lattice constants in the $`\widehat{x}`$ and $`\widehat{y}`$ directions equal to those computed theoretically for the corresponding bulk material (3.89 Å, 3.95 Å, and 3.86 Å for PbTiO<sub>3</sub>, BaTiO<sub>3</sub>, and SrTiO<sub>3</sub>, respectively). The symmetry-allowed displacements of the atoms in the $`z`$ (surface-normal) direction were then fully relaxed.

Each of the three materials studied displays a different sequence of structural phase transitions from the cubic paraelectric phase as the temperature is lowered. PbTiO<sub>3</sub> undergoes a single transition into a tetragonal ferroelectric (FE) phase at 763 K and then remains in this structure down to zero temperature. BaTiO<sub>3</sub> displays a series of three transitions to tetragonal, orthorhombic, and rhombohedral FE phases at 403 K, 278 K, and 183 K, respectively. SrTiO<sub>3</sub> remains cubic down to 105 K, at which point it undergoes an antiferrodistortive transition involving rotation of the oxygen octahedra and doubling of the unit cell. The material nearly goes ferroelectric at about $`T=30`$ K, but is evidently prevented from doing so by quantum zero-point fluctuations.

Because we are primarily interested in the room-temperature structures of these materials and their surfaces, we have chosen to focus on the tetragonal FE phases of PbTiO<sub>3</sub> and BaTiO<sub>3</sub> for our surface studies. We consider only the case of the tetragonal $`c`$ axis (i.e., polarization) lying parallel to the surface, since polarization normal to the surface is strongly suppressed by the depolarization fields that would arise from the accumulated charge at the surfaces. We take the tetragonal axis to lie along $`\widehat{x}`$, and relax the $`M_x`$ symmetry while retaining the $`M_y`$ and $`M_z`$ symmetries with respect to the center of the slab. For PbTiO<sub>3</sub>, which is tetragonal at $`T=0`$, this will indeed be the ground-state structure of the slab. For BaTiO<sub>3</sub>, on the other hand, the $`M_y`$ symmetry is artificially imposed so that the theoretical $`T=0`$ calculation will mimic the experimental room-temperature surface structure. In both cases, the slab lattice constants in the $`\widehat{x}`$ and $`\widehat{y}`$ directions were set equal to the corresponding theoretical equilibrium lattice constants computed for the bulk tetragonal phase: $`c`$=4.04 Å and $`a`$=3.86 Å for PbTiO<sub>3</sub>, and $`c`$=3.99 Å and $`a`$=3.94 Å for BaTiO<sub>3</sub>.

### C Surface Energies

A comparison of the relative stability of the AO and TiO<sub>2</sub> surface terminations is problematic because the corresponding surface slabs contain different numbers of AO and TiO<sub>2</sub> formula subunits. We treat this problem by introducing chemical potentials $`\mu _{\mathrm{AO}}`$ and $`\mu _{\mathrm{TiO}_2}`$ for these subunits, defined in such a way that $`\mu _{\mathrm{AO}}=0`$ and $`\mu _{\mathrm{TiO}_2}=0`$ correspond to thermal equilibrium with bulk crystalline AO and TiO<sub>2</sub>, respectively. We have computed the cohesive energies $`E_{\mathrm{AO}}`$ and $`E_{\mathrm{TiO}_2}`$ of crystalline PbO, BaO, SrO, and TiO<sub>2</sub> using the same first-principles pseudopotential method in order to provide these reference values. The grand potential for a given surface structure can then be computed as

$$F_{\mathrm{surf}}=\frac{1}{2}[E_{\mathrm{slab}}N_{\mathrm{TiO}_2}(E_{\mathrm{TiO}_2}+\mu _{\mathrm{TiO}_2})N_{\mathrm{AO}}(E_{\mathrm{AO}}+\mu _{\mathrm{AO}})],$$
(1)   
where $`N`$ is the number of formula subunits contained in the slab, and the factor of $`1/2`$ accounts for the fact that each slab contains two surfaces. Assuming that the surface of the ATiO<sub>3</sub> is in equilibrium with its own bulk, it follows that

$$\mu _{\mathrm{AO}}+\mu _{\mathrm{TiO}_2}=E_\mathrm{f},$$
(2)   
where $`E_\mathrm{f}`$ is the heat of formation of bulk ATiO<sub>3</sub> from bulk AO and bulk TiO<sub>2</sub>. The two chemical potentials are thus not independent, and we choose to treat $`\mu _{\mathrm{TiO}_2}`$ as the independent variable when presenting our results. Accordingly, $`\mu _{\mathrm{TiO}_2}`$ is allowed to vary over the range

$$E_\mathrm{f}\mu _{\mathrm{TiO}_2}0,$$
(3)   
the lower and upper limit corresponding to the precipitation of particulates of AO and TiO<sub>2</sub> on the surface, respectively.

## III Results and Discussions

### A Structural relaxations

We begin by presenting our new results on the structural properties of the PbTiO<sub>3</sub> (001) surfaces. The equilibrium atomic positions for both surface terminations in the two phases were obtained by starting from the ideal structures of the surfaces and then relaxing the atomic positions while preserving the symmetries described in section II B. The results for the fully relaxed geometries are summarized in Table I and II. By symmetry, there are no forces along $`\widehat{x}`$ and $`\widehat{y}`$ for the cubic surface, and no forces along $`\widehat{y}`$ for the tetragonal surface.

Tables I and II show for both surfaces a substantial inward contraction towards the bulk for the uppermost surface layers, whereas for the second layers we find an outward relaxation of the atoms relative to the positions of the atoms on the ideal surface. Generally, the metal and the oxygen atoms move in the same direction, but the relaxations of the metal atoms are much larger, leading to a rumpling of the layers. The single exception is the surface layer of the tetragonal TiO<sub>2</sub>–terminated surface, where one of the two oxygen atoms moves in the opposite direction to the metal atom. Therefore we can see here a significant asymmetry between the O atoms with respect to their positions perpendicular to the surface. This asymmetry between the oxygen atoms in the topmost surface layer of the tetragonal TiO<sub>2</sub>–terminated surface was also found for BaTiO<sub>3</sub> but with a much smaller amplitude. As expected, we find the largest relaxations for the surface–layer atoms, but the displacement of the Pb atom in the second layer of the TiO<sub>2</sub>–terminated surface is of the same magnitude.

In order to compare these results with previous calculations for SrTiO<sub>3</sub> and BaTiO<sub>3</sub>, we have calculated the changes in the interlayer distances $`\mathrm{\Delta }d_{ij}`$ and the amplitudes of the rumpling $`\eta _i`$ of the layers in the surface slabs for all three perovskites. The results for both surface terminations and the different phases are given in the Tables III and IV. We denote the change in the $`z`$ position of a

metal atom relative to the ideal unrelaxed structure as $`\delta _z(\mathrm{M})`$, and $`\delta _z(\mathrm{O})`$ is the same for the oxygen atom in the same layer (defined as $`[\delta _z(\mathrm{O}_\mathrm{I})+\delta _z(\mathrm{O}_{\mathrm{II}})]/2`$ for a TiO<sub>2</sub> layer). We then define the change of the interlayer distance $`\mathrm{\Delta }d_{ij}`$ as the difference between the averaged atomic displacements $`[\delta _z(\mathrm{M})+\delta _z(\mathrm{O})]/2`$ of layer $`i`$ and $`j`$, and the rumpling $`\eta _i`$ is defined as the amplitude of these displacements $`|\delta _z(\mathrm{M})\delta _z(\mathrm{O})|`$. From Tables III and IV we can see that, for all three perovskites and for both terminations, the surfaces display a similar oscillating relaxation pattern with a reduction of the interlayer distance $`d_{12}`$, an expansion of $`d_{23}`$ and again a reduction for $`d_{34}`$. However, compared to BaTiO<sub>3</sub> and SrTiO<sub>3</sub>, the amplitudes of the relaxations in PbTiO<sub>3</sub> are significantly increased.

The second interesting feature of Tables III and IV is that for BaTiO<sub>3</sub>, there is almost no difference in the relaxations of the surface layers between the cubic and the tetragonal phase. The same is true for the TiO<sub>2</sub>–terminated surface of PbTiO<sub>3</sub>. For the PbO–terminated surface, in contrast, the changes in the interlayer distances and the layer rumplings are strongly reduced in the tetragonal phase. We will come back to this point at the end of the next subsection.

### B Influence of the surface upon ferroelectricity

We turn now to the question of whether the presence of the surface has a strong effect upon the near–surface ferroelectricity. To analyze whether the ferroelectric order is enhanced or suppressed near the surface, we introduce average ferroelectric distortions $`\delta _{\mathrm{FE}}`$ for each layer of the surface slabs:

$$\begin{array}{cccc}\hfill \delta _{\mathrm{FE}}& =& |\delta _x(\mathrm{A})\delta _x(\mathrm{O}_{\mathrm{III}})|\hfill & \text{for AO planes and}\hfill \\ \hfill \delta _{\mathrm{FE}}& =& |\delta _x(\mathrm{Ti})[\delta _x(\mathrm{O}_\mathrm{I})+\delta _x(\mathrm{O}_{\mathrm{II}})]/2|\hfill & \text{for TiO}\text{2}\text{ planes.}\hfill \end{array}$$
(4)   
The calculated values of $`\delta _{\mathrm{FE}}`$ for BaTiO<sub>3</sub> and PbTiO<sub>3</sub> are given in Table V; the last row of the table gives the bulk values for reference.

For the PbO–terminated surface of PbTiO<sub>3</sub>, one can see a clear increase in the average ferroelectric distortions $`\delta _{\mathrm{FE}}`$ when going from the bulk values to the surface layer. On the other hand, for the TiO<sub>2</sub>–terminated surface, the average distortions are slightly decreased at the surface. Surprisingly, this is just the opposite of what one observes for BaTiO<sub>3</sub>, where one sees a reduction of the ferroelectric distortions for the BaO–terminated surface and a moderate enhancement for the TiO<sub>2</sub>–terminated surface. (Of course, the distortions are also much smaller for BaTiO<sub>3</sub> surfaces, as they are in the bulk, compared to PbTiO<sub>3</sub>.) These results tend to confirm that Pb is a much more active constituent in PbTiO<sub>3</sub> than is Ba in BaTiO<sub>3</sub>, presumably because of the partially covalent nature of the Pb–O bonds as discussed in Sec. I.

In any case, the present results again confirm that the presence of the surface does not lead to any drastic suppression of the ferroelectric order near the surface, supporting the view that extrinsic effects must be responsible for degradation of ferroelectricity in thin-film geometries.

Finally, we note that there are interesting signs of interplay between the relaxations parallel and perpendicular to the surface for PbTiO<sub>3</sub>. In particular, the relaxations perpendicular to the surface are substantially reduced (by a factor of $``$3) on the PbO-terminated surface when going from the cubic to the tetragonal case. This can be rationalized as follows. Because of the partial covalency of the Pb–O bonds, there is a tendency to reduce the Pb–O bond length (this length is 2.75, 2.51, and 2.30 Å in cubic PbTiO<sub>3</sub>, tetragonal PbTiO<sub>3</sub>, and PbO, respectively). For the cubic surface, by symmetry, the only possibility to shorten this bond length is by a strong movement of the Pb atom towards the bulk and a strong movement upwards of the O atoms in the second layer. This leads to the strong rumpling and the decrease of $`d_{12}`$. But in the tetragonal phase there is also the possibility to enlarge the ferroelectric distortion in order to shorten the Pb–O bond length. Evidently, the enlargement of the ferroelectric distortion is preferred to the relaxation perpendicular to the surface.

### C Surface energies

In this section we discuss the surface energetics of the three perovskite compounds. In order to compare the relative stability of the AO– and TiO<sub>2</sub>–terminated surfaces, we have calculated the grand thermodynamic potential $`F_{\mathrm{surf}}`$ (as introduced in Sec. II C) for the different surfaces as a function of the chemical potential $`\mu _{\mathrm{TiO}_2}`$. The results for the tetragonal surfaces of BaTiO<sub>3</sub> and PbTiO<sub>3</sub> are shown in Fig. 3. The graphs of the grand thermodynamic potentials for the SrTiO<sub>3</sub> surfaces are very similar to those of BaTiO<sub>3</sub> and are therefore not shown separately.

Figure 3 shows a very different behavior for the BaTiO<sub>3</sub> and PbTiO<sub>3</sub> surfaces. First of all, the formation energy $`E_\mathrm{f}`$ of PbTiO<sub>3</sub> (when formed from bulk PbO and TiO<sub>2</sub>) is 0.36 eV, much lower than the formation energies of SrTiO<sub>3</sub> and BaTiO<sub>3</sub> which are about 3.2 eV. This leads to a much smaller range for the chemical potential $`\mu _{\mathrm{TiO}_2}`$ for which PbTiO<sub>3</sub> surfaces can grow in thermodynamic equilibrium. Second, for BaTiO<sub>3</sub> the two different surfaces have a comparable range of thermodynamic stability, indicating that either BaO–terminated surfaces or TiO<sub>2</sub>–terminated surfaces could be formed depending on whether growth occurs in Ba–rich or Ti–rich conditions. In contrast, for PbTiO<sub>3</sub> only the PbO–terminated surface can be obtained in thermodynamic equilibrium.

To get a quantity describing the surface energetics that is independent of the chemical potential $`\mu _{\mathrm{TiO}_2}`$ and therefore allows a more direct comparison of the three compounds, we define the average surface energy per surface unit cell

$$E_{\mathrm{surf}}=\frac{1}{4}\left(E_{\mathrm{slab}}^{\mathrm{AO}}+E_{\mathrm{slab}}^{\mathrm{TiO}_2}7E_{\mathrm{bulk}}\right),$$
(5)   
which is equal to the average of the grand thermodynamic potential $`F_{\mathrm{surf}}`$ for the two kinds of surfaces. Again, the results for $`E_{\mathrm{surf}}`$ shown in Table VI are very similar for SrTiO<sub>3</sub> and BaTiO<sub>3</sub>, whereas the value for PbTiO<sub>3</sub> is significantly lower.

Finally we have computed the average relaxation energy $`E_{\mathrm{relax}}`$ of the three perovskite compounds. $`E_{\mathrm{relax}}`$ is defined as the difference between the average surface energy $`E_{\mathrm{surf}}`$ of the ideal surface without relaxation of the atoms, and the fully relaxed surfaces. The largest and smallest value for $`E_{\mathrm{relax}}`$ (see Table VI) were found for PbTiO<sub>3</sub> and BaTiO<sub>3</sub>, respectively, which is in agreement with the observation that the atomic relaxations are largest in PbTiO<sub>3</sub> and smallest in BaTiO<sub>3</sub>.

For all three compounds the average relaxation energy $`E_{\mathrm{relax}}`$ is many times larger than a typical bulk ferroelectric well depth, which is approximately 0.03 eV for BaTiO<sub>3</sub> and 0.05 eV for PbTiO<sub>3</sub>. This would indicate that the surface is capable of acting as a strong perturbation on the ferroelectric order. As we have shown in Sec. III B, this is not the case for BaTiO<sub>3</sub> and PbTiO<sub>3</sub>. One reason why the ferroelectric order is not as strongly affected by the surface as one might have thought has been pointed out in Ref. : the soft phonon eigenmode, which is responsible for the ferroelectric distortion, is only one of three zone center modes having the same symmetry. By looking at how strongly the surface relaxations are related to each of these zone center modes it has turned out that the distortions induced by the presence of the surface are to a large extent of non–ferroelectric character.

### D Surface band structure

For all three perovskite compounds we have carried out LDA calculations of the bulk and the surface electronic structure for our various surface slabs. It is well known that the LDA is quantitatively unreliable regarding excitation properties such as band gaps. Since we are in the following only looking at differences between band structures, we think that our conclusions drawn from the LDA results are nevertheless qualitatively correct.

As has already been shown in Ref. , the bulk band structures of SrTiO<sub>3</sub> and BaTiO<sub>3</sub> are very similar, whereas PbTiO<sub>3</sub> shows some significant differences. In SrTiO<sub>3</sub> and BaTiO<sub>3</sub> the upper edge of the valence band is very flat throughout the Brillouin zone. On the other hand, in PbTiO<sub>3</sub> the shallow $`6s`$ semicore states of the Pb atoms hybridize with the $`2p`$ states of the O atoms, leading to a lifting of the upper valence band states near the X point of the Brillouin zone.

This fact is responsible for a different behavior of the PbTiO<sub>3</sub> surface band structure compared to SrTiO<sub>3</sub> and BaTiO<sub>3</sub>. If we look at the calculated band gaps in Table VII, we see that for TiO<sub>2</sub>–terminated surfaces the band gap is significantly reduced for SrTiO<sub>3</sub> and BaTiO<sub>3</sub>, whereas for PbTiO<sub>3</sub> the band gap is almost unchanged. The reduction of the band gap in SrTiO<sub>3</sub> and BaTiO<sub>3</sub> is mainly due to an upward intrusion of the upper valence band states near the M point into the lower part of the band gap (as pointed out in Ref. , this is caused by the suppression of the hybridization of certain O $`2p`$ and Ti $`3d`$ orbitals in the surface layer). In PbTiO<sub>3</sub> we find the same upward movement of the upper valence band states near the M point, but these states stay just below the highest valence states at the X point, and so the band gap is almost unchanged.

On the other hand, for the AO–terminated surfaces we see no reduction of the band gap for any of the three perovskite compounds. Even here, however, there is a subtle difference between PbTiO<sub>3</sub> and the other materials, this time concerning the conduction band edge. According to our calculations, the Pb $`6p`$ states overlap the Ti $`3d`$ states to some degree in bulk PbTiO<sub>3</sub>, and this effect is accentuated at the $`\mathrm{\Gamma }`$ point of the surface Brillouin zone on the Pb–O terminated surface, where the lowest Pb $`6p`$ state falls just below the lowest Ti $`3d`$ state. We thus suggest that the conduction band minimum may actually have Pb $`6p`$ character at this surface, although the effect is too small to affect the band gaps in Table VII substantially. This might be an interesting target of investigation for future spectroscopic experimental studies.

## IV Summary

In summary, we have calculated structural and electronic properties of PbTiO<sub>3</sub> (001) surfaces using a first-principles density-functional approach. The results are compared and contrasted with corresponding previous calculations on BaTiO<sub>3</sub> and SrTiO<sub>3</sub> surfaces. We observe qualitatively different behavior of the PbTiO<sub>3</sub> surfaces in several respects. First, within the narrow range of PbO and TiO<sub>2</sub> chemical potentials permitted by bulk thermodynamics, we find that the TiO<sub>2</sub>-terminated surface is never thermodynamically stable. Thus, the PbO-terminated surface is expected to be the one observed experimentally. Second, the interaction between the ferroelectric distortion and the presence of the surface is quite different for PbTiO<sub>3</sub>, compared to BaTiO<sub>3</sub>. In particular, the ferroelectricity is strongly enhanced at the AO-terminated surface and suppressed at the TiO<sub>2</sub>-terminated surface, just the opposite of the behavior found for BaTiO<sub>3</sub>. Moreover, the ferroelectric distortion at the surface allows for a drastic reduction of the rumpling of the surface layer on the PbO-terminated surface, an effect which is not seen on the BaO-terminated of BaTiO<sub>3</sub>. Third, the surface electronic band structure is qualitatively modified in the case of PbTiO<sub>3</sub> by the presence of Pb $`6s`$ and $`6p`$ states in the upper valence and lower conduction regions.

###### Acknowledgements.

This work was supported by the ONR grant N00014-97-1-0048.