Brillouin Scattering and First-Principles Studies of BaMO3 (M = Ti, Zr, and Cu) Perovskites

Perovskite oxides with the general formula ABO3 comprise a large number of families among the structures of oxide-based materials, and currently, several perovskite structures have been identified. From a variety of compositions and structures, various functions are observed in perovskite compounds, and therefore, they became very useful for various applications in the electronic and medical industries. One of the most puzzling issues for perovskite compounds is the understanding of the vibration and relaxation dynamics in the gigahertz range. In that sense, the micro-Brillouin scattering system is a very effective tool to probe the gigahertz dynamics, and also, first-principles calculations can be used to describe the phonon structure with different atomic contributions. The micro-Brillouin scattering system and first-principles calculations provide the fundamental information on a variety of vibration and relaxation processes related to structural phase transitions under different external conditions such as temperature, electric field, and pressure. This review article summarizes the Brillouin scattering and first-principles studies on BaMO3 (M = Ti, Zr, and Cu). Through a detailed analysis of the existing results, we summarize the existing limitations and future perspectives in these research areas, which may propel the development of different perovskite ferroelectrics and extend their practical application areas.


Introduction
In solid-state physics, crystals with the ABO 3 -type (A: Large cation with different valence state and B: Transition metal) perovskite family structure are being widely studied because they put up most of the metal ions in the periodic table owing to its substantial number of different ions [1]. These types of oxides are significant in the field of various technological applications in modern electronics, non-linear optics, catalysis, refractories, geophysics, particle accelerators, etc. [2][3][4][5][6]. One of the major benefits of perovskite-structured oxides is that they can accept considerable substitutions in one or both cationic sites (A and B sites) while retaining their original crystal structures [7]. Currently, these types of compounds have several outstanding applications in the electronic industries such as piezoelectric devices, transducers, random access memories, tunable microwave device displays, sensors, capacitors, and wireless communications [8][9][10][11][12].
The ideal structure of perovskite is cubic, which is shown in Figure 1. Although few compounds have this pure cubic phase, many perovskites have a slightly distorted phase with lower symmetry (e.g., tetragonal, orthorhombic, or rhombohedral). The lattice distortion in a crystal results from the off-center and/or displacement of the A or B site cations, which cause a ferroelectric phase transition owing to the polar phonon mode softening at a Γ point of the Brillouin zone. Meanwhile, the softening of the phonon mode can also appear due to the rotation of the oxygen octahedron at the R or M point of the Brillouin zone boundary [13]. Furthermore, many oxygen-or cation-deficient structures have been observed. The different categories of distortion in the perovskite materials Since its discovery, barium titanate (BaTiO 3 (BTO)) has been studied extensively due to its remarkable ferroelectric properties, and it has been the subject of many studies for more than 50 years. Barium zirconate (BaZrO 3 (BZO)) has attracted much attention from the scientific community due to its possible applications in hydrogen sensors and fuel cells. Barium cuprate (BaCuO 3 (BCO)) is another promising perovskite material, which can be used to improve catalytic activity. The process of phase transition is different in these three crystals.
BTO is stable in a paraelectric cubic phase (space group Pm3m) at high temperatures with five atoms per unit cell. Upon cooling, BTO exhibits three successive phase transitions from the paraelectric cubic phase to ferroelectric tetragonal (space group P4mm), orthorhombic (space group Amm2), and rhombohedral (R3m) phases at 393, 278, and 183 K, respectively. There exist numerous experimental and theoretical studies explaining the ferroelectric behavior of BTO single crystals. In spite of numerous reports, still, the true nature of the ferroelectric phase transition in BTO is controversial. The ferroelectricity and concomitant cubic to tetragonal distortion were initially discussed on the basis of displacive-type phase transition from the condensation of the soft optic mode [14,15]. However, with the lowering of temperature, the subsequent phase transitions of the orthorhombic and rhombohedral phases are inconsistent with the displacive nature [16,17].
This anomaly together with the results of low-temperature diffuse X-ray scattering [18] results in the development of a phenomenological order-disorder model [19]. In addition, the deviation of the refractive index from high-temperature linear behavior [20], first-order Raman spectra [21], the Curie-Weiss law for the permittivity [22], and significant softening of the longitudinal acoustic phonon [23,24] in the paraelectric cubic phase were observed in BTO crystals, which were ascribed to the existence of local polar clusters. Although BTO single crystals were studied more than 50 years, there is a significant lack of some decisive properties concerning the existence of the local polar cluster and the effect of external stress on this cluster, which must be attained to unlock the full potential of BTO for possible future technological applications.
The lattice dynamics and structural stability of BZO are found to be more complex, which has drawn great attention owing to its low-temperature structure and properties. As other perovskites, so far, no phase transition has been noticed in BZO at ambient pressure at any temperature. BZO has a cubic structure up to 2K according to neutron and X-ray diffraction studies [25]. However, these structural data disagree with the first-principles calculations [26,27], where the instability of the R or M phonon mode is observed due to the rotations of oxygen octahedra. In addition, a new transverse optic mode at about 2 THz was observed by the terahertz time-domain spectroscopy technique, which was forbidden in the consideration of cubic symmetry according to the infrared (IR) selection rule [28]. The origin of this peak was attributed to an IR active mode of the P1 symmetry caused by local distortion. In BZO crystals, the first-order Raman active modes are observed, which is unfavorable for the cubic Pm3m symmetry, indicating the lowering of the cubic symmetry [29,30]. However, the observed Raman spectra in BZO are consistent with the local tetragonal and orthorhombic symmetry [30][31][32].
One of the most interesting properties of perovskite materials is superconductivity. The discovery of superconductivity in La-Ba-Cu-O perovskite oxide with the K 2 NiF 4 -type layered structure containing the CuO 2 plane was observed by Bednorz and Müller in 1984 [33]. After their report, significant attention has been paid to new types of hightemperature superconductors, mainly Cu-based perovskite oxides due to their excellent physical properties. It has been found that the presence of Cu on the B site cation is needed for superconductivity to occur in such perovskites. Due to the tremendous effort, a high value of transition temperatures (T c ) was found in the YBa 2 Cu 3 O 7 [34] and Bi 2 Sr 2 Ca 2 Cu 3 O 10 [35] systems in 1987 and 1988, respectively, and currently, T c has been further increased to 130-155 K in the HgBa 2 Ca 2 Cu 3 O 8+δ system [36]. So far, all the hightemperature superconductors are Cu-based oxides, and T c is clearly associated with the Cu-O chain. In addition to superconductivity, Cu-based perovskites exhibit high electrical conductivity and catalytic activity, which are close to those of metals such as Cu. Very recently, another typical Cu-based superconducting BaCuO 3 (BCO) perovskite has been reported [37,38], and T c was found to be around 10 k.
In conventional dielectric spectroscopy measurement, it is very difficult to measure the broadband spectra between 1 and 1000 GHz. However, Brillouin scattering techniques enable the acquisition of such broad spectra by a single scan. To observe the dynamics of local polarization fluctuation in perovskite-type ferroelectrics, Brillouin scattering has been used as a very significant tool due to its sensitivity to local symmetry breaking related to the heterogeneity among a few methods, which can be used to directly identify the local structure. Another advantage of Brillouin scattering is that the spatial resolution is comparable to the wavelength of an excited laser, and it is possible to measure the acoustic phonon and optic soft mode in different types of perovskites.
In several applications, perovskites have been proven to be of great interest because of their meaningful physical properties for different device applications. In the present paper, we review the recent systematic studies on BaMO 3 (M = Ti, Zr, and Cu) perovskites and discuss the existing limitations and future perspectives and research in these areas.

Experimental Section
In solid-state physics, the determination of the velocity of an acoustic wave is given by the elastic constant and density. There are many methods to determine the elastic constants. Among them, the ultrasonic pulse-echo method and the mechanical resonance method have been well established [39,40]. However, as a non-contact and non-destructive method, Brillouin scattering spectroscopy is a powerful tool to measure the elastic constants with different external conditions, for example pressure, temperature, and electric field. Such a spectroscopic technique covers a large frequency range between 1 and 1000 GHz, and recently, it has been used to study not only the acoustic phonons, but also the optical soft mode [41], magnetic excitation [42], and localized excitation, such as fracton [43].
The frequency shift υ B is proportional to the phase velocity V of the sound wave, which is measured by the following equation: where λ 0 , n m , and θ are the wavelength of incident light, the refractive index of the material to be studied, and the scattering angle, respectively. The optical arrangement of a micro-Brillouin scattering apparatus is shown in Figure 2 [44]. The tandem multi-pass Fabray-Pérot interferometer, (FP1 and FP2), was designed by Sandercock (JRS, Zwillikon, Switzerland). The laser is focused on the sample with the help of a small mirror in the optical microscope arrangement. The scattered light enters FP1 via pinhole P1 and aperture A1 and is collected by the objective lens. The laser light is then directed via mirror M3 to FP2. The light strikes the 90 • -oriented prism PR1 after the transmission through FP2, then is reflected backwards and returns to FP2. It continues via aperture A2 to FP1. The light is then passed via L1 and is focused onto mirror M4 with mirror M1 underneath, after the transmission through FP1. Mirror M4 then returns the laser light back to lens L1, and the light is projected to FP1 through M2 after collimation. The light is directed to PR2 through M5 after the final pass through FP2. A bandpass filter is formed by the combination of lens L2, prism PR2, and output pinhole P2. The light is projected onto the photomultiplier (PM) from mirror M6 via aperture A3, lens L2, and output pinhole P2. Finally, a conventional photon counter records the scattered signal.

Theoretical Section
The density functional theory (DFT) formalism is the basis of first-principles calculations upon which the different properties of a material can be derived. Using DFT, without solving the multibody Schrödinger equation, the problem can be described with an equivalent non-interacting problem. The ground state structure is obtained after solving the Kohn-Sham equations [46]. The ultrasoft pseudopotential, given by Vanderbilt, To measure the micro-Brillouin scattering spectra, a 3 + 3 pass Sandercock-type tandem Fabray-Pérot interferometer (JRS TEP-1, JRS, Zwillikon, Switzerland), was used with a photon counting system. All the Brillouin scattering spectra were collected in a backwards scattering geometry [44,45]. The objective lens of the microscope focuses an incident laser beam with a wavelength of 532 nm and a power of about 100 mW onto the samples, and simultaneously, the scattered light is collected. The mirror space was set at 2 mm with a free spectral range of 75 GHz. The temperature of the sample was controlled by a heating and/or a cooling stage (THMS600, Linkam, Salfords, UK) with a stability of ±0.1 K. Two surfaces of the sample were coated with silver paste, and gold contact wires were attached to the electrodes to apply the electric field.

Theoretical Section
The density functional theory (DFT) formalism is the basis of first-principles calculations upon which the different properties of a material can be derived. Using DFT, without solving the multibody Schrödinger equation, the problem can be described with an equivalent non-interacting problem. The ground state structure is obtained after solving the Kohn-Sham equations [46]. The ultrasoft pseudopotential, given by Vanderbilt, was used to optimize the geometry, because it allows for a lower cut-off energy, without sacrificing transferability too much [47]. The ultrasoft pseudopotential method performs the calculations of the material by considering the electron density in the valance as the soft part and the core region as the hard part, where the cut-off energy is dramatically reduced and the computational efficiency increased appreciably. The bulk modulus, independent elastic constants, and shear modulus were estimated using the stress-strain method contained in the Cambridge Serial Total Energy Package (CASTEP) code. On the other hand, the normconserving pseudopotential was used for the phonon calculations. The system geometry was optimized under the energy convergence at a total energy of 0.5 × 10 −5 eV/atom, the maximum interatomic force 0.01 eV/Å, the maximum stress 0.02 GPa, and an ionic displacement of 5 × 10 −4 Å. Understanding the energy-dependent optical parameters is very crucial to explore the potential application of a new compound. The optical properties are very important from the viewpoint of a possible optoelectronic device applications. The optical constants' spectra were extracted from the frequency-dependent dielectric function, ε(ω) = ε 1 (ω) + iε 2 (ω). The real part of the dielectric constant, ε 1 (ω), is calculated from the imaginary part, ε 2 (ω), using Kramers-Kronig transformations. The imaginary dielectric function, ε 2 (ω), can be calculated through the following equation: Here, ȓ is the electron radius vector, û is the polarization field, e is the charge of the electron,  is the frequency conduction and valence band wave functions at k, respec

Precursor Dynamics and Domain Switching in BaTiO3 Sin
BaTiO3 (BTO) is the first produced and most studi superior dielectric, ferroelectric, and piezoelectric prop constants of BTO results from its crystal structure. In B quencies were obtained from the FTIR measurements [49 482 cm −1 . The lowest frequency phonon mode is called th in a paraelectric phase. From the measured frequency, L dielectric constant using the Lyddane-Sachs-Teller (LST) dynamics, the dielectric constant increases markedly in tion temperature upon cooling. The obtained dielectric c by approaching the Curie temperature, while the soft m constant. To explain this discrepancy, two possible origin comes from the contribution of impurities, which gives dielectric constant due to the coupling with the soft mod istence of an additional relaxational mechanism in the mi had no direct conclusion about this point. Therefore, it is ational mechanism by considering another approach. In ion is too large compared with the Ti ion, which is quite Here, Here, ȓ is the electron radius vector, û is the polarization vector of the incident electric field, e is the charge of the electron,  is the frequency of light, and Ψ k c and Ψ k v are the conduction and valence band wave functions at k, respectively.

Precursor Dynamics and Domain Switching in BaTiO3 Single Crystal
BaTiO3 (BTO) is the first produced and most studied perovskite material with its superior dielectric, ferroelectric, and piezoelectric properties [48]. The high dielectric constants of BTO results from its crystal structure. In BTO crystals, three phonon frequencies were obtained from the FTIR measurements [49], located at around 40, 180, and 482 cm −1 . The lowest frequency phonon mode is called the soft mode and is overdamped in a paraelectric phase. From the measured frequency, Luspin et al. [49] estimated the dielectric constant using the Lyddane-Sachs-Teller (LST) relation. Due to the soft mode dynamics, the dielectric constant increases markedly in the vicinity of the phase transition temperature upon cooling. The obtained dielectric constant gives a diverging value by approaching the Curie temperature, while the soft mode frequency remains nearly constant. To explain this discrepancy, two possible origins can be invoked. The first one comes from the contribution of impurities, which gives rise to the rapid increase of the dielectric constant due to the coupling with the soft mode. Another possibility is the existence of an additional relaxational mechanism in the microwave region. However, they is the electron radius vector, û is the polarization vector of the incident electric field, e is the charge of the electron, ω is the frequency of light, and Ψ c k and Ψ v k are the conduction and valence band wave functions at k, respectively.

Precursor Dynamics and Domain Switching in BaTiO 3 Single Crystal
BaTiO 3 (BTO) is the first produced and most studied perovskite material with its superior dielectric, ferroelectric, and piezoelectric properties [48]. The high dielectric constants of BTO results from its crystal structure. In BTO crystals, three phonon frequencies were obtained from the FTIR measurements [49], located at around 40, 180, and 482 cm −1 . The lowest frequency phonon mode is called the soft mode and is overdamped in a paraelectric phase. From the measured frequency, Luspin et al. [49] estimated the dielectric constant using the Lyddane-Sachs-Teller (LST) relation. Due to the soft mode dynamics, the dielectric constant increases markedly in the vicinity of the phase transition temperature upon cooling. The obtained dielectric constant gives a diverging value by approaching the Curie temperature, while the soft mode frequency remains nearly constant. To explain this discrepancy, two possible origins can be invoked. The first one comes from the contribution of impurities, which gives rise to the rapid increase of the dielectric constant due to the coupling with the soft mode. Another possibility is the existence of an additional relaxational mechanism in the microwave region. However, they had no direct conclusion about this point. Therefore, it is necessary to observe this relaxational mechanism by considering another approach. In BTO crystals, the size of the Ba ion is too large compared with the Ti ion, which is quite small and remains stable in the octahedral position. The Ti ion having a charge of +4 has the minimum energy position outside the center, in the direction of each of the O 6 ions surrounding it. Due to the off-centering position of the Ti ions, the mechanism of the phase transition, in particular, the order-disorder-and/or displacive-type, is still under discussion. It was shown by Takahasi [50] and Comes et al. [18] that Ti ions go to the off-center position along eight equivalent [111] directions inside the cubic matrix with a non-zero local dipole moment. Therefore, this leads to a local polar nanoregion (PNR) in a cubic phase, which appears a few hundred degrees above the cubic to tetragonal phase transition temperature (T c ), called the Burns temperature (T B ). It is believed that PNRs play an important role in the appearance of the outstanding piezoelectric performance in electrical devices. In order to identify such types of nano-metric local polar regions, a low-frequency inelastic light scattering spectrum can be a very effective tool, because it can be used to directly probe the local symmetry and structure. Figure 3 shows a group of Brillouin scattering spectra at 398 and 403 K without any applied electric fields. dynamics, the dielectric constant increases markedly in the vicinity of the phase tr tion temperature upon cooling. The obtained dielectric constant gives a diverging v by approaching the Curie temperature, while the soft mode frequency remains n constant. To explain this discrepancy, two possible origins can be invoked. The firs comes from the contribution of impurities, which gives rise to the rapid increase o dielectric constant due to the coupling with the soft mode. Another possibility is th istence of an additional relaxational mechanism in the microwave region. However, had no direct conclusion about this point. Therefore, it is necessary to observe this r ational mechanism by considering another approach. In BTO crystals, the size of th ion is too large compared with the Ti ion, which is quite small and remains stable i octahedral position. The Ti ion having a charge of +4 has the minimum energy pos outside the center, in the direction of each of the O6 ions surrounding it. Due to th centering position of the Ti ions, the mechanism of the phase transition, in particula order-disorder-and/or displacive-type, is still under discussion. It was show Takahasi [50] and Comes et al. [18] that Ti ions go to the off-center position along equivalent [111] directions inside the cubic matrix with a non-zero local dipole mom Therefore, this leads to a local polar nanoregion (PNR) in a cubic phase, which appe few hundred degrees above the cubic to tetragonal phase transition temperature called the Burns temperature (TB). It is believed that PNRs play an important role i appearance of the outstanding piezoelectric performance in electrical devices. In ord identify such types of nano-metric local polar regions, a low-frequency inelastic scattering spectrum can be a very effective tool, because it can be used to directly p the local symmetry and structure. Figure 3 shows a group of Brillouin scattering sp at 398 and 403 K without any applied electric fields.  It is pertinent to note that a longitudinal acoustic (LA) mode and a weak central peak (CP) appear in the Brillouin scattering spectrum above T c (402 K). It is now believed that the CP appears in low-frequency light scattering spectra owing to the appearance of random fields, which originate from the displacement of Ti ions inside the oxygen octahedra [23,51]. It was also suggested from different experiments [52,53] that the motions of off-center Ti ions are correlated and create PNRs in which all the Ti ions are oriented in the same direction and fluctuate together. Above T c , previous theoretical reports also supported the creation of PNRs with local symmetry breaking [54,55]. Therefore, the presence of the CP in BTO crystals can be attributed to the polarization fluctuation of PNRs, which has also been observed in BTO-based ferroelectrics. However, to make further conclusions, it is very important to measure the broad CP by means of different experimental techniques. Ko et al. [56] measured the temperature dependence of the broad CP and discussed the critical slowing down behavior of the CP. They explained their result by the existence of polar clusters and the growth of these clusters with longer correlation lengths a few ten degrees above T c . To make still further conclusions, besides the temperature dependence, the electric field dependence of the broad CP is very important, which can reveal some unanswered questions. So far, there is no available result on the electric field dependence of the broad CP. Therefore, future experiments will test the growth of PNRs by measuring the broad CP by applying an electric field in different directions.
Another approach to the effect of PNRs is to observe the behavior of LA phonons in the zero field cooling (ZFC) and field cooling (FC) process. At first, the Brillouin spectra are measured without applying the electric field, called the ZFC process. To apply the electric field, the sample is heated above T B and short circuited to remove any kind of memory effects. Then, an electric field of 1 kV/cm at a high temperature is applied, and the Brillouin spectra upon cooling are collected, called the FC process. After measuring the Brillouin scattering spectra in ZFC and FC, they are fit by a Lorentzian function convoluted by a Gaussian instrumental function from which the LA frequency shift (ϑ LA ) and the FWHM (Γ LA ) can be obtained. Figure 4a and b shows the ϑ LA and Γ LA in the ZFC and FC processes, respectively. A sudden change of both ϑ LA and Γ LA at around 402 K reflects the transition from the cubic to tetragonal phase.
CP and discussed the critical slowing down behavior of the CP. They explained their re sult by the existence of polar clusters and the growth of these clusters with longer corr lation lengths a few ten degrees above Tc. To make still further conclusions, besides th temperature dependence, the electric field dependence of the broad CP is very im portant, which can reveal some unanswered questions. So far, there is no available resu on the electric field dependence of the broad CP. Therefore, future experiments will te the growth of PNRs by measuring the broad CP by applying an electric field in differen directions.
Another approach to the effect of PNRs is to observe the behavior of LA phonons i the zero field cooling (ZFC) and field cooling (FC) process. At first, the Brillouin spectr are measured without applying the electric field, called the ZFC process. To apply th electric field, the sample is heated above TB and short circuited to remove any kind o memory effects. Then, an electric field of 1 kV/cm at a high temperature is applied, an the Brillouin spectra upon cooling are collected, called the FC process. After measurin the Brillouin scattering spectra in ZFC and FC, they are fit by a Lorentzian function con voluted by a Gaussian instrumental function from which the LA frequency shift ( and the FWHM (Γ LA ) can be obtained. Figure 4a and b shows the and Γ in the ZF and FC processes, respectively. A sudden change of both and Γ at around 402 K re flects the transition from the cubic to tetragonal phase. Here, the interesting feature is that the LA frequency shift begins to deviate from its linear temperature dependence at around 533 K. It has been reported that the local strain fluctuation is created by the local polarization fluctuation inside a PNR. The macroscopic elastic anomaly and the piezoelectric coupling in a PNR occur due to the coupling between the local strain fluctuation in a PNR and acoustic phonons [57]. Therefore, we assigned this temperature as the Burns temperature, T B , where the dynamic PNRs begin to appear. On the other hand, the sound attenuation, related to Γ LA , begins to increase abruptly towards the transition temperature at around 503 K, where the permanent PNRs arise, and the temperature is known as the intermediate temperature, T*. The observed T B and T* in this study were somewhat higher than previously found by the Brillouin scattering results [56] carried out on BTO single crystals. The discrepancies of these results may appear because of the quality of the sample, because the authors noticed T c at 373 K, which is far behind the actual T c in BTO crystals. Interestingly, these two temperatures were also measured by the acoustic emission measurement [58] carried out on a single crystal of BTO, and the reported values of T B were approximately between 530-570 K, while T * was 506 K, which is in good agreement with the present Brillouin scattering results. Very recently, T B was also measured by the resonant ultrasound spectroscopy (RUS) technique [59]. The onset of the frequency vs. temperature curve in the RUS method shows the deviation from the exponential softening towards a weaker temperature dependence near 560 K, which they denoted as T B . These two characteristic temperatures have been noticed in many relaxor ferroelectrics [18,60]. However, so far, this temperature has not been clearly identified [58].
The electric field has a serious effect on the crystal structure and polarization properties of BTO. Under the applied electric field, Ti ions shift from a random position to aligned positions, obtaining high polarization and a high dielectric constant [61]. To analyze the effect of an electric field, we applied a 1 kVcm −1 electric field at a higher temperature (above T B ), and then, the Brillouin spectra were taken in the field cooling process. A significant difference was noticed between ϑ LA and Γ LA in both the cubic and tetragonal phases. It has already been discussed and reported that in BTO crystals, the Ti ions fluctuate along eight equivalent [111] directions inside a PNR [53]. Owing to the application of an electric field along the [100] direction, the density of dynamic/static PNRs can be decreased, which is expected to weaken the LA phonons' scattering by PNRs, resulting in a higher frequency shift and lowering the damping of the LA phonons. However, it is significant to note that T B is not influenced by the applied electric field, and T* and T c change after applying electric field. This observation is in good agreement with acoustic emission results [58]. According to acoustic emission results, the electric field dependence is stronger in T* than T c due to the possible field ordering of PNR dipoles, which could favor their freezing. In the tetragonal phase, a significant change between ϑ LA and Γ LA has been observed. Due to the sufficiently strong electric field, the correlation among the PNRs becomes sluggish, and the alignment of nanodomains/static PNRs is facilitated and enables complete switching into the macro/single-domain state at T c , giving rise to the formation of the ferroelectric state. In addition, the large differences in ϑ LA and Γ LA between ZFC and FC indicate the transformation from the ferroelectric [100] direction to the [001] direction in the tetragonal phase. Upon the subsequent heating (heating after field cooling), the crystals transform back to their initial state, which results in hysteresis below T c .
In ferroelectrics, to understand the domain switching process, many research works have been performed and remarkable progress has been achieved [62][63][64][65][66]. Previously, it was reported that under an applied electric field antiparallel to the polar direction, only 180 • domain switching can be possible in ferroelectrics. However, non-180 • domain switching was also suggested by both the theoretical and experimental observations [67,68]. Under the application of an external electric field, now, it has been established that two types of domains can nucleate: 90 • domains with polarization perpendicular to each other. which minimizes strain via the formation of twin walls, and 180 • domains having antiparallel spontaneous polarization, which diminishes the depolarization fields [69]. To observe the domain switching effect, a set of typical Brillouin scattering spectra in the tetragonal phase at 303 K is shown in Figure 5.
When there is no electric field, the transverse acoustic (TA) mode is denoted as the a-domain, since the sample is cut along the [100]-oriented plate. When the field is increased up to 2.94 kV/cm, another TA mode at around 30 GHz begins to appear and persists up to 3.9 kV/cm. The additional TA mode in BTO crystals means the appearance of the c-domain. Within 2.94 and 3.9 kV/cm, both domains co-exist. Above 3.9 kV/cm, the intensity of the additional TA mode goes to the maximum and saturates, suggesting that the 90 • domain switching has completed. The 90 • domain switching was also noticed by Li et al. [70] and Salje et al. [58]. Moreover, whether the enhanced polar properties are related to the engineered domain structures or whether the increased density of the domain boundaries between domains plays an important role in BTO crystals is an open question. switching was also suggested by both the theoretical and experimental observation [67,68]. Under the application of an external electric field, now, it has been establishe that two types of domains can nucleate: 90° domains with polarization perpendicular t each other. which minimizes strain via the formation of twin walls, and 180° domain having anti-parallel spontaneous polarization, which diminishes the depolarizatio fields [69]. To observe the domain switching effect, a set of typical Brillouin scatterin spectra in the tetragonal phase at 303 K is shown in Figure 5. When there is no electric field, the transverse acoustic (TA) mode is denoted as th a-domain, since the sample is cut along the [100]-oriented plate. When the field is in creased up to 2.94 kV/cm, another TA mode at around 30 GHz begins to appear and per sists up to 3.9 kV/cm. The additional TA mode in BTO crystals means the appearance o the c-domain. Within 2.94 and 3.9 kV/cm, both domains co-exist. Above 3.9 kV/cm, th intensity of the additional TA mode goes to the maximum and saturates, suggesting tha the 90° domain switching has completed. The 90° domain switching was also noticed b Li et al. [70] and Salje et al. [58]. Moreover, whether the enhanced polar properties are re lated to the engineered domain structures or whether the increased density of the do main boundaries between domains plays an important role in BTO crystals is an ope question.

Local Lattice Disorder in BaZrO 3 Single Crystals
The perovskite structure of BZO single crystals, long known to maintain their cubic symmetry at room temperature without any structural phase transition with wide temperature variations, has recently been exhibited to have various ground state structures with low symmetries involving octahedral rotation. However, the ground state structure in BZO is questionable. Several first-principles calculations reported that the true structural state of BZO may be tetragonal I4/mcm because of the unstable antiphase antiferrodistortive phonons of oxygen octahedra [69,71]. These local tetragonal distortions can be averaged out and may not be possible to detect using X-ray or neutron diffraction experiments, which results in the cubic symmetry. From the quantum mechanical viewpoint, Brillouin scattering is assumed as one of the interactions of light photons with acoustic or vibrational quanta such as phonons. No coupling between the ultrasonic transducer and other devices with samples is needed to produce acoustic excitation and detection in Brillouin scattering, which leads to accurately measuring the acoustic properties. Figure 6a shows a set of Brillouin scattering spectra at some selected temperatures, and a contour color map is displayed in Figure 6b. iments, which results in the cubic symmetry. From the quantum mechanical viewpoint, Brillouin scattering is assumed as one of the interactions of light photons with acoustic or vibrational quanta such as phonons. No coupling between the ultrasonic transducer and other devices with samples is needed to produce acoustic excitation and detection in Brillouin scattering, which leads to accurately measuring the acoustic properties. Figure  6a shows a set of Brillouin scattering spectra at some selected temperatures, and a contour color map is displayed in Figure 6b. We observed clearly an LA and a TA phonon mode. From the temperature dependence of these phonon modes, no sharp discontinuity related to a phase transition was observed. Therefore, macroscopically, it remains as the ground state structure at room temperature. However, the most striking feature of this Brillouin spectrum is the appearance of the TA mode, which exists up to 1273 K. The Brillouin spectra were recorded We observed clearly an LA and a TA phonon mode. From the temperature dependence of these phonon modes, no sharp discontinuity related to a phase transition was observed. Therefore, macroscopically, it remains as the ground state structure at room temperature. However, the most striking feature of this Brillouin spectrum is the appearance of the TA mode, which exists up to 1273 K. The Brillouin spectra were recorded in a backward scattering geometry, and only the LA mode should propagate along the [100] direction based on the selection rules in the cubic symmetry [72]. Due to this, several previous experiments noticed only the LA phonon mode in the cubic symmetry [73]. Therefore, the strong TA mode, not allowed in the cubic symmetry, reflects the presence of local lattice distortion, which can destroy the center of the symmetry in a cubic matrix. This local lattice distortion is also evident from terahertz time-domain spectra (THz-TDS). Figure 7a In cubic perovskite with the space group Pm3m, three IR active transverse optic (TO) modes are allowed according to the selection rule. Two former measurements on the IR phonon spectra have been reported [74,75]. In addition, the IR phonon frequencies have been reported by first-principles calculations [76,77]. In comparison with those published in BZO ceramics by IR and first-principles calculations, the observed TO mode was observed in THz-TDS at a very low frequency, which reflects the additional new mode in a BZO crystal. The FTIR spectra of BZO ceramics [74] show three narrow phonon peaks, as expected for a simple cubic perovskite structure. Among three phonon modes, the TO2 mode was the strongest one and was assigned as the Slater mode (Zr-O 6 stretching vibrations). The eigenvector of the lowest-frequency TO1 mode was close to the last mode, playing the role od the ferroelectric soft mode. The highest-frequency mode above 15 THz was due to the O 6 deformation and called the Axe mode. In addition, to fit the FTIR spectra, Nuzhnyy et al. [74] used the nine oscillator models and concluded that not all the oscillators have a clear meaning. However, the appearance of a new mode at 2 THz in a BZO crystal and an additional mode at 9 THz in BZO ceramics [74] indicated that BZO has a lower symmetry than cubic PmPm3m and/or these modes come from the local lattice distorted structure. Very recently, D.H. Gim et al. measured the Raman and IR spectra of both the ceramics and single crystals of BZO [75]. They also observed the additional phonon modes in both samples and reported that these additional phonon modes are due to the presence of local distortion or disorder. Therefore, the additional phonon modes in Brillouin scattering, THz-TDS, Raman spectra [75], as well as IR spectra [74] can be considered evidence of the presence of the local lattice distorted structure in BZO. in a backward scattering geometry, and only the LA mode should propagate along the [100] direction based on the selection rules in the cubic symmetry [72]. Due to this, several previous experiments noticed only the LA phonon mode in the cubic symmetry [73]. Therefore, the strong TA mode, not allowed in the cubic symmetry, reflects the presence of local lattice distortion, which can destroy the center of the symmetry in a cubic matrix. This local lattice distortion is also evident from terahertz time-domain spectra (THz-TDS). Figure 7a  In cubic perovskite with the space group Pm3 ̅ m, three IR active transverse optic (TO) modes are allowed according to the selection rule. Two former measurements on the IR phonon spectra have been reported [74,75]. In addition, the IR phonon frequencies have been reported by first-principles calculations [76,77]. In comparison with those published in BZO ceramics by IR and first-principles calculations, the observed TO mode was observed in THz-TDS at a very low frequency, which reflects the additional new mode in a BZO crystal. The FTIR spectra of BZO ceramics [74] show three narrow phonon peaks, as expected for a simple cubic perovskite structure. Among three phonon modes, the TO2 mode was the strongest one and was assigned as the Slater mode (Zr-O6 stretching vibrations). The eigenvector of the lowest-frequency TO1 mode was close to the last mode, playing the role od the ferroelectric soft mode. The highest-frequency

Lattice Instability and IR Spectra of BaZrO 3 Crystals
Tolerance factor t is a well-established parameter, which is used to identify the stable structure. The formula for t of ABO 3 -type perovskite is given by Here, R A−O is the sum of A and O ionic radii and R B−O is the sum of B and O ionic radii. Generally, t < 1 leads to the rotation and expansion of B-O 6 octahedra. A lowtemperature anti-ferroelectric phase is generated due to the octahedral rotations. If t > 1, the B-O 6 octahedra are stretched from their preferred B-O bond lengths, resulting in the distortions of the B cation by creating room for the B cations to move to the off-center position. In BZO crystals, the values of the ionic radii of Ba and Zr atoms are 1.61 Å and 0.72 Å, respectively, which leads to t = 1, so the cubic structure should have no driving force to deform to a lower symmetry. However, the calculation of the phonon spectrum of BZO using DFT calculations is unfavorable to achieve the stable structure of BZO. The phonon dispersion relation for the five-atom unit cell of BZO is shown in Figure 8. The very important characteristics of the phonon spectra are related to the lowest-frequency acoustic phonon branch, which goes to negative at the R-and M-points in the Brillouin zone. However, theses unstable phonon modes were absent in. [78,79], but appeared in [69], where all the calculations were performed using the DFT method. The negative phonon frequencies correspond to the imaginary frequencies, which implies the change of the crystal symmetry. Due to the anharmonic phonon interaction, the negative phonon frequencies may also appear, which are responsible for a structural phase transition of BZO with respect to some external conditions. Therefore, BZO is dynamically unstable, which can change in phase under pressure or an electric field. This observation is consistent with the present Brillouin scattering and THz-TDS results. Indeed, under the applied pressure, very recently, the cubic to rhombohedral R3c and rhombohedral to tetragonal I4/mcm phase transitions at 8.4 and 19.2 GPa, respectively, were observed [75]. A few years before, Yang et al. also observed the pressure dependence of the phase transition from the cubic to tetragonal I4/mcm at 17.2 GPa [70]. We also observed the pressure dependence of the phase transition by first-principles calculations [13]. The atomic contribution in the phonon spectra can be seen from the partial phonon density of states (PPDOS), which is shown in Figure 9. It can be noticed from the PPDOS that the negative phonon frequency appears due to the vibration of Zn and O atoms. Therefore, it can be concluded that within the frequency range from 0.5 to 3.0 THz, an unstable phonon mode at the R-and/or M-points of the Brillouin zone is responsible for the phase transition and the superior dielectric and piezoelectric responses of BZO, which correspond to the tilts of the ZrO6 octahedra. Among the perovskite ferroelectric oxides such as PZT [80], this type of unstable phonon mode is observed commonly at the R-or M-points of the Brillouin zone. The atomic contribution in the phonon spectra can be seen from the partial phonon density of states (PPDOS), which is shown in Figure 9. It can be noticed from the PPDOS that the negative phonon frequency appears due to the vibration of Zn and O atoms. Therefore, it can be concluded that within the frequency range from 0.5 to 3.0 THz, an unstable phonon mode at the Rand/or M-points of the Brillouin zone is responsible for the phase transition and the superior dielectric and piezoelectric responses of BZO, which correspond to the tilts of the ZrO 6 octahedra. Among the perovskite ferroelectric oxides such as PZT [80], this type of unstable phonon mode is observed commonly at the Ror M-points of the Brillouin zone.
Finally, we also estimated the IR spectra of BZO using first-principles calculations. According to the group theoretical analysis, due to the cubic symmetry, three phonon branches of the T 1u symmetry is IR active, and one branch of the T 2u symmetry is not IR, nor Raman active. The estimated IR spectra of the BZO single crystal are displayed in Figure 10 where the frequencies of the three TO modes are 2.4, 5.3, and 14.8 THz, respectively. The atomic contribution in the phonon spectra can be seen from the partial phonon density of states (PPDOS), which is shown in Figure 9. It can be noticed from the PPDOS that the negative phonon frequency appears due to the vibration of Zn and O atoms. Therefore, it can be concluded that within the frequency range from 0.5 to 3.0 THz, an unstable phonon mode at the R-and/or M-points of the Brillouin zone is responsible for the phase transition and the superior dielectric and piezoelectric responses of BZO, which correspond to the tilts of the ZrO6 octahedra. Among the perovskite ferroelectric oxides such as PZT [80], this type of unstable phonon mode is observed commonly at the R-or M-points of the Brillouin zone. Finally, we also estimated the IR spectra of BZO using first-principles calculations. According to the group theoretical analysis, due to the cubic symmetry, three phonon branches of the T1u symmetry is IR active, and one branch of the T2u symmetry is not IR, nor Raman active. The estimated IR spectra of the BZO single crystal are displayed in Figure 10 where the frequencies of the three TO modes are 2.4, 5.3, and 14.8 THz, respectively. These values of IR active modes are in good agreement with previous theoretical results [72,76]. However, Nuzhnyy et al. [74] obtained the values of these IR active modes to be 3.48, 6.4, and 15.6 THz, respectively, which are higher than the present calculated values. It is well known that the mode frequencies obtained from DFT calculations may vary due to the choice of different functionals and pseudopotentials, as shown for BZO, as well as for BTO. The origin of the additional modes observed by Nuzhnyy et al. and our THz-TDS has not been discussed anywhere using first-principles calculations. J. W. Bennett et al. estimated the values of IR modes by considering the P1 ̅ symmetry [71]. However, their reported values also cannot explain the experimentally observed additional modes in BZO. Therefore, further theoretical calculations are not only necessary, but will also be interesting.

First-Principles Calculations of BaCuO3 Crystals
Very recently, another important type of perovskite, BCO, was discovered [36]. BCO belongs to the family of chemically simple perovskite oxides. At a fundamental level, to know the dynamical stability of BCO over the whole Brillouin zone boundary, the phonon dispersion relation is very important. As discussed above, the observation of negative phonon frequencies corresponds to the instability of the crystal symmetry; otherwise, they are considered as dynamically stable. The calculated phonon spectra and These values of IR active modes are in good agreement with previous theoretical results [72,76]. However, Nuzhnyy et al. [74] obtained the values of these IR active modes to be 3.48, 6.4, and 15.6 THz, respectively, which are higher than the present calculated values. It is well known that the mode frequencies obtained from DFT calculations may vary due to the choice of different functionals and pseudopotentials, as shown for BZO, as well as for BTO. The origin of the additional modes observed by Nuzhnyy et al. and our THz-TDS has not been discussed anywhere using first-principles calculations. J. W. Bennett et al. estimated the values of IR modes by considering the P1 symmetry [71]. However, their reported values also cannot explain the experimentally observed additional modes in BZO. Therefore, further theoretical calculations are not only necessary, but will also be interesting.

First-Principles Calculations of BaCuO 3 Crystals
Very recently, another important type of perovskite, BCO, was discovered [36]. BCO belongs to the family of chemically simple perovskite oxides. At a fundamental level, to know the dynamical stability of BCO over the whole Brillouin zone boundary, the phonon dispersion relation is very important. As discussed above, the observation of negative phonon frequencies corresponds to the instability of the crystal symmetry; otherwise, they are considered as dynamically stable. The calculated phonon spectra and PPDOS of BCO at absolute zero temperature along the principal symmetry direction of the Brillouin zone are shown in Figures 11 and 12, respectively. As for BZO, negative phonon frequencies are absent in BCO. On the other hand, at the G-point, the positive value of the low-frequency acoustic phonon (indicated by red color) clearly reflects that BCO single crystals are dynamically stable. So far, there have been no available data, neither experimental, nor theoretical, to compare the present results. From the PPDOS, it is seen that the low-frequency region (1)(2)(3)(4)   The elastic properties of a material are essential to study because of its response under different external loads. this also helps predict the numerous applications of a solid in different sectors by exploring the outstanding physical properties, such as elastic stiffness, bulk, modulus, compressibility, Young's modulus, ductility, etc. The inde-  The elastic properties of a material are essential to study because of its response under different external loads. this also helps predict the numerous applications of a solid in different sectors by exploring the outstanding physical properties, such as elastic stiffness, bulk, modulus, compressibility, Young's modulus, ductility, etc. The inde- The elastic properties of a material are essential to study because of its response under different external loads. this also helps predict the numerous applications of a solid in different sectors by exploring the outstanding physical properties, such as elastic stiffness, bulk, modulus, compressibility, Young's modulus, ductility, etc. The independent elastic constants for cubic symmetry, namely C 11 , C 12 , and C 44 , of BCO have been calculated using the stress-strain method [81,82]. All the estimated elastic constants are positive and satisfy the mechanical stability criteria (C 11 − C 12 ) > 0, (C 11 + 2C 12 ) > 0, and C 44 > 0. Therefore, BCO crystals are stable from the viewpoint of both mechanics and thermodynamics.
For the material containing the cubic phase, the equations for the calculation of B and G are where B V and G V represent the Voigt approximation, B R and G R indicate the Reuss approximation, and B and G have been calculated following the Voigt-Reuss-Hill approximation. The hardness of a material strongly depends on C 44 and G. In order to check the hardness of the other well-known perovskites, we compared our results with BTO, BZO and SrZrO 3 , which is shown in Figure 12. From Figure 13, it can be noticed that BCO has a smaller value of not only B and C 44 , but also the other elastic moduli, which should result in BCO being a soft material. Therefore, in comparison with those perovskites, BCO is more suitable for use as a thermal barrier coating substance.
where BV and GV represent the Voigt approximation, BR and GR indicate the Reuss approximation, and B and G have been calculated following the Voigt-Reuss-Hill approximation. The hardness of a material strongly depends on C44 and G. In order to check the hardness of the other well-known perovskites, we compared our results with BTO, BZO and SrZrO3, which is shown in Figure 12. From Figure 13, it can be noticed that BCO has a smaller value of not only B and C44, but also the other elastic moduli, which should result in BCO being a soft material. Therefore, in comparison with those perovskites, BCO is more suitable for use as a thermal barrier coating substance. The most interesting phenomena observed in BCO is that it exhibits superconductivity at low temperatures. Recently, for Ba2CuO4-y, Tc was found to be around 73 K [83]. We also estimated Tc by using the Allen-Dynes modified McMillan equation [ where is the electron-phonon coupling constant, defined as = ( ) − ℎ , with − ℎ denoting the electron-phonon interaction energy responsible for the Fermi surface instability and Cooper pairing, and * is the repulsive Coulomb pseudopotential. The values of , , and * were obtained as 121, 0.84, and 0.01, respectively, using the calculation by the quantum ESPRESSO code. These values were used to calculate Tc using the The most interesting phenomena observed in BCO is that it exhibits superconductivity at low temperatures. Recently, for Ba 2 CuO 4-y , T c was found to be around 73 K [83]. We also estimated T c by using the Allen-Dynes modified McMillan equation [84] as where λ is the electron-phonon coupling constant, defined as λ = N(E F )V e−ph , with V e−ph denoting the electron-phonon interaction energy responsible for the Fermi surface instability and Cooper pairing, and µ * is the repulsive Coulomb pseudopotential. The values of ω log , λ, and µ * were obtained as 121, 0.84, and 0.01, respectively, using the calculation by the quantum ESPRESSO code. These values were used to calculate T c using the above equation, and we obtained T c = 9.96 K. The large difference arose because of the oxygen defect. It has already been mentioned that the oxygen defect can increase the T c of a material [85,86]. However, here, we report the T c for the pure BCO structure. At some selected pressures, the real (ε 1 (ω)) and imaginary (ε 2 (ω)) parts of the dielectric constants are displayed in Figure 14a,b, respectively. The low-energy optical response is dominated by ε 2 (ω) for metallic compounds. From Figure 14a, it is seen that ε 1 (ω) crosses zero from the bottom level at 37.9 eV, and at the same energy, ε 2 (ω) approaches zero, which reflects the metallic nature of BCO. It also suggests that the material becomes transparent when the energy of the incident photon exceeds 37.9 eV. Due to the intraband transition, the prominent peaks are observed in both ε 1 (ω) and ε 2 (ω) at around 27 and 31 eV. The peak at around 27 eV almost diminishes, when the applied pressure is 18 GPa.  To explain these observations, we calculated the pressure dependence of the partial density of electronic states, which is shown in Figure 15. It can be noticed from the pressure dependence of the partial density of states that there is a prominent hybridization between the neighboring Cu-dx 2 -y 2 and O-2px orbitals at 18 GPa, which is supposed to screen the Cu core holes, resulting in the decrease of the peak intensity. To explain these observations, we calculated the pressure dependence of the partial density of electronic states, which is shown in Figure 15. It can be noticed from the pressure dependence of the partial density of states that there is a prominent hybridization between the neighboring Cu-d x 2 −y 2 and O-2p x orbitals at 18 GPa, which is supposed to screen the Cu core holes, resulting in the decrease of the peak intensity.
Significant absorption at zero photon energy was observed in BCO, as shown in Figure 16a, indicating the metallic nature. Many peaks were observed in the absorption spectrum. However, the intensities of these peaks in the IR and visible regions were small compared to the UV region. At around 27 and 31 eV, two intense peaks were observed. As a result, the solar spectrum in UV regions can be absorbed by BCO. The real part of the photoconductivity (σ) of the BCO crystal is displayed in Figure 16b. It is seen that σ has a finite value at zero photon energy, which reflects the metallic conductivity of BCO, totally consistent with the results of the absorption and dielectric constant. The finite value of σ can be increased owing to the absorption of the photon with sufficient energy. The maximum value of σ was obtained in the UV region. A sudden jump and decrease of intensity of a peak around 27 eV can be noticed from Figure 16a,b with increasing pressure, which might be associated with the transition from the initial state 2p 6 3d 9 to the final state 2p 5 3d 10 . Overall, there are no studies on pure BCO to compare the result. Many studies are required for the superconducting BCO to understand the basic mechanism of its superconductivity. To explain these observations, we calculated the pressure dependence of the partial density of electronic states, which is shown in Figure 15. It can be noticed from the pressure dependence of the partial density of states that there is a prominent hybridization between the neighboring Cu-dx 2 -y 2 and O-2px orbitals at 18 GPa, which is supposed to screen the Cu core holes, resulting in the decrease of the peak intensity.   Significant absorption at zero photon energy was observed in BCO, as shown in Figure 16a, indicating the metallic nature. Many peaks were observed in the absorption spectrum. However, the intensities of these peaks in the IR and visible regions were small compared to the UV region. At around 27 and 31 eV, two intense peaks were observed. As a result, the solar spectrum in UV regions can be absorbed by BCO. The real part of the photoconductivity (σ) of the BCO crystal is displayed in Figure 16b. It is seen that σ has a finite value at zero photon energy, which reflects the metallic conductivity of BCO, totally consistent with the results of the absorption and dielectric constant. The finite value of σ can be increased owing to the absorption of the photon with sufficient energy. The maximum value of σ was obtained in the UV region. A sudden jump and decrease of intensity of a peak around 27 eV can be noticed from Figure 16a,b with increasing pressure, which might be associated with the transition from the initial state 2p 6 3d 9 to the final state 2p 5 3d 10 . Overall, there are no studies on pure BCO to compare the result. Many studies are required for the superconducting BCO to understand the basic mechanism of its superconductivity.

Summary
The crystal structure of perovskite oxide has a variety of compositions and component elements, which give rise to a wide range of applications due to its unique properties. Here, three different types of perovskites were reviewed by means of Brillouin scattering and first-principles studies.
In BTO single crystals, a significant softening of the frequency shift and FWHM of the LA phonon is observed in the vicinity of a cubic-tetragonal phase transition temperature. The Burns temperature, TB = 533 K, and intermediate temperature, T* = 503 K,

Summary
The crystal structure of perovskite oxide has a variety of compositions and component elements, which give rise to a wide range of applications due to its unique properties. Here, three different types of perovskites were reviewed by means of Brillouin scattering and first-principles studies.
In BTO single crystals, a significant softening of the frequency shift and FWHM of the LA phonon is observed in the vicinity of a cubic-tetragonal phase transition temperature. The Burns temperature, T B = 533 K, and intermediate temperature, T* = 503 K, were identified, which show very close agreement with recent acoustic emission measurement result. When the sample was cooled down with an electric field, an increase of the frequency shift and a decrease of the FWHM in both the cubic and tetragonal phase were observed due to the reduced density of the twin domain walls and static/dynamic PNRs. A complete 90 • domain switching process was noticed under the application of an electric filed in the ferroelectric tetragonal phase. In the BZO single crystal, no phase transition was observed from 8 K to 1273 K, reflecting that the BZO system is macroscopically cubic. The presence of a TA phonon in the Brillouin spectra suggests local lattice distortion, which destroys the center of symmetry in the cubic matrix. This observation is also justified with the observation of an additional absorption peak near 2 THz in the THz-TDS experiment. The imaginary part of the frequency in the phonon dispersion spectra suggests structural instability with respect to oxygen octahedral rotations. In BCO, the dynamical and mechanical stabilities were confirmed from the calculation of the phonon spectra and elastic constants using first-principles calculations. The magnitude of different elastic constants of BCO compared with other related perovskites shows that BCO is a mechanically stable soft perovskite with improved ductility. The superconducting transition temperature, T c = 9.96 K, was obtained from phonon calculations. The change of the peak position and the reduction of the peak intensity in various optical constants with pressure were in good agreement with the observation of a pressure-dependent partial DOS. The shift of the peak position at 18 GPa in different optical parameters was attributed to the transition from the initial state 2p 6 3d 9 to the final state 2p 5 3d 10 of the Cu atom. Since there are no experimental reports on cubic BCO available, the results discussed here will be very helpful, not only for future experimental and theoretical research, but also for the assessment of the suitability of BCO in different optoelectronic devices.