Resonant Dampers for Parametric Instabilities in Gravitational Wave Detectors

Resonant Dampers for Parametric Instabilities in Gravitational Wave Detectors

S. Gras    P. Fritschel    L. Barsotti    M. Evans Massachusetts Institute of Technology, 185 Albany St. NW22-295, 02139 MA, USA

Advanced gravitational wave interferometric detectors will operate at their design sensitivity with nearly 1MW of laser power stored in the arm cavities. Such large power may lead to the uncontrolled growth of acoustic modes in the test masses due to the transfer of optical energy to the mechanical modes of the arm cavity mirrors. These parametric instabilities have the potential of significantly compromising the detector performance and control. Here we present the design of “acoustic mode dampers” that use the piezoelectric effect to reduce the coupling of optical to mechanical energy. Experimental measurements carried on an Advanced LIGO-like test mass shown a 10-fold reduction in the amplitude of several mechanical modes, thus suggesting that this technique can greatly mitigate the impact of parametric instabilities in advanced detectors.


I Introduction

The network of advanced gravitational wave detectors currently under construction (two LIGO Harry and the LIGO Scientific Collaboration (2010) detectors in the USA, the Advanced Virgo Bitossi (2013) detector in Italy, and Kagra Somiya (2012) in Japan) promises to open the new window of gravitational wave astronomy within this decade.

Sensitivity to gravitational wave strains of order requires high optical power circulating in the arm cavities of these detectors. For instance, up to 750 kW of optical power will be sustained in the steady-state regime inside the Advanced LIGO arm cavities.

It has been experimentally observed that the stored energy in a resonant cavity can leak from the optical modes to the mechanical modes of the cavity mirrors via a 3-mode interaction Zhao et al. (2008). Given sufficiently high circulating optical power, and mirror materials with very low mechanical loss as required to avoid thermal noise, the uncontrolled growth of test mass acoustic modes can lead to Parametric Instabilities (PI) Evans et al. (2010); Gras et al. (2010). If left unaddressed PI will prevent high power operation, and thus limit the astrophysical output of gravitational wave detectors. While Advanced LIGO will serve as the primary example in this paper, all advanced gravitational wave detectors are susceptible to these instabilities.

Several schemes have been proposed to damp PI Zhao et al. (2005); Zhang et al. (2010). In particular, solutions directly applicable to Advanced LIGO involve active damping of acoustic modes by means of the test mass electro-static drive actuators Miller et al. (2011), and thermal tuning of the optical modes using the test mass ring heaters Gras et al. (2009). A significant constraint on any technique is that it must preserve the inherently low mechanical loss of the test mass in the gravitational wave frequency band to maintain a low level of thermal noise.

Here we present a novel method to passively control PI by reducing the Q-factor of the test mass acoustic modes with small resonant dampers. These “acoustic mode dampers” (AMD) dissipate the strain energy of the test mass mode through a resistive element after converting it to electrical energy via the piezoelectric effect (see figure 2).

The resonant nature of the AMD allows it to effectively damp test mass acoustic modes without introducing significant mechanical loss at lower frequencies where thermal noise can limit detector performance. With respect to other proposed solutions, this approach has the advantage of being simple, self-contained, and completely passive. Models indicate that AMDs can provide a broadband reduction in the Q of mechanical modes relevant to PI, and are therefore particularly beneficial in the presence of a large number of unstable modes.

The structure of this paper is as follows. In section II we set the stage by giving a brief overview of parametric instabilities, including equations of particular relevance to evaluating AMD performance. Section III presents a simple 1-dimensional model of the AMD which highlights the principles of AMD operation. This is followed by a description of the detailed finite element model (FEM) used to analyze the AMDs ability to suppress PI when attached to an Advanced LIGO test mass. The FEM predictions are compared with experimental results from a full-scale prototype in section IV. Finally, in section V we discuss an AMD design that will provide Advanced LIGO with protection from instability, without significantly increasing test mass thermal noise.

Ii Parametric Instabilities (PI)

The acousto-optic interactions responsible for parametric instabilities have been extensively studied Braginsky et al. (2001); Zhao et al. (2005); Evans et al. (2010). They consist of a scattering process and radiation pressure operating together in an optical cavity in a closed-loop manner. The graphical representation of this process is shown in Fig.1.

Figure 1: PI described as a positive feedback process. A steady state cavity field inside the interferometer arm cavity is distorted by scattering off a vibrating mirror surface. Two transverse optical sidebands are created. Both sidebands exert force on the test mass via radiation pressure. When the energy dissipation of the acoustic mode and the rate of work done by the radiation pressure is unbalanced, one of the sidebands excites the exponential growth of the acoustic mode amplitude. The dynamic of this process is commonly described in terms of the parametric gain , with in the case of instability.

The e-folding growth time, or “ring-up time”, of an acoustic mode in the presence of a 3-mode interaction is Ju et al. (2006), where and are the Q-factor and angular frequency of the mode, respectively, and is the parametric gain. When the amplitude increases exponentially until a saturation point is reached Vyatchanin (2013). The parametric gain can be defined as


where is the optical power circulating in the arm cavity, is the optical frequency of the light, and and are the cavity linewidths (full width, half maximum) for the fundamental and the higher-order optical mode, respectively. The parameter measures the spatial overlap between the acoustic and the higher-order optical modes; is representative of the interferometer optical configuration and is a function of the 3-mode interaction tuning , where is the frequency difference between the fundamental and higher-order optical modes. For , the parameter reaches its maximum (see Strigin and Vyatchanin (2007) for a more detailed description).

Unstable acoustic modes with parametric gain up to may arise in Advanced LIGO in the 10-90 kHz band Evans et al. (2010). To prevent these instabilities, a damping mechanism must be introduced to reduce the Q-factor of all unstable acoustic modes in this frequency band without introducing excess thermal fluctuation in the detection band of 10 Hz to 1 kHz.

Iii Model of the Acoustic Mode Damper (AMD)

In order to reduce the Q test mass acoustic modes we designed a resonant AMD which can be attached to the test mass and provide dissipation via the piezo-electric effect.

Figure 2: Overview of the Acoustic Mode Damper (AMD). The AMD can be described as a small spring-mass system attached to a larger mass . represents a modal mass of an acoustic test mass mode which vibrates due to an external force (radiation pressure). Such induced vibration causes deformation of the lossy spring of the damper. The non-zero loss angle of the damping spring (piezoelectric material) results in the lag of strain with respect to stress and thus energy dissipation. As a consequence of this dissipation process the amplitude of the acoustic mode vibration is reduced.

In this section we first describe the interaction between the AMD and the test mass with a simple 1-D model, then we present a complete Finite Element Model of the entire system.

iii.1 Simplified 1-D Model

The AMD and test-mass system can be described as a pair of coupled oscillators with a large mass ratio. The AMD mass is attached to the much more massive test mass via piezo electric shear plates, which are modeled as a lossy spring with complex spring constant of magnitude and loss angle .

The test mass acoustic mode for which we would like to estimate the impact of the AMD is simplified in this model to a mass , equal to the modal mass of the acoustic mode, attached to a fixed reference by a lossless spring . The coupled systems is then excited by the radiation pressure force applied to the TM mode, as shown in Fig. 2.

At frequencies near the resonance of the AMD, the lossy spring produced by the piezoelectric material and resistive load will dissipate the energy of the excited acoustic mode, as seen in Fig. 3.

Figure 3: A viscous-like behavior of the loss angle for shunted piezoelectric material. The total loss factor is a combination of the shunting loss and structural loss of the piezoelectric material. The maximum peak position is located at  40 kHz which corresponds to the shunt with 10 k resistor. The larger the resistor the lower is frequency of the peak loss. The peak height is proportional to electromechanical coupling coefficient . Note, the loss angle at low frequency section is mainly dominated by the structural loss angle of piezoelectric material.

For this system of coupled oscillators, the amplitude of the acoustic mode driven by force at angular frequency is


The resulting effective Q-factor is


where , and we assume .

When the acoustic mode resonance is near that of the AMD, , the acoustic mode Q is reduced to . When the acoustic mode resonance is well above the AMD resonance, , and when it is well below the AMD resonance, , assuming .

To suppress PIs, the test mass acoustic mode Q-factors only need to be reduced from to . Using this simple model we can estimate that a AMD with on a test mass, can give for resonances with .

However, this model ignores a number of important factors. One of these is the location of the AMD relative to the nodes and anti-nodes of each test mass acoustic mode. Quantitatively speaking, the modal mass of a given mode should be increased in this model by the ratio of the displacement at the AMD location to that of the mode’s anti-node squared . Thus, an AMD located near a node will have a reduced value of , and will provide little damping.

Other important factors include the multiple coupled degrees of freedom of the AMD and the directional nature of the piezo damping material, both of which are covered in the follow section.

. Component Material Dimenssions Loss angle Resistor Ceramic, surface mount mm 10 k Reaction mass Tungsten mm 4e-5111Assumed loss value. Epoxy conductive TruDuct 2902 25 0.15111Assumed loss value. Piezo shear plate TRS 200HD mm 0.014 Base fused silica + gold coating mm 7.6f 222Loss obtained from Penn et al. (2006). Epoxy non-conductive (EP30, MasterBond) 25 m 0.1 111Assumed loss value.

Table 1: List of the components for acoustic mode damper prototype. The dimensions and loss angle value are used in finite element modeling described in Sec. III.2

iii.2 Finite Element Model (FEM)

The simple 1-D model introduced in the previous section is useful to provide an intuitive understanding of the AMD damping mechanism. However, it is inadequate to represent the details of the interaction between the AMD and the test mass acoustic modes (TMAMs). Each AMD has at least six resonant modes and hundreds of TMAMs are present in the frequency band of interest; a Finite Element Model is required in order to properly reproduce these modes and calculate the mode overlap between them.

A FEM of the Advanced LIGO test mass with two attached AMDs was constructed with the ANSYS program ANS (). The AMD model corresponds to the parameters of our prototype AMD (see Table 1), and the test mass model parameters are reported in Table 6. All dissipation mechanisms in the test mass substrate, coating and bonds were included in the FEM, along with all the losses related to the acoustic mode damper structure (see table 5 for a full list).

A piezoelectric material (PZT), for which energy dissipation can be easily controlled, is ideal for AMD construction. The AMD design modeled here has 2 PZTs sandwiched between a reaction mass and the interferometer test mass (see figure 2). The PZTs operate in shear, and are poled orthogonally to ensure that all but 1 of the 6 lowest frequency AMD modes are damped. See table for a list of modes; the compression mode is not damped by the shear plates.

For the shear resonant damper, the spring constant can be associated with the shear deformation of the piezomaterial


where is the shear stiffness matrix element, is the area, and is the height of the shear plate, respectively.
The loss factor is induced by shunting the shear plate with a resistor. The active stiffness component becomes a complex quantity with a nonzero imaginary stiffness due to the shunt. The magnitude of is strictly related to the impedance of the shear plate-resistor circuit. The loss factor of the shunted shear plate can be defined as


Because the impedance of any PZT is capacitive, the loss factor is frequency dependent. As such, careful selection of the shunting resistor and PZT dimensions can be used to maximize loss in the band of interest, as shown in figure 3. A more detailed discussion of the piezoelectric loss angle can be found in Appendix A.

We validated the FEM of the AMD by computing its principle resonances, and comparing them with direct measurements performed on a prototype AMD. A total of five principle resonances were identified, with three types of modes: two flag, two anti-flag, and one rotation mode (see Table 2). All these modes are characterized by large shear deformation for the double piezo configuration in the AMD. The sixth compression mode was not measured as it does not involve shear of the piezo plate. Table 2 shows the good agreement between the output of the model and measurements on the AMD prototype.

Mode type FEM [kHz] Measured [kHz]
Fy-mode 9.96 9.77
Fx-mode 12.87 12.61
R-mode 23.28 24.13
aFy-mode 38.36 37.39
aFx-mode 50.51 48.86
Table 2: Verification of the finite element model for AMD. Five principle resonances (two flag, two anti-flag and one rotational mode) obtained with the model are compared to the measured values. The x, y suffix corresponds to the shear plate polarization direction. The principle resonances of the AMD were measured with a capacitive bridge circuit where one of the matching capacitors was the AMD prototype. The difference between flag and anti-flag modes corresponds to the location of the rotation axis about which reaction mass rocks. For the flag pole the rotation axis is at the bottom surface of the shear plate whereas for the anti-flag mode the rotation axis is at the height of mass center of the reaction mass.

Harmonic analysis (finite element analysis with an excitation at a fixed frequency) was conducted to estimate the Q-factor of each of the TMAMs. This approach allows us to include frequency dependent variables such as the shear plate stiffness , and its loss angle (see Appendix A).

The modal Q-factor of each resonant mode of the system was calculated as


where is the frequency of the acoustic mode (see figure 4). The loss associated with each structural component is treated separately; is the modulus of the strain energy in the component and is its loss factor (see tables 1 and 6).

The modeled resonant frequencies of 12 TMAMs were then compared to the measured resonance frequencies (see table 3), obtained according to methods described in section IV. The front face displacement of the modes are shown in Fig. 4.

Note that the agreement between calculated and measured Q-values increases with frequency. This may indicate an additional dissipation process which is missing in the model. Nevertheless, a small relative frequency offset below 1% in Table 3 indicates a good agreement of the FEM with measured values.

Mode # Freq. [Hz] Freq. [Hz] f % Q-factor
FEM Measured FEM Measured Ratio
1 8128.3 8150.9 0.3 37M 1.9M 19.5
2 10391.1 10418.1 0.3 63M 14M 4.5
3 12999.1 12984.7 0.1 29M 15M 1.9
4 15101.4 15047.2 0.4 56M 17M 3.3
5 15151.0 15539.1 2.5 55M 16M 3.4
6 19487.0 19544.9 0.3 30M 7.0M 4.3
7 20113.6 20185.5 0.3 27M 13M 2.1
8 24824.0 24901.8 0.3 32M 16M 2.0
9 26504.3 26681.2 0.7 48M 18M 2.7
10 29767.4 29699.5 1.0 18M 15M 1.1
11 30912.1 31003.3 0.3 18M 12M 1.5
12 32664.1 32743.2 0.2 14M 13M 1.1
Table 3: The Q-factor of the test mass modes without AMDs installed. The 4th and 7th column correspond to the relative frequency noise f% and Q-factor ratio between calculated and measured value, respectively.
Figure 4: Test mass drumhead modes for which Q-factor was measured. The color-code corresponds to the test mass front-face displacement amplitude. This figure was obtained from FEM analysis of the test mass model without AMDs.

Iv Experimental Results

Several AMD prototypes were constructed, each consisting of six components (see Table 1):

Reaction Mass: A 10 g tungsten cylinder which tunes the AMD principle resonances to frequencies above 10 kHz is located on top of the shear plates.

Shear Plate: Two piezoelectric shear plates, oriented with perpendicular polarizations are glued to the reaction mass and base with conductive epoxy Tra-Duct2902. The epoxy serves to electrically connect the PZT electrodes to the reaction mass and base.

Base: The interface between flat shear plates and curved barrel of the test mass. The top flat surface is gold coated with separate sections to which shear plates are glued. The bottom surface is curved and match to the test mass barrel curvature. The base is made from fused silica and glued to the test mass with nonconductive epoxy.

Shunting circuit: One 10 k resistor for each shear plate is glued to the reaction mass with conductive epoxy. The circuit is closed with 100 m diameter silver coated copper wire, which is soldered to the resistor and gold coated surface of the base.

To measure Q-factors of the test mass modes, both with and without AMDs attached, a 16 m optical cavity was used, with the test mass forming the end mirror of the cavity. The optical cavity supports a 2 mm diameter Gaussian mode; the resonant beam probes the motion in the center of the test mass face, so that modes with an anti-node in the center of the mass are easily measured. A total of 12 modes were identified (Fig. 4) and measured, with the results shown in Table 3.

In order to measure the impact of the AMD on the Q of the TMAMs, several modes of the test mass were excited using electrostatic actuators Carbone et al. (2012) and their ring-down times observed. This measurement was repeated in three configurations: with no AMD, with 2 AMDs, with 2 AMDs which lacked resistive damping (shunt wires cut).

Each acoustic mode was detected in the cavity locking error signal as a peak in the Fourier domain. After being excited, each mode amplitude was recorded as a function of time to estimate the decay time . The Q-factor was determined from the decay time , according to


where is the resonant frequency in Hz and is the exponential decay time constant.

Figure 5: Two acoustic mode dampers (AMDs) attached to the barrel of a suspended aLIGO test mass. The AMDs are on opposite sides of the test mass, down from the midline.

Measurements of the test mass mode Q-factor performed after installing two AMDs clearly show the substantial damping capability of the AMD prototype, as reported in Table 3 (see also Table 7 in Appendix for additional details).

The results indicate that out of 12 modes, 11 are suppressed by at least factor of 10 and in some cases by more than two orders of magnitude. The relatively large discrepancy between model and measurement for mode #2 can be explained by the AMD and TM interaction condition for this particular mode. The FEM predicts that mode #2 will be within 500Hz of the AMD Fy-mode at 10 kHz, while the measured values give a separation of 800Hz (mostly due to the AMD resonance being off), reducing the interaction between the AMD and TM mode.

In the off-resonance interactions, which are more common and set the lower limit to AMD performance, the discrepancy is generally less than a factor of a few.

Surprisingly, modes #2 and #6 show a counterintuitive behavior; opening the resistive circuit of the AMDs to decreases the TMAM Q-factor. However, since electrical circuit of AMD affects mechanical stiffness of the shear plates it is expected that principle resonances also changed when the circuit is opened. Also, if the AMD resonance is close to the TMAM frequency, equation 4 indicates that the mode Q can be decreased by lowering the AMD loss.

The large Q reduction for modes #3 and #4 is due the on-resonance interaction between TMAM and AMD, whereas for modes #5 and #12 the large Q reduction is due to the anti-node AMD location on the TMAM.

Mode # Freq. [Hz] Damping Resistive
factor contribution [%]
1 8150.9 32.2 12
2 10418.1 31.8 -57
3 12984.7 441.2 21
4 15047.2 81.0 9
5 15539.1 320 71
6 19544.9 15.6 -2
7 20185.5 13.8 41
8 24901.8 5.2 4
9 26681.2 360 0
10 29699.5 23.8 43
11 31003.3 307.7 68
12 32743.2 260 4
Table 4: Test mass mode suppression obtained with two AMDs.

V Analysis of the AMD thermal noise

The AMD is designed to increase the mechanical damping of the test mass acoustic modes above 10 kHz. At the same time, the AMD must introduce minimal additional mechanical loss in the 10-1000 Hz band, where low mechanical loss is required to keep test mass thermal noise small Callen and Greene (1952). Thus it is critical to calculate the thermal noise resulting from the AMDs in our overall evaluation of this PI mitigation technique.

We used the FEM described above to calculate the thermal noise resulting from our experimental test of two AMDs attached to a test mass. The AMD thermal noise was calculated numerically at 100 Hz, the most sensitive part of the detection band, using Levin’s approach Levin (1998). The energy dissipation per cycle was computed using Eqn. 23 for a pressure profile corresponding to the Advanced LIGO geometry (laser beam radius cm, incident on the front face of the TM). Results are shown in Table 5.

The test mass thermal noise level of corresponds to the design level for Advanced LIGO, and is dominated by optical coating loss Harry et al. (2002); Evans et al. (2008). As Table 5 shows, while our prototype AMD would contribute significantly to thermal noise it is not orders of magnitude above the more fundamental sources of thermal noise.

Component Thermal noise @ 100 Hz
Test Mass:
Substrate 0.8
Optical coating 5.1
Suspension ears 0.0
Ears bond 0.6
Total per TM 5.2
Reaction mass (RM) 0.3
Epoxy (RM-PZT) 4.8
Shear plate (PZT-X)111A product of a structural+resistive loss angle. Note that a major contribution to that value comes from the structural loss, see Fig. 3. 5.5
Shear plate (PZT-Y)111A product of a structural+resistive loss angle. Note that a major contribution to that value comes from the structural loss, see Fig. 3. 5.5
Epoxy Base-PZT 7.2
Base 0.0
Epoxy Base-TM 6.0
Total per AMD 13
Table 5: Thermal noise budget for test mass and acoustic mode damper prototype when attached to the test mass. The thermal noise was calculated for the laser beam spot size of = 5.5 cm.

Relative to the prototype device, our model points to several design and material improvements that can be made to significantly reduce the thermal noise impact. The major AMD thermal noise contributors are the epoxies used to bond the AMD elements, and structural loss in the piezo shear plates. The former can be improved with lower loss epoxy and thinner bond layers. The latter can be improved with a more suitable choice of piezo material.

Other design elements can also be altered. The mass of the reaction mass can be reduced to lower the thermal noise without significantly affecting the acoustic mode damping performance. Another modification would be to avoid alignment of the piezo shear plate polarization with the laser beam axis, to minimize the contribution of resistive loss to the thermal noise. These and other design optimizations will be explored in a future paper.

Vi Conclusion

Acoustic mode dampers represent a simple yet effective approach to damping parametric instabilities. The great advantage of this approach over active damping Miller et al. (2011) is that many test mass acoustic modes are effected simultaneously, and no further intervention is necessary. This is likely to be a critical feature in instruments that would otherwise suffer from multiple acoustic modes simultaneously excited by parametric instabilities.

The investigation presented here involved modeling and construction of a prototype AMD, which was shown to effectively damp test mass acoustic modes. The thermal noise associated with this prototype AMD was also computed. Though the prototype AMD does not meet the stringent thermal noise requirements of gravitational wave detectors, several design elements were identified for improvement, making this approach a viable solution to parametric instabilities.

The authors gratefully acknowledge the support of the National Science Foundation and the LIGO Laboratory, operating under cooperative Agreement No. PHY-0757058 This paper has been assigned LIGO Document No. LIGO-P1400257.


Appendix A Piezoelectric material

Any piezoelectric material is strictly characterized by electromechanical properties. The fact that both electrical and mechanical properties are inseparable allows us to represent the dissipation process in a shunted piezoelectric material either as a Jonson heat or mechanical loss. For the purpose of this work we focus on disspation process in terms of the mechanical loss.
From the stress-charge of the piezoelectric constitutive equation it is straightforward to derive the total induced stress in the shunted piezoelectric material Hagood and von Flotow (1991). Assuming no external current plugged to the piezoelement electrode we get


where is the mechanical stiffness matrix under constant electric field, and are the piezoelectric stress constant. The upper script refers to transpose operator 265 (1988). Electrode area is represented by matrix , whereas thickness of the piezomaterial between electrodes by matrix . The strain is a product of the acting stress on piezo-element and charge accumulation in the piezo-element. Each bracket corresponds to the matrix dimenssion. It is assumed that piezoelectric element has a brick shape thus matrices ,, and are diagonal.
The total impedance in Eqn.9 is inverse sum of the piezoelement admittance under constant electric field and the admittance of the external circuit connected to the piezoelement electrodes, thus


The admittance of piezelement is assumed to be exclusively capacitive. Hence,


where , is the capacitance and the dielectric constant matrices under constant strain and constant stress , respectively. It is convenient to operate with since this quantity can be easily measured at stress free conditions and no shortened piezelement electrodes.
According to Eqn.9, if the piezoelement is integrated with a nonzero impedance electric circuit, the imaginary part stiffness tensor arises. The imaginary part can be interpreted as a dissipative component of the stiffness tensor. The total shunted stiffness matrix is


Because matrices , , and are diagonal and knowing that Eqn. 12 can be written as


where is the nondemensional impedance


Note, for zero impedance shunted circuit sets to be a unity matrix whereas for the nonzero impedance shunting is a complex quantity. Using indexing Eqn.12 can be written as


where and corresponds to the electrode position, is the dielectric permittivity of the piezoelectric material.

From Eqn.11 and 15 becomes clear that loss can be control with the shunting circuit. Moreover, such loss is in fact a frequency dependent quantity with arbitrary shape of the loss curve dependent on shunting circuit.

For the shear plate with the single pair of electrodes only a stiffness matrix element is affected by shunting resistor, thus


and the nondemensional complex impedance has a single nonunity element


Since the stiffness matrix is a complex quantity we can define the loss factor as


and thus


where . Index is associated with the electrodes orientation, whereas indices and correspond to stress-strain directions in the stiffness matrix . The shunting loss factor of the shear plate becomes


This is the main loss mechanism based on which shear AMD operates. Note that according to Eqn.16, the material stiffness (the real part of ) is also reduced and should be included in analysis.

a.1 Total mechanical loss angle

An assumption of the stiffness matrix being real in the absence of shunt is not sufficient for accurate estimation of energy dissipation. It becomes especially important in the thermal noise analyisis, see Sec. V.
For known material loss angle of the nonshunted piezoelement, the total loss factor of piezoelectric material can be defined as


where is the material loss factor matrix of piezo-element. The loss factor of can be easily computed using Eqn. 18.

All piezoelectric material has anisotropic structure what implies that strain energy dissipation in such material must depend on the piezoelement geometric shape. It is more convenient to use effective loss angle rather than the loss angle for each stiffness matrix component. We define the energy dissipation per cycle in the piezo-element


where is the stored strain energy such that , where S, T are the complex strain and stress, respectively and V is the volume of the piezo-element.
Since loss factor of the piezoelement is inverse of its Q-factor or equally a ratio of dissipated energy per cycle to energy stored in the piezoelement, we can defined the effective noise fuctor


where , are the complex strain and stress, respectively and is the volume of the piezo-element. This is a key equation in the finite element analysis which properly estimate contribution of piezoelemnt in the strain energy dissipation in the test mass.

In our analysis we assumed constant material intrinsic loss factor for all elements equal to =0.014 111The information about the imaginary part of the stiffness matrix is not available. We assumed that the material loss factor corresponds to the Q-factor of the piezomatrial provided by the manufacturer. In our opinion low Q PZT should have fairly constant loss for all stiffness components, what leads to


Note, the remaining elements of shunting induced loss are equal to the material structural loss factor. For the shear plate configuration is simply a function of electromechanical coupling coefficient and is equal to . The total loss angle for the stiffness element is shown in Fig.3.

Appendix B Additional Tables

Optical parameters:
End test mass transmissivity 5.0 ppm
1 layer each of TaO/SiO with thickness of 2.32 m, and 3.49 m111For the purpose of numerical analysis, a multilayer coating was reduced to double layers with a total thickness which corresponds to the sum of all 18/19 layers of TaO/SiO for the ETM., respectively
Mechanical properties of optical coating and ear bond:
TaO SiO                             Ear bond
Young’s modulus 140 GPa 70 GPa                         7.2 GPa
Poisson ratio 0.23 0.17                               0.17
Density 8300 kg/m 2201 kg/m                   2202 kg/m
Refractive index 2.06539 1.45                               -
Loss angle222f-frequency. Note, the difference in loss angles for the substrate, and for the optical coating made from fused silica. Losses were obtained from Harry et al. (2007). Additionally, coating loss angles were revised to the current measured values. 2.4+f1.8 0.4+f1.1      0.1
Test mass dimensions:
Radius 0.17 m
Thickness 0.2 m
Flats width 0.095 m
Wedge angle 0.07 deg
Mass 40 kg
Loss angle333Both test mass and suspension ears. 7.6f
Material333Both test mass and suspension ears. fused silica

Table 6: Values for the aLIGO end test mass parameters used in this paper.
Mode # Resistive Open Circuit
FEM Measured FEM Measured
1 52k 59k 70k 67k
2 7.9k 440k 9.9k 280k
3 23k 34k 23k 43k
4 420k 210k 510k 230k
5 2.4M 50k 2.2M 170k
6 1.9M 450k 4.8M 440k
7 1.1M 940k 1.6M 1.6M
8 6.7M 3.1M 9.5M 3.2M
9 49k 50k 64k 50k
10 1.9M 630k 3.3M 1.1M
11 61k 39k 116k 120k
12 25k 50k 260k 52k
Table 7: Results of the Q-factor measurement with attached AMDs for shunted and non-shunted cases.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description