^{1}

^{*}

^{1}

^{2}

^{3}

^{1}

First principle calculation within the Density Functional Theory (DFT) and Density Functional Perturbation Theory (DFPT) using Local Density Approximation as implemented in Quantum ESPRESSO has been significantly used to investigate the structural and Piezoelectric, properties of Perovskite ZrTi(PbO
_{3}
)
_{2}
. From structural properties calculation, the ground state total energy of -
2417.12 eV has been obtained which led to an equilibrium lattice constant of a= 5.620Å for ZrTi(PbO_{3})_{2}. Our obtained optimized atomic positions and atomic effective charge shows that the optimized ZrTi(PbO_{3})_{2} is stable and the Piezoelectric stress tensor is calculated using Berry-phase approach within density functional perturbation theory (DFPT). From our calculation
,
we have obtained the stress tensor elements with values of d_{1,5}
= 6.81, d_{3,1} = 1.69, and d_{3,3}
= 6.18, which is in agreement with the values obtained for tetragonal PbTiO_{3}.

Perovskite oxides have wide range of structural, electrical and mechanical properties, making them vital materials for many nanotechnological applications, which controls production particles about nanometric dimensions. The task of nanotechnology is the development and implementation of methods of research and modeling of nanoparticles [

The benefits of these materials are connected with new properties of materials that are better in nanoscale [

In 1915 to 1940, there was little progress in development of sonar devices. In 1940, the US sonars typically consisted of a magnetostrictive transducer and an array of nickel tubes connected to a 1-foot-diameter steel plate attached back-to-back to a Rochelle salt crystal in a spherical housing. This assembly penetrated the ship hull and was manually rotated to the desired angle. The piezoelectric Rochelle salt crystal had better parameters, but the magnetostrictive unit wasn’t much reliable. With this, Early World War II losses prompted rapid research in the field, pursuing both improvements in magnetostrictive transducer parameters and Rochelle salt reliability. a superior alternative, an Ammonium dehydrogenate phosphate (ADP), was found as a replacement for Rochelle salt; the first application was a replacement of the 24 kHz Rochelle-salt transducers. Within nine months, Rochelle salt was obsolete. The ADP manufacturing facility grew from few dozen personnel in early 1940 to several thousands in 1942 [

One of the earliest applications of ADP crystals was hydrophones for acoustic mines which crystals were specified for low-frequency cutoff at 5 Hz, withstanding mechanical shock for deployment of aircraft from 3000 m (10,000 ft), and the ability to survive neighboring mines explosions. One of the key features of ADP reliability is that it can keep its parameters even over prolonged storage, and it has a zero aging characteristics. Another application was for acoustic homing torpedoes where two pairs of directional hydrophones were mounted on the torpedo nose, in the horizontal and vertical plane; the different signals from the pairs were used to steer the torpedo left-right and up-down. A countermeasure was developed that targeted the submarine discharged as an effervescent chemical, that made the torpedo goes after the noisier fizzy decoy. The countermeasure gives a torpedo with active sonar. Since active sonar reveals the presence and position of the operator, and does not allow exact classification of targets, it is used by fast (planes, helicopters) and by noisy platforms (most surface ships) but rarely by submarines. It is typically activated very briefly at intermittent periods to minimize the risk of detection. Consequently, active sonar is normally considered a backup to passive sonar. In aircraft, active sonar is used in the form of disposable sonbuoys that are dropped in the aircraft’s patrol area or in the vicinity of possible enemy sonar contacts [

Also Passive sonar arrays for submarines can be developed from ADP crystals. The standard US Navy scanning sonar at the end of World War II operated at 18 kHz using an array of ADP crystals that desired longer range and required the use of lower frequencies. The required dimensions were too big for ADP crystals, so in the early 1950s magnetostrictive and barium titanate piezoelectric systems were developed, but these had problems achieving uniform impedance characteristics, and the beam pattern suffered. Barium titanate was then replaced with more stable lead zirconate titanate (PZT), and the frequency was lowered to 5 kHz. The US fleet used this material in the AN/SQS-23 sonar for several decades. The SQS-23 sonar first used magnetostrictive nickel transducers, but these weighed several tons, and nickel was expensive and considered a critical material. After World War II Magnetostrictive transducers were pursued as an alternative to piezoelectric materials. Nickel scroll-wound ring transducers were used for high-power low-frequency operations, with size up to 13 feet (4.0 m) in diameter, probably the largest individual sonar transducers ever. The advantage of metals is their high tensile strength and low input electrical impedance, but they have electrical losses and lower coupling coefficient than PZT, whose tensile strength can be increased [

But it has been observed that most of these Perovskite oxides are complex systems with some degree of disorder, making them challenging to study experimentally and theoretically. However, their complexity gives them their favorable properties, highly accurate modeling which captures the essential features of the disordered structure necessary to explain the behavior of current materials and predict favorable compositions for new materials.

Current synthetic aperture sonar (SAS) systems are exploring device technology with broad bandwidth capabilities for improved system performance. Bandwidth for these systems is not only a function of the transducer but transmit system capabilities, such as available power, maximum pulse length, duty cycle, and waveform. Because these systems are designed for UUVs as small as 9 inches in diameter, transmitted systems must be housed with limited interior volume. This places increased performance constraints on the transducer technology. Bandwidth and power capabilities of the transducer are a function of the electromechanical coupling coefficient, mechanical quality factor, and the drive level limits of the material. At nanoscale dimensions, ferroelectrics exhibit almost all the above mentioned behaviors. High-performance nanoscale ferroelectric perovskites are also required for miniaturization of devices and high-density nonvolatile memory storage. The processes governing the stability of polarization are not entirely understood and are now under intense experimental and theoretical investigation. Recently, a combination of methodological improvements and rise in computer speed has made first-principles calculations a viable tool for understanding these type of complex systems. In particular, the density functional theory (DFT) approach [

As a result of the above stated problems, the main aim of this research work is to use density functional theory (DFT) simulations as implemented in Quantum ESPRESSO code to model and investigate the properties of new Perovskite oxide ZrTi(PbO_{3})_{2} and also to find out if the material can served as potential material for naval sonar applications.

Density functional theory (DFT) as implemented in Quantum-ESPRESSO package was used to carry out the calculations on Perovskite oxide ZrTi(PbO_{3})_{2}. Total energy calculations were performed within the linear density approximation (LDA). The Perdew-Zunger (LDA) exchange-correlation parameterization was used in the calculations. Zr (4s, 4p, 4d), Pb (5d, 6s, 6p), Ti (3s, 3p, 3d, 4s), and (O 2s and 2p) orbitals were treated as valence orbitals. In all case, the atomic core states were calculated relativistically, spin-orbit coupling was ignored, while the valence states were treated semi-relativistically. We have used a special point’s for Brillouin-zone samplings with a 12 × 12 × 12 k-point mesh. While for the Berry’s phase calculations, we have used the k-space integrations of 8 × 8 × 10 k-point mesh. The results of the calculations were checked for convergence with respect to the number of k points and the plane wave cutoff energy.

In order to obtain the polarization of our compound Perovskite oxide ZrTi(PbO_{3})_{2}, we first carried out self-consistent total-energy calculations as mentioned above. Then, the polarization was computed through the Berry-phase approach within DFPT. Different strains were applied in the z direction to compute d_{33}, and different biaxial strains were also applied simultaneously in the x and y directions to compute d_{31} and d_{15} respectively.

_{3})_{2} while _{3})_{2}. Since it is known that the energetic of the system supports the most stable crystal parameter with the lowest energy It is confirmed that once the minimum energy and the corresponding most stable parameters is obtained one can also have the most stable structure as shown in _{3})_{2}.

The analysis of the convergences with respect to plane-wave energy cutoff, and with respect to k-point sampling grids density are the key starting points of any plane-wave self-consisted field (PWSCF) calculations, this is because, a poor convergence test results always resulted in obtaining incorrect converged total energy. _{3})_{2}, while

However, in

Atom | x | y | z |
---|---|---|---|

Zr | 0.000000000 | 0.000000000 | 0.000000000 |

Ti | 0.500000000 | 0.500000000 | 0.500000000 |

Pb | 0.729799147 | 0.829605131 | 0.729799147 |

Pb | 0.270200853 | 0.170394869 | 0.270200853 |

O | 0.756926987 | 0.249398040 | 0.339505879 |

O | 0.339505879 | 0.249398040 | 0.756926987 |

O | 0.839916652 | 0.421491919 | 0.839916652 |

O | 0.243073013 | 0.750601960 | 0.660494121 |

O | 0.660494121 | 0.750601960 | 0.243073013 |

S/N | E_{CUT} (Ry) | E_{TOT} (Ry) |
---|---|---|

1 | 12 | −2391.11794851 |

2 | 16 | −2405.62513615 |

3 | 20 | −2412.08433421 |

4 | 24 | −2415.06686281 |

5 | 28 | −2416.41157927 |

6 | 32 | −2416.43214121 |

7 | 36 | −2417.28052871 |

S/N | K-point | E_{TOT} (Ry) |
---|---|---|

1 | 2× 2×2 | −2391.11794851 |

2 | 4× 4×4 | −2405.62513615 |

3 | 6× 6×6 | −2412.08433421 |

4 | 8× 8×8 | −2415.06686281 |

5 | 10× 10×10 | −2416.41157927 |

6 | 12× 12×12 | −2416.43214121 |

used in most DFT calculations is used in this research [_{3})_{2}, since the supercell approach used have the same length along each lattice vectors, which is also assume to be the same with the reciprocal lattice vector, and therefore this is the reason why 8 × 8 × 8 k-points are used in each direction. However, by considering the idea of numerical integration, it is expected that using the higher k-points for example, 12 × 12 × 12 k-points will give a more accurate result than a calculation using 8 × 8 × 8 k-points but it will be time consuming.

The best way of choosing specific k-points points to be use in practice is by computing the calculated total energies together with the selected number of k-points on convergence graph. In this situation, it can be seen that with the increasing number of k points, for instance, from 2 × 2 × 2 k points to 8 × 8 × 8 k points, the total energy of the system is (almost) continuously dependent on the number of k-points. But, when the k-points are increased from 8 × 8 × 8 k-points to 12 × 12 × 12 k-points, the total energy is seen to be independent of the number of k-points.

Also as an initial step of understanding the structural properties of any material, the total energy was calculated as a function of the lattice constant for Perovskite ZrTi(PbO_{3})_{2}. The energy versus lattice constant data is plotted in

To extract the value of _{0} from any DFT calculations, it is important to think about the functional form of E_{tot}(a). One of the simplest methods is to express the total energy using a truncated Taylor expansion:

E T ( a ) ≈ E T ( a 0 ) + α ( a , a 0 ) + β ( a , a 0 ) (1)

with α = d E / d a | a 0 and β = 1 2 d 2 E / d a 2 | a 0 .

By definition, if a 0 is the lattice parameter corresponding to the minimum energy, that can fit our numerical data to

E T ( a ) ≅ E 0 + β ( a , a 0 ) 2 (2)

where E 0 , β , a 0 are treated as fitting parameters. Even though we have treated β as a fitting parameter in equation 1, the graph of _{0} = 5.6 Å, this indicate that the value is greater than the equilibrium lattice parameter of cubic Perovskite PbTiO_{3} [_{3})_{2} or using LDA functional in our calculations as GGA-PBE functional adopted in their calculations has proved to always overestimate the lattice parameters [

Therefore, for our Perovskite ZrTi(PbO_{3})_{2} one can clearly see that the bulk modulus at the fully relaxed state is comparable with the other results of similar materials obtained by [

_{3})_{2}. Since it is known that the energetic of the system supports the most stable crystal parameter with the lowest energy.

Using the fully-relaxed structures, the piezoelectric tensor is determined by the density-functional perturbation theory (DFPT). In the present study, we calculate the piezoelectric strain constants which connects the induced polarization P_{i} and strain tensor element as

d i j k = ∂ P i ∂ ε j k (5)

where we computed P i using Berry phase method as a function of strain as

Piezoelectric materials | Lattice constants (Å) | Bulk modulus (GPa) | ||||
---|---|---|---|---|---|---|

This work | Theory | Expt. | This work | Theory | Expt. | |

ZrTi(PbO_{3})_{2} ( a = b = c ). | 5.620 | 5.544 | --- | 189.56 | 170.55 | -- |

BaTiO_{3} ( a = b = c ) | --- | 3.977 | 4.000 | - | 178.65 | 195.00 |

PbTiO_{3} ( a = b = c ) | --- | 3.966 | 3.930 | - | 173.21 | 195.00 |

Materials | Method | d 1 , 5 | d 3 , 1 | d 33 |
---|---|---|---|---|

ZrTi(PbO_{3})_{2} | DFT-LDA | 6.81 | 1.69 | 6.18 |

PbTiO_{3} | DFT-LDA | 6.73 | 1.68 | 6.14 |

PbTiO_{3} | GGA-PBEsol | 3.68 | 1.49 | 3.94 |

implemented in QE code, Wu and Cohen applied the same method for piezoelectricity of PbTiO_{3}. For our Perovskite ZrTi(PbO_{3})_{2}, we have noticed that the polarization is along (001) axis, which is in agreement with the results obtained in case of PbTiO_{3} therefore to describe the piezoelectric effects in our compound we computed three piezoelectric constants ( d 3 , 1 = d 3 , 2 ) , d 3 , 3 and ( d 1 , 5 = d 2 , 4 ) .

Additionally, ( d 3 , 1 = d 3 , 2 ) and d 3 , 3 describe the zero field polarization induced along the z-axis when crystal is uniformly strained in the basal xy-plane or along the z-axis, respectively. _{3} for comparison [

As shown in _{3}. However, there is consistently larger difference with those of GGA as expected.

We have successfully used First principle theory to contribute significantly to our understanding and systematically presented a theoretical study of the structural, piezoelectric properties of the Perovskite ZrTi(PbO_{3})_{2} materials based on Density Functional Theory (DFT) and Density Functional Perturbation Theory (DFPT) within Local Density Approximation.

The ground state total energy of −2417.12 eV is obtained which gives the stable structure of the Perovskite ZrTi(PbO_{3})_{2} material with an equilibrium lattice constant of a = 5.620Å. The optimized atomic positions of Perovskite ZrTi(PbO_{3})_{2} shows that our optimization is within the linear density approximation (LDA). Also our calculated Piezoelectric stress tensor elements of ferroelectric ZrTi(PbO_{3})_{2}, gives d 1 , 5 = 6.81, d 3 , 1 = 1.69, and d 3 , 3 = 6.18, which is in agreement with the values obtained in tetragonal PbTiO_{3}.

The authors declare no conflicts of interest regarding the publication of this paper.

Adeyemi, O.J., Mawashi, G.M., Shuaibu, A., Akinade, B.J. and Ushiekpan, U.R. (2019) First Principle Study of the Structural and Piezoelectric Properties of Perovskite ZrTi(PbO_{3})_{2} Materials for Naval SONAR. Journal of Applied Mathematics and Physics, 7, 271-280. https://doi.org/10.4236/jamp.2019.72022