## Tuning photostriction in (PbTiO ) /(SrTiO ) 3 n 3 m superlattices via chemical composition: An ab-initio study

Carmel Dansou , Charles Paillard 1 1,2 ,Laurent Bellaiche 1,3

1 Smart Ferroic Materials Center, Institute for Nanoscience &amp; Engineering and Department of Physics, University of Arkansas, Fayetteville, Arkansas 72701, USA 2 Universit´ e Paris-Saclay, CentraleSup´lec, e CNRS, Laboratoire SPMS, 91190 Gif-sur-Yvette, France. and 3 Department of Materials Science and Engineering, Tel Aviv University, Ramat Aviv, Tel Aviv 6997801, Israel. (Dated: August 2, 2024)

Light-induced mechanical deformations in single-domain (PbTiO ) 3 n /(SrTiO ) 3 m superlattices were simulated using first-principle calculations. By varying the chemical fraction of PbTiO , we 3 discover that these heterostructures' photostrictive behavior can be tuned quantitatively and qualitatively. Additionally, we present simple analytical models to explain the calculated deformations and predict a critical PbTiO 3 fraction with no photostriction. In addition to the report in [1], our results present another way for tuning the photostrictive behavior of (PbTiO ) 3 n /(SrTiO ) 3 m superlattices, which could be utilized for innovative optomechanical applications.

## I. INTRODUCTION

In the last decade, ferroelectric/dielectric superlattices (FE/DE SLs) have gained significant interest because of their potential in next-generation multifunctional devices [2]. Among these FE/DE SLs, (PbTiO ) 3 n /(SrTiO ) 3 m have been widely studied as a model system [3-13]. Many peculiar properties have been evidenced in such nanosystems, such as polar vortices and skyrmions or negative capacitance effects [4-13]. Recently, significant photostriction -the light-induced mechanical deformation in a material -has been evidenced in these superlattices [14-16], adding yet another dimension to the rich number of their functionalities. Yet, there is little microscopic understanding of the origin of this effect in FE/DE SLs. There exists thus an urgent need to develop models and understanding of the photo-induced deformation in FE/DE SLs.

In a parallel report [1], we show how the electrostatic coupling strength between the layers of a PTO/STO SL with fixed chemical composition can strongly change the photostrictive response from a compression to an extension by varying the overall SLs period from 2 to 10 perovskite unit cells. In the present report, we demonstrate how the chemical composition, at fixed SLs period, can also strongly impact the amplitude and sign of the photostriction. Specifically, we present a first-principle study of the photostrictive behavior of PbTiO ) 3 n /(SrTiO ) 3 m SLs ( m + n = 10). Here, we generalize and extend the discussion of the results presented in [1], changing the chemical period of the SLs. Our aim is to provide, in addition to the results in [1], a general theoretical understanding into the mechanisms and ways to tune photostriction in these SLs with any chemical composition. In the present report, we aim at uncovering the dependence of photostriction in PTO/STO SLs on the PTO volume fraction. Our report is organized as follows. Construction of the simulated PTO/STO SLs and the computational methods are detailed in Sec II. In Sec III, we briefly summarize the results from [1] on ( n n | ) SLs. In Sec IV, we present and analyze the photostrictive behavior of the ( n m | ) SLs as function of PTO ratio. Section V is concerned with phenomenological models to quantify lightinduced strains in the SLs. And we close this report by a summary of our findings and future perspectives. We thus hope that the photostrictive tunability of PTO/STO SLs explored in this article may be of significant interest to the application of FE/DE SLs for optomechanics application and may open new avenues for fundamental explorations.

## II. SYSTEMS AND METHODS

## A. Computational details

Calculations in this report are based on density functional theory (DFT) using the projector augmented wave method [17] as implemented in the Abinit package [1820]. The Perdew-Burke-Ernzerhof generalized gradient approximation revised for solids is used for the exchangecorrelation functional [21]. The plane-wave energy cutoff is set to 20 Hartree. Integration over the Brillouin zone is done using a Monkhorst-Pack [22] scheme equivalent to 6 × 6 × 6 k-point mesh in a five-atom perovskite cell. To simulate the c phase, the Ti atoms are displaced up along the z axis from their centrosymmetric positions, and the out-of-plane cell dimension, atomic positions and angles are relaxed until forces on the supercell are less than 10 -8 Hartree.Bohr -1 and the force on ions less than 2 × 10 -6 Hartree.Bohr -1 . At the end of the relaxation, the obtained structures have the P4mm space group symmetry with polarization along the +ˆ direction. Starting z from the optimized structure in the dark, photoexcitation is simulated using the two quasi-Fermi levels method to mimic thermalized electrons and holes recently implemented in the Abinit [23] code by constraining a fixed number n ph of electrons/holes to the conduction/valence

bands. Under the constraint of n ph electrons/holes in the conduction/valence bands, the structures were relaxed again with the same convergence criteria as in the dark ( n ph = 0).

## B. Systems

To simulate the PTO/STO SLs, we adopt a supercell approach used in many previous studies on FE/DE SLs [5-7, 24, 25]. Stoichiometric SLs are built by stacking alternatively n unit cells of PTO on top of m unit cells of STO with n + m = 10 along the [001] direction, corresponding to a total of 50 atoms in our simulation box. To construct the SLs, we first optimized the lattice parameter of both PTO and STO. For cubic STO, the optimized lattice value is 3 897 ˚ in good agreement with existing . A experimental value 3 905 ˚ [26-28] and previously com-. A puted values [29]. For cubic PTO, we found a = 3 9248 . ˚ which also A is in good agreement with experimental [30-32] and previously reported theoretical values [3335]. For tetragonal PTO, we found a = 3 885 ˚ . A and c/a = 1 0712 well . comparable to experimental values 3 902 ˚ ( . A c/a = 1 0650) [30-32] and previous computa-. tional studies 3 888 ˚ and . A c/a = 1 0739 [33, 34]. . Before structural relaxation, atoms are placed in their ideal cubic positions, and the out-of-plane lattice constant of the supercell is set to ( ma STO + na PTO ), with a STO and a PTO being the theoretical optimized values of cubic PTO and STO, respectively. A schematic of the SLs is shown in Fig 1. Only SLs with polarization along the [001] direction commonly denoted the c phase in literature are considered in this study [36]. Structural, electronic and physical properties of the c phase of PTO/STO SLs are well documented [5, 8, 9, 11, 13, 25, 37, 38]. Additionally, we build reference paraelectric phases for the calculation of the layer-by-layer polarization. We aim at simulating PTO/STO SLs grown on an STO substrate as commonly done in experiment. We therefore fix the in-plane lattice constant of the SLs to the theoretical computed value for STO. We did not consider oxygen octahedral tilts modes in our simulation as they would require doubling the cell along the in-plane directions making our simulations computationally involved especially during photoexcitation. Also only single domain SLs are considered in this study.

## III. RESULTS FROM THE ( n n | ) PTO/STO SUPERLATTICES

To give context to the results presented in this report, we first summarize the findings presented in the companion paper [1]. In [1], we investigated the light-induced deformation of single domain PTO/STO SLs with equal PTO and STO chemical fraction denoted ( n n | ) SLs, n varying from 1 to 5. We show that, for very thin SLs ( n = 1 -2), free carriers are mostly delocalized and screen existing dipoles, which induces contractions in the outof-plane lattice of the SLs. On the other hand, when the film is thick enough ( n = 3 -5), the free charges are strongly localized in the interface regions. They produce a field that opposes existing fields (polarizing field in STO and depolarizing field in PTO) in the SLs and thereby induces a contraction in the STO layer and oppositely an expansion in the PTO layer. In the SLs, we found the scaling P z ∝ η α T = [( c -c p ) /a ] α relation ( c and c p are the out-of-plane lattice constants of the polar SLs and the corresponding paraelectric SLs, and a the inplane lattice parameter of the SLs), consistent with earlier reports [9, 15]. This relation enables us to relate the light-induced strains to the relative light-induced change in the polarization as follows:

$$\eta _ { 3 3 } = \frac { \eta _ { d } } { \alpha } \left ( \frac { \Delta P _ { z } } { P _ { z } } \right ) ;$$

with η d = ( c -c p ) /c . As shown in [1], the estimated light induced strain using Eq. 1 agrees well with the DFT data.

## IV. PHOTOSTRICTION AS FUNCTION OF PTO VOLUME FRACTION

We start by discussing the photostrictive behavior of the ( n m | ) SLs as a function of the PTO fraction. As shown in Fig. 2, the photostrictive behavior of the ( n m | ) SLs is dependent on the PTO fraction. Like for the ( n n | ) SLs [1], two behaviors can be distinguished when we look at the sign of the light-induced deformation in the out-of plane lattice of the SLs that we denote as η 33 = ( c -c 0 ) /c 0 , where c is the out-of-plane lattice constant of the SLs under illumination and c 0 its outof-plane lattice constant in the absence of illumination. As shown in Fig. 2, for PTO fractions less than a critical fraction x c ≈ 0 34 . (estimated from the graphs on Fig. 2; this value is also independent of n ph ), the SLs shows negative photostriction ( η 33 &lt; 0) while for fraction between x c and 0 9, the SLs shows positive photo-. striction ( η 33 &gt; 0). Note also that beyond x c , the lightinduced deformation in the SLs, for each considered n ph , as function of PTO fraction, reaches a maximal value at x = 0 6 and subsequently decreases back to negative . values in pure bulk PTO, consistent with previous calculations [35, 39]. The layer-by-layer decomposition of the photo-induced average tetragonality ( c/a ) shows that, for PTO fraction x ≤ 0 2, optical excitation penalizes distor-. tions in both PTO and STO layers (see Fig 3a where we show data for the (2 8) SLs), similar to (1 1) and (2 2) | | | SLs [1]. On the other hand, for PTO fraction x ≥ 0 3, . the free carriers decrease distortion in STO layer while they enhance distortion in the PTO layer (see Figs 3b, c&amp; d for selected SLs, see [40] for the other SLs), similar to our observation in the ( n n | ) SLs with n ≥ 3

We estimate the polarization in the SLs in a layer-bylayer fashion as done in the companion paper [1]. On

aSTO

FIG. 1. Schematic of the simulated PTO/STO superlattices. The atoms are shown in their perfect centrosymmetric positions.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-9dd74d9c-97a8-4f72-b056-c1f2040226bf.jpg)

0.4

Average c/a

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-1f09cfcc-9000-4d64-a438-25106fdd5bbe.jpg)

Average c/a

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-6b08c30b-fa0f-4ad6-a262-e973a9149ff0.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-78f67790-b564-42e8-8f74-58dde70d2b37.jpg)

FIG. 2. Light-induced strains as function of PTO fraction in the ( n m | ) SLs.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-fbd333db-fde7-4bfa-8c88-767c67ccb018.jpg)

Fig 4, we show the estimated total polarization as well as polarization in each layer for selected SLs (See [40] for the rest of the SLs) as function of n ph . The decomposed polarization per local five atoms unit cell of the SLs are also shown in Fig 5 (we only show for selected SLs, see the rest in the supplemental materials). For PTO fraction x ≤ 0 2, . optical excitation reduces polarization in both PTO and STO layer and thus the reducing of the total polarization in the SLs, see Figs 4a and Figs 5a. When the fraction of PTO is increased above 0 2 . ( x ≥ 0 3), . polarization in STO is decreased by optical excitation whereas polarization in PTO is enhanced, see Figs 4 and Figs 5 (b,c,d). The total polarization of the SLs decreases for x ≤ 0 3 and increases for . x ≥ 0 4 PTO . fraction. These light-induced change in the polarization in the SLs are in line with the observations from the lattice distortion in the SLs discussed in the precedent section. The overall behaviors (polarization and lattice distortion) of these SLs is thus the resultants of two actions

FIG. 3. a) Average c/a for PTO and STO layer in the (2 8) SLs, b) Average c/a for PTO and STO layer in the (3 7) | | SLs, c) Average c/a for PTO and STO layer in the (7 3) SLs, | d) Average c/a for PTO and STO layer in the (8 2) SLs as | function of n ph = 0 0 . -0 2 e/f.u. .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-32a6e571-4092-4135-8c4b-cb4fd11cd4b5.jpg)

(the total strain in the SLs is the resultant of the lightinduced strain in PTO layer and STO layer as shown in the piezoelectric model in section VB). This seems to be a general feature of FE/DE SLs, since a recent experiment on (2 4) BaTiO /CaTiO | 3 3 SLs showed that optical excitation induced compression in CaTiO 3 and oppositely an expansion in BaTiO 3 [41].

## V. MODELS AND DISCUSSIONS

## A. First Model: P z ∝ η α T

To describe the light-induced mechanical response within a phenomenological method, we applied the P z ∝ η α T model, described in our companion paper [1], to each

a)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-9a53f27d-0eb6-4187-8581-ef92b50f990e.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-80274760-99a7-4c15-9c25-50edffd8bbac.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-b6403a7c-5844-47a8-b1b7-103f5220cb39.jpg)

FIG. 4. Polarizations as function of number of excited carriers for the a) (2 8) SLs, b) (3 7) SLs, c) (7 3) SLs, d) (8 2) | | | | SLs.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-053d5a1c-1ee8-4b6a-86f7-0709db934160.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-484212b0-e1fc-466a-b2b9-537a49466085.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-c816d540-449f-4784-a30c-d1a04bedabdd.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-ac520d6e-d33c-4cee-9aab-ede90c067ab3.jpg)

FIG. 5. Layer-by-layer decomposition of the polarizations as function of number of excited carriers for the a) (2 8) SLs, | b) (3 7) SLs, c) (7 3) SLs, d) (8 2) SLs. | | | are to separate local cells corresponding to PTO and STO layer in each SLs.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-338bc124-db3b-487c-b961-a5b85919a3bc.jpg)

of the SLs. In this phenomenological model, the lightinduced strain is related to the relative change in the polarization via Eq. 1. For each of the SLs, we obtained the value of α by a linear fit of ln( P z ( n ph )) vs ln( η T ( n ph )) as described in [1] (see [40] for the values of α ). On figure 6, we show the results (for selected SLs, the results for other SLs can be found in the Supplemental Material [40]) of our fit along with the estimation of the strain using α = 0 5 and . α = 1 1 for comparaison. . The model with estimated α from the scaling relation fits better the DFT data for the corresponding SLs. As discussed in [1] and by Lee et al. [15], we are unable to attach a physical meaning to α and its changes as function of photoexcited charges and the type of SLs. We leave this as an

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-e5ba5afc-a8dc-4f01-8dbb-85206f822bb8.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-0a2d71be-329f-4f49-b014-e15c62bad2e5.jpg)

FIG. 6. Different fits of the DFT data ( η 33 vs n ph ) using Eq 1 for a) (2 8) SLs, b) (3 7) SLs, c) (7 3) SLs. | | |

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-e0c13b2e-f30c-4929-9743-a4ab5f16aa2e.jpg)

interesting future exploration.

## B. Converse piezoelectric effect

Here, we aim at finding another simple model to describe and reproduce our DFT data. The total lightinduced strain can be approximated as the weighted average of the deformation in STO and PTO layers in the SLs (see Supplemental Material [40] for details and discussion of the derivation):

$$\eta _ { 3 3 } \approx x \delta \eta _ { 3 3 } ^ { p } + \left ( 1 - x \right ) \delta \eta _ { 3 3 } ^ { s } \text{ \quad \ \ } \text{(2)}$$

̸

̸

where x is the volume fraction of PTO in the SLs. δη p 33 is the light-induced strain in the PTO layer and δη s 33 the light-induced strain in the STO layer. The volume fraction x changes as function of n ph . Introducing the changes in the polarization ( δP 1 = P 1 ( n ph = 0) -P 1 ( n ph = 0)) for PTO and ( δP 2 = P 2 ( n ph = 0) -P 2 ( n ph = 0)) for STO), and writing the following converse piezoelectric relation between strain and change in polarization [35]:

$$\delta \eta _ { 3 3 } ^ { p } = \epsilon _ { 0 } ^ { - 1 } \frac { d _ { 3 3 } ^ { 1 } } { \chi _ { 3 3 } ^ { 1 } } \delta P _ { 1 } ; \ \delta \eta _ { 3 3 } ^ { s } = \epsilon _ { 0 } ^ { - 1 } \frac { d _ { 3 3 } ^ { 2 } } { \chi _ { 3 3 } ^ { 2 } } \delta P _ { 2 }, \quad ( 3 )$$

we have:

$$\underset { \alpha 1 } { \overset { \infty } { \text{al} } } \quad \eta _ { 3 3 } \approx \epsilon _ { 0 } ^ { - 1 } \left [ x \left ( \frac { d _ { 3 3 } ^ { 1 } } { \chi _ { 3 3 } ^ { 1 } } \delta P _ { 1 } - \frac { d _ { 3 3 } ^ { 2 } } { \chi _ { 3 3 } ^ { 2 } } \delta P _ { 2 } \right ) + \frac { d _ { 3 3 } ^ { 2 } } { \chi _ { 3 3 } ^ { 2 } } \delta P _ { 2 } \right ] ; \quad ( 4 )$$

where ( d 1 33 , χ 1 33 ) and ( d 2 33 , χ 2 33 ) are the coupled piezoelectric constants, dielectric susceptibility of the PTO layer and the STO layer respectively. Note that although STO is cubic in its bulk structure, it becomes tetragonal, with a finite internal field, in the SLs due to its proximity

with the PTO layer and thus has a nonvanishing piezoelectric constant d 33 . The quantities d i 33 , and χ i 33 change with illumination and their values can be computed, in principle, within density functional perturbation theory (DFPT). However, the current implementation of DFPT in Abinit does not account for the presence of free charge in the conduction band. As such, we define the ratio r i = d i 33 /χ i 33 and estimate r i from the DFT data. The estimation of r i are detailed in the Supplemental Materials [40]. Eq. 4 can therefore be rewritten as:

$$\eta _ { 3 3 } \approx \epsilon _ { 0 } ^ { - 1 } \left [ x \left ( r ^ { 1 } \delta P _ { 1 } - r ^ { 2 } \delta P _ { 2 } \right ) + \delta P _ { 2 } r ^ { 2 } \right ] \quad ( 5 ) \quad \text{mar}$$

As observed in the DFT data, this model predicts a critical PTO fraction for a critical SLs that is at the boundary of negative and positive photostrictive SLs:

$$x _ { c } = \frac { 1 } { 1 - \frac { \delta P _ { 1 } } { \delta P _ { 2 } } \frac { r ^ { 1 } } { r ^ { 2 } } } \quad \quad \ \ ( 6 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \end{cases}$$

Interestingly, the DFT data suggests that this critical volume fraction is independent of n ph . On Fig 7, we show the results from this model and the DFT data for n ph = 0 00 . -0 20 e/f.u for few selected SLs (the rest of the . SLs is shown in the Supplemental Material [40]). For such a simple model, it agrees rather well and captures the underlying physics. The model, however, fails to quantitatively capture the light induced strains at higher n ph for some of the SLs. This is because at high n ph , it can be noticed that the photostrictive behavior of some of the SLs becomes very nonlinear and a rigorous calculation of the r i ratio in Eq. 5 with DFPT might improve the model at large n ph . We also would like to relate this

(%)

FIG. 7. Different fits of the DFT data ( η 33 vs n ph ) using Eq 4 for a) (2 8) SLs, b) (3 7) SLs, c) (7 3) SLs, d) (8 2) SLs. | | | |

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-ac0f7105-680b-4f47-b4ac-a29cdf0b147c.jpg)

study to the two only existing experimental studies we are aware of on light-induced strain in PTO/STO SLs. In Refs. [14, 15], a film of 8 units of PTO and 3 units of STO

was used. This corresponds to roughly a 0 72 PTO frac-. tion. They all reported an overall positive photostrictive behavior, qualitatively consistent with the observation in our study.

## VI. CONCLUSION AND PERSPECTIVES

In summary, using DFT within a constrained occupation scheme, we investigate photostriction in single domain PTO/STO SLs with varying chemical composition. We find that, for PTO fraction x ≤ 0 2, both layers of the . SLs contract upon optical excitation, leading an overall contraction of the SLs. On the other hand, for x ≥ 0 3 , . optical excitation induces a contraction in the STO layers and an expansion in the PTO layers. Depending on the fraction of PTO in the SLs, an overall contraction or expansion is observed. This study, in conjunction with the report in [1], show that the photostrictive behavior of these SLs can be tuned not only through the SLs composition but also its thickness. The possibility of tuning the light-induced strain by simply changing the SLs composition or thickness is appealing as it presents potential for novel applications using light as modulator.

In the present study, as well as in the one in [1], we restrict ourselves, because of the computational cost of simulating optical excitation in large system, to single domain PTO/STO SLs with no octahedral tilts (oxygen octahedral tiltings only appear in STO below 105K [42] and does not exist in pure PTO. Additionally existing experimental studies[14-16] were done at much higher temperature and did not consider any octahedral tiltings, thus neglecting their effect is a realistic consideration), consisting of a total of 10 five-atom perovskite cells. Consequently, this choice limits the predictive power of our studies and a quantitative comparison with reported experimental values in Refs.[14, 15] can be misleading. However, note that our results are qualitatively in line with these reports. We believe that, through these combined investigations, we have provided first-principle insights into the photostrictive behavior of PTO/STO SLs. Additionally, our studies can be used by experimentalists as a starting point to optimize the photostrictive response of PTO/STO SLs. A straightforward future development is to extend these studies to thicker SLs (i.e., SLs with more than 10 five-atoms perovskite cells) and investigate how the dependence on PTO chemical fraction of the photostrictive behavior (shown on Fig 2) scales with the SLs thickness. Another approach would be to extend this study to SLs with domain structures, octahedral tilts or super-orders [43, 44]. Expanding on the ideas presented in our reports in these directions offers an interesting avenue for future explorations. However, pursuing these ideas may require treating the 'constrained occupation scheme' in an effective way through large scale methods such as effective Hamiltonians or second-principles approaches, as simulating photoexcitation in these large systems using DFT becomes computationally cumber-

## ACKNOWLEDGMENTS

The work is supported by the ARO Grants No. W911NF-21-2-0162 (ETHOS) and W911NF-21-1-0113,

- [1] C. Dansou, C. Paillard, and L. Bellaiche, Tailoring photostriction via superlattices engineering, Submitted to Physical Review Letters .
- [2] S. Das, Z. Hong, M. McCarter, P. Shafer, Y. T. Shao, D. Muller, L. Martin, and R. Ramesh, A new era in ferroelectrics, APL Mater. 8 , 120902 (2020).
- [3] J. Jiang, X. Pan, W. Tian, C. Theis, and D. Schlom, Abrupt PbTiO /SrTiO 3 3 superlattices grown by reactive molecular beam epitaxy, Appl. Phys. Lett. 74 , 2851 (1999).
- [4] A. Y. Abid, Y. Sun, X. Hou, C. Tan, X. Zhong, R. Zhu, H. Chen, K. Qu, Y. Li, M. Wu, J. Zhang, J. Wang, K. Liu, X. Bai, D. Yu, X. Ouyang, J. Wang, J. Li, and P. Gao, Creating polar antivortex in PbTiO /SrTiO 3 3 superlattice, Nature Commun. 12 , 1 (2021).
- [5] P. Aguado-Puente, P. Garcia-Fernandez, and J. Junquera, Interplay of couplings between antiferrodistortive, ferroelectric, and strain degrees of freedom in monodomain pbtio /srtio 3 3 superlattices, Phys. Rev. Lett. 107 , 217601 (2011).
- [6] P. Aguado-Puente and J. Junquera, Structural and energetic properties of domains in PbTiO /SrTiO 3 3 superlattices from first principles, Phys. Rev. B 85 , 184105 (2012).
- [7] E. Bousquet, M. Dawber, N. Stucki, C. Lichtensteiger, P. Hermet, J. M. T. S. Gariglio, and P. Ghosez, Improper ferroelectricity in perovskite oxide artificial superlattices, Nature 452 , 732 (2008).
- [8] M. Dawber, C. Lichtensteiger, M. Cantoni, M. Veithen, P. Ghosez, K. R. K. Johnston, and J. M. Triscone, Unusual behavior of the ferroelectric polarization in PbTiO /SrTiO 3 3 superlattices, Phys. Rev. Lett. 95 , 177601 (2005).
- [9] M. Dawber, N. Stucki, C. Lichtensteiger, S. Gariglio, P. Ghosez, and J. M. Triscone, Tailoring the properties of artificially layered ferroelectric superlattices, Adv. Mater. 19 , 4153 (2007).
- [10] S. Das, Y. Tang, Z. Hong, M. Gon¸ calves, M. McCarter, C. Klewe, K. Nguyen, F. G´mez-Ortiz, P. Shafer, o E. Arenholz, V. A. Stoica, S. L. Hsu, B. Wang, C. Ophus, J. F. Liu, C. T. Nelson, S. Saremi, B. Prasad, A. B. Mei, D. G. Schlom, J. ´ I˜iguez, n P. Garc´ ıa-Fern´ndez, a D. A. Muller, L. Q. Chen, J. Junquera, L. W. Martin, and R. Ramesh, Observation of room-temperature polar skyrmions, Nature 568 , 368 (2019).
- [11] P. Zubko, N. Jecklin, N. Stucki, C. Lichtensteiger, G. Rispens, and J. M. Triscone, Ferroelectric domains in PbTiO /SrTiO 3 3 superlattices, Ferroelectrics 433 , 127 (2012).
- [12] P. Zubko, J. C. Wojdeglyph[suppress] l, M. Hadjimichael, S. FernandezPena, I. L. A. Sen´, J. M. Triscone, and J. ´ ˜iguez, Nege In

and the Vannevar Bush Faculty Fellowship (VBFF) grant no. N00014-20-1-2834 from the Department of Defense. C.P. acknowledges the support from a public grant overseen by the French Agence Nationale de la Recherche under grant agreement no. ANR-21-CE24-0032 (SUPERSPIN). We also acknowledge the computational support from the Arkansas High Performance Computing Center (AHPCC) for computational resources.

ative capacitance in multidomain ferroelectric superlattices, Nature 534 , 524 (2016).

- [13] P. Zubko, N. Jecklin, A. Torres-Pardo, P. AguadoPuente, A. Gloter, C. Lichtensteiger, J. Junquera, O. St´ ephan, and J.-M. Triscone, Electrostatic coupling and local structural distortions at interfaces in ferroelectric/paraelectric superlattices, Nano letters 12 , 2846 (2012).
- [14] Y. Ahn, J. Park, A. Pateras, M. B. Rich, Q. Zhang, P. Chen, M. H. Yusuf, H. Wen, M. Dawber, and P. G. Evans, Photoinduced domain pattern transformation in ferroelectric-dielectric superlattices, Phys. Rev. Lett. 119 , 057601 (2017).
- [15] H. J. Lee, Y. Ahn, S. D. Marks, E. C. Landahl, S. Zhuang, M. H. Yusuf, M. Dawber, J. Y. Lee, T. Y. Kim, S. Unithrattil, S. H. Chun, S. Kim, I. Eom, S. Y. Park, K. S. Kim, S. Lee, J. Y. Jo, J. Hu, and P. G. Evans, Structural evidence for ultrafast polarization rotation in ferroelectric/dielectric superlattice nanodomains, Phys. Rev. X 11 , 031031 (2021).
- [16] D. Daranciang, M. J. Highland, H. Wen, S. M. Young, N. C. Brandt, H. Y. Hwang, M. Vattilana, M. Nicoul, F. Quirin, J. Goodfellow, et al. , Ultrafast photovoltaic response in ferroelectric nanolayers, Phys. Rev. Lett. 108 , 087601 (2012).
- [17] P. E. Bl¨chl, Projector augmented-wave method, Physio cal review B 50 , 17953 (1994).
- [18] X. Gonze, B. Amadon, G. Antonius, F. Arnardi, L. Baguet, J.-M. Beuken, J. Bieder, F. Bottin, J. Bouchet, E. Bousquet, et al. , The abinit project: Impact, environment and recent developments, Computer Physics Communications 248 , 107042 (2020).
- [19] F. Bottin, S. Leroux, A. Knyazev, and G. Z´rah, Largee scale ab initio calculations based on three levels of parallelization, Computational Materials Science 42 , 329 (2008).
- [20] A. H. Romero, D. C. Allan, B. Amadon, G. Antonius, T. Applencourt, J. B. L. Baguet, F. Bottin, J. Bouchet, E. Bousquet, F. Bruneval, G. Brunin, M. C. D. Caliste, J. Denier, C. Dreyer, P. Ghosez, M. Giantomassi, Y. Gillet, O. Gingras, D. R. Hamann, G. Hautier, F. Jollet, G. Jonard, A. Martin, H. P. C. Miranda, F. Naccarato, G. Petretto, N. A. Pike, V. Planes, S. Prokhorenko, T. Rangel, F. Ricci, G. M. Rignanese, M. Royo, M. Stengel, M. Torrent, M. J. van Setten, B. V. Troeye, M. J. Verstraete, J. Wiktor, J. W. Zwanziger, and X. Gonze, Abinit: Overview and focus on selected capabilities, J. Chem. Phys. 152 , 124102 (2020).
- [21] J. P. Perdew, A. Ruzsinszky, G. I. Csonka, O. A. Vydrov, G. E. Scuseria, L. A. Constantin, X. Zhou, and K. Burke, Restoring the density-gradient expansion for exchange in

- solids and surfaces, Phys. Rev. Lett. 100 , 136406 (2008).
- [22] H. J. Monkhorst and J. D. Pack, Special points for brillouin-zone integrations, Physical review B 13 , 5188 (1976).
- [23] C. Paillard, E. Torun, L. Wirtz, J. ´ ˜iguez, and L. BelIn laiche, Photoinduced phase transitions in ferroelectrics, Phys. Rev. Lett. 123 , 087601 (2019).
- [24] N. Bristowe, E. Artacho, and P. Littlewood, Oxide superlattices with alternating p and n interfaces, Phys. Rev. B 80 , 045425 (2009).
- [25] M. Gu, J. Wang, Q. Xie, and X. Wu, Structural and electronic properties of PbTiO /SrTiO 3 3 superlattices from first principles, Phys. Rev. B 82 , 134102 (2010).
- [26] V. M. Longo, M. d. G. S. Costa, A. Z. Simoes, I. L. V. Rosa, C. O. P. Santos, J. Andres, E. Longo, and J. A. Varela, On the photoluminescence behavior of samariumdoped strontium titanate nanostructures under uv light. a structural and electronic understanding, Phys. Chem. Chem. Phys. 12 , 7566 (2010).
- [27] N. F. Muhamad, R. A. M. Osman, M. S. Idris, and M. N. M. Yasin, Physical and electrical properties of SrTiO 3 and SrZrO , in 3 EPJ Web of Conferences , Vol. 162 (EDP Sciences, 2017) p. 01052.
- [28] O. Nakagawara, M. Kobayashi, Y. Yoshino, Y. Katayama, H. Tabata, and T. Kawai, Effects of buffer layers in epitaxial growth of SrTiO 3 thin film on si (100), Journal of applied physics 78 , 7226 (1995).
- [29] R. Eglitis and D. Vanderbilt, First-principles calculations of atomic and electronic structure of SrTiO (001) and 3 (011) surfaces, Phys. Rev. B 77 , 195408 (2008).
- [30] R. Nelmes and W. Kuhs, The crystal structure of tetragonal PbTiO 3 at room temperature and at 700 k, Solid state communications 54 , 721 (1985).
- [31] G. Shirane, R. Pepinsky, and B. Frazer, X-ray and neutron diffraction study of ferroelectric PbTiO 3 , Phys. Rev. 97 , 1179 (1955).
- [32] A. Glazer and S. Mabud, Powder profile refinement of lead zirconate titanate at several temperatures. ii. pure PbTiO , Acta Crystallographica 3 Section B: Structural Crystallography and Crystal Chemistry 34 , 1065 (1978).
- [33] R. A. Ribeiro and S. R. d. L´zaro, Theoretical investigaa tions of the bulk modulus in the tetra-cubic transition of PbTiO 3 material, Qu´ ımica Nova 37 , 1165 (2014).
- [34] Y. Umeno, B. Meyer, C. Els¨sser, and P. Gumbsch, Ab a initio study of the critical thickness for ferroelectricity in ultrathin Pt/PbTiO /Pt films, Phys. Rev. B 3 74 , 060101 (2006).
- [35] C. Paillard, S. Prosandeev, and L. Bellaiche, Ab initio approach to photostriction in classical ferroelectric materials, Phys. Rev. B 96 , 045205 (2017).
- [36] N. Pertsev, A. Zembilgotov, and A. Tagantsev, Equilibrium states and phase transitions in epitaxial ferroelectric thin films, Ferroelectrics 223 , 79 (1999).
- [37] Z. Zhen-Ye, W. Biao, W. Hai, Z. Yue, and L. Qing-Kun, The first-principles study of ferroelectric behaviours of PbTiO /SrTiO 3 3 and BaTiO /SrTiO 3 3 superlattices, Chinese Physics 16 , 1780 (2007).
- [38] R. Eglitis, S. Piskunov, and Y. F. Zhukovskii, Ab initio calculations of PbTiO /SrTiO 3 3 (001) heterostructures, physica status solidi (c) 13 , 913 (2016).
- [39] C. Paillard, B. Xu, B. Dkhil, G. Geneste, and L. Bellaiche, Photostriction in ferroelectrics from density functional theory, Phys. Rev. Lett. 116 , 247401 (2016).
- [40] See supplemetal material at [url] for the data on the other superlattices., .
- [41] D. Sri Gyan, H. J. Lee, Y. Ahn, R. B. Carson, J. Carnis, T. Y. Kim, S. Unithrattil, J. Y. Lee, S. H. Chun, S. Kim, et al. , Optically induced picosecond lattice compression in the dielectric component of a strongly coupled ferroelectric/dielectric superlattice, Advanced Electronic Materials 8 , 2101051 (2022).
- [42] J. Scott, Soft-mode spectroscopy: Experimental studies of structural phase transitions, Reviews of Modern Physics 46 , 83 (1974).
- [43] R. Xu, F. Delodovici, B. Dkhil, and C. Paillard, Oxygen tilt driven polar superorders in BiFeO -based superlat3 tices, Phys. Rev. B 109 , L220101 (2024).
- [44] R. Gu, R. Xu, F. Delodovici, B. Carcan, M. Khiari, G. Vaudel, V. Juv´ e, M. Weber, A. Poirier, P. Nandi, et al. , Superorders and acoustic modes folding in BiFeO /LaFeO 3 3 superlattices, arXiv preprint arXiv:2401.06647 (2024).

## Supplemental Materials for 'Tuning photostriction in (PbTiO ) /(SrTiO ) 3 n 3 m superlattices via chemical composition: An ab-initio study'

## VII. AVERAGE C/A IN THE (1 9) | , (4 6) | , (6 4) | AND (9 1) | SLS

(PTO) /(STO)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-9b3a0ddc-afe3-4f20-b125-394f6b2f872e.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-201a98c7-6dea-472c-a12a-876b42615347.jpg)

(PTO) /(STO)

FIG. 8. a) Average c/a for PTO and STO layer in the (1 9) SLs, b) Average c/a for PTO and STO layer in the (4 6) SLs, c) | | Average c/a for PTO and STO layer in the (6 4) SLs, d) Average c/a for PTO and STO layer in the (9 1) SLs as function of | | n ph = 0 0 . -0 2 e/f.u. .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-8bc30a23-db2a-40c0-876e-7d4cf5d8bb1e.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-a3dcf437-6853-42b9-98ad-3b16a171d322.jpg)

(PTO) /(STO)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-3a91f5d7-01f3-4eb7-b584-0a35524eaf8c.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-15f05d61-86ed-495b-b121-a2233b1f9842.jpg)

c)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-a60a3842-7f4b-4e04-b5aa-ea3209ded61c.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-75051457-a16c-4610-a155-2096fbb6e4fe.jpg)

d)

FIG. 9. Polarizations as function of number of excited carriers for the a) (1 9) SLs, b) (4 6) SLs, c) (6 4) SLs, d) (9 1) SLs. | | | |

/s945

FIG. 10. a) Values of α as function of the ( n n | ) SLs; b) Values of α as function of PTO fraction in the ( n m | ) SLs.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-ba4e4266-88c5-4ddc-9217-328330ece09b.jpg)

## X. RESULTS OF THE MODEL USING P z ∝ η α T FOR THE (1 9) | , (4 6) | , (6 4) | AND (9 1) | SLS

(%)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-febb7422-4b1e-4d6d-8f2a-726ede5af509.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-822f180a-7a55-400a-a0cf-1193abb763ba.jpg)

c)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-3a7b87df-d6ce-4aba-8a56-a8afa15912a7.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-1ddedf24-2135-4c5a-859d-9855d46fde99.jpg)

d)

FIG. 11. Different fits of the DFT data ( η 33 vs n ph ) using the ' P z ∝ η α T model' in the main text for a) (1 9) SLs, b) (4 6) SLs, | | c) (6 4) SLs, d) (9 1) SLs. | |

a)

(%)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-cc1bb189-4c22-4be2-837d-f7ddcd635835.jpg)

c)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-59f7a3cd-ced8-429e-98b2-dc21ee8e3d1e.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-cd8cddd7-5cee-4760-919b-605d390d9892.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-fcba7e09-b5f6-48ff-9dea-8380397e6d67.jpg)

d)

(PTO) /(STO)

(PTO) /(STO)

(%)

FIG. 12. Different fits of the DFT data ( η 33 vs n ph ) using the converse piezoelectric model in the main text for a) (1 9) SLs, | b) (4 6) SLs, c) (6 4) SLs, d) (9 1) SLs. | | |

## XII. EXPRESSION OF THE TOTAL STRAIN IN THE SLS

̸

Let l p and l s be the lengths of the PTO and STO layer in the ground state structures ( n ph = 0). Also let l l p and l l s be their lengths after illumination ( n ph = 0). The light-induced strain in the SLs is as follows:

$$\eta _ { 3 3 } = \frac { l _ { s } ^ { l } + l _ { p } ^ { l } - l _ { s } - l _ { p } } { l _ { s } + l _ { p } } = \frac { l _ { s } ^ { l } - l _ { s } } { l _ { s } + l _ { p } } + \frac { l _ { p } ^ { l } - l _ { p } } { l _ { s } + l _ { p } } ;$$

now we write:

$$\frac { l _ { s } ^ { l } - l _ { s } } { l _ { s } + l _ { p } } = \frac { l _ { s } } { l _ { s } + l _ { p } } \frac { l _ { s } ^ { l } - l _ { s } } { l _ { s } } ; \ \frac { l _ { p } ^ { l } - l _ { p } } { l _ { s } + l _ { p } } = \frac { l _ { p } } { l _ { s } + l _ { p } } \frac { l _ { p } ^ { l } - l _ { p } } { l _ { p } }.$$

In the optimized ground state structure, the volume fraction x of PTO and (1 -x ) of STO in the SLs can be calculated by the followings:

$$x = \frac { l _ { p } } { l _ { s } + l _ { p } } \text{ and } 1 - x = \frac { l _ { s } } { l _ { s } + l _ { p } }.$$

Defining the following induced-strains in PTO and STO layers:

$$\eta _ { 3 3 } ^ { p } = \frac { l _ { p } ^ { l } - l _ { p } } { l _ { p } } \text{ and } \eta _ { 3 3 } ^ { s } = \frac { l _ { s } ^ { l } - l _ { s } } { l _ { s } },$$

we arrive at the formula for Eq. 2 in the main text:

$$\eta _ { 3 3 } & \approx x \eta _ { 3 3 } ^ { p } + ( 1 - x ) \eta _ { 3 3 } ^ { s }$$

Note that in this derivation, x as defined in Eq. 3 is independent of n ph . However, the lengths of the PTO and STO layer of the superlattices obtained from DFT change with n ph and thus the variation of x as function of n nph as shown on Figs 13. That is why we carefully use the approximation sign in written the final expression in Eq. 11. To have the converse piezoelectric model to better agree with the DFT data, as presented in this paper, we use the values of x as shown on Figs 13.

## XIII. DETAILS ON THE ESTIMATION OF THE r i RATIO FROM THE DFT DATA

Following [41], for each SLs, the polarization in each layer can be related to its tetragonal distortion via the following relation:

$$P _ { i } ^ { 2 } \in \frac { c _ { i } } { a } - 1$$

with c i being the average out-of-plane lattice parameter of the layer of the SLs ( i =1 for PTO and 2 for STO). Note that c i is function of n ph . Taking the logarithmic derivative of Eq.(12) gives:

$$2 \frac { \delta P _ { i } } { P _ { i } } = \frac { \delta c _ { i } } { c _ { i } - a },$$

which can be rewritten as:

$$\delta P _ { i } & = \frac { c _ { i } P _ { i } } { 2 } \frac { \eta _ { i } } { c _ { i } - a },$$

writting δP i = ϵ 0 χ i 33 δE i and η i = d i 33 δE i , we have:

$$\epsilon _ { 0 } \chi _ { 3 3 } ^ { i } = \frac { c _ { i } P _ { i } } { 2 } \frac { d _ { 3 3 } ^ { i } } { c _ { i } - a }$$

and thus we get:

$$r ^ { i } = \frac { d _ { 3 3 } ^ { i } } { \chi _ { 3 3 } ^ { i } } = 2 \epsilon _ { 0 } \frac { c _ { i } - a } { c _ { i } P _ { i } }$$

a)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-71813f45-494a-4873-9526-6ef728d76161.jpg)

PTO fraction

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-cf181ea7-f20f-4994-9daa-2e2e106a4d4f.jpg)

g)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-f65ecf2b-536f-4992-a7bf-fcc6664c09ba.jpg)

e)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-5433e9c7-a2ef-44b2-91eb-fb58641ff600.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-a3f32c02-8705-4f77-b8fc-cc8304c7372b.jpg)

PTO fraction

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-476c7dcb-a06c-4f50-b282-e15a34c2321e.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-70fee1f5-0b7e-4386-adf8-8417aa739e78.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-cbcd2f1e-796a-4c91-a70e-75ddab264e6c.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00105v1-figure-8ed72722-82cb-4bec-9878-8b5f9dfb8816.jpg)

h)

PTO fraction

FIG. 13. Variation of the volume fraction x as function of n ph for each of the SLs. a) to i) correspond to the (1 9), | (2 8),(3 7),(4 6),(5 5),(6 4),(7 3),(8 2),(9 1) respectively. | | | | | | | |
Carmel Dansou Charles Paillard Laurent Bellaiche
Tuning photostriction in (PbTiO$_{3}$)$_{n}$/(SrTiO$_{3}$)$_{m}$ superlattices via chemical composition: An $\textit{ab-initio}$ study
Light-induced mechanical deformations in single-domain
(PbTiO$_{3}$)$_{n}$/(SrTiO$_{3}$)$_{m}$ superlattices were simulated using
first-principle calculations. By varying the chemical fraction PbTiO$_{3}$, we
discover that these heterostructures' photostrictive behavior can be tuned
quantitatively and qualitatively. Additionally, we present simple analytical
models to explain the calculated deformations and predict a critical
PbTiO$_{3}$ fraction with no photostriction. In addition to the report in [1],
our results present another way for tuning the photostrictive behavior of
(PbTiO$_{3}$)$_{n}$/(SrTiO$_{3}$)$_{m}$ superlattices, which could be utilized
for innovative optomechanical applications.