Saltar al contenido principal

# 2.2.2: Solución para la carga ligada

$$\newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} }$$

$$\newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}}$$

$$\newcommand{\id}{\mathrm{id}}$$ $$\newcommand{\Span}{\mathrm{span}}$$

( \newcommand{\kernel}{\mathrm{null}\,}\) $$\newcommand{\range}{\mathrm{range}\,}$$

$$\newcommand{\RealPart}{\mathrm{Re}}$$ $$\newcommand{\ImaginaryPart}{\mathrm{Im}}$$

$$\newcommand{\Argument}{\mathrm{Arg}}$$ $$\newcommand{\norm}[1]{\| #1 \|}$$

$$\newcommand{\inner}[2]{\langle #1, #2 \rangle}$$

$$\newcommand{\Span}{\mathrm{span}}$$

$$\newcommand{\id}{\mathrm{id}}$$

$$\newcommand{\Span}{\mathrm{span}}$$

$$\newcommand{\kernel}{\mathrm{null}\,}$$

$$\newcommand{\range}{\mathrm{range}\,}$$

$$\newcommand{\RealPart}{\mathrm{Re}}$$

$$\newcommand{\ImaginaryPart}{\mathrm{Im}}$$

$$\newcommand{\Argument}{\mathrm{Arg}}$$

$$\newcommand{\norm}[1]{\| #1 \|}$$

$$\newcommand{\inner}[2]{\langle #1, #2 \rangle}$$

$$\newcommand{\Span}{\mathrm{span}}$$ $$\newcommand{\AA}{\unicode[.8,0]{x212B}}$$

$$\newcommand{\vectorA}[1]{\vec{#1}} % arrow$$

$$\newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow$$

$$\newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} }$$

$$\newcommand{\vectorC}[1]{\textbf{#1}}$$

$$\newcommand{\vectorD}[1]{\overrightarrow{#1}}$$

$$\newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}}$$

$$\newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}}$$

$$\newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} }$$

$$\newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}}$$

$$\newcommand{\avec}{\mathbf a}$$ $$\newcommand{\bvec}{\mathbf b}$$ $$\newcommand{\cvec}{\mathbf c}$$ $$\newcommand{\dvec}{\mathbf d}$$ $$\newcommand{\dtil}{\widetilde{\mathbf d}}$$ $$\newcommand{\evec}{\mathbf e}$$ $$\newcommand{\fvec}{\mathbf f}$$ $$\newcommand{\nvec}{\mathbf n}$$ $$\newcommand{\pvec}{\mathbf p}$$ $$\newcommand{\qvec}{\mathbf q}$$ $$\newcommand{\svec}{\mathbf s}$$ $$\newcommand{\tvec}{\mathbf t}$$ $$\newcommand{\uvec}{\mathbf u}$$ $$\newcommand{\vvec}{\mathbf v}$$ $$\newcommand{\wvec}{\mathbf w}$$ $$\newcommand{\xvec}{\mathbf x}$$ $$\newcommand{\yvec}{\mathbf y}$$ $$\newcommand{\zvec}{\mathbf z}$$ $$\newcommand{\rvec}{\mathbf r}$$ $$\newcommand{\mvec}{\mathbf m}$$ $$\newcommand{\zerovec}{\mathbf 0}$$ $$\newcommand{\onevec}{\mathbf 1}$$ $$\newcommand{\real}{\mathbb R}$$ $$\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}$$ $$\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}$$ $$\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}$$ $$\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}$$ $$\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}$$ $$\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}$$ $$\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}$$ $$\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}$$ $$\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}$$ $$\newcommand{\laspan}[1]{\text{Span}\{#1\}}$$ $$\newcommand{\bcal}{\cal B}$$ $$\newcommand{\ccal}{\cal C}$$ $$\newcommand{\scal}{\cal S}$$ $$\newcommand{\wcal}{\cal W}$$ $$\newcommand{\ecal}{\cal E}$$ $$\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}$$ $$\newcommand{\gray}[1]{\color{gray}{#1}}$$ $$\newcommand{\lgray}[1]{\color{lightgray}{#1}}$$ $$\newcommand{\rank}{\operatorname{rank}}$$ $$\newcommand{\row}{\text{Row}}$$ $$\newcommand{\col}{\text{Col}}$$ $$\renewcommand{\row}{\text{Row}}$$ $$\newcommand{\nul}{\text{Nul}}$$ $$\newcommand{\var}{\text{Var}}$$ $$\newcommand{\corr}{\text{corr}}$$ $$\newcommand{\len}[1]{\left|#1\right|}$$ $$\newcommand{\bbar}{\overline{\bvec}}$$ $$\newcommand{\bhat}{\widehat{\bvec}}$$ $$\newcommand{\bperp}{\bvec^\perp}$$ $$\newcommand{\xhat}{\widehat{\xvec}}$$ $$\newcommand{\vhat}{\widehat{\vvec}}$$ $$\newcommand{\uhat}{\widehat{\uvec}}$$ $$\newcommand{\what}{\widehat{\wvec}}$$ $$\newcommand{\Sighat}{\widehat{\Sigma}}$$ $$\newcommand{\lt}{<}$$ $$\newcommand{\gt}{>}$$ $$\newcommand{\amp}{&}$$ $$\definecolor{fillinmathshade}{gray}{0.9}$$

Resolveremos la ecuación de movimiento para un caso particular de onda incidente muy significativo, que es el que nos interesará en el resto del curso: la onda armónica plana

\begin{aligned} &\mathbf{E}=\mathbf{E}_{0} e^{i(\mathbf{k} \cdot \mathbf{r}-\omega t)} \\ &\mathbf{B}=\frac{1}{\omega} \mathbf{k} \wedge \mathbf{E} \end{aligned}
no hay más que utilizar estos campos en la expresión de la ecuación 1.1.5. La ecuación resultante es muy complicada,
$m \ddot{\mathbf{r}}=q(\mathbf{E}+\dot{\mathbf{r}} \wedge \mathbf{B})-m \omega_{0}^{2} \mathbf{r}-\gamma m \dot{\mathbf{r}} \notag$
de modo que vamos a introducir dos aproximaciones:

1. Para una oap sabemos que
$\left|\mathbf{B}_{0}\right|=\frac{1}{c}\left|\mathbf{E}_{0}\right| \notag$
si esto se lleva a la ecuación
$\mathbf{F}=q(\mathbf{E}+\dot{\mathbf{r}} \wedge \mathbf{B}) \notag$
el segundo término va en $$\frac{\dot{r}}{c}$$. Así, podemos aproximar la fuerza sólo por la fuerza eléctrica, $$\mathbf{F} \simeq q \mathbf{E}$$, lo que es como suponer que $$\dot{r} \ll c$$. Si el campo em fuera suficientemente fuerte (por ejemplo, el debido a un laser) esto no sería cierto. En todo caso, dejamos pendiente de comprobación la sensatez de la hipótesis. La ecuación queda
\begin{aligned} m \ddot{\mathbf{r}} &=q \mathbf{E}-m \omega_{0}^{2} \mathbf{r}-\gamma m \dot{\mathbf{r}} \\ &=q \mathbf{E}_{0} e^{i(\mathbf{k} \cdot \mathbf{r}-\omega t)}-m \omega_{0}^{2} \mathbf{r}-\gamma m \dot{\mathbf{r}} \end{aligned}

2. La ecuación $$2.1$$ es no lineal en la medida en que la incógnita $$\mathbf{r}$$ aparece en un exponente. Admitiremos que en la solución la separación entre cargas positivas y negativas cumple:
$|\mathbf{r}(t)-\mathbf{r}(0)| \ll \lambda \notag$
El tamaño típico de los átomos es del orden del nm, y $$\mathbf{r}(t)$$ cabe esperar que será como máximo el tamaño del átomo. La longitud de onda del visible es del orden de 500 veces mayor. Así,
$\mathbf{k} \cdot \mathbf{r}(t) \simeq \mathbf{k} \cdot \mathbf{r}(0)=c t e \notag$
Entonces podemos incluir este factor constante en la amplitud $$\mathbf{E}_{0}$$. Con ello queda la siguiente expresión para la fuerza de LORENTZ
$\mathbf{F} \simeq q \mathbf{E}_{0} \mathbf{e}^{-i \omega t} \notag$
La ecuación que tenemos que resolver se reduce una vez admitidas ambas hipótesis a:
$\ddot{\mathbf{r}}+\gamma \dot{\mathbf{r}}+\omega_{0}^{2} \mathbf{r}=\frac{q}{m} \mathbf{E}_{0} e^{-i \omega t} \notag$

que es la ecuación de un oscilador armónico amortiguado y forzado. La teoría de edos nos aporta
$\mathbf{r}_{g i}=\mathbf{r}_{p i}+\mathbf{r}_{g h}$
La solución general de la homogénea $$\left(\mathbf{r}_{g h}\right)$$ la tomamos con $$\mathbf{E}_{0}=0$$, ya que dependerá de algo parecido a $$e^{-\gamma t}$$. Siendo $$\gamma$$ tan grande $$\left(\simeq 10^{8} s^{-1}\right)$$ el tiempo del transitorio que es la solución de la homogénea es del orden del 10ns. Aquí no nos interesaremos por esos transitorios rápidos. Nos quedaremos pues con una solución particular de la inhomogénea $$\left(\mathbf{r}_{p i}\right)$$, que dependerá del campo eléctrico.
$\mathbf{r}=\mathbf{r}_{0} e^{-i \omega t} \notag$
$$\mathbf{r}_{0}$$ se obtiene insertando esta solución en la ecuación 2.2. Se calculan las derivadas $$\dot{\mathbf{r}}=$$ $$-i \omega \mathbf{r}_{0} e^{-i \omega t}$$ y $$\ddot{\mathbf{r}}=-\omega^{2} \mathbf{r}_{0} e^{-i \omega t}$$ y la ecuación queda
$\left(-\omega^{2}-i \omega \gamma+\omega_{0}^{2}\right) \mathbf{r}_{0} e^{-i \omega t}=\frac{q}{m} \mathbf{E}_{0} e^{-i \omega t} \notag$
despejando, la amplitud y la ecuación de la trayectoria de la carga ligada resultan ser
\begin{align} \mathbf{r}_{0} &=\frac{\frac{q}{m}}{\omega_{0}^{2}-\omega^{2}-i \gamma \omega} \mathbf{E}_{0} \notag \\ \mathbf{r}(t) &=\frac{\frac{q}{m}}{\omega_{0}^{2}-\omega^{2}-i \gamma \omega} \mathbf{E}_{0} e^{-i \omega t} \end{align}
En el caso de un átomo $$m$$ y $$q$$ son la masa y la carga del electrón, respectivamente.

## Comentarios a la ecuación de la trayectoria $$2.3$$

1. El desplazamiento de las cargas es lineal en el campo.
2. La constante de proporcionalidad es compleja, lo que deriva del hecho de que hay pérdidas (profundizaremos en esto más adelante).
3. La constante de proporcionalidad depende de la frecuencia de forzamiento. Sabemos que si la frecuencia de forzamiento es muy próxima a la de resonancia, la respuesta será muy fuerte.

## Cumplimiento de las hipótesis

En cuanto al cumplimiento de las hipótesis que nos han conducido a 2.3, se darán dos situaciones:

• $$\omega \neq \omega_{0}$$ (aproximación muy buena) y
• $$\omega=\omega_{0}$$ (resonancia). En este segundo caso la aproximación de estos cálculos clásicos es peor, aunque suficiente. Para este caso se podría pasar al modelo cuántico de la materia.
##### Example $$\PageIndex{1}$$

verificar numéricamente estas dos afirmaciones sobre las aproximaciones.

##### Exercise $$\PageIndex{1}$$

Evaluar el cumplimiento de las hipótesis $$\dot{r} \ll c$$ y $$|\mathbf{r}(t)-\mathbf{r}(0)| \ll \lambda$$ para luz solar y electrones como carga ligada.

Para la luz solar
$\left|\mathbf{E}_{0}\right|_{\text {solar }} \simeq 10^{3} \frac{V}{m}$
una frecuencia típica del visible es
$\omega \simeq 3 \times 10^{15} \frac{\mathrm{rad}}{\mathrm{s}}$
utilizamos el valor de $$\gamma \simeq 10^{8} \mathrm{~s}^{-1}$$ y los datos de carga y masa del electrón. Puesta en forma real, la solución para una carga libre es una oscilación de amplitud
$\left|\mathbf{r}_{0}\right| \simeq 2 \times 10^{-17} \mathrm{~m}=2 \times 10^{-8} \mathrm{~nm}$
y longitud de onda $$\lambda \simeq 600 \mathrm{~nm}$$. Vemos que la diferencia de orden de magnitud entre el tamaño del átomo ( $$\simeq$$ separación de las cargas) y la longitud de onda incidente justifica la prumera aproximación tomada
$|\mathbf{r}(t)-\mathbf{r}(0)| \leq 2 \times 10^{-8} n m \ll \lambda=600 \mathrm{~nm}$
En cuanto a la velocidad máxima, es del orden de
$\left|\mathrm{r}_{0}\right| \omega \simeq 6 \times 10^{-2} \frac{\mathrm{m}}{\mathrm{s}}$
lo cual queda muy lejos de la velocidad de la luz, lo que justifica (son 10 órdenes de magnitud...) el haber despreciado el efecto del campo magnético.

2.2.2: Solución para la carga ligada is shared under a CC BY-SA 1.0 license and was authored, remixed, and/or curated by Alvaro Tejero Cantero.