Optimization of Dimensions of Cylindrical Piezoceramics as Radio-Clean Low Frequency Acoustic Sensors

Circular piezoelectric transducers with axial polarization are proposed as low frequency acoustic sensors for dark matter bubble chamber detectors. The axial vibration behaviour of the transducer is studied by three different methods: analytical models, FEM simulation, and experimental setup. To optimize disk geometry for this application, the dependence of the vibrational modes in function of the diameter-to-thickness ratio from 0.5 (a tall cylinder) to 20.0 (a thin disk) has been studied. Resonant and antiresonant frequencies for each of the lowest modes are determined and electromechanical coupling coefficients are calculated. From this analysis, due to the requirements of radiopurity and little volume, optimal diameter-to-thickness ratios for good transducer performance are discussed.


Introduction
The first discovery of the piezoelectric effect was found in 1880 by Jacques and Pierre Curie.They discovered that a quartz crystal submitted to mechanical action will produce an electrical potential that is proportional to the force applied.This effect is called direct piezoelectric effect.In the next year, Gabriel Lippmann described an analytical method for the reverse piezoelectric effect and the Curie brothers demonstrated this theory in the laboratory.Since its discovery, the piezoelectricity has been used in different fields such as defence, medical diagnosis, particle detection, ultrasonics, motors, and echolocation [1,2].Due to the different vibrations modes in a piezoelectric material, the applications of this phenomenon can be different in each field.The first piezoelectric material used was the quartz; however, this material is not very good to transform the electrical energy into mechanic energy and its sensitivity is small.Later, in the middle of XIX century, lead zirconate titanate crystals (PZT) were developed.This material has a better sensitivity and frequency response.This material is produced using a very intense electrical field that polarizes the ceramic in a specific direction.
Piezoceramic transducers have usually a regular geometry and the knowledge of the behaviour of piezoceramic disks is important for transducer design and applications [3,4].With these transducers, it is possible to generate vibrations from a few to several hundreds of kHz, demonstrating its feasibility to be used as ultrasonic sensors and actuators in this frequency range.
For example, recent applications use piezoelectric sensors in order to detect the acoustic signal emitted by the interaction of elementary particles in a fluid target.In this sense, the dark matter detectors PICO bubble chambers use a superheated fluid target filled in a glass vessel [5,6].Under some thermodynamic conditions, the particle interaction produces a bubble nucleation within the fluid and acoustic waves are emitted during the bubble growth.These sensors are glued to the external walls of the vessel that contains the metastable fluid.The acoustic discernment of the different acoustic events depends forcefully on the complete sensor properties and there are constraints as well in the little size and radiopurity of the piezoelectric ceramics [7].These detectors use circular piezoceramics (both cylindrical and disk types) in the transducer design for a frequency bandwidth up to 150 kHz.These show the influence of bonding the ceramic to the vessel in its acoustic response.Among others, it is extracted that there is an increase of the sensors sensitivity, mainly in low frequencies.This is due to a better adaptation of acoustic impedances between the medium and the ceramic, through the glass of the vessel [7].However, the first step in the design of the final transducer is the choice of the type and size of the ceramic taking into account its final use.For this, and since low radioactivity is a must, the amount of ceramic used in the transducer is a factor of great importance.This is because the ceramic material contains lead, which is usually accompanied with heavy radioisotopes that are alpha emitters and thus a source of background through alphaneutron reactions.
In this paper, we explain a proven methodology to study the optimization of this type of sensors.For this, several circular section (disks or cylinders) PIC255 piezoceramics with different width and height are studied with analytical and numerical methods and the results obtained are contrasted by experimental measurements.In previous works with circular PZT, some authors studied the natural vibration modes of axial symmetric piezoceramics using Finite Element Methods (FEM) [3] focusing on the efficiency and transducer design.Theoretical analysis contrasted with experimental measurements of resonant vibration with interferometry and laser Doppler vibrometer have been used [8].Additional characterization methods and theoretical approach formulas can also be found.Here, the main parameters of piezoceramics will be compared: the resonance,   , and antiresonance,   , frequencies of the electrical impedance and its amplitude, the product of each frequency resonance and the length (thickness, , or diameter, ) associated with the mode of vibration (  ⋅     ⋅ ), and the piezoelectric coupling factor,  2 =  2  −  2  / 2  .Moreover, the ratio of coupling factors ( 1 / 2 ) of the lowest piezoelectric modes gives us a quantified estimation of the energy distribution in this frequency range.The conclusions of these studies lead us to have design principles to select a specific piezoceramic circular geometry with a radio-clean piezoelectric material that can be used in the next generation of dark matter bubble chamber detectors (PICO 500 L) [9].

Theoretical Analysis
As discussed, the optimization of the volume is based on a free bounded circular ceramic such as that shown in Figure 1 that shows the geometrical scheme of a piezoceramic disk with thickness t and radius .The piezoelectric disk is polarized along the thickness (-axis), and the two opposite plane faces are covered with complete electrodes.

Analytical Method.
The vibrational behaviour of a circular piece of the characteristics shown in Figure 1 can be obtained from the constitutive equations for a piezoelectric material.See [10][11][12] for a comprehensive review in this regard.For the purposes of this text, it is important to note that there are three different modes of vibration of thin piezoceramic disks: tangential, transverse, and radial extensional modes.However, theoretical and experimental analysis of these modes [8] shows that only the resonance frequencies of the radial vibrations modes can be measured in an impedance analysis.For this reason, we only analyse the vibrational characteristics of the extensional modes.
Assuming that there is a harmonic radial vibration in time () with a known frequency (), the radial displacement in the median plane along the radius  can be expressed as The piezoelectric equations lead to the following general solution for the radial component of the vibration: where  1 is the Bessel function of the first kind, for order 1, and the parameters  and  are given by where V  =   12 /  11 is Poisson's ratio,  is the radius,  is the density of the material,  is the electrical potential difference between the electrodes, and    11 ,    12 , and  31 are mechanical and piezoelectric coefficients that can be seen in the Appendix.
The resonance and antiresonance frequencies can be found where the current becomes infinite and becomes zero, respectively.With this, the characteristic equations for both resonance frequencies can be written as where  0 and  1 are the Bessel function of the first kind, for orders 0 and 1, respectively, and  = .Equation ( 4) is the well-known result for two-dimensional analysis of the radial modes [8].From ( 4) and ( 5), the resulting resonant and antiresonant frequencies of the radial extensional mode for a circular piezoelectric piece with free boundary conditions can be expressed as As mentioned, the basis for comparing these analytical results with experimental measurements or numerical simulations of the electrical impedance is that it reaches a local minimum when the sample vibrates at the resonance frequency and a maximum at antiresonance.

Numerical Simulation Method.
To simulate the behaviour of piezoelectric ceramics, the COMSOL Multiphysics software has been used.With this, it has obtained the impedance curves that have been used in this study.In the sizing of tetrahedral meshing elements, it was taken into account that the smaller wavelength (maximum frequency) was discretized in almost ten parts.A mechanical free boundary condition was set in all the contours of the transducer and uniformly electrical open-circuit voltages were set to the electrodes.The samples were polarized in the length direction (Figure 1) and the applied voltage was set to 500 mV.
In Figure 2, an example of the type of mesh used to discretize the solutions is shown.The number of meshing elements with tetrahedral structure was 6698.Since the electrical impedance measurements were made in air, a surface boundary impedance of 450 kRayl was used.
Table 1 shows the thickness () and diameter () of the different cylindrical ceramics that have been simulated in order to study different vibrational modes.The selected dimensions are the typical ones used in PICO detectors.In addition, it has been decided to use a larger set of ceramics of  = 2 mm to study the variation with the diameter and a set of ceramics with  = 10 mm to study the variation with the thickness.
The geometries with a / ratio equal to or more than 5 can be considered as disks [3].The smaller the / ratio, the more the ceramic behaves as a cylinder.With this in mind, the FEM method can be used to study the ceramics with different radius and length in order to obtain the resulting impedance response and the frequencies and amplitudes of each peak.Subsequently, the resulting behaviour will be studied taking into account the volume of each testing piece.

Experimental Setup
In order to contrast theoretical predictions of models with experimental results, we have measured the bandwidth frequency electrical impedance response of several cylindrical PIC 255 piezoceramics with different diameter and thickness (see Table 2).The measurement was done through the resonance method by using a Wayner Kerr Electronics 6500P LF impedance analyser.The samples are held in air only by three metal wires at the bottom and another one at the top to avoid loads and constraints for the vibration.Figure 3 shows the experimental setup to measure a piezoceramic with  = 25 mm and  = 2 mm.
When the ceramic is excited in a resonance frequency of a radial vibration, the value of the impedance reaches a minimum and the admittance is maximum.

4.1.
Width Frequency Band Impedance.Figure 4 shows both numerical and measured results for one of the studied transducers ( = 25 mm and  = 2 mm).The frequency was swept from 100 Hz to 1.5 MHz with an increment step of 100 Hz.The modulus of the electrical impedance obtained is shown.We can distinguish between the low frequency radial modes and the high frequency thickness modes.The local minima and maxima appearing in the impedance curve correspond to resonance and antiresonance frequencies, respectively [13].
For the first three modes of Figure 4 that correspond to radial modes (left square), a good agreement, either in the frequency position of the peaks or in their amplitude, can be seen with a difference smaller than 0.01% in frequency and smaller than 0.05% in amplitude.For the thickness mode (right square), it should be pointed out that there is a large relative deformation and the effects of the losses are larger than those in radial modes.Anyway, we are not interested in very high frequencies, and thus the high frequency thickness modes are not of application in this work.Even so, when the / ratio is smaller and smaller, this mode becomes of lower frequency, and experimental results and the simulations are in better agreement.

Modes
Behaviour.With these experimental, numerical, and analytical methods, the relationship between the resonance frequencies and geometry parameters of the disk have to be studied.All these approaches will be taken into account and the results will be compared and discussed.For instance, Figure 5 shows the relationship between the frequency and diameter  of the transducer for the first four modes when the thickness is  = 2 mm.The solid line represents the analytical solutions, the circle-symbols represent the numerical simulation results, and the dashed line with squares denotes the situation of  = 25, 20, and 10 mm related to the experimental measurements.
Generally, the analytical and numerical results are in good agreement except for higher modes than the third one.In this range, the frequencies are close to the thickness vibration mode frequencies of the piezoelectric cylinder and strong coupling vibration of thickness and radial direction occurs.In addition, in these regions, the analytical solutions may have some errors due to the small diameter-to-thickness (/) ratio; also the analytical method is based on a membrane theory [8].We can see that the resonance frequency increases with decreasing diameter in a nonlinear relationship, especially for small diameters.
Figure 6 shows the relationship between the frequency  and thickness of the transducer  for the first four modes in the calculation for  = 10 mm.The dashed line denotes the situation of  = 2 and 5 mm.Both experimental and simulated values show that frequency decreases with an increase in thickness and it converges to the analytical solution for low thickness.Since the analytical solution is valid only for ideal disks, it is independent of the thickness.For this, the analytical solutions are constant and close to the FEM solutions for larger diameter-to-thickness ratio, that is, for small thickness.

Frequency Constant.
The frequency coefficient  describes the relationship between the geometrical dimensions of a body and the corresponding (series) resonant frequency.When mechanical free piezoelectric ceramic element is excited with an alternating electric field and there is an impedance minimum, the radial resonance frequency matches with the series resonance frequency.The relationship between the radial mode resonance frequency constant,   , and the diameter of the ceramic element, , is expressed by   =   ⋅ .At higher frequencies, there is another impedance minimum, the thickness resonance frequency.Thickness mode frequency constant,   , is related to the thickness of the ceramic element, , by   =   ⋅ . Figure 7 shows the relation   ⋅  to different / ratio obtained both in FEM methods (dashed line) and experimentally (squares).It can be seen that when / > 10 the frequency coefficient decreases significantly.This behaviour offers a geometrical distinction between disks and cylinders.In accordance with [1], the measures of resonant and antiresonant frequencies of a piezoelectric disk should be made better in terms of / ratio than the conventional standards of measurement of piezoelectric coefficients [14,15].The black solid line indicates the frequency constant   in the theoretical model, which supports the above statement.

Thickness Mode Behaviour.
Although high modes present a larger deviation between experimental and predicted results, it is interesting to show the distinction between antiresonant frequency constant in low and high modes.Figure 8 shows the value of   =   ⋅  for the first four low modes (radial modes).As expected, for lower / ratio,   increases.Both numerical and analytical results agree  with experimental values above   = 1.However, there is a clear trend upwards for lower / ratios.Additionally,   for thickness mode is shown.Since there is no thickness variation of these measures, the frequency constant remains quite constant and it is quite different to the low modes pattern.This behaviour shows clearly the difference between radial and thickness modes in piezoelectric disks.With this, we can predict that, for low / ratios, both radial and thickness modes will be at similar frequencies, being possible that the thickness mode frequency is lower than the radial one.

Electromechanical Coupling Coefficients.
The electromechanical coupling coefficient is an important characteristic of piezoelectric ceramics used to measure the energy conversion efficiency.The dynamic electromechanical coupling coefficient, proposed by Mason et al., can be obtained by measuring the resonant and antiresonant frequencies through the following expression [16,17]: This coefficient is of particular interest in designing piezoelectric elements.From Figure 9, we can find the value of  for each of the extensional modes in function of the / ratio.
As mentioned above, the theoretical approach used in this study is based on the assumption that the piezoelectric element is thin so that the theoretical predictions of resonant frequencies are accurate only if the diameter-to-thickness ratio / is large.Figures 7 and 8 show the discrepancy between the analytical model and numerical simulations of the resonant frequencies for extensional vibrations in the performed values of / ranging from 1 to 20.It can be seen that there is a good agreement for / > 6 with discrepancies within 2% for the first three extensional vibration modes.For the specific case of the first mode, the theoretical prediction is excellent, with discrepancy for / = 10 of less than 0.5%.
In general, it can be seen that the best efficiency for lowest radial mode occurs when / < 1 but there is a minimum in the second radial mode that suggests that there is a weak interaction between radial modes [18].For intermediate / ratio, from 1 to 10, a complicated vibrational spectrum appears, increasing the coupling of higher modes.In ratios / > 10, it can be seen that the coupling factors of the lowest mode remain constants and those of higher modes decrease smoothly [19].These conclusions about the coupling coefficient  could be of great importance for the transducer design.

Discussion and Conclusions
The above studies have been done in order to develop tools to optimize the design of piezoelectric sensors for dark matter bubble chamber detectors.These studies started investigating the vibration of radial extensional modes in piezoelectric disks with mechanical free boundary by an analytical model, numerical simulations, and experimental measurements.Results with the diameter-to-thickness / ratios, ranging from 0.1 (cylindrical shape) to 20 (disk shape), are presented.It is shown that the results agree fairly well for the four lower resonant frequencies.In addition, electromechanical coupling  indicates a different interaction between low and high modes for various / values.The studies and results presented here are especially relevant for high frequencies, a complicated vibrational spectrum appears, and the transducer design becomes more complicated.The results show that a simple length expander mode, with minimum coupling to extraneous vibrational modes, can be achieved for / < 1.These conclusions can be seen as well in Figure 10 where the ratio between the electromechanical coupling factor of the first mode ( 1st,mode ) and that of the second mode ( 2nd,mode ) is shown.
From Figures 8 and 10, it is possible to understand the different behaviour according to / ratio.In situations with / < 1, the thickness mode frequency is below or quite close to the radial modes.Then, a high variation on the ratio between the electromechanical coupling coefficients  can happen due to the different nature of the modes observing values of  1 / 2 up to 4.5, which results in a significant contrast in the efficiency between the two lower modes of the ceramic.For geometries 1 < / < 6, the lowest resonance frequencies already correspond only to the radial mode, with relatively high values but with a lower coefficient  1 / 2 (minimum is about 1.4), that is, with a greater interaction between modes, resulting in a more homogeneous response in this frequency range.For the case / > 6, the resonance frequency decreases and the  1 / 2 ratio increases but stabilizing at a value of approximately 2.2.
In the design of piezoelectric ceramics for low background particle detectors, we must take into account the amount of material used in the construction of the sensors.In this sense, an additional parameter to take into account is the mass (or, equivalently, the volume) of the ceramic.Thus, since it depends on the square of the diameter and it is linear with respect to the thickness, we must reinterpret the previous results in terms of the volumetric ratio for each / value.
In order to have a first approximation to this optimization problem, we compared the behaviour of three of the ceramics we have measured due to their different / and volume values.Table 3 summarizes the results seen in this article.
Among these ceramics, for a wideband sensor application, no geometry meets the desired requirements of a small volume, low frequencies, and reduced  1 / 2 ratio.However, considering that the electromechanical factor strongly  influences the acoustic behaviour of a piezoceramic used as acoustic receiver, we would choose the case / = 2 because it presents a more flat acoustic response, with no prominent peaks [7].
Keeping this in mind, we can generalize the study to the set of simulated ceramics.Figure 11 shows the resulting volume of the studied geometries, distinguishing those corresponding to  = 10 mm (red line) and  = 2 mm (blue line).Hence, from radioactivity contamination conditions, we can establish a volume threshold.As an example, a threshold of 0.5 cm 3 was chosen, so that, all ceramics below this threshold could be used matching the radio-clean limitation of the sensor.
Analogous to Figure 10, Figure 12 shows the ratio between electromechanical coupling factors of the lowest modes with respect to the volume of each sample.In this case, the curves of both geometries have been adjusted to the results by least squares, obtaining a good agreement with a relationship  1 / 2 ∼ − −1 .It is observed that, for ceramics with  = 2 mm, we must use volumes smaller than 0.1 cm 3 to have a coupling between the lowest modes,  1 / 2 < 2. On the other hand, for geometries with  = 10 mm the situation is more relaxed and the ceramic measured is below the threshold of 0.5 cm 3 , matching therefore both criteria simultaneously.
Once the volumetric thresholds and coupling between modes have been established, it is necessary that the ceramic responds to the lowest possible frequency.We must be cautious because, in general, when the volume is smaller, so is the coupling, but the fundamental frequency increases.Figure 13 shows the fundamental resonance frequency for each studied sample, observing the geometry that minimizes the frequency response.From these three figures, we conclude that within 1.5 < / < 9 we can find geometries that are below the volumetric threshold (Figure 11).Among them, those that fulfil 1.5 < / < 4.5 satisfy a ratio  1 / 2 < 2 (Figure 12), that is, a smaller coupling between lowest modes.Within the resulting geometries, the one with a lower fundamental frequency (Figure 13) is  = 10 mm and  = 7 mm.We see that this relation is similar to the experimentally measured ceramic with  = 10 mm and  = 5 mm, and therefore, it could also be a good option.
The compromise between low  1 / 2 ratio and low fundamental frequency can also be studied from the measurement of the acoustic response (received voltage response) of each sample [7].However, with the method presented in this text, we can limit the number of geometries that are adapted to the radio-clean requirements of the piezoelectric acoustic sensor

Figure 2 :
Figure 2: Modelling with FEM.(a) Photograph of a sample of circular piezoelectric transducer and (b) 3D meshing model used in numerical simulations.

Figure 3 :
Figure 3: Experimental setup to measure the electrical impedance of the ceramics.

Figure 4 :
Figure 4: Comparison of impedance modulus obtained from FEM simulation and measurements.

Figure 5 :
Figure 5: Dependence of the resonance frequencies with the diameter.

Figure 6 :
Figure 6: Resonant frequencies as a function of the thickness.

Figure 9 :
Figure 9: Behaviour of the electromechanical coupling coefficients K obtained by theory, FEM simulations, and experimental measures for the lower extensional vibration modes.

Figure 10 :
Figure 10: Ratio between the electromechanical coupling coefficients  of the first and second mode.

Figure 11 :
Figure 11: Relationship between volume and / ratio for the studied geometries.

Figure 12 :Figure 13 :
Figure 12: Ratio between the electromechanical coupling coefficients  of the first and second lowest modes as a function of the volume.