17.01.2014 Views

Vibrational resonance in the Morse oscillator - Indian Academy of ...

Vibrational resonance in the Morse oscillator - Indian Academy of ...

Vibrational resonance in the Morse oscillator - Indian Academy of ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

<strong>Vibrational</strong> <strong>resonance</strong> <strong>in</strong> <strong>the</strong> <strong>Morse</strong> <strong>oscillator</strong><br />

2<br />

1<br />

2<br />

1<br />

0<br />

0<br />

250<br />

Figure 4. Three-dimensional plot <strong>of</strong> <strong>the</strong> <strong>the</strong>oretically computed response amplitude<br />

Q as a function <strong>of</strong> β and g.<br />

500<br />

This is shown <strong>in</strong> figure 3a forβ = 1.5 and 2 with ω = 1. Figure 4 presents <strong>the</strong><br />

dependence <strong>of</strong> Q on β and g for ω = 1. Resonance is seen only for β>ω 2 (= 1).<br />

• At <strong>resonance</strong> ωr 2 = ω 2 and hence Q max = 1/dω and it depends only on d and ω.<br />

The response amplitude at <strong>resonance</strong> is <strong>in</strong>dependent <strong>of</strong> <strong>the</strong> parameters β, g and .<br />

• An analytical expression for g VR (at which Q becomes maximum) is difficult to<br />

obta<strong>in</strong> because ωr 2 is a complicated function <strong>of</strong> g. However, g VR can be calculated<br />

from <strong>the</strong> <strong>resonance</strong> curve. It depends on , ω and β and <strong>in</strong>dependent <strong>of</strong> d and f .<br />

Q decreases with <strong>in</strong>crease <strong>in</strong> <strong>the</strong> value <strong>of</strong> d. In figure 5 we plot <strong>the</strong> <strong>the</strong>oretically<br />

predicted g VR and numerically computed g VR vs. <strong>the</strong> parameter β for a few fixed<br />

values <strong>of</strong> ω. g VR <strong>in</strong>creases with <strong>in</strong>crease <strong>in</strong> <strong>the</strong> value <strong>of</strong> β.<br />

• For very large values <strong>of</strong> g, ωr 2 → 0 and Q approaches <strong>the</strong> limit<strong>in</strong>g value Q L given by<br />

1<br />

Q L (g →∞) =<br />

ω √ ω 2 + d . (14)<br />

2<br />

That is, Q does not decay to zero but approaches <strong>the</strong> above limit<strong>in</strong>g value (see figure<br />

3a). This limit<strong>in</strong>g value depends only on <strong>the</strong> parameters ω and d (note that<br />

Q max = 1/dω). The po<strong>in</strong>t is that when ωr<br />

2 → 0, eq. (9) becomes <strong>the</strong> damped free<br />

particle driven by <strong>the</strong> periodic force whose solution is<br />

Y (t) = Q L f cos(ωt + ), = tan −1 (d/ω). (15)<br />

250<br />

150<br />

50<br />

0<br />

1<br />

2<br />

Figure 5. Variation <strong>of</strong> <strong>the</strong>oretically predicted (cont<strong>in</strong>uous curve) and numerically<br />

computed (solid circles) g VR with <strong>the</strong> parameter β for a few fixed values <strong>of</strong> <strong>the</strong><br />

parameter ω. The value <strong>of</strong> is 10ω.<br />

Pramana – J. Phys., Vol. 81, No. 1, July 2013 133

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!