Ruido cuántico

En comunicaciones ópticas y submilimétricas, se llama ruido cuántico al proceso estocástico (aleatorio) subyacente que, al combinarse con la señal, produce el ruido de disparo fotónico, característico de la naturaleza discreta de los fotones. O dicho de otra forma, el ruido fotónico de disparo puede atribuirse, indistintamente, a la naturaleza discreta de los fotones, o al efecto del ruido cuántico sobre la señal.

El ruido cuántico se debe a los fenómenos de decoherencia presentes en el mundo cuántico, es uno de los principales desafíos a superar en el desarrollo experimental, industrial y comercial de los dispositivos de procesamiento de información cuántica.

En electrodinámica cuántica se identifica con las fluctuaciones del punto cero (o fluctuaciones del vacío).

Densidad espectral de potencia

El ruido cuántico no puede detectarse directamente (no aparece en la medida de un Analizador de Espectro Óptico para fibra, por ejemplo), por ser una fluctuación de punto cero.[1]

No obstante, para explicar los efectos que induce, hay que asignarle una densidad espectral de potencia de h ν / 4 {\displaystyle h\nu /4} por cada grado de libertad ( k {\displaystyle k} = dimensiones del espacio de fase) asociado a la señal que se esté midiendo. Así, para una señal monomodo en una única polarización, k = 2 fases (I, Q):

n q = h ν 2 B o {\displaystyle n_{\mathrm {q} }={\frac {\mathrm {h\nu } }{\mathrm {2} }}B_{\mathrm {o} }}

donde nq es la potencia de ruido cuántico en vatios, h es la constante de Planck en J·s, ν {\displaystyle \nu } es la frecuencia en Hz, y Bo es en ancho de banda óptico (FWHM) en Hz.

Para una señal monomodo en dos polarizaciones, k = 4 fases (HI, HQ, VI, VQ), n q = h ν B o {\displaystyle n_{\mathrm {q} }=h\nu B_{\mathrm {o} }} .

Aunque la densidad espectral de potencia y la energía del fotón se miden en las mismas unidades, conceptualmente son diferentes. La densidad espectral de potencia debe interpretarse como w/Hz, en vez de como energía (J).

Ruido de disparo como manifestación del ruido cuántico

Se puede llegar a la fórmula de Schottky para la varianza del ruido de disparo fotónico, suponiendo que el voltaje de ruido cuántico v ¯ n {\displaystyle {\bar {v}}_{n}} es una variable aleatoria que se suma al voltaje de señal (vs, determinista). De esta forma, la potencia instantánea ( w = | v s + v ¯ n | 2 {\displaystyle w=\left\vert v_{s}+{\bar {v}}_{n}\right\vert ^{2}} ) es una nueva variable aleatoria, y podemos identificar su varianza con la del ruido de disparo fotónico.[2]

Si v ¯ n {\displaystyle {\bar {v}}_{n}} es una variable aleatoria Gaussiana k-dimensional de media cero y V a r [ v ¯ n ] σ 2 k = h ν B o k / 4 {\displaystyle Var[{\bar {v}}_{n}]\equiv \sigma ^{2}\cdot k=h\nu B_{o}\cdot k/4} , entonces w / σ 2 {\displaystyle w/\sigma ^{2}} es una variable aleatoria Chi-cuadrado no central de parámetros k, v s 2 / σ 2 {\displaystyle v_{s}^{2}/\sigma ^{2}} . Lo que lleva a que V a r [ w ] = 4 σ 2 v s 2 + 2 k σ 4 {\displaystyle Var[w]=4\sigma ^{2}v_{s}^{2}+2k\sigma ^{4}} ,

Suponiendo que la potencia de señal s v s 2 σ 2 {\displaystyle s\equiv v_{s}^{2}\gg \sigma ^{2}} , y usando el ancho de banda HWHM (Be = Bo/2), en vez de Bo (FWHM), llegamos a la fórmula de Schottky para fotones:

V a r [ w ] = 2 h ν B e s {\displaystyle Var[w]=2h\nu B_{\mathrm {e} }s}

Es el mismo resultado que obtendríamos suponiendo que el número de fotones, l {\displaystyle l} , detectados en un intervalo de duración T, es una variable aleatoria de Poisson de media y varianza iguales a s T / h ν {\displaystyle sT/h\nu } . En este caso, la potencia instantánea sería w = h ν l / T {\displaystyle w=h\nu \cdot l/T} , y su varianza:

V a r [ w ] = ( h ν ) 2 T 2 V a r [ l ] = h ν T s = 2 h ν B e s {\displaystyle Var[w]={\frac {(h\nu )^{2}}{T^{2}}}Var[l]={\frac {h\nu }{T}}s=2h\nu B_{\mathrm {e} }s}

donde hemos tenido en cuenta que B e = 1 / 2 T {\displaystyle B_{e}=1/2T} (ancho de banda de Nyquist).

Relación con el ruido térmico

La densidad espectral de potencia de ruido térmico u ( ν , T ) {\displaystyle u(\nu ,T)} que genera una fuente pasiva monomodo, en una única polarización, a la temperatura física T y a la frecuencia ν {\displaystyle \nu \,} es:

u T ( ν , T ) = h ν e h ν k T 1 {\displaystyle u_{T}(\nu ,T)={\frac {h\nu }{e^{\frac {h\nu }{kT}}-1}}}

Corresponde a la densidad espectral de potencia térmica radiada por un cuerpo negro en un solo modo, y una sola polarización, según la ley de Planck.

Algunos autores sostienen que el ruido cuántico es indisociable del ruido térmico y hay que sumárselo, aunque no pueda medirse directamente.[3]​ De esta forma quedaría:

n T ( ν , T ) = h ν e h ν k T 1 + h ν 2 = h ν 2 coth h ν / 2 k T {\displaystyle n_{T}(\nu ,T)={\frac {h\nu }{e^{\frac {h\nu }{kT}}-1}}+{\frac {h\nu }{2}}={\frac {h\nu }{2}}\coth {\frac {h\nu /2}{kT}}}

Para bajas frecuencias o altas temperaturas ( h ν k T {\displaystyle h\nu \ll kT} ), la densidad espectral de potencia de ruido térmico se reduce a la fórmula de Johnson (ley de Rayleigh-Jeans):

n T ( ν , T ) k T {\displaystyle n_{T}(\nu ,T)\simeq kT}

y para altas frecuencias o bajas temperaturas ( h ν k T {\displaystyle h\nu \gg kT} ), por ejemplo, en transmisión por fibra óptica:

n T ( ν , T ) h ν 2 {\displaystyle n_{T}(\nu ,T)\simeq {\frac {h\nu }{2}}}

es decir, en este último caso, el ruido térmico propiamente dicho, u T ( ν , T ) {\displaystyle u_{T}(\nu ,T)} , es despreciable frente al ruido cuántico.

En telecomunicaciones, u T ( ν , T ) {\displaystyle u_{T}(\nu ,T)} es lo que se mide con un Analizador de Espectro a la salida de una línea de transmisión óptica de longitud semi-infinita, que se encuentra en equilibrio térmico a la temperatura T, en todos sus puntos z. En esta situación el coeficiente de atenuación, α [ N p / K m ] {\displaystyle \alpha [Np/Km]} , es constante y mayor que el coeficiente de ganancia, γ [ N p / K m ] {\displaystyle \gamma [Np/Km]} , también independiente de z, y su cociente es el factor de Boltzmann:

γ α = N 2 N 1 = e x p ( h ν k T ) {\displaystyle {\frac {\gamma }{\alpha }}={\frac {N_{2}}{N_{1}}}=exp(-{\frac {h\nu }{kT}})}

donde estamos suponiendo un sistema cuántico de dos niveles para explicar la atenuación y ganancia de la línea, y que α {\displaystyle \alpha } y γ {\displaystyle \gamma } son proporcionales a las poblaciones de nivel inferior ( N 1 {\displaystyle N_{1}} ) y superior ( N 2 {\displaystyle N_{2}} ), respectivamente.

En este ejemplo, el ruido térmico puede tratarse como Emisión Espontánea Amplificada (ASE, por sus siglas en inglés) y como tal responde a la ecuación diferencial:

d u d z = ( γ α ) u + γ h ν {\displaystyle {du \over dz}=(\gamma -\alpha )u+\gamma h\nu }

La solución para la línea de transmisión descrita ( < z 0 {\displaystyle -\infty <z\leq 0} ) es:

u T ( z ) = h ν γ α γ = h ν e h ν k T 1 {\displaystyle u_{T}(z)=h\nu {\frac {\gamma }{\alpha -\gamma }}={\frac {h\nu }{e^{\frac {h\nu }{kT}}-1}}}

Si en vez de u {\displaystyle u} consideramos n = u + h ν / 2 {\displaystyle n=u+h\nu /2} , la ecuación diferencial quedaría:

d n d z = ( γ α ) n + ( γ + α ) h ν 2 {\displaystyle {dn \over dz}=(\gamma -\alpha )n+(\gamma +\alpha ){\frac {h\nu }{2}}}

y su solución:

n T ( z ) = h ν 2 α + γ α γ = h ν 2 ( 2 α γ 1 + 1 ) = h ν e h ν k T 1 + h ν 2 {\displaystyle n_{T}(z)={\frac {h\nu }{2}}\cdot {\frac {\alpha +\gamma }{\alpha -\gamma }}={\frac {h\nu }{2}}({\frac {2}{{\frac {\alpha }{\gamma }}-1}}+1)={\frac {h\nu }{e^{\frac {h\nu }{kT}}-1}}+{\frac {h\nu }{2}}}

conforme a lo considerado.

El caso γ > α {\displaystyle \gamma >\alpha } (inversión de población) tendría la misma solución si admitimos energías negativas para representar temperaturas negativas.[4]

Relación con el ruido de emisión espontánea amplificada

La emisión espontánea puede explicarse como una emisión estimulada por el ruido cuántico, en vez de por fotones reales como ocurre en la emisión estimulada propiamente dicha.

En una línea de transmisión óptica, la evolución de la densidad espectral de potencia de emisión espontánea amplificada (o ruido ASE, por sus siglas en inglés) se rige por una ecuación diferencial ordinaria lineal de primer orden. Según la referencia que se tome para definir la d.e.p. de ruido ASE, la ecuación adquiere una de las siguientes formas:

d r d z = ( γ α ) r + γ h ν {\displaystyle {dr \over dz}=(\gamma -\alpha )r+\gamma h\nu }
d n d z = ( γ α ) n + ( γ + α ) h ν 2 {\displaystyle {dn \over dz}=(\gamma -\alpha )n+(\gamma +\alpha ){\frac {h\nu }{2}}}
d m d z = ( γ α ) m + α h ν {\displaystyle {dm \over dz}=(\gamma -\alpha )m+\alpha h\nu }

donde m ( z ) = n ( z ) + h ν / 2 = r ( z ) + h ν {\displaystyle m(z)=n(z)+h\nu /2=r(z)+h\nu } , y γ ( z ) {\displaystyle \gamma (z)} , α ( z ) {\displaystyle \alpha (z)} son los coeficientes de ganancia y atenuación en cada punto z {\displaystyle z} de la línea de transmisión.

Estas tres definiciones de la d.e.p. de ruido ASE sólo difieren en la forma de incorporar el ruido cuántico h ν / 2 {\displaystyle h\nu /2} , y provienen de la mecánica cuántica. Cada una ellas se relaciona con una de las tres representaciones más habituales de estados en el espacio de fase:[5]

r ( z ) {\displaystyle r(z)} : Representación - P de Sudarshan, asociada al orden normal de operadores.
n ( z ) {\displaystyle n(z)} : Representación de Wigner, asociada al orden simétrico de operadores.
m ( z ) {\displaystyle m(z)} : Representación - Q de Husimi, asociada al orden anti-normal de operadores.

La solución a la ecuación diferencial en cada caso es:

r ( z ) = g r i n + g r e {\displaystyle r(z)=g\cdot r_{in}+g\cdot r_{e}} , con r e ( z ) = h ν 0 z γ ( z ) g ( z ) d z {\displaystyle r_{e}(z)=h\nu \int _{0}^{z}{\frac {\gamma (z')}{g(z')}}dz'} :
n ( z ) = g n i n + g n e {\displaystyle n(z)=g\cdot n_{in}+g\cdot n_{e}} , con n e ( z ) = h ν 2 0 z γ ( z ) + α ( z ) g ( z ) d z {\displaystyle n_{e}(z)={\frac {h\nu }{2}}\int _{0}^{z}{\frac {\gamma (z')+\alpha (z')}{g(z')}}dz'}
m ( z ) = g m i n + g m e {\displaystyle m(z)=g\cdot m_{in}+g\cdot m_{e}} , con m e ( z ) = h ν 0 z α ( z ) g ( z ) d z {\displaystyle m_{e}(z)=h\nu \int _{0}^{z}{\frac {\alpha (z')}{g(z')}}dz'}

donde m i n = n i n + h ν / 2 = r i n + h ν {\displaystyle m_{in}=n_{in}+h\nu /2=r_{in}+h\nu } , son el ruido ASE a la entrada de la línea de transmisión ( z = 0 {\displaystyle z=0} ), y g ( z ) {\displaystyle g(z)} es la ganancia de potencia de señal entre la entrada y el punto z {\displaystyle z} , que cumple:

d g d z = ( γ α ) g {\displaystyle {dg \over dz}=(\gamma -\alpha )g} , y por lo tanto: g ( z ) = e x p 0 z [ γ ( z ) α ( z ) ] d z {\displaystyle g(z)=exp{\int _{0}^{z}[\gamma (z')-\alpha (z')]dz'}}

Aunque son equivalentes, cada representación tiene sus ventajas. Por ejemplo r ( z ) {\displaystyle r(z)} se corresponde con la lectura de densidad espectral de potencia que nos daría un analizador de espectro óptico, por lo que es la definición preferida en la mayoría de textos de ingeniería.[6]​ Sin embargo, la fórmula del factor de ruido ( f {\displaystyle f} ), se simplifica mucho cuando se usa la definición de n ( z ) {\displaystyle n(z)} :

f = h ν / 2 n r 1 g + r e n r + 1 h ν / 2 n r = 1 + n e n r {\displaystyle f={\frac {h\nu /2}{n_{r}}}\cdot {\frac {1}{g}}+{\frac {r_{e}}{n_{r}}}+1-{\frac {h\nu /2}{n_{r}}}=1+{\frac {n_{e}}{n_{r}}}}

por lo que también está muy extendida,[7]​ sobre todo entre los físicos.[8]

n T 1 {\displaystyle n_{T1}} es la ordenada en el origen de la recta tangente a n ( s ) {\displaystyle n(s)} , en el punto ( s 1 , n 1 ) {\displaystyle (s_{1},n_{1})}
Ejemplo 1: Línea de transmisión compuesta por tramos puros de amplificación ( α = 0 {\displaystyle \alpha =0} ) alternados con tramos puros de atenuación ( γ = 0 {\displaystyle \gamma =0} ).
Ejemplo 2: Línea de transmisión genérica (en cada punto hay amplificación y atenuación simultáneamente). Puntos 1 y 3 con γ > α > 0 {\displaystyle \gamma >\alpha >0} , puntos 2 y 5 con α > γ > 0 {\displaystyle \alpha >\gamma >0} , y punto 4 con α = γ > 0 {\displaystyle \alpha =\gamma >0} .

Es interesante advertir que, usando r ( z ) {\displaystyle r(z)} estamos considerando implícitamente que todo el ruido de la línea de transmisión se genera allí donde existe amplificación, ya que los atenuadores puros ( γ = 0 {\displaystyle \gamma =0} ) no introducen ruido ( r e = 0 {\displaystyle r_{e}=0} ). En el extremo opuesto, al usar m ( z ) {\displaystyle m(z)} estamos atribuyendo todo el ruido generado en la línea de transmisión a los puntos donde existe atenuación, ya que los amplificadores puros ( α = 0 {\displaystyle \alpha =0} ) no introducen ruido ( m e = 0 {\displaystyle m_{e}=0} ). Y si usamos n ( z ) {\displaystyle n(z)} estamos considerando que amplificadores y atenuadores contribuyen en la generación de ruido por igual.

La equivalencia de las tres definiciones se aprecia mejor en una gráfica de evolución del ruido ( r {\displaystyle r} , n {\displaystyle n} , ó m {\displaystyle m} ) frente a la señal ( s {\displaystyle s} ), siendo s = g s i n {\displaystyle s=g\cdot s_{in}} . De esta forma se puede escribir la derivada de n {\displaystyle n} respecto a s {\displaystyle s} como:

d n d s = ( γ α ) n + ( γ + α ) h ν 2 ( γ α ) s = n n T s {\displaystyle {dn \over ds}={\frac {(\gamma -\alpha )n+(\gamma +\alpha ){\frac {h\nu }{2}}}{(\gamma -\alpha )s}}={\frac {n-n_{T}}{s}}}

donde:

n T ( z ) = h ν 2 α + γ α γ {\displaystyle n_{T}(z)={\frac {h\nu }{2}}\cdot {\frac {\alpha +\gamma }{\alpha -\gamma }}}

es la d.e.p. de potencia del ruido térmico en cada punto z {\displaystyle z} de la línea de transmisión, y en la gráfica queda representado como la ordenada en el origen de la recta tangente a la curva n ( s ) {\displaystyle n(s)} en el punto z {\displaystyle z} .

Se cumple que n T h ν / 2 {\displaystyle n_{T}\geq h\nu /2} donde la línea actúa como atenuador neto ( α > γ {\displaystyle \alpha >\gamma } ), con n T = h ν / 2 {\displaystyle n_{T}=h\nu /2} para los atenuadores puros ( γ = 0 {\displaystyle \gamma =0} ). Y que n T h ν / 2 {\displaystyle n_{T}\leq -h\nu /2} donde la línea de transmisión actúa como amplificador neto ( γ α {\displaystyle \gamma \geq \alpha } ), con n T = h ν / 2 {\displaystyle n_{T}=-h\nu /2} para los amplificadores puros ( α = 0 {\displaystyle \alpha =0} ). Es decir, n T {\displaystyle n_{T}} puede tomar cualquier valor excepto los comprendidos entre h ν / 2 {\displaystyle -h\nu /2} y h ν / 2 {\displaystyle h\nu /2} .

La relación señal a ruido decrece con z {\displaystyle z} , en todo z {\displaystyle z} , para las tres definiciones de ruido ( r {\displaystyle r} , n {\displaystyle n} , ó m {\displaystyle m} ), como puede verse en las gráficas de los dos ejemplos (las flechas indican el sentido de las z {\displaystyle z} crecientes).

Véase también

Referencias

  1. «A. R. Kerr, J. Randa, Thermal Noise and Noise Measurements, IEEE Microwave Magazine.». 
  2. «B. Bristiel et al. Intrinsic Noise Figure Derivation for Fiber Raman Amplifiers From Equivalent Noise Figure Measurement, IEEE LTIMC 2004». 
  3. «A. R. Kerr et al, Receiver Noise Temperature, the Quantum Noise Limit, and the Role of the Zero-Point Fluctuations, 8th Int Symp. on Space Terahertz Tech, 1997». 
  4. Smorodinski, Ya. (1983). La Temperatura. Mir Moscú. p. 112. ISBN 8000085909. 
  5. Glauber, Roy J. (2006). «Chapter 10: Density Operators and Quasiprobability Distributions». Quantum Theory of Optical Coherence: Selected Papers and Lectures (en inglés). Wiley-VCH Verlag. p. 390. ISBN 9783527406876. 
  6. Desurvire, E. (202). «Chapter 2: Fundamentals of noise in optical fibere amplifiers». Erbium Doped Fiber Amplifiers. Principles and Applications (en inglés). John Wiley & Sons. p. 76. ISBN 0471589772. 
  7. Bristiel, B. (2004). «Intrinsic noise figure derivation for fiber Raman amplifiers from equivalent noise figure measurement». Proceedings of the Lightwave Technologies in Instrumentation and Measurement Conference. doi:10.1109/LTIMC.2004.1371010. 
  8. Caves, C. M. (2012). «Quantum limits on phase-preserving linear amplifiers». Physical Review A. doi:10.1103/PhysRevA.86.063802. 

Enlaces externos

  • Receiver Noise Temperature, the Quantum Noise Limit, and the Role of the Zero-Point Fluctuations
  • Optical Amplification
  • Intrinsic Noise Figure derivation for fiber Raman amplifiers
Control de autoridades
  • Proyectos Wikimedia
  • Wd Datos: Q7269075
  • Wd Datos: Q7269075