Oscil·lador harmònic

De Viquipèdia
Dreceres ràpides: navegació, cerca

Es diu que un sistema qualsevol, mecànic, elèctric, pneumàtic, etc. és un 'oscil·lador harmònic' si quan es deixa en llibertat, fora de la seva posició d'equilibri, torna cap a ella descrivint oscil·lacions sinusoïdals, o sinusoïdals esmorteïdes entorn d'aquesta posició estable.

La massa penjada del ressort forma un oscil·lador harmònic.

L'exemple és el d'una massa penjada a un ressort. Quan s'allunya la massa de la seva posició de repòs, el ressort exerceix sobre la massa una força que és proporcional al desequilibri (distància a la posició de repòs) i que està dirigida cap a la posició de equilibri. Si es deixa anar la massa, la força del ressort accelera la massa cap a la posició d'equilibri. A mesura que la massa s'apropa a la posició d'equilibri i que augmenta la seva velocitat, la energia potencial elàstica del ressort es transforma en energia cinètica de la massa. Quan la massa arriba a la seva posició d'equilibri, la força serà zero, però com la massa està en moviment, continuarà i passarà de l'altre costat. La força s'inverteix i comença a frenar la massa. L'energia cinètica de la massa va transformant ara en energia potencial del ressort fins que la massa s'atura. Llavors aquest procés torna a produir-se en direcció oposada completant una oscil·lació.

Si tota l'energia cinètica es transformés en energia potencial i viceversa, l'oscil·lació seguiria eternament amb la mateixa amplitud. En la realitat, sempre hi ha una part de l'energia que es transforma en una altra forma, a causa de la viscositat de l'aire o perquè el ressort no és perfectament elàstic. Així doncs, l'amplitud del moviment disminuirà més o menys lentament amb el pas del temps. Es començarà tractant el cas ideal, en el qual no hi ha pèrdues. S'analitzarà el cas unidimensional d'un únic oscil·lador (per la situació amb diversos oscil·ladors, vegeu moviment harmònic complex).

Casos[modifica | modifica el codi]

Oscil·lador harmònic sense pèrdues[modifica | modifica el codi]

Article principal: Moviment harmònic simple

S'anomenarà \scriptstyle {m} a la massa i \scriptstyle {i} a la distància entre la posició de la massa i la posició d'equilibri. Se suposarà que la força del ressort és estrictament proporcional al desequilibri: \scriptstyle {F =-ky} (llei de Hooke). \scriptstyle {F} és la força i \scriptstyle {k} la constant elàstica del ressort. El signe negatiu indica que quan \scriptstyle {y} és positiva la força està dirigida cap a les \scriptstyle {y} negatives. La segona llei de Newton ens diu:


F=ma=m{dv\over dt}=m{d^2y\over dt^2}

substituint la força obtenim:


 m{d^2y\over dt^2}= -ky

La solució d'aquesta equació diferencial ordinària és immediata: les úniques funcions reals (no complexos) la segona derivada és la mateixa funció amb el signe invertit són sinus i cosinus. Les dues funcions corresponen al mateix moviment. Escollim arbitràriament "cosinus". La solució s'escriu:

La corba de dalt dóna la posició de l'oscil·lador en funció del temps. La del mig dóna la velocitat. A continuació hi ha les corbes de les energies. En blau està l'energia cinètica \scriptstyle{{1\over 2}mv^2} i en vermell l'energia potencial del ressort \scriptstyle{{1\over 2}ky^2}
y = A\cos(\omega t + \phi)\,
  • y\, és l'elongació o diferència respecte a l'estat d'equilibri, les seves unitats són les de \scriptstyle{A}.
  • \scriptstyle{A} és l'amplitud, màxima diferència respecte a la posició d'equilibri.
  • \scriptstyle{\omega=2\pi f} és la pulsació (o freqüència angular) i \scriptstyle{f} la freqüència.
  • \scriptstyle{t} és el temps.
  • \scriptstyle{\phi} és la fase inicial (per a \scriptstyle{t= 0}).

És fàcil comprovar que el valor d' \scriptstyle{\omega} és:

\omega=\sqrt{{k\over m}}

El període d'oscil·lació és:

T=2\pi\sqrt{{m\over k}}

Com ja hem dit, durant un quart d'una oscil·lació l'energia potencial es transforma en energia cinètica. Durant un altre quart, l'energia cinètica es transforma en energia potencial. A la figura de la dreta s'ha traçat la posició en funció del temps (corba de dalt), la velocitat en funció del temps (al mig) i les energies potencials i cinètiques (a sota).

Oscil·lador harmònic esmorteït[modifica | modifica el codi]

Oscil·lador harmònic amb amortidor. La força viscosa és proporcional a la velocitat.

Afegint pèrdues d'energia, s'aconsegueix modelar una situació més propera a la realitat. Així, cal notar que l'oscil·lació descrita en l'apartat anterior es perllongaria indefinidament en el temps (la sinusoide que descriu la posició no convergeix a zero en cap moment). Una situació més versemblant es correspon amb la presència d'una força addicional que frena el moviment. Aquesta força pot ser constant (però sempre amb signe tal que freni el moviment). És el cas de fregaments secs: la força no depèn ni de la velocitat ni de la posició. Una altra situació que es produeix en la realitat és que la força sigui proporcional a la velocitat elevada a una potència, sencera o no. Així succeeix quan la força que frena prové de la viscositat o de les pèrdues aerodinàmiques. Es tractarà únicament el cas més simple, és a dir, quan la força sigui proporcional a la velocitat. En aquest cas la força serà:

F_f=-bv=-b{dy\over dt}

On \scriptstyle{b} és un coeficient que mesura l'amortiment causa de la viscositat. Si \scriptstyle{b} és petit, el sistema està poc esmorteït. Cal notar el signe negatiu que indica, com abans, que si la velocitat és positiva, la força té la direcció oposada a la velocitat. Amb aquest terme complementari l'equació diferencial del sistema és:

 m{d^2y\over dt^2}=-ky-b{dy\over dt}

Es tracta d'una equació diferencial ordinària, lineal, de segon ordre[1] (conté derivades segones) i homogènia (no hi ha terme independent de y). Té tres tipus de solucions segons el valor de \scriptstyle{b^2-4km}:

  • Si \scriptstyle{b^2-4km > 0} el sistema està sobreesmorteït (esmorteïment fort o supercrític)
  • Si \scriptstyle{b^2-4km = 0} el sistema té esmorteïment crític.
  • Si \scriptstyle{b^2-4km < 0} el sistema oscil·la amb amplitud decreixent (esmorteïment dèbil o subcrític)

Oscil·lador sobreesmorteït[modifica | modifica el codi]

Posició en funció del temps d’un oscil·lador harmònic esmorteït.
corba blava: esmorteïment crític.
corba vermella: esmorteïment doble que el crític.
corba verda: esmorteïment igual a 90% de l'esmorteïment crític.

En aquest cas el sistema no és realment un oscil·lador, ja que no oscil·la. La solució és de la forma:

y = A_1e^{\lambda_1 t}+ A_2e^{\lambda_2 t}

on els coeficients de les exponencials són menors que zero i reals (pel que no hi ha oscil·lació):

\lambda_1={-b -\sqrt{b^2-4km} \over 2 m}

y

\lambda_2={-b +\sqrt{b^2-4km} \over 2 m}

\scriptstyle{A_1} y \scriptstyle{A_2} depenen de les condicions inicials (és a dir, de la situació del sistema per a \scriptstyle{t = 0}). La posició no és oscil·lant i tendeix cap a la posició d'equilibri de manera asimptòtica. Les dues exponencials decreixents de les solucions tenen constants de temps diferents. Una és petita\scriptstyle{1/\lambda_1} i correspon a la ràpida cancel·lació de l'efecte de la velocitat inicial. La segona \scriptstyle{1/\lambda_2} és més gran i descriu la lenta tendència cap a la posició d'equilibri.

Oscil·lador amb esmorteïment crític[modifica | modifica el codi]

Aquest cas és el límit entre un sistema oscil·lant i un no oscil·lant. Passa quan

b^2=4km\,

La solució única és:

y = A_1e^{-{b\over 2m}t}+ A_2te^{-{b\over 2m}t}

com abans, \scriptstyle{A_1} i \scriptstyle{A_2} són constants que depenen de les condicions inicials.

L'esmorteïment crític correspon a la tendència més ràpida cap a la situació d'equilibri quan no sobrepassa aquesta posició. Si es disminueix una mica l'esmorteïment el sistema s'apropa més ràpidament a la posició d'equilibri, però sobrepassant la posició oscil·la entorn d'aquest punt (prenent valors positius i negatius).

Oscil·lador amb esmorteïment dèbil[modifica | modifica el codi]

Oscil·lacions esmorteïdes. L'amplitud de la sinusoide està controlada per l'exponencial.

En aquest cas, que és més interessant, tenim un oscil·lador que oscil·la al voltant de la posició d'equilibri amb amplitud decreixent. Passa quan:

b^2<4km\,

La solució és:

y = Ae^{-{b\over 2m}t}\cos(\omega t + \phi)

com abans, \scriptstyle{A} i \scriptstyle{\phi} són constants que depenen de les condicions inicials. La pulsació és:

\omega= \sqrt{{k\over m}-\left({b\over 2m}\right)^2}

La pulsació del sistema esmorteït és una mica menor que la pulsació del sistema no esmorteït \scriptstyle{\omega_\circ= \sqrt{{k\over m}}} perquè la força que l'esmorteeix, frena la massa i la retarda.

L'oscil·lació del sistema està descrita per una sinusoide de freqüència  \scriptstyle {f = {1 \ over 2 \ pi} \ sqrt {{k \ over m} - \ left ({b \ over 2m} \ right) ^ 2}} l'amplitud està multiplicada per una exponencial decreixent la constant de temps és \scriptstyle{\tau={2m\over b}}.

Factor de qualitat Q[modifica | modifica el codi]

En un sistema poc esmorteït és interessant de definir el factor de qualitat (Quality factor en anglès) o simplement Q com a:

Q={\sqrt{km}\over b}

aquesta quantitat és igual a  \scriptstyle {2 \ pi} vegades l'invers de les pèrdues relatives d'energia per període. Així, un sistema que perd 1% d'energia a cada cicle, tindrà un 'Q' de 628. Més interessant, Q és també  \textstyle {\ pi} vegades el nombre d'oscil·lacions que el sistema fa mentre la seva amplitud es divideix per un factor  \textstyle {e}. Si es pot acceptar una aproximació més grollera, Q és 3 vegades el nombre d'oscil·lacions que un sistema fa mentre la seva amplitud cau a 1/3 de l'amplitud inicial. Com a exemples, el Q d'un vehicle amb els amortidors en bon estat és una mica més gran que 1. El Q d'una corda de guitarra és de diversos milers. El Q dels cristalls de quars utilitzats en electrònica com a referència de freqüència és l'ordre d'1 milió. Una copa de vidre ordinari té un Q molt més petit que una copa de vidre de plom (vidre).

Oscil·lacions forçades[modifica | modifica el codi]

Podem iniciar el moviment un oscil·lador harmònic desplaçant de la seva posició d'equilibri i abandonant a la seva oscil·lació lliure (vegeu paràgrafs precedents).

Alternativament, podem aplicar-li una força la intensitat variï de manera sinusoïdal amb el temps. En aquesta situació, l'equació diferencial lineal és inhomogènia. La solució a aquest tipus d'equació està formada per dos termes: la solució general del sistema homogeni més una solució particular del cas homogeni.[2] Per tant, la solució està formada per dues parts, una part transitòria (que s'anul·la passat cert temps), similar a les que vàrem veure en els paràgrafs precedents, més una part estacionària. La solució de la part transitòria és la mateixa la que ja hem vist (equació homogènia). Les úniques diferències són les condicions inicials i finals, que no són idèntiques. Anem a interessar-nos a la solució estacionària. En l'equació diferencial del sistema cal afegir la força sinusoïdal:

 m{d^2y\over dt^2} + b{dy\over dt}+ky=F_m\cos(\omega t)

Per resoldre aquesta equació és més interessant utilitzar el mateix mètode que en electricitat i electrònica. Per a això, s'afegeix a la força real una força imaginària  \scriptstyle {jF_m \ text {sen} (\ omega t)}. Com en electrònica, s'utilitza  \scriptstyle {j = \ sqrt {-1}} en lloc de i. Ara l'equació a resoldre és:

 m{d^2y\over dt^2} + b{dy\over dt}+ky=F_me^{j\omega t}

Però per descomptat, com en electricitat, només la part real de y serà d'interès. La solució és immediata:

 y = Ae ^ {j \ omega t} \,

Si es deriva aquesta expressió i se substitueix en l'equació diferencial, es troba el valor de A:

A={F_m \over k-m\omega^2+jb\omega}

Però A es pot escriure com a \scriptstyle{A=\rho e^{j\phi}} i la solució de \scriptstyle{y} complexa és:

y=\rho e^{j\phi} e^{j\omega t}=\rho e^{j\left(\omega t+\phi\right)}

El valor de \scriptstyle{y} real és la part real de l'expressió precedent:

y=\rho \cos\left(\omega t+\phi\right)

on \scriptstyle{\rho} és el mòdul de \scriptstyle{A} i \scriptstyle{\phi} el seu argument:

\rho=\mbox{Re}(A)=\mbox{Re}\left( {F_m \over k-m\omega^2+jb\omega} \right) =
{F_m (k-m\omega^2)\over (k-m\omega^2)^2 +b^2m^{-2}\omega^2}
\phi=\arg\left(A\right)=\arg\left( {F_m \over k-m\omega^2+jb\omega} \right)

Com en electricitat, l'angle  \scriptstyle {\ phi} dóna el desfasament del moviment pel que fa a la força externa. Si  \scriptstyle {\phi} és positiu, el moviment està en avanç de fase i si  \scriptstyle {\phi} és negatiu el moviment està en retard de fase. En aquest cas el desfasament serà sempre negatiu.

Resposta en freqüència[modifica | modifica el codi]

L'amplitud de les oscil·lacions forçades dependrà, per descomptat, de l'amplitud de la força externa. Però per a una mateixa amplitud de la força, l'amplitud de l'oscil·lació dependrà també de la freqüència. Vegem com varia l'amplitud \scriptstyle {\rho} amb  \scriptstyle {\omega}. Utilitzant la definició de freqüència pròpia del sistema (sense esmorteïment ni força externa):

Resposta en freqüència d'un oscil·lador harmònic. A la freqüència de resonància, l'amplitud és Q vegades més gran que a molta baixa freqüència.
\omega_\circ=\sqrt{{k\over m}}

es pot escriure-hi:

A={F_m\over k}{1\over 1-\left({\omega\over\omega_o}\right)^2 +j{\omega\over\omega_o}\sqrt{{b^2\over km}}}

Si a més s'utilitza la definició de \scriptstyle{Q={\sqrt{km}\over b}}, s'obté:

A={F_m\over k}{1\over 1-\left({\omega\over\omega_o}\right)^2+j{\omega\over\omega_o}{1\over Q}}

En el dibuix de dreta s'ha representat l'amplitud de l'oscil·lació forçada en funció de la freqüència per diversos valors del factor de qualitat Q. A molt baixa freqüència l'amplitud és la mateixa que si la força fos estàtica  \ scriptstyle {F_m = kA}, i el sistema oscil·larà entre les posicions  \scriptstyle {{F_m \ over k}} i  \scriptstyle {- {F_m \ over k}}. Quan la freqüència augmenta, l'amplitud també, aconseguint un màxim quan la freqüència d'excitació és igual a la freqüència pròpia del sistema. A aquesta freqüència pròpia també se l'anomena freqüència de ressonància. També es diu que un sistema excitat a una freqüència propera a la freqüència de ressonància "ressona" o "entra en ressonància". A la freqüència de ressonància, l'amplitud de les oscil·lacions serà Q vegades més gran que la que s'obté en baixa freqüència.

L'ample del pic de ressonància a mitjana alçada, és a dir quan l'amplitud és igual a la meitat del màxim, és igual a la freqüència de ressonància dividida per Q. Aquest ample també es diu banda passant.

Oscil·lador forçat i caos[modifica | modifica el codi]

L'oscil·lador harmònic no pertorbat en una dimensió és un exemple de sistema integrable, amb comportament regular. Tanmateix, l'oscil·lador harmònic pertorbat pot presentar un comportament caòtic[3] caracteritzat per un atractor estrany. Per exemple en el cas d'una pertorbació de tipus x^3\; l'equació de moviment és:


\frac{d^2x}{dt^2} + a\frac{dx}{dt} + \omega^2 x + \varepsilon\omega^2 x^3 = b \cos(\omega t)

Aquest sistema és no integrable i el moviment tendeix ràpidament cap a l'anomenat atractor de Duffing.[4]

Oscil·lador de van der Pol[modifica | modifica el codi]

Comportament caòtic en l'oscil·lador de van der Pol amb excitació sinusoïdal. μ = 8.53, mentre que l’excitació externa té amplitud A = 1.2 i freqüència angular ω = 2π / 10.

L’oscil·lador de van der Pol és un cas especial d’oscil·lador amb esmorteïment no lineal, que respon a l'equació:

{d^2x \over dt^2}-\mu(1-x^2){dx \over dt}+x= 0

Va ser descrit per primera vegada el 1935 per Balthasar van der Pol[5] i presenta comportament caòtic.

Oscil·lador harmònic torsional[modifica | modifica el codi]

Article principal: Ressort de torsió

Importància en Física[modifica | modifica el codi]

Considerem el cas d'un cos sotmès a una força unidimensional:  F (y) . Desenvolupant aquesta força en sèrie de Taylor al voltant del punt d'equilibri (y =0):

F(y)=F(y=0)+y\left ( \frac{dF}{dy} \right )_{y=0}+{1\over 2}y^2\left ( \frac{d^2F}{dy^2} \right )_{y=0}+...

Com l'origen és el punt d'equilibri, el primer terme del desenvolupament és nul. Si les oscil·lacions al voltant de  y = 0 són prou petites, un es pot quedar amb l'aproximació lineal i menysprear els termes d'ordre superior:

F(y)\simeq \ y\left ( \frac{dF}{dy} \right )_{y=0}

Anomenant-lo k a la derivada de la força, s'obté de nou la força recuperadora de Hooke. Aquí rau la importància del oscil·lador harmònic: suposa una primera aproximació per a l'estudi d'un sistema quan es produeixen petites oscil·lacions entorn de la seva posició (o estat) d'equilibri.[6]

Exemples[modifica | modifica el codi]

Circuit LC[modifica | modifica el codi]

Circuit LC sense pèrdues[modifica | modifica el codi]

Circuit LC sense pèrdues.

En la figura de la dreta s'ha dibuixat un circuit oscil·lant LC ideal, és a dir sense pèrdues.

Suposem que, en la situació inicial, el condensador està carregat a una tensió V i que en aquest moment es connecta la inductància. La tensió present en les extremitats de la inductància va a fer aparèixer una corrent de sentit invers a la de la fletxa del dibuix, que augmentarà amb el temps. A mesura que el condensador subministra corrent a la inductància, es descarrega i la tensió disminueix. La disminució de la tensió fa que el corrent augmenti menys ràpidament. La situació continua així, amb la tensió del condensador que disminueix cada vegada més ràpidament (perquè el corrent augmenta) i el corrent que augmenta més lentament (perquè la tensió disminueix). Arriba un moment en el qual el condensador està completament descarregat i el corrent ha arribat a un màxim. Ara el corrent continua circulant perquè la inductància l'hi imposa. El condensador comença a carregar en l'altre sentit i fa aparèixer una tensió en els borns de la inductància que fa disminuir el corrent. La situació continua de la següent manera: el condensador es va carregant cada vegada més lentament (perquè el corrent disminueix), mentre que el corrent va disminuint cada vegada més ràpidament (perquè la tensió inversa augmenta). Així, s'arriba a la situació en la qual el corrent s'anul·la i la tensió del condensador és màxima i del mateix valor que la tensió inicial, però amb sentit oposat. La situació és anàloga a la d'una massa sostinguda per un ressort. La inductància juga el paper de la massa. La massa té inèrcia i impedeix que el moviment canviï bruscament. La inductància impedeix que el corrent canviï bruscament. Vegem les equacions.

El comportament elèctric del condensador està descrit per l'equació: \scriptstyle{I=C{dV\over dt}}. El de la inductància està descrit per \scriptstyle{V=L{dI\over dt}}. Com en el esquema \scriptstyle{I} és positiu quan surt del costat positiu de la inductància, cal afegir un signe negatiu: \scriptstyle{V=-L{dI\over dt}}. Es té, doncs, aquest sistema d'equacions diferencials:

 I=C{dV\over dt}
 V=-L{dI\over dt}

Per eliminar \scriptstyle{I}, n'hi ha prou derivar la primera equació, per reemplaçar la derivada de I en la segona:

 V=-LC{d^2V\over dt^2}

que es pot escriure com a:

 L{d^2V\over dt^2}=-{1\over C}V

Aquesta equació és la mateixa que la de la massa amb un ressort.  \scriptstyle {V} és equivalent a la posició  \scriptstyle {y} .  \scriptstyle {L} és equivalent a la massa  \scriptstyle {m} i  \scriptstyle {{1 \ over C}} és equivalent a la constant del ressort\scriptstyle{k}.

La solució és:

V=V_\circ\cos(\omega t + \phi)

amb

\omega = {1\over \sqrt{LC}}

Com de costum, \scriptstyle{V_\circ} i \scriptstyle{\phi} depenen de les condicions inicials.

Circuit LC amb pèrdues[modifica | modifica el codi]

Circuit LC amb pèrdues. La resistència s'adona de totes les pèrdues possibles.

L'esquema de la dreta representa un circuit oscil·lant LC amb pèrdues. Les pèrdues estan representades per les pèrdues en una resistència. En un circuit real, les pèrdues provenen de resistències en sèrie com la dibuixada. Aquestes resistències poden estar en l'exterior de la inductància o del condensador, però també poden ser resistències internes d'aquests components. També hi pot haver resistències en paral·lel, perdudes en el dielèctric del condensador o al nucli de la bobina (si és ferromagnètic). També hi pot haver pèrdues per radiació d'ones electromagnètiques. La resistència farà que la tensió sobre la bobina sigui diferent de la tensió sobre el condensador. El corrent creada serà menor que si no hi hagués hagut pèrdues i quan el corrent carregui de nou el condensador, la tensió a la qual arribarà serà menor. Per la seva banda, l'amplitud disminuirà i tendirà cap a zero. L'equació del nou sistema és:

 L{d^2V\over dt^2}+R{dV\over dt}+{1\over C}V = 0

L'equació és la mateixa que la d'una massa amb un ressort i amb un esmorteïdor. Aquesta vegada  \scriptstyle {R} és l'equivalent del coeficient de fregament \scriptstyle{b}. La solució és:

V=V_\circ e^{-{R\over 2L}t}\cos(\omega t + \phi)

amb

\omega = \sqrt{{1\over LC}-\left({R\over 2L}\right)^2}

i

Q = \sqrt{L \over RC}={\omega_\circ L \over R}= {R \over \omega_\circ C}

on \scriptstyle{\omega_\circ={1\over \sqrt{LC}}} és la freqüència pròpia del circuit (sense pèrdues).

Oscil·lacions forçades d’un circuit LC amb pèrdues[modifica | modifica el codi]

Circuit LRC atacat per un generador sinusoïdal.

L'esquema de la dreta mostra un generador connectat a un circuit LC en sèrie. Si la tensió del generador és \scriptstyle{V_f\cos(\omega t)}, l'equació és:

 LC{d^2V \over dt^2}+RC{dV \over dt} + V = V_f\cos(\omega t)

L’expressió es pot reescriure, donant-li un aspecte similar a les formes precedents:

 L{d^2V \over dt^2}+R{dV \over dt} + {1\over C}V = {V_f \over C}\cos(\omega t)

Com en l'exemple mecànic, en règim estacionari la solució és:

 V=V_\circ\cos(\omega t + \phi)

on

 V_\circ = V_f\left|{1 \over 1 - \left({\omega\over \omega_\circ}\right)^2 + j{1\over Q}{\omega\over \omega_\circ}}\right|

i

 \phi = Arg\left({1 \over 1 - \left({\omega\over \omega_\circ}\right)^2 + j{1\over Q}{\omega\over \omega_\circ}}\right)

\scriptstyle{\omega_\circ} i \scriptstyle{Q} són els mateixos que en el paràgraf precedent. L’amplitud de la tensió de sortida és màxima a la ressonància (quan \scriptstyle{\omega=\omega_\circ}) i val \scriptstyle{Q} vegades la tensió d'entrada.

Oscil·lador harmònic quàntic[modifica | modifica el codi]

Funcions d'ona per als primers sis autoestats, n = 0 a 5. L’eix horitzontal mostra la posició y en unitats (h/2πmω)1/2. Les gràfiques estan sense normalitzar.
Densitats de probabilitat dels primers autoestats (dimensió vertical,amb els de menor energia en la part inferior) per a les diferents localitzacions espacials (dimensió horitzontal).

En l'esmentada, l'oscil·lador harmònic es pot emprar per estudiar sistemes que realitzin petites oscil·lacions al voltant d'una posició d'equilibri. En particular, l' oscil·lador harmònic quàntic es pot emprar per estudiar les oscil·lacions dels àtoms d'una molècula diatòmica, com la d'hidrogen, H 2, o la de clorur d'hidrogen, HCl.[7] L'oscil·lador harmònic és un dels casos en què es pot obtenir una solució analítica senzilla de l'equació de Schrödinger. En aquesta situació, l'hamiltonià de la partícula considerada estarà descrit per:

\hat H = \frac{\hat p^2}{2m} + \frac{1}{2} m \omega^2 y^2

Cal notar que per al cas de molècules diatòmiques, la massa m seria, en realitat, la massa reduïda del sistema. Es veu clarament que el primer sumand és un terme cinètic, mentre que el segon és l'harmònic. Com el hamiltonià no depèn del temps, només resta resoldre l'equació de Schrödinger independent del temps, per tal de trobar els autoestats de l'energia E:

 \hat H |\psi(y)\rang = E |\psi(y)\rang

Es pot demostrar que les funcions d'ona, \psi(y), el mòdul al quadrat descriu la densitat de probabilitat que la partícula tingui una determinada posició y, són el producte d'exponencials pels polinomis d'Hermite. La figura de la dreta mostra la forma d'aquestes funcions pels sis autoestats amb energia més baixa (l'estat de menor energia és el que figura a la part superior de la mateixa). En particular, l'energia del nivell n-èsim serà:

 E_n = \hbar \omega \left(n + {1\over 2}\right) \qquad n=0, 1, 2, ...

on \scriptstyle{\hbar} és la constant de Planck.

És important assenyalar un parell de fets:

  • Els nivells d'energia es troben quantitzats, és a dir, només poden prendre una sèrie de valors discrets.
  • El nivell mínim d'energia no és zero, sinó  \scriptstyle {(1/2) \ hbar \omega}. Noteu que la funció d'ona d'aquest estat mostra que la partícula no es troba en tot moment en la posició d'equilibri  y = 0 .

A la segona figura, es mostren les densitats de probabilitat espacial de la partícula per als diferents autoestats. Cal notar que a mesura que creix l'energia de l'autoestat considerat (és a dir, l'ordre  n), les distribucions de probabilitat tendeixen a concentrar-se en els punts de retorn, o màxima amplitud. Aquesta situació és la que es dóna en el cas clàssic, si es defineix per a ell una densitat de probabilitat inversament proporcional a la velocitat de la partícula en cada punt.[8] Per tant, es compleix el principi de correspondència (és a dir, es poden predir els resultats que s'obtindrien en el límit clàssic).

Vegeu també[modifica | modifica el codi]

Referències[modifica | modifica el codi]

  1. Simmons, capítol 3
  2. Simmons, pàgines 84-87
  3. T.N. Palmer (1995): "A local deterministic model of Quantum Spin Measurement", Proceedings: Mathematical and Physical Sciences, Volume 451, Issue 1943, pp. 585-608
  4. Atractor de Duffing
  5. Cartwright, M.L., "Balthazar van der Pol", J. London Math. Soc., 35, 367-376, (1960).
  6. Marion, pàgines 103 i 104
  7. Tipler, pàgina 1190
  8. Guillén, pàgines 114 i 115

Bibliografia[modifica | modifica el codi]

  • Feynman, Leighton and Sands. Lectures on physics. Addison-Wesley. ISBN 0-8053-9045-6. 
  • R. Resnick and D. Halliday. Physics. John Wiley & Sons, 1996. ISBN 0-471-83202-2. 
  • Marion, Jerry B.. Dinámica clásica de las partículas y sistemas. Barcelona: Ed. Reverté, 1996. ISBN 84-291-4094-8. 
  • Simmons, George F.. Ecuaciones Diferenciales. Con aplicaciones y notas históricas. Aravaca (Madrid): McGraw-Hill, 1999. ISBN 84-481-0045-X. 
  • Tipler, Paul A.. Física para la ciencia y la tecnología (volumen 2). Barcelona: Ed. Reverté, 2000. ISBN 84-291-4382-3. 
  • Sánchez Guillén, Joaquín. Braun, Mijail A.. Física cuántica. Madrid: Alianza Editorial, 1993. ISBN 84-206-8145-8. 

Enllaços externs[modifica | modifica el codi]

A Wikimedia Commons hi ha contingut multimèdia relatiu a: Oscil·lador harmònic Modifica l'enllaç a Wikidata