Access the full text.
Sign up today, get DeepDyve free for 14 days.
Loading next page...
/lp/arxiv-cornell-university/modeling-temperature-frequency-and-strain-effects-on-the-linear-QKyoTR4XA5
References
References for this paper are not available at this time. We will be adding them shortly, thank you for your patience.
- eISSN
- ARCH-3341
- Publisher site
- See Article on Publisher Site
Abstract
An electro-optic modulator offers the function of modulating the propagation of light in a material with electric field and enables seamless connection between electronics-based computing and photonics-based communication. The search for materials with large electro-optic coefficients and low optical loss is critical to increase the efficiency and minimize the size of electro-optic devices. We present a semi-empirical method to compute the electro-optic coefficients of ferroelectric materials by combining first-principles density-functional theory calculations with Landau- Devonshire phenomenological modeling. We apply the method to study the electro-optic constants, also called Pockels coefficients, of three paradigmatic ferroelectric oxides: BaTiO , LiNbO , and 3 3 LiTaO . We present their temperature-, frequency- and strain-dependent electro-optic tensors calculated using our method. The predicted electro-optic constants agree with the experimental results, where available, and provide benchmarks for experimental verification. Page 1 of 31 Introduction AMO -type ferroelectric oxides offer strong coupling between electrical, thermal, and optical properties, and enable novel applications that leverage the coupled phenomena. They are currently used in nonvolatile memories, actuators, transducers, and electro-optic (EO) devices, owing to their excellent dielectric, piezoelectric and pyroelectric properties, and optical response. [1–4] For optical applications, ferroelectric oxide perovskites exhibit large EO coefficients with low optical loss, and are the materials of choice for low-power electro-optic devices. Since the 1970s, EO modulators based on LiNbO have been used widely in fiber-optic systems due to its good linear EO or Pockels effect (𝑟𝑟 : 32 pm/V) and high transparency over a large range of wavelength. [5] Thin film deposition of EO oxides, characterization of their optical response, [6] and fabrication of optical devices have undergone significant refinement since the 1990s. [7–10] In line with this development, there is a growing interest in achieving epitaxially grown ferroelectric thin films integrated on silicon-based chips for optical waveguide modulators. [11–16] Among different ferroelectric oxides, BaTiO , LiNbO , and LiTaO have been investigated 3 3 3 intensively for on-chip EO applications due to their sizable linear EO effect (tetragonal BaTiO on Si 𝑟𝑟 : 105 pm/V, 𝑟𝑟 : 148 pm/V [17] and LiNbO on Si 𝑟𝑟 : 17.6 pm/V [18]). However, 42 𝑒𝑒𝑒𝑒𝑒𝑒 33 modeling methods for the EO response of these ferroelectric materials as a function of temperature, frequency, strain and electric dipole orderings has not been well-established. [19–23] In fact, EO effects are shown to be sensitive to the microstructure, and an accurate assessment of this intrinsic property requires single crystals or high-quality thin films, which are not easily accessible or prepared. Therefore, theoretical prediction of the nonlinear optical properties of crystalline materials along with the effect of various experimental conditions, such as strain and temperature, can help to establish performance limits for subsequent experimental verification. In the past decade, sustained efforts on theoretical investigations of nonlinear optical phenomena in oxide perovskites have resulted in accurate methods for predicting these properties. DiDomenico and Wemple revealed the importance of oxygen octahedra in perovskites on their optical properties. [24,25] Ghosez and co-workers calculated the optical susceptibilities, Raman efficiencies, and electro-optic tensors based on density functional perturbation theory. [26–29] More recently, Hamze et al., Qiu et al., and Paillard et al. studied the effect of strain on the electro- optic tensor. [19–23,30] Furthermore, with the ability to prepare atomically precise Page 2 of 31 heterostructures and superlattices, it is of both scientific and practical importance to understand the mechanism of the EO effect in these complex systems and predict the EO coefficients reliably. [31,32] While first-principles calculation methods used in previous studies are effective in predicting the EO effect for single crystals, modeling EO effects in superlattices and multilayers presents a formidable challenge. The periodicity of the superlattices and multilayers, which span few nm to few 10s of nm, and breadth of phase space in terms of materials and periodicities needed to model EO effects and identify high efficiency structures make first-principles methods computationally expensive and impractical. Phenomenological models, such as those based on Landau-Devonshire theory [33] enable fast, accurate, and highly scalable calculations of the functional properties of complex structures. It is important to note that Landau-Devonshire model uses input from experimental results or first- principles calculations to fit the coefficients used in the model. Hence, the accuracy of Landau- Devonshire expansion coefficients in subsequent estimation of functional properties is determined by these inputs. For a multicomponent system, such as superlattices and multilayers, one can simulate their physical properties by summing up the thermodynamic free energies of each component as a function of strain, electric fields, and their gradients. [34,35] This approach has been extensively applied for the simulation of dielectric and piezoelectric responses of ferroelectric materials and multilayer heterostructures. [36–38] The objective of this study is to establish a semi-empirical model to simulate the EO behavior of perovskite ferroelectrics. This model uses the phenomenological Landau-Devonshire model with parameters obtained from first-principles calculations to improve the scalability of EO calculations for complex structures without compromising on speed and accuracy. We show that this model can be applied to prototypical ferroelectric oxides such as LiNbO , LiTaO , and BaTiO . We 3 3 3 obtained the free-energy landscape associated with the transition between ferroelectric and paraelectric phases using density-functional theory (DFT) calculations. We extracted Landau- Devonshire coefficients using a polynomial fitting to the energy landscape and calculated the EO coefficients at room temperature for these three prototypical ferroelectric oxides. We find that our model can predict EO coefficients that have good agreement with experimental results, wherever available. Using this model, we have calculated the temperature-dependence of the EO coefficients Page 3 of 31 for LiNbO and BaTiO and find them to be within 30% of experimental results for most cases. 3 3 Moreover, the strain effect on the EO coefficient is discussed in the range of -5 to 5% misfit strain for BaTiO . Our model is able to capture the ferroelectric to paraelectric phase transition, which is associated with a divergence of the EO tensor. Methods Density-Functional Theory Calculations The landscapes of free energy for the different AMO oxides were computed using DFT as implemented in Vienna Ab-initio Simulation Package (VASP). [39] We used projector augmented-wave (PAW) potentials. [40] In general, the accuracy in the estimation of ferroelectric properties is sensitive to the adopted exchange-correlation functionals such as the local density approximation (LDA), [41] and the semi-local generalized gradient approximation (GGA) in the standard from of Perdew-Burke-Ernzerhof (PBE) [42]. GGA is known to suffer from the so-called super-tetragonality error, which significantly overestimates the structural distortion in conventional perovskite ferroelectrics [43]. For the most-studied oxide ferroelectrics, BaTiO and PbTiO , the lattice distortion, spontaneous polarization, and lattice dynamics predicted by LDA functional agree well with the experimental results. [44] Therefore, we chose LDA to describe the electronic exchange-correlation interactions. We have considered three paradigmatic AMO oxides, BaTiO (P4mm, Amm2, R3m), LiNbO (R3c), and LiTaO (R3c), for determining structural 3 3 3 transition and ferroelectric polarization. A cutoff energy of 700 eV was used to determine the number of planewave basis sets in the calculations. We used Γ-centered 10×10×10 k-points mesh for sampling the Brillouin zone of BaTiO and 10×10×4 k-points mesh for LiNbO and LiTaO . 3 3 3 −3 The crystal structures were fully optimized until residual forces were less than 10 eV/Å. The spontaneous polarization induced by the polar soft-phonon modes was calculated based on the modern theory of polarization [45], which is a sum over the contribution from the ionic and electronic charges. Symmetry and distortion-mode analyses were conducted using programs from the Bilbao crystallographic server. [46] The intermediate images corresponding to soft-phonon distortion were interpolated using the ISOTROPY software suite. [47] Density-Functional Perturbation Theory Calculations Page 4 of 31 We also calculated the EO tensor of the ground-state R3m phase of BaTiO purely from first- principles as a comparison to that obtained using the Landau-Devonshire model. The theoretical framework developed by Veithen et al. [26–29] for the computation of EO response under a static or low-frequency electric field perturbation has been implemented in the ABINIT software package. [48,49] Teter extended norm-conserving pseudopotentials [19–23] for BaTiO system were used for these calculations and the exchange-correlation interactions were described within LDA. [50] We used 12×12×12 k-points mesh and 55 Hartree cutoff energy for all the calculations. To study the effect of strain on EO response, we adopted the same strategy described by Fredrickson et al. [21] Varying epitaxial strains between –2 to +2 % with the negative values denoting compressive strain were applied to a and b lattice constants. The optimal c lattice constant for a given epitaxial strain was calculated using the elastic constants of tetragonal BaTiO (C = 3 11 222 GPa, C = 108 GPa, C = 111 GPa, C = 151 GPa). [51] The ionic positions in the strained 12 13 33 −5 lattices were optimized until the forces were less than 1×10 eV/A. Landau-Devonshire Model Landau phenomenological theory is widely used to describe phase transitions and temperature dependence of physical properties of ferroelectrics. [52] Here, we use Helmholtz free energy to describe the thermodynamics due to the convenience in choosing the internal variables: polarization (P) and strain (S) as independent variables, whereas the electric field (E) and the stress are external applied variables. The free energy and free energy density in this article refer to Helmholtz free energy and Helmholtz free energy density, unless noted otherwise. The Helmholtz free energy density (f ) of a ferroelectric system under no external field can be written as an expansion of the order parameter - the polarization (P), as: [53] 2 2 2 4 4 4 2 2 2 2 2 2 𝑓𝑓 = 𝑎𝑎 (𝑃𝑃 +𝑃𝑃 +𝑃𝑃 ) +𝑎𝑎 (𝑃𝑃 +𝑃𝑃 +𝑃𝑃 ) +𝑎𝑎 (𝑃𝑃 𝑃𝑃 +𝑃𝑃 𝑃𝑃 +𝑃𝑃 𝑃𝑃 ) + 0 1 1 2 3 11 1 2 3 12 1 2 1 3 2 3 6 6 6 4 2 2 4 2 2 4 2 2 ( ) ( ) ( ) 𝑎𝑎 (𝑃𝑃 +𝑃𝑃 +𝑃𝑃 ) +𝑎𝑎 (𝑃𝑃 𝑃𝑃 +𝑃𝑃 +𝑃𝑃 𝑃𝑃 +𝑃𝑃 +𝑃𝑃 𝑃𝑃 +𝑃𝑃 ) + 111 1 2 3 112 1 2 3 2 1 3 3 1 2 2 2 2 𝑎𝑎 𝑃𝑃 𝑃𝑃 𝑃𝑃 , (1) 1 2 3 where the subscripts 1,2,3 refer to [100], [010], and [001] directions in the crystal, 𝑎𝑎 ,𝑎𝑎 ,𝑎𝑎 are 𝑖𝑖 𝑖𝑖𝑖𝑖 𝑖𝑖 the phenomenological Landau-Devonshire coefficients, and 𝑃𝑃 is the polarization along direction 𝑖𝑖 i. The temperature dependence of ferroelectricity is governed by the coefficient 𝑎𝑎 and it is defined as (2) 𝑎𝑎 = (𝑇𝑇 − T )/2𝜀𝜀 C. 1 0 0 Page 5 of 31 𝑖𝑖𝑖𝑖 The other coefficients are all assumed to be temperature independent. Here T and C are the Curie- Weiss temperature and constant above which the system transitions to a paraelectric state, and 𝜀𝜀 is the dielectric constant of free space, respectively. We set T to be 388 K in the entire simulation for BaTiO and 1480 K and 950 K for LiNbO and LiTaO , which were observed from 3 3 3 experiments. [5,54–56] Classical Landau theory ignores the temperature effect on the higher-order coefficients in the expansion. Nevertheless, it is shown that the higher order terms are actually temperature dependent. [57] We include temperature effects in our calculations, and for simplicity, we only consider the temperature dependent 𝑎𝑎 in this work. The effect of temperature-dependent high-order terms, such as 𝑎𝑎 , on EO responses will be the target of future work. Model Fitting and Parameters The ferroelectric transition from a centrosymmetric reference can be expressed as the result of ionic displacements along a specific direction with charge separation leading to a net electrical dipole moment. [58] By interpolating the ionic displacements from a centrosymmetric structure to a polar phase, the energy as a function of ionic displacements can be mapped using DFT calculations. As has been shown recently by Paoletta and Demkov [19], phonons causing the ionic displacements will in turn alter the electronic energy of the system, and this is the origin of electron-phonon interactions under the adiabatic approximation. That is to say, our DFT calculations for free energy landscape of each distortion mode also reflects the electron-phonon interactions. For the subsequent Landau-Devonshire fittings, we have converted the ionic displacements into spontaneous polarization based on the modern theory of polarization. [59] The landscape of the change in free energy density (J/m ) for the three ferroelectric phase transitions from paraelectric BaTiO (P4/mmm) as a function of the electric polarization are shown in Figure 1. By fitting the Landau-Devonshire expansion to the change in energy density with polarization, quadratic and higher-order coefficients of the polynomial can be derived for ferroelectric transition along [001], [011], and [111] direction for tetragonal (P4mm), orthorhombic (Amm2), and rhombohedral (R3m) structures, respectively. The free energy density with respect to the 2 2 2 2 2 � � polarization 𝑃𝑃 = 𝑃𝑃 , 𝑃𝑃 = 𝑃𝑃 +𝑃𝑃 , and 𝑃𝑃 = 𝑃𝑃 +𝑃𝑃 +𝑃𝑃 can be described by the 001 3 011 111 2 3 1 2 3 following equations, respectively: [60] 2 4 6 𝑓𝑓 = 𝑎𝑎 𝑃𝑃 +𝑎𝑎 𝑃𝑃 +𝑎𝑎 𝑃𝑃 , (3.a) 001 1 001 11 001 111 001 2 𝑂𝑂 4 𝑂𝑂 6 𝑓𝑓 = 𝑎𝑎 𝑃𝑃 +𝑎𝑎 𝑃𝑃 +𝑎𝑎 𝑃𝑃 , (3.b) 011 1 011 11 011 111 011 Page 6 of 31 2 𝑅𝑅 4 𝑅𝑅 6 𝑓𝑓 = 𝑎𝑎 𝑃𝑃 +𝑎𝑎 𝑃𝑃 +𝑎𝑎 𝑃𝑃 , (3.c) 111 1 111 11 111 111 111 where the superscripts O and R indicate the orthorhombic and rhombohedral phase for BaTiO , 1 1 1 1 1 𝑂𝑂 𝑂𝑂 𝑅𝑅 𝑅𝑅 𝑎𝑎 = 𝑎𝑎 + 𝑎𝑎 , 𝑎𝑎 = (𝑎𝑎 +𝑎𝑎 ), 𝑎𝑎 = (𝑎𝑎 +𝑎𝑎 ) , and 𝑎𝑎 = (3𝑎𝑎 + 11 11 12 111 111 112 11 11 12 111 111 2 4 4 3 27 6𝑎𝑎 +𝑎𝑎 ). We used the “Curve Fitting Toolbox” in MATLAB to fit the free energy density 112 123 curves obtained from DFT calculations. We fitted the energy density landscape of tetragonal BaTiO with the eqn. (3.a) to obtain the 𝑎𝑎 , 𝑎𝑎 , and 𝑎𝑎 . To get the 𝑎𝑎 and 𝑎𝑎 , the 1 11 111 12 112 orthorhombic energy density is fitted to the eqn. (3.b). 𝑎𝑎 is derived by fitting the energy density of rhombohedral phase using the eqn. (3.c) using all the other parameters obtained from the previous steps. Then all the parameters are manually tuned to minimize the coefficients of determination (R ) of three equations (3, a-c) by slightly changing only one parameter at a time while fixing all remaining parameters. Thus, the whole sets of the Landau-Devonshire coefficients can be derived. For the case of LiNbO and LiTaO , we applied the same procedure as BaTiO but simplified it to 3 3 3 [001] direction since we are only interested in the most intense EO tensor component -- 𝑟𝑟 . [61] Therefore, we calculated the free energy density curve for R3c LiNbO and LiTaO as a function 3 3 of polarization along [001] direction. The Landau-Devonshire coefficients 𝑎𝑎 and 𝑎𝑎 were 1 11 obtained by fitting the energy density curve to the eqn. (3.a). Then, we applied the strain and electrostrictive energy terms to the Landau-Devonshire model to investigate the strain-induced phase changes in BaTiO . The free energy density f of the thin film as a function of polarization and misfit strain 𝑆𝑆 = (𝑎𝑎 −𝑎𝑎 )/𝑎𝑎 , [36] where 𝑎𝑎 is the substrate 𝑚𝑚 𝑠𝑠 𝑒𝑒 𝑠𝑠 𝑠𝑠 lattice parameter and 𝑎𝑎 is the lattice constant of the film in its bulk form, is given by: [62] 𝑒𝑒 ∗ 2 2 ∗ 2 ∗ 4 4 ∗ 4 ∗ 2 2 ∗ 2 2 𝑓𝑓 = 𝑎𝑎 (𝑃𝑃 +𝑃𝑃 ) +𝑎𝑎 𝑃𝑃 +𝑎𝑎 (𝑃𝑃 +𝑃𝑃 ) +𝑎𝑎 P +𝑎𝑎 𝑃𝑃 𝑃𝑃 +𝑎𝑎 (𝑃𝑃 𝑃𝑃 + 1 1 2 3 3 11 1 2 33 3 12 1 2 13 1 3 2 2 6 6 6 4 2 2 4 2 2 ( ) ( ) 𝑃𝑃 𝑃𝑃 ) +𝑎𝑎 (𝑃𝑃 +𝑃𝑃 +𝑃𝑃 ) +𝑎𝑎 (𝑃𝑃 𝑃𝑃 +𝑃𝑃 +𝑃𝑃 𝑃𝑃 +𝑃𝑃 + 2 3 111 1 2 3 112 1 2 3 2 1 3 𝑆𝑆 𝑚𝑚 4 2 2 2 2 2 𝑃𝑃 (𝑃𝑃 +𝑃𝑃 )) +𝑎𝑎 𝑃𝑃 𝑃𝑃 𝑃𝑃 + , 3 1 2 123 1 2 3 𝑠𝑠 +𝑠𝑠 11 12 (4) 𝑄𝑄 +𝑄𝑄 11 12 where 𝑎𝑎 = 𝑎𝑎 − 𝑆𝑆 , (4.a) 1 1 𝑚𝑚 𝑠𝑠 +𝑠𝑠 11 12 2𝑄𝑄 𝑎𝑎 = 𝑎𝑎 − 𝑆𝑆 , (4.b) 3 1 𝑚𝑚 𝑠𝑠 +𝑠𝑠 11 12 Page 7 of 31 2 2 1�𝑄𝑄 +𝑄𝑄 �𝑠𝑠 −2𝑄𝑄 𝑄𝑄 𝑠𝑠 11 12 11 11 12 12 𝑎𝑎 = 𝑎𝑎 + , (4.c) 11 11 2 2 2 𝑠𝑠 −𝑠𝑠 11 12 𝑄𝑄 𝑎𝑎 = 𝑎𝑎 − , (4.d) 𝑠𝑠 +𝑠𝑠 11 12 2 2 2 �𝑄𝑄 +𝑄𝑄 �𝑠𝑠 −2𝑄𝑄 𝑄𝑄 𝑠𝑠 𝑄𝑄 11 12 12 11 12 11 44 𝑎𝑎 = 𝑎𝑎 − + , (4.e) 12 12 2 2 𝑠𝑠 −𝑠𝑠 2𝑠𝑠 11 12 𝑄𝑄 (𝑄𝑄 +𝑄𝑄 ) 12 11 12 𝑎𝑎 = 𝑎𝑎 + , (4.f) 13 12 𝑠𝑠 +𝑠𝑠 11 12 where 𝑄𝑄 are the electrostriction coefficients and 𝑠𝑠 are the elastic compliances. The 𝑄𝑄 and 𝑠𝑠 𝑖𝑖𝑖𝑖 𝑖𝑖𝑖𝑖 𝑖𝑖𝑖𝑖 𝑖𝑖𝑖𝑖 values in Table 1 for BaTiO are taken from Ref. 63. [63] Table 1. Elastic compliance (s ) and ij electrostrictive coefficients (Q ) of ij BaTiO taken from Ref. 63. -12 2 s (10 m /N) 8.33 -12 2 s (10 m /N) -2.68 -12 2 s (10 m /N) 9.24 4 2 Q (m /c ) 0.10 4 2 Q (m /c ) -0.034 4 2 Q (m /c ) 0.029 The variation of the free energy density under external electric field is written as: ∆𝑓𝑓 = 𝑓𝑓 −𝐸𝐸 𝑃𝑃 −𝐸𝐸 𝑃𝑃 −𝐸𝐸 𝑃𝑃 , (5) 0 1 1 2 2 3 3 where E , E , and E is the applied electric field along x, y, and z principal crystallographic 1 2 3 directions, respectively. The equilibrium configuration is determined by finding the minima of ∆𝑓𝑓 , 𝜕𝜕 ∆𝑒𝑒 where we shall have = 0 . Then, the electric field E as a function of polarization can be 𝜕𝜕𝜕𝜕 determined by: 𝜕𝜕 𝑒𝑒 𝐸𝐸 = , (6.a) 𝜕𝜕 𝜕𝜕 𝜕𝜕 𝑒𝑒 𝐸𝐸 = , (6.b) 𝜕𝜕 𝜕𝜕 𝜕𝜕 𝑒𝑒 𝐸𝐸 = . (6.c) 𝜕𝜕 𝜕𝜕 In this article, every time an external electric field is applied, we solve eqns. (6.a-c) to deduce the field-induced polarizations. Then the obtained polarizations are applied to solve the corresponding Page 8 of 31 eqns. (8-9) in the following paragraphs. In this way, it will always maintain the thermodynamic 𝜕𝜕 ∆𝑒𝑒 equilibrium: = 0. 𝜕𝜕𝜕𝜕 The dielectric tensor 𝜀𝜀 can be defined in terms of the first-order derivative of polarization with 𝑖𝑖𝑖𝑖 respect to the external electric field. Here, we summarize the derived dielectric constants for tetragonal (𝑃𝑃 = 𝑃𝑃 = 0,𝑃𝑃 ≠ 0) and orthorhombic (orthorhombic, 𝑃𝑃 = 𝑃𝑃 ≠ 0,𝑃𝑃 = 0) BaTiO 1 2 3 1 2 3 phases in the box I and II. We don’t include the low-temperature rhombohedral phase since the rhombohedral phase is not accessible in the experiments though strain engineering at room temperature. 𝜀𝜀 = 𝜀𝜀 = , 11 22 ∗ ∗ 2 4 2𝑎𝑎 +2𝑎𝑎 𝜕𝜕 +2𝑎𝑎 𝜕𝜕 1 13 3 112 3 𝜀𝜀 = . 33 ∗ ∗ 2 4 2𝑎𝑎 +12𝑎𝑎 𝜕𝜕 +30𝑎𝑎 𝜕𝜕 3 33 3 3 Box I. Expression of spontaneous polarization and dielectric constants for tetragonal BaTiO . ∗ ∗ 2 ∗ 2 4 2 2 4 ( ) 𝑋𝑋 = 2𝑎𝑎 + 12𝑎𝑎 𝑃𝑃 + 2𝑎𝑎 𝑃𝑃 + 30𝑎𝑎 𝑃𝑃 +𝑎𝑎 12𝑃𝑃 𝑃𝑃 + 2𝑃𝑃 , 11 1 12 1 12 2 111 1 112 1 2 2 ∗ ∗ 2 ∗ 2 4 2 2 4 ( ) 𝑋𝑋 = 2𝑎𝑎 + 12𝑎𝑎 𝑃𝑃 + 2𝑎𝑎 𝑃𝑃 + 30𝑎𝑎 𝑃𝑃 +𝑎𝑎 12𝑃𝑃 𝑃𝑃 + 2𝑃𝑃 , 22 111 112 1 12 2 12 1 2 1 2 1 ∗ 3 3 𝑋𝑋 = 4𝑎𝑎 𝑃𝑃 𝑃𝑃 + 8𝑎𝑎 (𝑃𝑃 𝑃𝑃 +𝑃𝑃 𝑃𝑃 ), 12 12 1 2 112 1 2 1 2 ∗ ∗ 2 2 3 3 ( ) 𝑋𝑋 = 2𝑎𝑎 + 2𝑎𝑎 (𝑃𝑃 +𝑃𝑃 ) + 8𝑎𝑎 𝑃𝑃 𝑃𝑃 +𝑃𝑃 𝑃𝑃 , 33 3 13 1 2 112 1 2 1 2 𝑋𝑋 𝑋𝑋 1 22 11 𝜀𝜀 = , 𝜀𝜀 = , 𝜀𝜀 = . 11 2 22 2 33 𝑋𝑋 𝑋𝑋 −𝑋𝑋 𝑋𝑋 𝑋𝑋 −𝑋𝑋 𝑋𝑋 11 22 11 22 33 12 12 Box II. Expression of spontaneous polarization and dielectric constants for orthorhombic BaTiO . The propagation of light in a crystal is determined by the refractive index 𝑛𝑛 . The relation between 𝑖𝑖𝑖𝑖 the dielectric constant and the refractive index is 𝑛𝑛 = 𝜀𝜀 /𝜀𝜀 . The linear EO tensor 𝑟𝑟 describes 𝑖𝑖𝑖𝑖 0 𝑖𝑖 𝑖𝑖𝑖𝑖 the change of refractive index of a crystal in response to the applied electric field. Therefore, we write the linear EO tensor 𝑟𝑟 as first-order dependence of the inverse of refractive index square 𝑖𝑖 when a static or low-frequency modulating electric field 𝐸𝐸 is applied: 𝑖𝑖 −2 Δ�𝑛𝑛 �= 𝑟𝑟 𝐸𝐸 (7) 𝑖𝑖𝑖𝑖 𝑖𝑖 𝑖𝑖 The index ijk refers to the ij component of the refractive index and the dielectric tensor, for an Page 9 of 31 𝑖𝑖𝑖𝑖 𝑖𝑖𝑖𝑖 𝑖𝑖𝑖𝑖 applied electric field along the k direction. [20,28,64] For the following paragraphs, we denote the index ij with Voigt notations, i.e. 11 → 1, 22 → 2, 33 → 3, 23 → 4, 13 → 5, and 12 → 6. From eqn. (6), we have the electric field as a function of polarization, and equations in the box I and II give the dielectric constant as a function of polarization by substituting them into eqn. (7). Thus, given all the Landau-Devonshire coefficients obtained using polynomial fitting, the EO coefficients can be obtained. Here, we consider the case of tetragonal and orthorhombic phases of BaTiO , as examples. The EO tensors in the ferroelectric tetragonal P4mm phase of BaTiO have 3 3 three independent elements (Voigt notations), 𝑟𝑟 , 𝑟𝑟 , and 𝑟𝑟 . [64] 13 33 42 𝜀𝜀 (4𝑎𝑎 𝜕𝜕 +8𝑎𝑎 𝜕𝜕 ) 0 12 3 112 3 𝑟𝑟 = , (8.a) ∗ ∗ 2 4 2𝑎𝑎 +12𝑎𝑎 𝜕𝜕 +30𝑎𝑎 𝜕𝜕 1 11 3 111 3 ∗ 3 𝜀𝜀 �24𝑎𝑎 𝜕𝜕 +120𝑎𝑎 𝜕𝜕 � 0 11 3 111 3 𝑟𝑟 = , (8.b) 33 ∗ ∗ 2 4 2𝑎𝑎 +12𝑎𝑎 𝜕𝜕 +30𝑎𝑎 𝜕𝜕 1 11 3 3 ∗ 3 ∗ 2 8𝑎𝑎 𝜕𝜕 4𝑎𝑎 𝜕𝜕 +8𝑎𝑎 𝜕𝜕 4𝑎𝑎 +24𝑎𝑎 𝜕𝜕 123 3 13 3 112 3 13 112 3 𝑟𝑟 = 𝜀𝜀 ( + + ). (8.c) 42 0 ∗ 2 ∗ ∗ 2 4 ∗ 3 4𝑎𝑎 +4𝑎𝑎 𝜕𝜕 2𝑎𝑎 +2𝑎𝑎 𝜕𝜕 +2𝑎𝑎 𝜕𝜕 4𝑎𝑎 𝜕𝜕 +8𝑎𝑎 𝜕𝜕 12 123 3 1 13 3 112 3 13 3 112 For the derivation details, please see the appendix. The orthorhombic phase of BaTiO is not a thermodynamically stable phase at room temperature. However, it could be stabilized under tensile strain, such as epitaxially grown orthorhombic BaTiO films on MgO. [7] The EO tensors of orthorhombic BaTiO are 3 3 ∗ 3 2 ∗ 2 3 𝜀𝜀 (24𝑎𝑎 𝜕𝜕 +120𝑎𝑎 𝜕𝜕 +24𝑎𝑎 𝜕𝜕 𝜕𝜕 +4𝑎𝑎 𝜕𝜕 +24𝑎𝑎 𝜕𝜕 𝜕𝜕 +8𝑎𝑎 𝜕𝜕 ) 0 11 1 111 1 112 1 2 12 2 112 2 1 112 2 𝑟𝑟 = , (9.a) ∗ ∗ 2 2 4 4 2 2 (2𝑎𝑎 +2𝑎𝑎 �𝜕𝜕 +𝜕𝜕 �+2𝑎𝑎 (𝜕𝜕 +𝜕𝜕 )+2𝑎𝑎 �𝜕𝜕 𝜕𝜕 � ) 3 13 1 2 112 1 2 123 1 2 ∗ 3 3 𝜀𝜀 (4𝑎𝑎 (𝜕𝜕 +𝜕𝜕 )+8𝑎𝑎 (𝜕𝜕 +𝜕𝜕 )) 0 13 1 2 112 1 2 𝑟𝑟 = , (9.b) 33 ∗ ∗ 2 2 4 4 2 2 �2𝑎𝑎 +2𝑎𝑎 �𝜕𝜕 +𝜕𝜕 �+2𝑎𝑎 �𝜕𝜕 +𝜕𝜕 �+2𝑎𝑎 𝜕𝜕 𝜕𝜕 � 3 13 1 2 112 1 2 123 1 2 ∗ 3 ∗ 3 4𝑎𝑎 𝜕𝜕 +8𝑎𝑎 𝜕𝜕 4𝑎𝑎 𝜕𝜕 +8𝑎𝑎 𝜕𝜕 13 1 112 1 13 2 112 2 𝑟𝑟 = 𝜀𝜀 ( + ). (9.c) 42 0 ∗ ∗ 2 ∗ 2 4 2 2 4 ∗ 3 3 2𝑎𝑎 +12𝑎𝑎 𝜕𝜕 +2𝑎𝑎 𝜕𝜕 +30𝑎𝑎 𝜕𝜕 +12𝑎𝑎 𝜕𝜕 𝜕𝜕 +2𝑎𝑎 𝜕𝜕 4𝑎𝑎 𝜕𝜕 𝜕𝜕 +8𝑎𝑎 𝜕𝜕 𝜕𝜕 +8𝑎𝑎 𝜕𝜕 𝜕𝜕 111 112 112 1 2 112 2 112 1 1 11 1 12 2 1 1 2 2 12 1 2 It is worth noting that the 𝑎𝑎 coefficient is temperature-dependent, as shown in eqn. (2). Hence, the temperature-dependent EO responses could also be obtained using this method. For LiNbO3 and LiTaO3, the EO coefficient 𝑟𝑟 is 𝜀𝜀 24𝑎𝑎 𝜕𝜕 0 11 3 𝑟𝑟 = . (10) 33 2 2𝑎𝑎 +12𝑎𝑎 𝜕𝜕 1 11 Page 10 of 31 Microscopically, there are three contributions induced by a modulating electric field to the EO tensor: the electronic contribution from the valence electrons, the ionic contribution from the displacement of the ions, and the piezoelectric contribution from the distortion of the unit cell through the converse piezoelectric effect. [19-25] At moderately high modulating frequencies, that are low compared to the optical phonon modes, ionic contributions to the EO tensor dominate. At such frequencies, strain relaxations can be avoided. To investigate the frequency dispersion of the coefficients, we also applied the time-dependent Ginzburg-Landau (TDGL) equation [6,33]: 𝜕𝜕 𝜕𝜕 (𝑡𝑡 ) 𝜕𝜕𝜕𝜕 (𝜕𝜕 ) 𝑖𝑖 𝑖𝑖 = −𝐿𝐿 (11) 𝜕𝜕𝑡𝑡 𝜕𝜕 𝜕𝜕 𝑖𝑖 where L is the kinetic coefficient (proportional to the dipole motion velocity) and t is time. The energy function F is Δ𝑓𝑓 in equation (5) except the fact that the applied electric field is static but here, dynamic electric field is used as a triangle wave function: −1 𝐸𝐸 (𝑡𝑡 ) = 𝐸𝐸 sin [sin(𝑡𝑡𝑓𝑓𝑓𝑓 )] (12) where E is the amplitude of the electric field and f the frequency. It is not easy to solve TDGL explicitly, as we have done for the static calculations. Therefore, we performed the calculations using finite element method to obtain the EO coefficients from 10 Hz to 1 GHz. In this work, we provide an example of the frequency dependent r for the tetragonal BaTiO . We have ignored optical frequencies where the electronic contributions dominate, as they are not included in our model. It is acceptable to ignore the electronic contributions at low frequencies, as they have been shown to be relatively small compared to the dipole contributions in perovskites [20,28]. In summary, we have extracted the Landau-Devonshire coefficients from the free-energy landscape calculated using first-principles DFT using equations (2.a-c). We have simulated the dependence of polarization on the applied electric field using equations (6.a-c) and calculated the dielectric constant using the relations shown in Box I and II. The electric field and dielectric constant as a function of polarization are then plugged into equation (7) to obtain the electro-optic tensors. In this case, we give two solutions for tetragonal and orthorhombic BaTiO in equations (8.a-c) and (9.a-c), respectively. Page 11 of 31 Results and Discussion We have used the LDA functional to calculate the Helmholtz free energy density as a function of the polarization for LiNbO , LiTaO , and BaTiO . In the case of BaTiO , the high-temperature 3 3 3 3 phase has a centrosymmetric cubic structure. However, as the temperature decreases, a sequence of phase transitions are observed experimentally as follows: 388𝐾𝐾 273𝐾𝐾 183𝐾𝐾 cubic�⎯⎯� tetragonal�⎯⎯� orthorhombic�⎯⎯� rhombohedral. [54] These three ferroelectric phase transitions result in a change in the direction of the spontaneous polarization from the [001] axis (tetragonal, 𝑃𝑃 = 𝑃𝑃 = 0,𝑃𝑃 ≠ 0), to the [110] axes (orthorhombic, 𝑃𝑃 = 𝑃𝑃 ≠ 0,𝑃𝑃 = 0), and to 1 2 3 1 2 3 the [111] axes (rhombohedral, 𝑃𝑃 = 𝑃𝑃 = 𝑃𝑃 ≠ 0) as the temperature decreases. In Figure 1, the 1 2 3 energy density changes for the phase transition from the cubic structure directly to tetragonal, orthorhombic, and rhombohedral phase of BaTiO are shown. The spontaneous polarization is corresponding to the polarization value, where the energy reaches the minima at the bottom of the double-well. The calculated spontaneous polarization values for tetragonal, orthorhombic, and 2 2 2 rhombohedral BaTiO are 0.24 C/m , 0.27 C/m , and 0.32 C/m , respectively. (a) (b) (c) Figure 1. Energy density as a function of polarization along (a) [001], (b) [011], and (c) [111] directions in BaTiO obtained using DFT calculations. Page 12 of 31 In the case of LiNbO and LiTaO , both undergo a transition from a high temperature 3 3 rhombohedral paraelectric 𝑅𝑅 3𝑐𝑐 phase to a low temperature ferroelectric R3c phase at 1480 K and 950 K, respectively. [5,55,56] At room-temperature, the spontaneous polarization points along the c-axis direction (𝑃𝑃 = 𝑃𝑃 = 0,𝑃𝑃 ≠ 0) for the ferroelectric rhombohedral phase of LiNbO and 1 2 3 LiTaO . We have compared the change in Helmholtz free energy density with respect to the polarization along [001] direction for BaTiO , LiNbO and LiTaO , as shown in Figure 2. The 3 3 3 double-well depth here is a quantitative indicator of the energetic stability of the ferroelectric phase with respect to the paraelectric phase. Tetragonal BaTiO yields a shallow double-well indicating a relatively easier transition from the ferroelectric to paraelectric phase. Table 2. Extracted Landau-Devonshire coefficients from the DFT calculated free energy curves. BaTiO LiNbO LiTaO Landau-Devonshire 3 3 3 Coefficient This work Ref. 67 This work Ref. 67 This work Ref. 67 2 2 7 9 8 9 9 𝑎𝑎 (Nm /C ) -4.74×10 -1.20×10 -6.28×10 -1.54×10 -1.006×10 1 -6.07×10 6 4 8 8 9 9 8 𝑎𝑎 (Nm /C ) -2.10×10 9.03×10 1.26×10 2.21×10 9.01×10 11 4.32×10 6 4 𝑎𝑎 (Nm /C ) 6.29×10 7.97×10 10 6 𝑎𝑎 (Nm /C ) 1.29×10 1.29×10 10 6 𝑎𝑎 (Nm /C ) -1.44×10 -1.95×10 10 6 𝑎𝑎 (Nm /C ) -1.67×10 -2.50×10 Page 13 of 31 Figure 2. Dependence of free energy density on the polarization along the [001] crystallographic direction in LiNbO , LiTaO , and BaTiO . The solid lines 3 3 indicated polynomial fitting to the polarization dependent energy curves for LiNbO , LiTaO , and BaTiO along the 3 3 3 [001] axis. The Landau-Devonshire coefficients are extracted by fitting the double-well energy curves obtained from the first-principles calculations to a polynomial expression (eqn. (1)). Under rigid th symmetry framework, the 6 order series expansions are commonly accepted as the basic free energy format describing the ferroelectric phase transitions in BaTiO . [60] A polynomial with higher-order expansion terms will yield better accuracy but the large number of fitting parameters can also lead to overfitting. Moreover, the higher-order terms require additional information at high polarization region and their physical meaning still remains unclear. [65] Hence, we th performed the polynomial fitting to the free energy density curves of BaTiO up to 6 order expansions to ensure that the critical aspects of the electro-optic phenomena can be sufficiently described without overfitting. We fit the polynomial of eqns. (3.a), (3.b), and (3.c) to free energy density curves of the three ferroelectric phases: tetragonal, orthorhombic, and rhombohedral, respectively, by simultaneously and manually adjusting the fitting parameters to find the smallest 2 2 𝑅𝑅 (coefficient of determination). The 𝑅𝑅 values are 0.998, 0.999, and 0.996 for tetragonal, orthorhombic, and rhombohedral phase of BaTiO , respectively. Order parameter just below the Page 14 of 31 Curie temperature has rich information for all the polynomial coefficients, including double P-E loops, [66] while the free energy curve for a system much below the phase transition temperature is a simple double well without much information about higher order polynomials to the expansion coefficients, which is the case for LiNbO and LiTaO [5,56] The fourth-order polynomial can 3 3. already provide a good fit to the first-principles data, where 𝑅𝑅 is already 0.999. Hence, the higher- order terms were omitted for LiNbO and LiTaO . The Landau-Devonshire coefficients extracted 3 3 for the noted materials are presented in Table 2. The coefficients from Long-Qing Chen [67] are listed in the table as well for comparison. We now apply the previously fitted parameters to calculate the electro-optic coefficients using equations (8.a-c). Figure 3 shows the comparison between the calculated values using the model and experimental results. [61] Overall, our method predicts the sign of the EO constants correctly and the values of the EO coefficients are in good agreement with the experimental values. The calculated 𝑟𝑟 for LiNbO3 is in close agreement (deviation is ~ 1.8%) with the experimentally reported results. [61] The calculated values of 𝑟𝑟 are ~30% larger than the experimental value of tetragonal BaTiO and LiTaO . Our theoretical 𝑟𝑟 value of BaTiO also matches well with the 3 3 3 experimental ones. [61] However, the calculated 𝑟𝑟 is 88% lower than the experimentally reported value. [61] The deviation from experimental values could be due to extrinsic factors such as stoichiometry and structural quality of the samples and domain structures. Furthermore, DFT calculations using the LDA functional underestimate the polar distortion from paraelectric phases, which results in relatively shallow double-well depths. [44] Therefore, the accuracy of Landau- Devonshire fitting parameters could be further improved by applying more reliable exchange- correlation functionals, such as the recently developed strongly constrained and appropriately normed (SCAN) meta-GGA functional, which has been shown to systematically improve over LDA for structural properties and ferroelectric transitions of diversely bonded materials. [68] These advances will be part of future work. The temperature dependence of the EO coefficients is another way to examine the accuracy of the model. We compare the calculated results with the experimental values for the LiNbO in Figure 4 (a). [69] The model captures the experimentally observed increase in 𝑟𝑟 ; however, it consistently over-estimates the value with a larger deviation on decreasing the temperature. As Page 15 of 31 (b) (a) Figure 3. Theoretical and experimental (Ref. [61]) electro-optic coefficients of (a) tetragonal BaTiO , (b) LiNbO and LiTaO . 3 3 3 temperature approaches 0 K, the deviation is < 20%, and, as the temperature reaches room temperature, the deviation narrows down to < 10%. To the best of our knowledge, the experimental results of temperature dependency of r for LiTaO and BaTiO have not been 33 3 3 reported in the literature, and hence, we could not make this comparison for these two compounds. Nevertheless, we compare our results (noted as green hollow circles with a trendline) with the two available first-principles results from Veithen et al. [27] and Pietro et al. [70] for tetragonal BaTiO , as depicted with red and yellow markers in Figure 3 (b), respectively. The electro-optic coefficient 𝑟𝑟 increases with temperature below the Curie temperature and shows a divergent Figure 4. Temperature dependence of electro-optic coefficients in (a) LiNbO and (b) BaTiO 3 3. Page 16 of 31 trend close to the Curie temperature, as the spontaneous polarization abruptly drops to zero at and above the Curie temperature. It is worth reemphasizing that in our model, we assume 𝑎𝑎 as the only temperature-dependent term and hence, the temperature dependence of EO response is essentially attributed to it. The contribution of other terms to EO coefficients in equation (8.a) and (9.a), such as 𝑎𝑎 𝑃𝑃 , are typically much smaller than 𝑎𝑎 . Hence, we resorted to this simplification. Any corrections to this simplification will be explored in the future depending on the availability of the experimental results. EO coefficients tend to be extremely large as the temperature approaches the Curie temperature, T . By definition, 𝑎𝑎 , given in eqn. (2), converges to 0 as the temperature reaches 𝑇𝑇 . 1 0 This hints that small 𝑎𝑎 is desirable to achieve a large EO coefficient. To obtain a small 𝑎𝑎 , the 1 1 energy barrier for switching the polarization from one energy well to the other has to be low. It indicates that the origin of this electro-optic enhancement is attributed to the ease of the ferroelectric switching as manifested in the free energy landscape. From Figure 2 and Table 1, BaTiO has the shallowest energy well and the smallest absolute value of 𝑎𝑎 , which leads to the largest 𝑟𝑟 among the three ferroelectric oxides. This also explains why the EO coefficients of relaxor ferroelectric oxide alloys are high, for which the corresponding ferroelectric switching energy barrier is relatively low. [31] Figure 5. (a) Dependence of the electro-optic coefficient, r , in tetragonal BaTiO on the 33 3 frequency of the applied electric field. (b) The frequency-dependent normalized dielectric constant for tetragonal BaTiO in comparison with experimental data. [80-82] 3, Page 17 of 31 The frequency dispersion of r and the normalized dielectric constant are shown in Figure 5. The EO coefficient 𝑟𝑟 turns out to be 100 pm/V at 1 kHz, which is the most common frequency for dielectric measurements. Compared to its experimental value of 105 pm/V, the deviation has dropped from 30 % for static-field calculations (𝑟𝑟 : 137 pm/V) to -5 % for dynamic calculations. The EO coefficient decreases with increasing frequency of the electric field. As the frequency approaches 10 MHz, the EO coefficient is 63.3 pm/V, as shown in Figure 5(a). We believe the magnitude of EO coefficients are limited by the polarization switching velocity and hence, show a decreasing trend with increasing frequency. Typically, the polarization switching includes three stages: nucleation of domains, domain growth, and sideways motion of newly formed domains. [71–74] Besides nucleation and domain motion, polarization switching can also occur with the simultaneous reversal of all local dipoles, as proposed by Landau and Lifshitz. [75] This type of polarization switching is called intrinsic polarization switching or “Landau switching”. In this work, we will follow the Landau approach. The switching time for Landau switching used in this work is sub nanoseconds, which is comparable to the fastest domain wall driven switching record (220 ps). [76,77] In comparison, the measured polarization switching time usually varies from μs to ns scale. [78,79] The intrinsic switching is believed to be much faster than domain nucleation and motion. To evaluate the switching time assumption, we also calculate the frequency-dependent dielectric constant across the 10 kHz to 1 GHz range and qualitatively compare it with the experimental results, as shown in Figure 5 (b). The dielectric dispersion from our method falls in between the data from a bulk single domain BaTiO single crystal [80] and single crystal thin films data, [81,82] which proves that the selected switching time in this work is reasonable. In this work, we demonstrated the dipole contribution to the EO coefficient, specifically 𝑟𝑟 for tetragonal BaTiO by varying the frequency of the electric field. Admittedly, we have not 3, considered the electronic contribution in our model, which dominates at even higher frequencies ~ 100 GHz level. This contribution is relatively small compared to the dipole contribution in perovskites, and is acceptable to ignore it. [20,28] In the following, we demonstrate the role of the piezoelectric term to the EO coefficient by applying strain and adding the electrostrictive terms to the free energy. Page 18 of 31 Figure 6 shows the EO tensor as a function of strain with misfit strains ranging from -5% to 5% along the in-plane a and b axes for BaTiO , which was obtained using eqns. 8(a-c) and 9(a-c). The thermodynamically stable phases are obtained by minimizing the total free energy F under a given misfit strain. The calculation was performed by setting the temperature in the free energy expansion coefficient, 𝑎𝑎 , as room temperature. We obtained two stable single-phase states. The stable strain condition for tetragonal phase is denoted as T with yellow shade (𝑃𝑃 = 𝑃𝑃 = 0,𝑃𝑃 ≠ 1 2 3 0) below -1% compressive strain, and for orthorhombic phase, it is denoted as O with blue shade (𝑃𝑃 = 𝑃𝑃 ≠ 0,𝑃𝑃 = 0) above ~1.9% tensile strain. For the phase region between the T and O phase, 1 2 3 it contains a two-phase mixture (T+O). At zero strain state, we obtain almost identical values of electro-optical tensors as bulk values summarized in Figure 4. The small deviation comes from that the electrostrictive energy term is not included in the previous calculations. The strain-induced polarization variations under compression and tensile strain generate a large contribution to the relevant EO constants. We predict a surprisingly high value, up to thousands of pm/V, for 𝑟𝑟 coefficient of BaTiO under ~ −1.3% compressive strain and 𝑟𝑟 ≈ 800 pm/V under 1.8% tensile strain, which is 1 – 2 orders of magnitude larger than 𝑟𝑟 = 32 pm/V for LiNbO . [61] In this specific case, at the transition region between T and T-O, the sign change of 𝑎𝑎 leads to a large value of the 𝑟𝑟 coefficient; similarly, 𝑟𝑟 is large at the T-O and O phase boundary due to the 42 33 vanishingly small 𝑎𝑎 values. Similar EO coefficient enhancement has been demonstrated as a function of temperature near a ferroelectric phase transition due to the divergence of the dielectric constant. [24,70] The polarization rotation at the phase boundary is easier to achieve than in a single phase region, as the corresponding switching barrier is small. This could also explain why ferroelectrics at the phase boundary usually have large dielectric, piezoelectric, and electro-optic responses. [83–85] Page 19 of 31 SrTiO GdScO (a) (b) (c) Figure 6. Evolution of EO coefficients (a) 𝑟𝑟 , (b) 𝑟𝑟 , (c) 𝑟𝑟 with misfit strain for BaTiO . 42 13 The EO coefficients calculated from LD and DFPT are indicated with solid lines and dotted lines with marker. The blue region and yellow region, respectively, represent orthorhombic (O) and tetragonal (T) phase. The misfit strain imposed on coherently strained BaTiO films grown on SrTiO and GdScO substrates are 3 3 marked to identify the potentially achievable EO coefficients in strained thin films. We have compared the change in EO response of BaTiO under strain using our model, with those obtained using computationally intensive first-principles calculations as discussed in the Methods section. We obtain qualitatively similar results, as depicted in Figure 5. The first-principles DFPT results also indicate the divergence of r and r at ~1.3% and −1.0% strain, respectively. This 33 42 clearly demonstrates that our model could qualitatively describe the physics of EO response under epitaxial strain conditions, although the absolute values should be taken with a grain of salt. We give the potential substrates in Figure 5 that one could use to achieve misfit compressive strain close to 1% to obtain large EO coefficients in coherently strained thin films of BaTiO using the established strain engineering approach. The realization of low loss thin films in these limits remains a challenge, where intrinsic and extrinsic defects may play an important role. Other factors, Page 20 of 31 such as the formation and distribution of ferroelectric domains at low fields, any inhomogeneity in application of electric fields will also influence dielectric losses. Overcoming these issues are critical to enable high EO coefficient materials and devices for future photonic technologies. Conclusion A methodology to predict the EO coefficients in ferroelectric oxides as a function of strain, modulation frequency, and temperature is demonstrated here. This method enables highly scalable calculation of EO coefficients by combining computationally expensive, but accurate, first- principles calculations with scalable phenomenological Landau-Devonshire theory. We applied our approach to three representative ferroelectric oxides, namely, LiNbO , LiTaO , and BaTiO . 3 3 3 The calculated EO coefficients are in good agreement with the experimental results. And the relevance of specific model parameters for EO effect are discussed. In the light of the previous discussion of temperature and strain effects on the EO responses, we conclude that small 𝑎𝑎 (𝑎𝑎 for 1 1 strained case) and shallow energy barrier are favorable for high-𝑟𝑟 EO materials. One way to reduce 𝑎𝑎 is fabrication of ferroelectric/dielectric (FE/DE) heterostructures. As we introduced earlier, Helmholtz free energy is simply the sum of all the energy components in Landau-Devonshire theory. The total 𝑎𝑎 is the sum of 𝑎𝑎 of the ferroelectric and dielectric phases weighted according 1 1 to their respective volume fractions. 𝑎𝑎 is negative for ferroelectrics and positive for dielectrics. By engineering the volume fraction of each component, one can achieve small 𝑎𝑎 coefficients with high EO response. Another aspect of EO performance enhancement in ferroelectrics or multilayer heterostructures is compositional engineering of the phase boundary involving the ferroelectric- paraelectric phases or two different ferroelectric phases such as the O- and T-phases in a mixed- phase BaTiO film. Low switching energy at phase boundary makes the polarization switching (between up and down) and/or polarization rotation (between two ferroelectric states) easier and hence, can generate the enhanced EO responses. [86] Our method can be extended to simulate the EO response of FE/DE heterostructures or film having phase boundary, once the complete set of the bulk, elastic, electrostrictive and gradient energy terms are available either from first-principles calculations or from experiments. We expect that this method will pave a way to discover new materials with high electro-optic performance. Page 21 of 31 Appendix The EO tensors, 𝑟𝑟 ,𝑟𝑟 , and 𝑟𝑟 , are defined in equation 7. Here we provide detail derivations to 13 33 42 them. From equation 6.b-c, we have E and E as 2 3 𝜕𝜕 𝑓𝑓 (12.a) ∗ ∗ 3 ∗ 2 ∗ 2 5 𝐸𝐸 = = 2𝑎𝑎 𝑃𝑃 + 4𝑎𝑎 𝑃𝑃 + 2𝑎𝑎 𝑃𝑃 𝑃𝑃 + 2𝑎𝑎 𝑃𝑃 𝑃𝑃 + 6𝑎𝑎 𝑃𝑃 2 1 2 11 2 12 1 2 13 2 3 111 2 𝜕𝜕 𝑃𝑃 4 4 3 2 2 2 2 +𝑎𝑎 � 2𝑃𝑃 (𝑃𝑃 +𝑃𝑃 ) + 4𝑃𝑃 (𝑃𝑃 +𝑃𝑃 )�+ 2𝑎𝑎 𝑃𝑃 𝑃𝑃 𝑃𝑃 112 2 3 1 2 1 3 123 1 2 3 𝜕𝜕 𝑓𝑓 (12.b) ∗ ∗ 3 ∗ 2 2 5 𝐸𝐸 = = 2𝑎𝑎 𝑃𝑃 + 4𝑎𝑎 𝑃𝑃 + 2𝑎𝑎 𝑃𝑃 (𝑃𝑃 +𝑃𝑃 ) + 6𝑎𝑎 𝑃𝑃 3 3 3 33 3 13 3 1 2 111 3 𝜕𝜕 𝑃𝑃 4 4 3 2 2 2 2 +𝑎𝑎 � 2𝑃𝑃 (𝑃𝑃 +𝑃𝑃 ) + 4𝑃𝑃 (𝑃𝑃 +𝑃𝑃 )�+ 2𝑎𝑎 𝑃𝑃 𝑃𝑃 𝑃𝑃 112 3 123 3 1 2 3 1 2 1 2 For the left hand side of the equation (7), we have: (13.a) 1 𝜕𝜕 𝑓𝑓 ∗ ∗ 2 ∗ 2 ∗ 2 4 = 𝜀𝜀 = 𝜀𝜀 (2𝑎𝑎 + 12𝑎𝑎 𝑃𝑃 + 2𝑎𝑎 𝑃𝑃 + 2𝑎𝑎 𝑃𝑃 + 30𝑎𝑎 𝑃𝑃 0 0 1 11 1 12 2 13 3 111 1 2 2 𝑛𝑛 𝜕𝜕 𝑃𝑃 1 1 2 2 2 4 4 2 2 ( ) +𝑎𝑎 12𝑃𝑃 (𝑃𝑃 +𝑃𝑃 + 2(𝑃𝑃 +𝑃𝑃 )) + 2𝑎𝑎 𝑃𝑃 𝑃𝑃 ) 112 1 2 3 2 3 123 2 3 1 𝜕𝜕 𝑓𝑓 (13.b) ∗ ∗ 2 ∗ 2 2 4 2 2 2 ( ) = 𝜀𝜀 = 𝜀𝜀 (2𝑎𝑎 + 12𝑎𝑎 𝑃𝑃 + 2𝑎𝑎 (𝑃𝑃 +𝑃𝑃 ) + 30𝑎𝑎 𝑃𝑃 +𝑎𝑎 12𝑃𝑃 (𝑃𝑃 +𝑃𝑃 0 0 3 33 3 13 1 2 111 3 112 3 1 2 2 2 𝑛𝑛 𝜕𝜕 𝑃𝑃 3 3 4 4 2 2 + 2(𝑃𝑃 +𝑃𝑃 )) + 2𝑎𝑎 𝑃𝑃 𝑃𝑃 ) 1 2 123 1 2 1 𝜕𝜕 𝑓𝑓 (13.c) ∗ 3 3 2 ( ) = 𝜀𝜀 = 𝜀𝜀 (4𝑎𝑎 𝑃𝑃 𝑃𝑃 +𝑎𝑎 8𝑃𝑃 𝑃𝑃 + 8𝑃𝑃 𝑃𝑃 + 4𝑎𝑎 𝑃𝑃 𝑃𝑃 𝑃𝑃 ) 0 0 13 2 3 112 2 3 2 3 123 1 2 3 𝑛𝑛 𝜕𝜕 𝑃𝑃 𝜕𝜕 𝑃𝑃 2 3 Then substituting the eqns. 12 and eqns. 13 to eqn.7 we obtain: (14.a) 𝜕𝜕 � � ∗ 3 𝜀𝜀 (4𝑎𝑎 𝑃𝑃 + 8𝑎𝑎 𝑃𝑃 ) 𝑛𝑛 0 12 3 112 3 𝑟𝑟 = = ∗ ∗ 2 4 𝜕𝜕 𝐸𝐸 2𝑎𝑎 + 12𝑎𝑎 𝑃𝑃 + 30𝑎𝑎 𝑃𝑃 3 111 1 11 3 3 (14.b) 𝜕𝜕 � � ∗ 3 ( ) 𝑛𝑛 𝜀𝜀 24𝑎𝑎 𝑃𝑃 + 120𝑎𝑎 𝑃𝑃 0 11 3 111 3 𝑟𝑟 = = ∗ ∗ 2 4 𝜕𝜕 𝐸𝐸 2𝑎𝑎 + 12𝑎𝑎 𝑃𝑃 + 30𝑎𝑎 𝑃𝑃 3 111 1 11 3 3 (14.c) 𝜕𝜕 � � 2 ∗ 8𝑎𝑎 𝑃𝑃 4𝑎𝑎 𝑃𝑃 + 8𝑎𝑎 𝑃𝑃 𝑛𝑛 123 3 3 112 13 3 𝑟𝑟 = = 𝜀𝜀 ( + 42 0 ∗ 2 ∗ ∗ 2 4 𝜕𝜕 𝐸𝐸 4𝑎𝑎 + 4𝑎𝑎 𝑃𝑃 2𝑎𝑎 + 2𝑎𝑎 𝑃𝑃 + 2𝑎𝑎 𝑃𝑃 2 123 112 12 3 1 13 3 3 ∗ 2 4𝑎𝑎 + 24𝑎𝑎 𝑃𝑃 13 112 3 + ) ∗ 3 4𝑎𝑎 𝑃𝑃 + 8𝑎𝑎 𝑃𝑃 3 112 13 3 We use the same procedure to derive the equations 9.a-c except the fact that the P and P are 1 2 non-zero term in the orthorhombic phase. Page 22 of 31 Acknowledgements This work was partially supported by the Army Research Office under award no. W911NF-21-1- 0327. The work at USC was also supported in part by ARO under award No. W911NF-19-1-0137, the National Science Foundation (NSF) through award DMR-2122071, and the Air Force Office of Scientific Research under contract FA9550-16-1-0335. Y. L. and J. R. acknowledge the Center for Advanced Research Computing (CARC) at the University of Southern California for providing computing resources that have contributed to the research results reported within this publication. URL: https://carc.usc.edu. The work at Washington University was also supported in part by NSF through awards DMR-1806147, DMR-2122070 and DMR-1931610. This work used the computational resources of the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by NSF grant ACI-1548562. The authors gratefully acknowledge Mr. Wente Li and Prof. Alex Demkov of UT Austin for feedback on the manuscript. Data Availability The data that support the findings of this study are available from the corresponding authors upon reasonable request. Page 23 of 31 References [1] P. K. Panda, Review: Environmental Friendly Lead-Free Piezoelectric Materials, J. Mater. Sci. 44, 5049 (2009). [2] Z. Sun, Z. Wang, Y. Tian, G. Wang, W. Wang, M. Yang, X. Wang, F. Zhang, and Y. Pu, Progress, Outlook, and Challenges in Lead-Free Energy-Storage Ferroelectrics, Adv. Electron. Mater. 6, 1 (2020). [3] J. Hao, W. Li, J. Zhai, and H. Chen, Progress in High-Strain Perovskite Piezoelectric Ceramics, Mater. Sci. Eng. R Reports 135, 1 (2019). [4] A. S. Bhalla, R. Guo, and R. Roy, The Perovskite Structure - A Review of Its Role in Ceramic Science and Technology, Mater. Res. Innov. 4, 3 (2000). [5] R. Weis and T. Gaylord, Lithium Niobate: Summary of Physical Properties and Crystal Structure R., Appl. Phys. A Mater. Sci. Process. 37, 191 (1985). [6] D. G. Schlom, L. Q. Chen, X. Pan, A. Schmehl, and M. A. Zurbuchen, A Thin Film Approach to Engineering Functionality into Oxides, J. Am. Ceram. Soc. 91, 2429 (2008). [7] D. H. Kim and H. S. Kwok, Pulsed Laser Deposition of BaTiO3 Thin Films and Their Optical Properties, Appl. Phys. Lett. 67, 1803 (1995). [8] D. M. Gill, C. W. Conrad, G. Ford, B. W. Wessels, and S. T. Ho, Thin-Film Channel Waveguide Electro-Optic Modulator in Epitaxial BaTiO3, Appl. Phys. Lett. 71, 1783 (1997). [9] G. Yi, Z. Wu, and M. Sayer, Preparation of Pb(Zr,Ti)O3 Thin Films by Sol Gel Processing: Electrical, Optical, and Electro-Optic Properties, J. Appl. Phys. 64, 2717 (1988). [10] D. H. Reitze, E. Haton, R. Ramesh, S. Etemad, D. E. Leaird, T. Sands, Z. Karim, and A. R. Tanguay, Electro-Optic Properties of Single Crystalline Ferroelectric Thin Films, Appl. Phys. Lett. 63, 596 (1993). [11] R. A. McKee, F. J. Walker, and M. F. Chisholm, Crystalline Oxides on Silicon: The First Five Monolayers, Phys. Rev. Lett. 81, 3014 (1998). [12] S. R. Bakaul, C. R. Serrao, O. Lee, Z. Lu, A. Yadav, C. Carraro, R. Maboudian, R. Ramesh, and S. Salahuddin, High Speed Epitaxial Perovskite Memory on Flexible Substrates, Adv. Mater. 29, 1605699 (2017). Page 24 of 31 [13] S. H. Baek and C. B. Eom, Epitaxial Integration of Perovskite-Based Multifunctional Oxides on Silicon, Acta Mater. 61, 2734 (2013). [14] S. Abel, F. Eltes, J. E. Ortmann, A. Messner, P. Castera, T. Wagner, D. Urbonas, A. Rosa, A. M. Gutierrez, D. Tulli, P. Ma, B. Baeuerle, A. Josten, W. Heni, D. Caimi, L. Czornomaz, A. A. Demkov, J. Leuthold, P. Sanchis, and J. Fompeyrine, Large Pockels Effect in Micro- and Nanostructured Barium Titanate Integrated on Silicon, Nat. Mater. 18, 42 (2019). [15] K. J. Kormondy, Y. Popoff, M. Sousa, F. Eltes, D. Caimi, M. D. Rossell, M. Fiebig, P. Hoffmann, C. Marchiori, M. Reinke, M. Trassin, A. A. Demkov, J. Fompeyrine, and S. Abel, Microstructure and Ferroelectricity of BaTiO Thin Films on Si for Integrated Photonics, Nanotechnology 28, 075706 (2017). [16] C. Wang, M. Zhang, X. Chen, M. Bertrand, A. Shams-Ansari, S. Chandrasekhar, P. Winzer, and M. Lončar, Integrated Lithium Niobate Electro-Optic Modulators Operating at CMOS-Compatible Voltages, Nature 562, 101 (2018). [17] S. Abel, T. Stöferle, C. Marchiori, C. Rossel, M. D. Rossell, R. Erni, D. Caimi, M. Sousa, A. Chelnokov, B. J. Offrein, J. Fompeyrine, S. Abel, T. Sto, D. Caimi, M. Sousa, A. Chelnokov, B. J. Offrein, and J. Fompeyrine, A Strong Electro-Optically Active Lead-Free Ferroelectric Integrated on Silicon, Nat. Commun. 4, 1671 (2013). [18] H. Akazawa and M. Shimada, Electro-Optic Properties of c-Axis Oriented LiNbO3 Films Grown on Si(100) Substrate, Mater. Sci. Eng. B 120, 50 (2005). [19] T. Paoletta and A. A. Demkov, Pockels Effect in Low-Temperature Rhombohedral BaTiO3, Phys. Rev. B 103, 14303 (2021). [20] A. K. Hamze, M. Reynaud, J. Geler-Kremer, and A. A. Demkov, Design Rules for Strong Electro-Optic Materials, Npj Comput. Mater. 6, (2020). [21] K. D. Fredrickson, V. V. Vogler-Neuling, K. J. Kormondy, D. Caimi, F. Eltes, M. Sousa, J. Fompeyrine, S. Abel, and A. A. Demkov, Strain Enhancement of the Electro-Optical Response in BaTi O3 Films Integrated on Si(001), Phys. Rev. B 98, 1 (2018). [22] A. K. Hamze and A. A. Demkov, First-Principles Study of the Linear Electro-Optical Response in Strained SrTi O3, Phys. Rev. Mater. 2, 1 (2018). [23] C. Paillard, S. Prokhorenko, and L. Bellaiche, Strain Engineering of Electro-Optic Constants in Ferroelectric Materials, Npj Comput. Mater. 5, 1 (2019). Page 25 of 31 [24] M. DiDomenico and S. H. Wemple, Oxygen-Octahedra Ferroelectrics. I. Theory of Electro-Optical and Nonlinear Optical Effects, J. Appl. Phys. 40, 720 (1969). [25] S. H. Wemple and M. DiDomenico, Oxygen-Octahedra Ferroelectrics. II. Electro-Optical and Nonlinear-Optical Device Applications, J. Appl. Phys. 40, 735 (1969). [26] M. Veithen and P. Ghosez, First-Principles Study of the Dielectric and Dynamical Properties of Lithium Niobate, Phys. Rev. B 65, 214302 (2002). [27] M. Veithen and P. Ghosez, Temperature Dependence of the Electro-Optic Tensor and Refractive Indices of BaTiO3 from First Principles, Phys. Rev. B 71, 132101 (2005). [28] M. Veithen, X. Gonze, and P. Ghosez, First-Principles Study of the Electro-Optic Effect in Ferroelectric Oxides, Phys. Rev. Lett. 93, 187401 (2004). [29] M. Veithen, X. Gonze, and P. Ghosez, Nonlinear Optical Susceptibilities, Raman Efficiencies, and Electro-Optic Tensors from First-Principles Density Functional Perturbation Theory, Phys. Rev. B 71, 125107Í‹ (2005). [30] J. H. Qiu, J. N. Ding, N. Y. Yuan, X. Q. Wang, and Y. Zhou, Film Thickness Dependence of Electro-Optic Effect in Epitaxial BaTiO 3 Thin Films, Solid State Commun. 151, 1344 (2011). [31] Y. Lu and R. J. Knize, Enhanced Dielectric and Electro-Optic Effects in Relaxor Oxide Heterostructured Superlattices, J. Phys. D. Appl. Phys. 37, 2432 (2004). [32] J. Hiltunen, D. Seneviratne, R. Sun, M. Stolfi, H. L. Tuller, J. Lappalainen, and V. Lantto, BaTiO3–SrTiO3 Multilayer Thin Film Electro-Optic Waveguide Modulator, Appl. Phys. Lett. 89, 242904 (2006). [33] P. Chandra and P. B. Littlewood, A Landau Primer for Ferroelectrics, in Physics of Ferroelectrics (Springer Berlin Heidelberg, Berlin, Heidelberg, 2007), pp. 69–116. [34] K. G. Lim, K. H. Chew, L. H. Ong, and M. Iwata, Recent Advances in Application of Landau-Ginzburg Theory for Ferroelectric Superlattices, Solid State Phenom. 232, 169 (2015). [35] K. H. Chew, Recent Applications of Landau-Ginzburg Theory to Ferroelectric Superlattices: A Review, Solid State Phenom. 189, 145 (2012). [36] L.-Q. Chen, Phase-Field Method of Phase Transitions/Domain Structures in Ferroelectric Thin Films: A Review, J. Am. Ceram. Soc. 91, 1835 (2008). Page 26 of 31 [37] A. I. Khan, K. Chatterjee, J. P. Duarte, Z. Lu, A. Sachid, S. Khandelwal, R. Ramesh, C. Hu, and S. Salahuddin, Negative Capacitance in Short-Channel FinFETs Externally Connected to an Epitaxial Ferroelectric Capacitor, IEEE Electron Device Lett. 37, 111 (2016). [38] P. Marton, I. Rychetsky, and J. Hlinka, Domain Walls of Ferroelectric BaTiO3 within the Ginzburg-Landau-Devonshire Phenomenological Model, Phys. Rev. B - Condens. Matter Mater. Phys. 81, (2010). [39] G. Kresse and J. Furthmüller, Efficient Iterative Schemes for Ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set, Phys. Rev. B 54, 11169 (1996). [40] D. Joubert, From Ultrasoft Pseudopotentials to the Projector Augmented-Wave Method, Phys. Rev. B - Condens. Matter Mater. Phys. 59, 1758 (1999). [41] R. O. Jones and O. Gunnarsson, The Density Functional Formalism, Its Applications and Prospects, Rev. Mod. Phys. 61, 689 (1989). [42] J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized Gradient Approximation Made Simple, Phys. Rev. Lett. 77, 3865 (1996). [43] D. I. Bilc, R. Orlando, R. Shaltaf, G. M. Rignanese, J. Íñiguez, and P. Ghosez, Hybrid Exchange-Correlation Functional for Accurate Prediction of the Electronic and Structural Properties of Ferroelectric Oxides, Phys. Rev. B - Condens. Matter Mater. Phys. 77, 1 (2008). [44] Y. Zhang, J. Sun, J. P. Perdew, and X. Wu, Comparative First-Principles Studies of Prototypical Ferroelectric Materials by LDA, GGA, and SCAN Meta-GGA, Phys. Rev. B 96, 35143 (2017). [45] R. D. King-Smith and D. Vanderbilt, Theory of Polarization of Crystalline Solids, Phys. Rev. B 47, 1651 (1993). [46] A. K. M. I. Aroyo, J. M. Perez-Mato, D. Orobengoa, E. Tasci, G. de la Flor, No Title, “Crystallography Online Bilbao Crystallogr. Server” Bulg. Chem. Commun. 43, 183 (2011). [47] B. J. C. H. T. Stokes, D. M. Hatch, ISOTROPY Software Suite, Iso.byu.edu. [48] X. Gonze, A Brief Introduction to the ABINIT Software Package, Zeitschrift Für Krist. - Cryst. Mater. 220, 558 (2005). Page 27 of 31 [49] X. Gonze, B. Amadon, P.-M. Anglade, J.-M. Beuken, F. Bottin, P. Boulanger, F. Bruneval, D. Caliste, R. Caracas, M. Côté, T. Deutsch, L. Genovese, P. Ghosez, M. Giantomassi, S. Goedecker, D. R. Hamann, P. Hermet, F. Jollet, G. Jomard, S. Leroux, M. Mancini, S. Mazevet, M. J. T. Oliveira, G. Onida, Y. Pouillon, T. Rangel, G.-M. Rignanese, D. Sangalli, R. Shaltaf, M. Torrent, M. J. Verstraete, G. Zerah, and J. W. Zwanziger, ABINIT: First-Principles Approach to Material and Nanosystem Properties, Comput. Phys. Commun. 180, 2582 (2009). [50] N. Troullier and J. L. Martins, Efficient Pseudopotentials for Plane-Wave Calculations. II. Operators for Fast Iterative Diagonalization, Phys. Rev. B 43, 8861 (1991). [51] M. Zgonik, P. Bernasconi, M. Duelli, R. Schlesser, P. Günter, M. H. Garrett, D. Rytz, Y. Zhu, and X. Wu, Dielectric, Elastic, Piezoelectric, Electro-Optic, and Elasto-Optic Tensors of BaTiO3 Crystals, Phys. Rev. B 50, 5941 (1994). [52] F. Jona and G. Shirane, Ferroelectric Crystals (Wiley Online Library, 1962). [53] Karin M. Rabe ; Charleshn ; Jean-Marc Triscone, Physics of Ferroelectrics: Topics in Applied Physics (2007). [54] G. A. Samara, Pressure and Temperature Dependence of the Dielectric Properties and Phase Transitions of the Ferroelectric Perovskites: Pbtio3 and Batio3, Ferroelectrics 2, 277 (1971). [55] R. C. Miller and A. Savage, TEMPERATURE DEPENDENCE OF THE OPTICAL PROPERTIES OF FERROELECTRIC LiNbO 3 AND LiTaO 3, Appl. Phys. Lett. 9, 169 (1966). [56] J. L. Serving and F. Gervais, Displacive-Type Phase Transition in LiTaO 3 and LiNbO 3, Ferroelectrics 25, 609 (1980). [57] Y. L. Wang, A. K. Tagantsev, D. Damjanovic, N. Setter, V. K. Yarmarkin, and A. I. Sokolov, Anharmonicity of BaTiO3 Single Crystals, Phys. Rev. B - Condens. Matter Mater. Phys. 73, 1 (2006). [58] R. E. Cohen, Origin of Ferroelectricity in Perovskite Oxides, Nature 358, 136 (1992). [59] R. Resta and D. Vanderbilt, Theory of Polarization: A Modern Approach, in Physics of Ferroelectrics: Topics in Applied Physics (Springer Berlin Heidelberg, Berlin, Heidelberg, 2007), pp. 31–68. Page 28 of 31 [60] X. Lu, H. Li, and W. Cao, Landau Expansion Parameters for BaTiO3, J. Appl. Phys. 114, (2013). [61] M. J. Weber, Handbook of Optical Materials, Vol. 41 (CRC press, 2002). [62] N. A. Pertsev, A. G. Zembilgotov, and A. K. Tagantsev, Effect of Mechanical Boundary Conditions on Phase Diagrams of Epitaxial Ferroelectric Thin Films, Phys. Rev. Lett. 80, 1988 (1998). [63] V. B. Shirokov, Y. I. Yuzyuk, B. Dkhil, and V. V. Lemanov, Phenomenological Theory of Phase Transitions in Epitaxial BaTi O3 Thin Films, Phys. Rev. B - Condens. Matter Mater. Phys. 75, 1 (2007). [64] R. W. Boyd, Nonlinear Optics, Third Edition, 3rd ed. (Academic Press, Inc., USA, 2008). [65] J. A. Gonzalo and J. M. Rivera, Ferroelectric Free Energy Expansion Coefficients from Double Hysteresis Loops, Ferroelectrics 2, 31 (1971). [66] W. J. Merz, Double Hysteresis Loop of BaTiO3 at the Curie Point, Phys. Rev. 91, 513 (1953). [67] L.-Q. Chen, Appendix A - Landau Free-Energy Coefficient, in Physics of Ferroelectrics: Topics in Applied Physics (2007), pp. 363–370. [68] J. Sun, R. C. Remsing, Y. Zhang, Z. Sun, A. Ruzsinszky, H. Peng, Z. Yang, A. Paul, U. Waghmare, X. Wu, M. L. Klein, and J. P. Perdew, Accurate First-Principles Structures and Energies of Diversely Bonded Systems from an Efficient Density Functional, Nat. Chem. 8, 831 (2016). [69] C. Herzog, G. Poberaj, and P. Günter, Electro-Optic Behavior of Lithium Niobate at Cryogenic Temperatures, Opt. Commun. 281, 793 (2008). [70] P. Bernasconi, M. Zgonik, and P. Gunter, Temperature Dependence and Dispersion of Electro-Optic and Elasto-Optic Effect in Perovskite Crystals, J. Appl. Phys. 78, 2651 (1995). [71] W. J. Merz, Domain Formation and Domain Wall Motions in Ferroelectric BaTiO3 Single Crystals, Phys. Rev. 95, 690 (1954). [72] E. A. Little, Dynamic Behavior of Domain Walls in Barium Titanate, Phys. Rev. 98, 978 (1955). Page 29 of 31 [73] S. V Kalinin, S. Jesse, B. J. Rodriguez, Y. H. Chu, R. Ramesh, E. A. Eliseev, and A. N. Morozovska, Probing the Role of Single Defects on the Thermodynamics of Electric-Field Induced Phase Transitions, Phys. Rev. Lett. 100, 155703 (2008). [74] S. Jesse, B. J. Rodriguez, S. Choudhury, A. P. Baddorf, I. Vrejoiu, D. Hesse, M. Alexe, E. A. Eliseev, A. N. Morozovska, J. Zhang, L.-Q. Chen, and S. V Kalinin, Direct Imaging of the Spatial and Energy Distribution of Nucleation Centres in Ferroelectric Materials, Nat. Mater. 7, 209 (2008). [75] L. D. Landau and E. M. Lifshitz, Electrodynamics of Continuous Media, Second Edi (1984). [76] J. Li, B. Nagaraj, H. Liang, W. Cao, C. H. Lee, R. Ramesh, J. Li, B. Nagaraj, H. Liang, W. Cao, C. H. Lee, and R. Ramesh, Ultrafast Polarization Switching in Thin-Film Ferroelectrics Ultrafast Polarization Switching in Thin-Film Ferroelectrics, 1174, 3 (2007). [77] A. K. Tagantsev, I. Stolichnov, N. Setter, J. S. Cross, and M. Tsukada, Non-Kolmogorov- Avrami Switching Kinetics in Ferroelectric Thin Films, 1 (2002). [78] J. F. Scott, A Review of Ferroelectric Switching, Ferroelectrics 503, 117 (2016). [79] Y. W. Li, J. F. Scott, D. N. Fang, F. X. Li, Y. W. Li, J. F. Scott, D. N. Fang, and F. X. Li, 90-Degree Polarization Switching in BaTiO3 Crystals without Domain Wall Motion 90- Degree Polarization Switching in BaTiO 3 Crystals without Domain Wall Motion, 232901, 0 (2014). [80] K. Laabidi, M. D. Fontana, M. Maglione, B. Jannot, and K. A. Müller, Indications of Two Separate Relaxators in the Subphonon Region of Tetragonal BaTiO 3, Europhys. Lett. 26, 309 (1994). [81] Y. K. Vayunandana Reddy and D. Mergel, Frequency and Temperature-Dependent Dielectric Properties of BaTiO3 Thin Film Capacitors Studied by Complex Impedance Spectroscopy, Phys. B Condens. Matter 391, 212 (2007). [82] M. Zhang and C. Deng, Orientation and Electrode Configuration Dependence on Ferroelectric, Dielectric Properties of BaTiO3 Thin Films, Ceram. Int. 45, 22716 (2019). [83] A. Piorra, A. Petraru, H. Kohlstedt, M. Wuttig, and E. Quandt, Piezoelectric Properties of 0.5(Ba0.7Ca0.3TiO3) - 0.5[Ba(Zr0.2Ti0.8)O3] Ferroelectric Lead-Free Laser Deposited Thin Films, J. Appl. Phys. 109, 3 (2011). Page 30 of 31 [84] D. Xue, Y. Zhou, H. Bao, C. Zhou, J. Gao, and X. Ren, Elastic, Piezoelectric, and Dielectric Properties of Ba(Zr Ti )O -50(Ba Ca )TiO Pb-Free Ceramic at 0.2 0.8 3 0.7 0.3 3 the Morphotropic Phase Boundary, J. Appl. Phys. 109, 054110 (2011). [85] A. D. Dupuy, Y. Kodera, and J. E. Garay, Unprecedented Electro-Optic Performance in Lead-Free Transparent Ceramics, Adv. Mater. 28, 7970 (2016). [86] G. Keiser, Optical Fiber Communications, in Wiley Encyclopedia of Telecommunications (John Wiley & Sons, Inc., Hoboken, NJ, USA, 2003). Page 31 of 31
Journal
Physics – arXiv (Cornell University)
Published: Jun 6, 2021