## Tailoring photostriction via superlattices engineering

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)

We report systematic first-principles investigation of light-induced mechanical deformations in monodomain (PbTiO ) 3 n /(SrTiO ) 3 n superlattices ( n = 1 -5). We reveal that photostriction in these heterostructures quantitatively and qualitatively depends on the chemical period n . Specifically, we show that by changing the chemical period, we can induce positive or negative photostriction. We also present a simple analytical model to account for the calculated deformations. Our findings indicate that superlattices architectures may be key to design novel optomechanical applications.

Photostriction is a physical phenomenon in which light induces mechanical strains in a material [1]. It therefore converts light into mechanical energy [1-3]. In recent years, materials such as ferroelectrics [4-8], multiferroics [9-14], organic polymers [15], inorganic semiconductors [16] and halide perovskites [17-19] were investigated for their photostrictive response. Photostrictive properties can be harnessed for practical applications such as optomechanical sensors [20], energy harvesting devices [21], and phostrictive actuators [22, 23], among others.

In the last decade, the field of ferroelectrics has also witnessed the emergence of ferroelectrics/dielectrics (FE/DE) heterostructures. They host exotic properties that make them prime candidates for the next generation of electronic devices [24-33]. Recent experiments have shown that FE/DE superlattices (SLs) also exhibit lattice deformations upon optical excitation [34, 35]. Due to the fine, atomic layer control achieved by nowadays stateof-the-art growth techniques, SLs offer a unique playground to engineer, from the bottom-up, the photostrictive response. However, existing experimental reports are limited to as-grown SLs, with no systematic exploration of the parameters involved in their photo-induced deformations. For instance, only light-induced lattice expansion has been reported in the prototype PbTiO /SrTiO 3 3 (PTO/STO) SLs [34, 36]. First-principles calculations, as conducted here, allow on the other hand to explore the space of parameters while saving time and avoiding material waste.

theory (DFT) calculations. We employ the constrained occupation number approach developed in Ref. [38], and thus enforce n ph electrons (resp. holes) to stay in the conduction (resp. valence) bands using a Fermi-Dirac distribution with its own quasi-Fermi level µ e (resp. µ h ). Here, we consider ( n n | ) SLs with polarization along the [001] direction, generally referred to as c phase in literature [39]. We vary the periodicity n from 1 to 5 (and thus the chemical period from 2 to 10 perovskite cells). We thus sample SLs with electrostatic coupling between the PTO and STO layers varying from strong ( n ≤ 3) to weak ( n ≥ 3) as defined in Ref. [40]. Polarization is induced by displacing Ti atoms in the PTO layers up along the z-axis from their centrosymmetric position at the beginning of our simulations. The technical details and convergence parameters are the same as presented in the companion paper [41].

In this Letter, we propose a first-principle study of photostrictive behavior of single domain (PbTiO ) 3 n /(SrTiO ) 3 n (( n n | ) PTO/STO) SLs and show that, by changing the chemical period n of the SLs, the magnitude and even sign of the photo-induced deformations can be tuned. We also develop a simple phenomenological model that reproduces the observed unusual photostrictive behaviors qualitatively and quantitatively.

We use the Abinit code [37] for our density functional

Let us first look at the photo-induced strain η 33 = ( c n ( ph ) -c (0)) /c (0) (where c n ( ph ) is the out-of-plane lattice constant at the concentration of photo-excited carriers n ph ) as a function of n in the ( n n | ) SLs. As shown in Fig. 1, one can distinguish two behaviors when focusing on the sign of the light-induced strain. For the ultrathin chemical period ( n = 1 2) SLs i.e., in the strong electro-, static coupling regime of the PTO and STO layers, light induces a compression, i.e. negative photostriction (see Fig 1(a) for different chosen values of n ph ). In contrast, photostriction is positive ( i.e., the SLs expands in the out-of-plane direction) for thicker SLs n ≥ 3 ( i.e., in the weak electrostatic coupling regime of the PTO and STO layers). To investigate in more details these two regimes, we examine the light-induced tetragonal distortion in the different layers. Consequently, in Fig 1(b) &amp; c) is shown the average tetragonal distortion in the PTO and STO layers of the SLs as a function of the concentration of photoexcited carriers for all values of n . STO is an incipient ferroelectric material, however, the proximity with PTO induces in DFT calculations, a polarization in the STO layer (via an internal polarizing field) and thus a

a)

c)

Average c/a

FIG. 1. a) Lattice deformation in the ( n n | ) SLs as function of n ph . b) Average c/a of the PTO layer as function of n ph , c) Average c/a of the STO layer as function of n ph .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-aba50195-cdc3-43ab-bde5-92b50bb97c8a.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-bcdb2378-e416-4ebf-9772-b669445663b9.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-793986b2-74a1-498c-a827-b6775af3f913.jpg)

significant tetragonal distortion. For utrathin SLs with n = 1 2, both the distortion in the PTO and STO layers , (black and red curves in Figs 1 b&amp;c) decrease under optical excitation leading to the overall observed compression in the out-of-plane lattice constant; a feature reminiscent of the behavior of bulk ferroelectric titanates [5, 38]. On the other hand, in the SLs n ≥ 3, upon optical excitation, the distortion in the STO layer decreases whereas the distortion in the PTO layer increases (see blue, green and purple curves in Fig. 1 b&amp;c). The positive distortion in the PTO layers dominates and thus the overall expansion is observed. The data depicted in Fig. 1 thus suggest that the photostrictive response (negative or positive) of the PTO layer and the SL as a whole is strongly dependent on the strong or weak electrostatic coupling between the layers.

average of the polarization in PTO and STO layers. As shown on Fig 2 for selected SLs (See Supplementary Material [44] for the other SLs), optical excitation reduces the total polarization in ultra-thin ( n n | ) SLs ( n = 1 and n = 2). Additionally, polarization in both the STO and PTO layers decreases upon optical excitation. This is consistent with the decrease of the lattice distortion shown on Fig 1 (b&amp; c) for the (1 1) and (2 2) SLs. When | | the SLs layers are thick enough, that is for n ≥ 3 , the overall polarization is enhanced by optical excitation as shown by the black curves on Figs 2 (b&amp;c). Additionally, the polarization in PTO layer increases (see red curves on Fig 2 b &amp; c) whereas the induced polarization in the STO layer is decreased (see blue curves on Fig 2 b &amp; c) - implying that STO becomes less polar under light.

We now examine the evolution of the electrical polarization of the SLs as a function of the concentration of photo-excited carriers. The polarization in the SLs is computed in a layer-by-layer fashion using the following equation:

$$P _ { i } = \frac { e } { \Omega _ { i } } \sum _ { j \in i } w _ { j } \delta _ { j } Z _ { j } ^ { * } \quad \ \ ( 1 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \cdot \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat{ \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat{\ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat {\ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat {
 } }$$

where Ω , i e , and δ j are the volume of the local unit cell i (local unit cell is used here to refer to an effective five atoms unit cell centered on a Ti atom), the elementary charge and the atomic displacement from a built reference nonpolar SLs (we built this reference structure for each of the SLs) respectively. The weight factor w j tells how many of such local cells share the j ions. Its values are 1 / 2, 1, 1 / 8 for the O, Ti and (Pb,Sr) atoms respectively. Z ∗ j is the z component of the Born effective charge of atom j , with j running through all the atoms in the local cell i . The Born effective charges Z ∗ j are those of the atoms in dark conditions and are taken from Ref. [42]. We do not expect a significant alteration of Z ∗ j under the concentration of photo-excited carriers investigated here, and thus do not employ the more refined approach devised in Ref. [43]. The total polarization in the SLs is the

To gain further insight on the role of the free charges in the above observed changes in the lattice constant and polarization inside the SLs, we computed the macroaverage [45] excited charges density and the electrostatic potential inside the SLs under illumination for the considered ( n n | ) SLs. As shown in Fig 3 (a)-(b), in the ultrathin ( n n | ) SLs ( n = 1 or n = 2), the photo-excited carriers are more delocalized through the SLs than for n ≥ 3. When the thickness of the films is increased ( n is equal or greater than 3), the free charges become strongly localized on TiO 2 planes forming the interface, as shown in Fig. 3 (a) to (c) for the (1 1), (3 3) and (5 5) | | | SLs respectively (see the Supplemental Materials [44] for (2 2) and (4 4) SLs). | | The spacial distribution of the free charges in the ( n n | ) SLs provide us with the following pictures: the ultrathin SLs, where there is a strong electrostratic coupling between PTO and STO layers, behave like bulk system [5, 13, 38] where photo-excited carriers are spread out within the SL and lead to an overall decrease of polarization and thus out-of-plane lattice compression in both PTO and STO layers. For thick enough SLs, the free charges strongly localize at the interfaces. Electrons and holes populate given interfaces and produce fields that oppose the existing field in the STO and PTO layers within the SLs. The light-induced field in

a)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-c24178e7-1b4d-407c-9013-39e7b3270f17.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-30e388eb-33cf-4252-a138-22dfd993ebbb.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-c0addeea-0857-4fac-9f6e-1531ead44511.jpg)

b)

c)

d)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-0c6f5712-b7d6-4425-be45-1da4cae4efa2.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-4059667e-1a5c-45bf-9e3b-eb4c5abab3f5.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-a89aeb11-5365-4cea-8202-f5f64caa23d2.jpg)

e)

FIG. 2. Top panel: polarization as a function of n ph in the ( n n | )SLs. a)-c) (1 1), (3 3), and (5 5) SLs respectively. | | | Bottom panel: projection to five-atom unit cell of the polarization in the d)-f) (1 1), (3 3), and (5 5) SLs respectively. | | | The vertical dashed lines in e) and f) are to separate local cells corresponding to PTO and STO layer in each SLs.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-f6583f60-2f2c-4f05-9a28-46d39915dcbc.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-c57d74d1-5726-44e9-b2b6-0060941a9e85.jpg)

by mapping out the electrostatic potential in the SLs before and after illumination (see Fig. 3d for the (4 4) SLs, | the other are shown in the supplemental materials). In dark, the electrostatic potential has a non-zero slope in both the STO and PTO layers (see blue curve in Figure 3d). This suggests the presence of a finite field as postulated above. Upon illumination, the electrostatic potential flattens leading to vanishingly small electric fields within layers of the SLs (see orange curve in Fig 3d in the case of the (4 4) SLs). | Essentially, in thick SLs, photo-excitation provides enough free carriers to screen the polarization bound charges and quench the resulting local electrostatic fields in each layer.

FIG. 3. Excited carriers density for a) (1 1) SLs, b) (3 3) SLs, | | c) (5 5) SLs, computed at | n ph = 0 4 e/f.u. . (d) Electrostatic potential in dark (blue) and under illumination (orange) for the (4 4) SLs. | Plots in a)- c) are the difference between the macro-average density under illumination and in dark.

STO thus opposes the polarizing field and leads to the decrease in the polarization in STO and to the observed compression of this specific layer. On the other hand, the light-induced field in PTO opposes the depolarizing field there too, which then results in an increase of the electrical polarization, therefore generating the observed expansion of the PTO layer (see the Supplemental Materials [44] for a schematic). We further elucidate this

Previous studies based on Landau-GuinzburgDevonshire (LGD) theory have shown that, when the direction of the polarization is locked along the z axis (the argument is also true for the other axes), the tetragonality of the SLs can serve as a probe of polarization of the SLs [24, 25]. It was established that the polarization is proportional to an exponent of the tetragonality η T = ( c -c p ) /a with c and c p being the out-of-plane lattice constant of the polar SLs and the corresponding paraelectric SLs, and a the in-plane lattice parameter of the SLs. Earlier works on single domain PTO/STO SLs in dark have shown that P z ∝ η 1 / 2 T [25], however, in a recent work, it was found that α differs from 1 2 / when free charges screening effect is considered. The study found that the values of

α depends on the SLs, its structure and on the pump used for the excitation. Typically, in presence of free charge, it was found that P z ∝ η 1 1 . T [35] for a polydomain (PbTiO ) /(SrTiO ) 3 8 3 3 SLs. In view of these previous works, we plot ln( P z ( n ph )) as function ln( η T ( n ph )) for the considered SLs. As shown on the top panel of Fig 4 (shown for (1 1), (3 3) and (5 5) SLs), there is a linear | | | relation between ln( P z ( n ph )) and ln( η T ( n ph )) - implying that P z ∝ η α T = [( c -c p ) /a ] α . We obtained α (0.65, 0.56, 0.40, 0.62 and 0.49 for the (1 1), (2 2), (3 3), (4 4) | | | | and (5 5) respectively) by taking the slope of a linear fit | into ln( P z ( n ph )) vs ln( η T ( n ph )) for each SLs.

$$P _ { z } \in \eta _ { T } ^ { \alpha } = \left [ ( c - c _ { p } ) / a \right ] ^ { \alpha }. \quad \quad ( 2 ) \quad \text{ an} \quad \text{no.}$$

Taking the logarithmic derivative of Eq 2, gives the following:

$$\frac { d P _ { z } } { P _ { z } } = \alpha \frac { d c } { ( c - c _ { p } ) }. \quad \begin{array} { c c c c } & & \stackrel { \circ } { \text{kn} } \\ & & \text{Hi} \\ & & \text{pu} \end{array}.$$

From Eq 3, it can be readily established that the strain dc/c in the out-of-plane lattice constant is given by

$$d c / c = \eta _ { 3 3 } = \frac { \eta _ { d } } { \alpha } \frac { d P _ { z } } { P _ { z } }. \quad \quad \ ( 4 )$$

where we define η d = ( c -c p ) /c . The scaling relation in Eq. 2 provides us with a simple analytical expression of the light induced strain η 33 as function of the light-induced relative change in the polarization. The estimated strains using Eq. 4, along with the DFT data, are shown on the right panel in Fig 4 for each of the SLs. It is in excellent agreement with the DFT data. Also we plot the estimated strain using reported values of α from Refs. [25, 35] for comparison.

Although not the focus of this study, we would like to comment on the obtained values of α . As found in Ref. [35], it would be expected that α depends on the SLs. Also one would expect the values of α to deviate from 1 / 2 as we find in our study. However, with the low excited carriers concentration considered in this study and the simple structure of the SLs, we do not expect a drastic deviation from 1 / 2, consistent with the range of values 0 4 . -0 65 found for the considered ( . n n | ) corresponding to a maximum deviation of 0 15 from 0 5. . . As discussed in Ref. [35], the physical significance of α is unclear to us. It would interesting to investigate in more details the effect of free charge on α and its behavior as function of the SLs.

In summary, we study, using first-principle calculations, light-induced deformation in single domain ( n n | ) PTO/STO SLs ( n = 1 -5). We show that, for very thin SLs ( n ≤ 2), free carriers are 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), the free charges are localized in the interfaces regions. They produce a field that opposes existing depolarizing fields in the SLs and thereby induces a contraction in the STO layers and oppositely an expansion in the PTO layers and leads to an overall expansion of the SLs. These show that the photostrictive behavior of PTO/STO SLs can be tuned through its thickness. Our findings show a promising route to controlling and tuning the light-induced strain in PTO/STO SLs by simply changing its thickness making them appealing for potential applications using light as a handle.

The work is supported by the ARO Grants No. W911NF-21-2-0162 (ETHOS) and W911NF-21-1-0113, and the Vannevar Bush Faculty Fellowship (VBFF) grant no. N00014-20-1-2834 from the Department of Defense. C.P. thanks the support from a public grant overseen by the French Agence Nationale de la Recherche under grant agreement no. ANR-21-CE24-0032. We also acknowledge the computational support from the Arkansas High Performance Computing Center (AHPCC) for computational resources.

- [1] B. Kundys, Photostrictive materials, Appl. Phys. Rev. 2 , 011301 (2015).
- [2] C. Chen and Z. Yi, Photostrictive effect: characterization techniques, materials, and applications, Advanced Functional Materials 31 , 2010706 (2021).
- [3] H. Finkelmann, E. Nishikawa, G. Pereira, and M. Warner, A new opto-mechanical effect in solids, Phys. Rev. Lett. 87 , 015501 (2001).
- [4] R. Haleoot, C. Paillard, T. P. Kaloni, M. Mehboudi, B. Xu, L. Bellaiche, and S. Barraza-Lopez, Photostrictive two-dimensional materials in the monochalcogenide family, Phys. Rev. Lett. 118 , 227401 (2017).
- [5] C. Paillard, S. Prosandeev, and L. Bellaiche, Ab initio approach to photostriction in classical ferroelectric materials, Phys. Rev. B 96 , 045205 (2017).
- [6] C. Paillard, X. Bai, I. C. Infante, M. Guennou, G. Geneste, M. Alexe, J. Kreisel, and B. Dkhil, Photovoltaics with ferroelectrics: Current status and beyond, Adv. Mater. 28 , 5153 (2016).
- [7] S. Matzen, L. Guillemot, T. Maroutian, S. K. Patel, H. Wen, A. D. DiChiara, G. Agnus, O. G. Shpyrko, E. E. Fullerton, D. Ravelosona, et al. , Tuning ultrafast photoinduced strain in ferroelectric-based devices, Advanced Electronic Materials 5 , 1800709 (2019).
- [8] L. Gao, C. Paillard, and L. Bellaiche, Photoinduced control of ferroelectricity in hybrid-improper ferroelectric superlattices, Phys. Rev. B 107 , 104109 (2023).
- [9] B. Kundys, M. Viret, C. Meny, V. Da Costa, D. Colson, and B. Doudin, Wavelength dependence of photoinduced deformation in BiFeO , Phys. Rev. B 3 85 , 092301 (2012).
- [10] B. Kundys, M. Viret, D. Colson, and D. O. Kundys, Light-induced size changes in BiFeO 3 crystals, Nature materials 9 , 803 (2010).
- [11] H. Wen, P. Chen, M. P. Cosgriff, D. A. Walko, J. H. Lee, C. Adamo, R. D. Schaller, J. F. Ihlefeld, E. M. Dufresne, D. G. Schlom, et al. , Electronic origin of ultrafast photoinduced strain in BiFeO 3 , Phys. Rev. Lett. 110 , 037601

a)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-098b10b6-35a7-4c89-bc29-839b2b82259e.jpg)

d)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-81c5d43c-48c1-49ff-a64a-b994e64f73f4.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-106f05e9-163e-4e96-a23d-05dabecc1786.jpg)

b)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-c375f697-3e45-4d52-951c-c7996ef7d643.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-e0a67293-4c6f-4cdd-9783-5c515c003654.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-2bfe6580-c870-49cb-8b69-812b0ff22589.jpg)

e)

(%)

FIG. 4. Top panel a) to c) plot of ln( P ) vs ln( η T ) for the (1 1), (3 3) and (5 5) SLs respectively. | | | Bottom panel d) to f) plot of the DFT data and fit using eq 4 for different values of α for the (1 1), (3 3) and (5 5) SLs respectively. | | |

(2013).

- [12] D. Schick, M. Herzog, H. Wen, P. Chen, C. Adamo, P. Gaal, D. G. Schlom, P. G. Evans, Y. Li, and M. Bargheer, Localized excited charge carriers generate ultrafast inhomogeneous strain in the multiferroic BiFeO , Phys. Rev. Lett. 3 112 , 097602 (2014).

on adsorption-induced and photoinduced stresses in microelectromechanical systems devices, Journal of Vacuum Science &amp; Technology B: Microelectronics and Nanometer Structures Processing, Measurement, and Phenomena 19 , 1173 (2001).

- [13] C. Paillard, B. Xu, B. Dkhil, G. Geneste, and L. Bellaiche, Photostriction in ferroelectrics from density functional theory, Phys. Rev. Lett. 116 , 247401 (2016).
- [14] R. Gu, V. Juv´ e, C. Laulh´ e, H. Bouyanfif, G. Vaudel, A. Poirier, B. Dkhil, P. Hollander, C. Paillard, M. C. Weber, et al. , Temporal and spatial tracking of ultrafast light-induced strain and polarization modulation in a ferroelectric thin film, Science Advances 9 , eadi1160 (2023).
- [15] M. Nakano and T. Ikeda, Photomechanics: directed bending of a polymer film by light, Nature 425 , 145 (2003).
- [16] J. R. Buschert and R. Colella, Photostriction effect in silicon observed by time-resolved x-ray diffraction, Solid state communications 80 , 419 (1991).
- [17] Y. Zhou, L. You, S. Wang, Z. Ku, H. Fan, D. Schmidt, A. Rusydi, L. Chang, L. Wang, P. Ren, et al. , Giant photostriction in organic-inorganic lead halide perovskites, Nature communications 7 , 11193 (2016).
- [18] C. Paillard and L. Bellaiche, Light: A new handle to control the structure of cesium lead iodide, Phys. Rev. B 107 , 054107 (2023).
- [19] B. Peng, D. Bennett, I. Bravi´, and B. Monserrat, Tunc able photostriction of halide perovskites through energy dependent photoexcitation, Phys. Rev. Mat. 6 , L082401 (2022).
- [20] P. Datskos, S. Rajic, M. Sepaniak, N. Lavrik, C. Tipple, L. Senesac, and I. Datskou, Chemical detection based
- [21] H.-P. Wang, D.-H. Lien, M.-L. Tsai, C.-A. Lin, H.-C. Chang, K.-Y. Lai, and J.-H. He, Photon management in nanostructured solar cells, Journal of Materials Chemistry C 2 , 3144 (2014).
- [22] P. Poosanaas, K. Tonooka, and K. Uchino, Photostrictive actuators, Mechatronics 10 , 467 (2000).
- [23] D. Sun and L. Tong, Modeling of wireless remote shape control for beams using nonlinear photostrictive actuators, International Journal of Solids and Structures 44 , 672 (2007).
- [24] 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).
- [25] 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).
- [26] 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).
- [27] 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).

- [28] V. R. Cooper, K. Johnston, and K. M. Rabe, Polarization enhancement in short period superlattices via interfacial intermixing, Phys. Rev. B 76 , 020103 (2007).
- [29] 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).
- [30] A. Yadav, C. Nelson, S. Hsu, Z. Hong, J. Clarkson, C. Schlep¨tz, u A. Damodaran, P. Shafer, E. Arenholz, L. Dedon, D. Chen, A. Vishwanath, A. M. Minor, L. Q. Chen, J. F. Scott, L. W. Martin, and R. Ramesh, Observation of polar vortices in oxide superlattices, Nature 530 , 198 (2016).
- [31] 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).
- [32] 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).
- [33] 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).
- [34] 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).
- [35] 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).
- [36] 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).
- [37] 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).
- [38] C. Paillard, E. Torun, L. Wirtz, J. ´ ˜iguez, and L. BelIn laiche, Photoinduced phase transitions in ferroelectrics, Phys. Rev. Lett. 123 , 087601 (2019).
- [39] N. Pertsev, A. Zembilgotov, and A. Tagantsev, Equilib-
- rium states and phase transitions in epitaxial ferroelectric thin films, Ferroelectrics 223 , 79 (1999).
- [40] 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).
- [41] C. Dansou, C. Paillard, and L. Bellaiche, Tuning photostriction in (PbTiO ) 3 n /(SrTiO ) 3 m superlattices via chemical composition: An ab-initio study, Submitted to Physical Review B .
- [42] W. Zhong, R. King-Smith, and D. Vanderbilt, Giant loto splittings in perovskite ferroelectrics, Phys. Rev. Lett. 72 , 3618 (1994).
- [43] C. E. Dreyer, S. Coh, and M. Stengel, Nonadiabatic born effective charges in metals and the drude weight, Physical review letters 128 , 095901 (2022).
- [44] See supplemetal material at [url] for the data on the other superlattices., .
- [45] J. Junquera, M. H. Cohen, and K. M. Rabe, Nanoscale smoothing and the analysis of interfacial charge and dipolar densities, J. Phys: Cond. Mat. 19 , 213203 (2007).

## Supplemental Materials for 'Tailoring photostriction via superlattices engineering'

POLARIZATION IN THE (2 2) | AND (4 4) | SLS AS FUNCTION OF n ph .

a)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-da9a4e4b-c6e6-469b-b298-78974204320e.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-66f9820a-e5a3-421b-93e0-322894f937e4.jpg)

b)

c)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-6a057859-1146-4c85-8f83-e8bbea72f034.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-319c383e-8f4d-471c-b172-f31086a8c488.jpg)

d)

FIG. 5. Polarization as a function of n ph in the ( n n | )SLs. a) (2 2), b) (4 4) SLs respectively. | | Projection to five-atom unit cell of the polarization in the c) (2 2), d) (4 4) SLs respectively. | |

FIG. 6. Schematic diagram of the action of light in the ( n n | ) SLs. ( /vector E , P p /vector p ) and ( /vector E , P s /vector s ) are the coupled field and polarization in PTO and STO layer in dark conditions, respectively. /vector E , E l p /vector l s are the light-induced field in PTO and STO layer, respectively.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-783cb55e-886b-4ba8-8235-4ee5536dce39.jpg)

## PHOTOEXCITED CHARGES DENSITY IN (2 2) | AND (4 4) | SLS COMPUTED AT n ph = 0 4 .

FIG. 7. Excited carriers density for a) (2 2) SLs, b) (4 4) SLs, computed at | | n ph = 0 4 e/f.u. .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-fd7e80f2-6dda-43c3-8e5d-5b80f18b5f53.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-7e7be7cf-11a0-4f48-a575-cced64be9070.jpg)

(%)

33

/s951

/s8722

/s8722

/s8722

/s8722

a)

0.00

/s48

/s46

/s48

/s50

/s48

/s48

/s48

/s46

/s46

/s46

/s48

/s48

/s48

/s52

/s54

/s56

0.00

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-dafdc4aa-0199-4ce0-8926-fb3a47e14752.jpg)

/s945

FIG. 8. Values of α as function of the ( n n | ) SLs

## MODEL RESULTS FOR THE (2 2) | AND (4 4) | SLS

b)

DFT data

Fit

α

= 0.56

Fit

α

= 0.50

Fit

α

= 1.10

0.05

0.10

0.15

0.20

n ph

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-629ae148-5e44-48c6-a833-7e0e16564952.jpg)

(e/f.u)

c)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-c099041f-f37f-4729-98c2-2ed6588111c8.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-76eee17a-119e-4b4b-98cf-2b4bcf0c2c5f.jpg)

d)

FIG. 9. a), b) plot of the DFT data and fit using Eq. 4 in the main text for the (2 2), (4 4) SLs respectively. | | c), d) plot of ln( P ) vs ln( η T ) for the (2 2), (4 4) SLs respectively. | |

(%)

1.2

## Supplemental Materials for 'Tailoring photostriction via superlattices engineering'

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.

## POLARIZATION IN THE (2 2) | AND (4 4) | SLS AS FUNCTION OF n ph .

a)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-ff98981a-b554-4069-b977-cb798beca9f2.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-bb2b0a29-f4dc-49a1-9b4f-79aeab3b21ba.jpg)

c)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-103166b3-cdf2-4162-9887-6fc1fd4bb282.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-816a1048-ab7d-4e8f-adae-f319e5243417.jpg)

d)

FIG. 1. Polarization as a function of n ph in the ( n n | )SLs. a) (2 2), b) (4 4) SLs respectively. | | Projection to five-atom unit cell of the polarization in the c) (2 2), d) (4 4) SLs respectively. | |

FIG. 2. Schematic diagram of the action of light in the ( n n | ) SLs. ( /vector E , P p /vector p ) and ( /vector E , P s /vector s ) are the coupled field and polarization in PTO and STO layer in dark conditions, respectively. /vector E , E l p /vector l s are the light-induced field in PTO and STO layer, respectively.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-783cb55e-886b-4ba8-8235-4ee5536dce39.jpg)

## PHOTOEXCITED CHARGES DENSITY IN (2 2) | AND (4 4) | SLS COMPUTED AT n ph = 0 4 .

FIG. 3. Excited carriers density for a) (2 2) SLs, b) (4 4) SLs, computed at | | n ph = 0 4 e/f.u. .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-fd7e80f2-6dda-43c3-8e5d-5b80f18b5f53.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-7e7be7cf-11a0-4f48-a575-cced64be9070.jpg)

(%)

33

/s951

/s8722

/s8722

/s8722

/s8722

a)

0.00

/s48

/s46

/s48

/s50

/s48

/s48

/s48

/s46

/s46

/s46

/s48

/s48

/s48

/s52

/s54

/s56

0.00

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-bb3e0967-ecc9-49b9-914c-f662df944a49.jpg)

/s945

FIG. 4. Values of α as function of the ( n n | ) SLs

## MODEL RESULTS FOR THE (2 2) | AND (4 4) | SLS

b)

DFT data

Fit

α

= 0.56

Fit

α

= 0.50

Fit

α

= 1.10

0.05

0.10

0.15

0.20

n ph

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-4895cddc-5353-419b-acd3-ee387b09b6eb.jpg)

(e/f.u)

c)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-c099041f-f37f-4729-98c2-2ed6588111c8.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00101v1-figure-28a539e9-cb5d-4c66-92d6-50f787699d45.jpg)

d)

FIG. 5. a), b) plot of the DFT data and fit using Eq. 4 in the main text for the (2 2), (4 4) SLs respectively. | | c), d) plot of ln( P ) vs ln( η T ) for the (2 2), (4 4) SLs respectively. | |

(%)

1.2
Carmel Dansou Charles Paillard Laurent Bellaiche
Tailoring photostriction via superlattices engineering
We report systematic first-principles investigation of light-induced
mechanical deformations in monodomain (PbTiO$_{3}$)$_{n}$/(SrTiO$_{3}$)$_{n}$
superlattices ($n=1-5$). We reveal that photostriction in these
heterostructures quantitatively and qualitatively depends on the chemical
period $n$. Specifically, we show that by changing the chemical period, we can
induce $\textit {positive}$ or $\textit {negative}$ photostriction. We also
present a simple analytical model to account for the calculated deformations.
Our findings indicate that superlattices architectures may be key to design
novel optomechanical applications.