Showing posts with label Álgebra. Show all posts
Showing posts with label Álgebra. Show all posts

El carácter del espectro del hamiltoniano cuántico

Considera la ecuación estacionaria de Schrödinger en una dimensión,
\begin{equation}{H}\psi=E\psi\hspace{0.25in}\text{con}\hspace{0.25in}{H}\equiv-\frac{\hbar^2}{2m}\frac{d^2}{dx^2}+V\end{equation} con el potencial V tal que $\displaystyle{\lim_{|x|\to\infty}V(x)=\text{cte}}$. Al tratar con problemas de juguete seguido asumimos que los eigenvalores o energía E de H están determinados por un valor entero positivo k y nos servimos denotar al k-ésimo nivel de energía como ${E_k}$. O bien, en el caso de la partícula libre por ejemplo, sabemos que aquél valor puede ser cualquier número real, o inclusive en el átomo de hidrógeno, que puede ser de ambos tipos.

Para sentar ideas, considera el siguiente esquema
Imagen (recarga la página)
uno puede asegurar que los valores ${E_1}$ y ${E_2}$ pertenecen a la parte discreta del espectro de H, sin embargo el valor ${E_3}$ (una mera etiqueta) pertenece en general a la parte continua del espectro de H.

En general se puede asegurar que el espectro de H será discreto siempre que
\begin{align}E<\lim_{|x|\to\infty}V(x)&\label{dagger}\\[0.1in]\lim_{|x|\to\infty}\psi(x)=0&\label{ddagger}\end{align} lo que a su vez no significa que, por ejemplo, si (\ref{dagger}) no se satisface, el espectro será siempre continuo. Es decir,
$$\begin{array}{lll}(\ref{dagger})\,\&\,(\ref{ddagger})&\Rightarrow&\left\{E_k\right\}_{k\in\mathbb{N}}\\[0.1in]\left\{E_k\right\}_{k\in\mathbb{N}}&\nRightarrow&(\ref{dagger})\,\&\,(\ref{ddagger})\end{array}$$ esto parece ser evidente de forma matemática, ya que ${\mathbb{Z}\subset\mathbb{R}}$, sin embargo físicamente es relevante distinguir entre los casos discreto y continuo del espectro. Cuando el espectro del hamiltoniano es discreto, se dice que los estados ${|\psi\rangle}$ son ligados (bounded states), mientras que si es continuo, se dice que los correspondientes estados son no ligados (unbounded o scattering states). Los casos en que no se satisface (\ref{ddagger}) y el espectro es discreto, se dice que se tienen estados ligados en el continuo (bounded states in the continuum), caso que seguido podría darse en nanotecnología u óptica cuántica, por ejemplo.

Una interpretación física intuitiva para la aseveración anterior puede ser una analogía con las órbitas cerradas y abiertas que se hallan en los sistemas clásicos. Supóngase un núcleo masivo y un electrón en interacción mediante un potencial coulombiano. Si se satisfacen las dos condiciones, el sistema será un átomo hidrogenoide y como es sabido, tendrá valores de energía múltiplos de un entero positivo. En cambio si alguien disparase el electrón contra el núcleo desde muy lejos y con suficiente energía cinética, la energía total sería positiva y el núcleo desviaría al electrón sin capturarlo y sin alterar su energía, que evidentemente podrá tomar cualquier valor real. El caso más simple para un espectro energético continuo es el de la partícula libre, esto es, ${V(\mathbf{r})=0,\,\forall\,\mathbf{r}\in\mathbb{R}^3}$ donde la condición de normalización ${\int\psi_m\psi_n\,d\mathbf{r}=\delta_m\delta_n}$ deja de ser válida y puede generalizarse mediante la delta de Dirac, además de que en general las eigenfunciones se salen del espacio de Hilbert ${\mathcal{L}^2}$.

Ahora bien, he asumido ambas condiciones, sin embargo la condición (\ref{ddagger}) es una puramente cuántica y se refiere a una de más características de lo que seguido llamamos funciones de onda físicamente aceptables. La condición en realidad puede aflojarse y pedir simplemente que la función de onda sea constante en ${|x|\to\infty}$. Una forma intuitiva de entender por qué es necesaria esta condición es la siguiente. Supón un sistema que satisface (\ref{dagger}) y que eliges una solución particular $\psi$ que decrece en ${x\to\infty}$. Para este propósito, sean ${x_1,\,x_2}$ con ${x_1<x_2}$ los puntos en los que E intersecta al potencial y que $\psi$ es positivo en ${x_2}$. Entonces $\psi$ debe ser apropiadamente convexa en ${x>x_2}$ de modo que (\ref{ddagger}) se satisfaga. Para que esto ocurra, de la ecuación de la forma ${\psi^{\prime\prime}=\zeta^2\psi}$, se sabe que la curvatura depende de la energía vía ${\zeta^2}$, y por tanto, debe haber sólo un cierto valor particular de la energía para esta $\psi$ tal que la curvatura sea la apropiada dado que debe satisfacer (\ref{dagger}).

Esta es básicamente la forma de identificar inmediatamente el carácter discreto del espectro energético. Ahora bien, las identificaciones físicas o la justificación como la anterior son suficientes de manera práctica, pero gracias a unas preguntas en Physics SE me he preguntado a detalle por la forma de mostrarlo matemáticamente. Lo que se me ha ocurrido de inmediato y hace ya algún tiempo, es utilizar la teoría de Sturm-Liouville.

El detalle es el siguiente, siempre que exista $\zeta$ en el conjunto resolvente de un operador simétrico T que satisface el problema ${T\varphi=\lambda\varphi}$ tal que el operador resolvente ${R(T,\zeta)}$ sea compacto y simétrico, el espectro de T será discreto ${\{\lambda_k\}_{k\in\mathbb{N}}}$. Para que esto ocurra, T debe satisfacer el problema general de Sturm-Liouville regular en forma autoadjunta ${(T-\zeta)\,v=f}$ en un intervalo ${[a,b]}$ con función de peso $\rho$, i.e. en ${\mathcal{L}^2_\rho(a,b)}$. Otras propiedades bastante importantes en este caso, son que los eigenvalores ${\lambda_k}$ son únicos y constituyen un conjunto infinito y creciente, y que el conjunto de eigenfunciones ${\{\varphi_{_k}\}_{k\in\mathbb{N}}}$ forma una base ortogonal completa del espacio ${\mathcal{L}^2_\rho(a,b)}$. Es un clavado un tanto profundo para llegar aquí y desarrollarlo podría llevar un buen espacio, de cualquier modo el lector puede acudir a cualquier libro estándar de análisis funcional o en específico de la teoría de Sturm-Liouville.

La solución explícita por método de función de Green del problema regular de Sturm-Liouville en forma autoadjunta ${(T-\zeta)\,v=f}$ es
\begin{equation}v=(T-\zeta)^{-1}f=\int_a^bG(x,s)\,f(s)\,ds\end{equation} donde la función de Green
\begin{equation}G(x,s)=\begin{cases}\frac{\varphi_b(x)\,\varphi_a(s)}{W(s)},&a\leq{s}\leq{x}\\[0.1in]\frac{\varphi_a(x)\,\varphi_b(s)}{W(s)},&x\leq{s}\leq{b}\end{cases}\end{equation} se construye con las soluciones ${\varphi_a,\,\varphi_b}$ de la ecuación homogénea ${T\varphi=0}$ que satisface las respectivas condiciones de frontera (Neumann, Dirichlet, mixtas) en ${x=a,\,b}$. De este modo el resolvente asociado a T, por definición, ${R\equiv(T-\zeta)^{-1}}$ se puede relacionar de manera explícita como
\begin{equation}R\,f=\int_a^bG(x,s)\,f(s)\,ds\end{equation} donde G se interpreta como el kernel o núcleo de R. Al mostrar que el kernel de R es de cuadrado integrable, i.e. que satisface
\begin{equation}\int_a^b\int_a^b\left|G(x,s)\right|^2ds\,dx<\infty\end{equation} uno puede mostrar que R es compacto. Para esto uno simplemente debe mostrar que G es continua en el intervalo finito ${[a,b]}$. Ahora bien, por definición del operador resolvente, uno puede hacer el álgebra y ver que
\begin{equation}R\varphi=\mu\varphi\,\Longrightarrow\,T\varphi=\left(\frac{1}{\mu}+\zeta\right)\varphi\end{equation} de donde se sigue que las eigenfunciones de T son las mismas que las de R y entonces que R también será simétrico en ${[a,b]}$, además se pueden seguir las propiedades mencionadas para los eigenvalores ${\lambda\equiv\frac{1}{\mu}+\zeta}$ tomando en cuenta que para un operador compacto y simétrico en un espacio de Hilbert $\mathcal{H}$ con producto interior ${\left\langle\cdot\,\middle|\,\cdot\right\rangle_\mathcal{H}}$ los eigenvalores forman un conjunto infinito, acotado y numerable que converge a cero, $\displaystyle{\lim_{k\to\infty}\mu_k=0}$. Así pues, si la ecuación ${H\psi=E\psi}$ puede reducirse a un problema regular de Sturm-Liouville, el espectro de H será discreto.

Al querer atacar el problema de este modo sin embargo, aparecen ciertas inconveniencias que aún no me parecen tan claras con este argumento (aunque las haya procurado defender en mis respuestas de Physics SE). La primera es que el intervalo (o dominio en general) en que la ec. de Schrödinger es válida no está restringido, como pide el problema de S-L regular. Aquí me parece crucial la restricción (\ref{dagger}), ya que en general para una energía que intersecta a V en ${x=x_1,\,x_2}$ con ${x_1<x_2}$, la función de onda tendrá la forma
\begin{equation}\psi=\begin{cases}\psi_1,&x<x_1\\\psi_2,&x_1\leq{x}\leq{x}_2\\\psi_3,&x>x_2\end{cases}\end{equation} de modo que presumo que puede emplearse ${\psi_2}$ para reducir el problema original a uno de S-L regular y obtener información del espectro de H. Para esto, H debe ser autoadjunto en ${[x_1,x_2]}$, lo que puede lograrse homogeneizando las condiciones de frontera para ${\psi_2}$, por ejemplo por interpolación de una nueva función a través de la original, que satisfaga las condiciones de frontera originales. Esto me parece suficiente así visto toscamente, sin embargo me sigue pareciendo que podría haber sutilezas o algo por el estilo; quizá pronto me decida a intentarlo con detalle. La segunda, es que originalmente H se ha establecido como un operador simétrico en general en todo el espacio en que actúa. Esta inconveniencia desaparece si las funciones de onda completas $\psi$ en general satisfacen (\ref{ddagger}). Nuevamente esto me parece dudoso, aunque a la vez, visto toscamente, cierto. De igual modo hay ciertas cuestiones quizá un tanto más complejas, como el carácter ya sea acotado o libre del operador H, pues por el teorema de Hellinger-Toeplitz un operador autoadjunto es acotado, sin embargo el hamiltoniano (la energía) en general no es acotado. En el enlace anterior a Wikipedia se puede leer que en ese caso H queda definido sólo en un subconjunto denso de ${\mathcal{L}^2}$, además de que el hamiltoniano presuntamente siempre está acodado por abajo. En fin, dejo estos detalles más sutiles al lector, por mi parte procuraré profundizar un tanto más en el tema conforme tenga oportunidad o vaya aprendiendo quizá sobre teorías un tanto más fundamentales como teoría cuántica de campos, por ejemplo.

El 211 aniversario del nacimiento de Niels Henrik Abel

Abel nació en Nedstrand, Noruega, el 5 de Agosto de 1802 y murió el 6 de Abril de 1829. Su vida fue un tanto trágica, y es otro personaje que en su corta vida hizo aportaciones vitales a las matemáticas. Su vida seguido se vio perseguida por los problemas económicos y de salud, e incluso falleció de tuberculosis esperando un mejor puesto de trabajo.

Acá comparto la biografía de Abel (en inglés) en la que resaltan nombres como Gauss, Jacobi, Legendre y el que se volviera su gran amigo, August Crelle. También comparto algunos tuits sobre el trabajo de Abel por la cuenta del Instituto de Matemáticas de la UNAM.





El teorema de Nicómaco

¿Habías notado que
\begin{equation}(1+2+3+\ldots+n)^2=1^3+2^3+3^3+\ldots+n^3\end{equation} es decir
\begin{equation}\left(\sum_{k=1}^n{k}\right)^2=\sum_{k=1}^n{k^3}\label{nico1}\end{equation} ? Bueno, esta identidad me ha resultado impresionante, y la verdad es que no la conocía, siendo que es bastante antigua. Se atribuye principalmente al trabajo de Nicómaco de Gerasa por ahí de los siglos I y II.

Nicómaco hizo la siguiente observación:

Teorema de Nicómaco:
La (suma) del primer número natural impar es igual al primer cubo, la suma de los dos siguientes números impares es igual al segundo cubo, la suma los tres siguientes impares es el tercer cubo, ...

Bueno, lo he dicho con mis palabras, pero entiéndase como
\begin{align}1&=1^3\label{nico2}\\3+5&=2^3\label{nico3}\\7+9+11&=3^3\label{nico4}\\13+15+17+19&=4^3\\21+23+25+27+29&=5^3\\&\vdots\nonumber\\(n^2-n+1)+(n^2-n+3)+\ldots+(n^2+n-3)+(n^2-n-1)&=n^3,\;\forall{n}\in\mathbb{N}\end{align} es decir,
\begin{equation}\sum_{k=1}^n\left[(n^2-n-1)+2k\right]=n^3\end{equation} de aquí, sumando los respectivos lados de cada ecuación se debería recuperar (\ref{nico1}), pero antes de mostrarlo mediante este teorema de Nicómaco original, probémoslo de -la que me parece- la forma más elemental de hacerlo: Inducción Matemática.

Sabemos que (\ref{nico1}) es cierta para ${n=1}$, entonces supóngase que es cierta para algún ${n=m}$. Así pues, averigüemos qué pasa para ${n=m+1}$; tenemos que
\begin{equation}\left(\sum_{k=1}^{m+1}k\right)^2=\sum_{k=1}^{m+1}k^3\end{equation} escribiendo explícitamente el término $m+1$,
\begin{equation}{\left(\sum_{k=1}^{m}k+(m+1)\right)^2}=\sum_{k=1}^{m}k^3+(m+1)^3\end{equation} expandiendo el cuadrado,
\begin{equation}{\left(\sum_{k=1}^m{k}\right)^2+2(m+1)\sum_{k=1}^m{k}}+(m+1)^2=\sum_{k=1}^{m}k^3+(m+1)^3\end{equation} empleando la serie aritmética, $\sum_{k=1}^m{k}=m(m+1)/2$,
\begin{equation}{\left(\sum_{k=1}^m{k}\right)^2+m(m+1)^2+(m+1)^2}=\sum_{k=1}^{m}k^3+(m+1)^3\end{equation} es decir,
\begin{equation}{\left(\sum_{k=1}^m{k}\right)^2}=\sum_{k=1}^m{k^3}\end{equation} por tanto, en efecto, (\ref{nico1}) es cierta ${\forall{n\in\mathbb{N}}}$.

Demostrémosla ahora partiendo del Teorema de Nicómaco. Notemos que sumando los respectivos lados de (\ref{nico2}) con (\ref{nico3}), y luego (\ref{nico2}) con (\ref{nico3}) con (\ref{nico4}) y así con las siguientes, el Teorema nos dice que la suma de cubos de números naturales hasta $n$ es igual a la suma de tantos números impares como la serie aritmética hasta $n$, es decir
\begin{align}1^3+2^3&=1+3+5\\1^3+2^3+3^3&=1+3+5+7+9+11\\&\vdots\nonumber\\\sum_{k=1}^n{k^3}&=\sum_{k=1}^{n(n+1)/2}(2k-1)\end{align} Ahora es necesario hacer uso del llamado Teorema del Número Impar, que me limitaré a presentar con esta demostración sin palabras extraída del blog Bill The Lizard:
Imagen (si ves este texto, recarga la página)
Lo que dice este teorema, y como muestra la imagen, es que
\begin{equation}n^2=\sum_{k=1}^{n}(2k-1)\end{equation} Así entonces concluímos que
\begin{equation}\sum_{k=1}^n{k^3}=\left(\frac{n(n+1)}{2}\right)^2\end{equation} esto es,
\begin{equation}\sum_{k=1}^n{k^3}=\left(\sum_{k=1}^n{k}\right)^2\end{equation} como se esperaba.

Otra forma bastante interesante de presentar esta identidad, es como sigue,
\begin{align}\sum_{k=1}^n{k^3}&=\left(\frac{(n+1)!}{2(n-1)!}\right)^2\nonumber\\[0.1in]&=\left(\frac{(n+1)!}{2!((n+1)-2)!}\right)^2\nonumber\\[0.1in]&=\binom{n+1}{2}^2\end{align} y digo interesante, porque finalmente lo que expresa son dos formas distintas de contar una cantidad dada, así que es de esperar que se pueda obtener una expresión que involucre un combinatorio. Quizá demostrar la identidad en esta última forma, utilizando únicamente combinatoria sería una tarea un poco más elaborada (o lo contrario) pero también ingeniosa y divertida.

Las siguientes demostraciones sin palabras podrían ayudar también con un poco de intuición. Da clic para ver la fuente original.

Imagen (si ves este texto, recarga la página)
y también,
Imagen (si ves este texto, recarga la página)

Yutaka Taniyama

The following is an excerpt of the article by J J O'Connor and E F Robertson.
Read the original @ www-history.mcs.st-andrews.ac.uk/Biographies/Taniyama

YUTAKA TANIYAMA
Born: 12 Nov 1927 in Kisai (north of Tokyo), Japan
Died: 17 Nov 1958 in Tokyo, Japan

Imagen (recarga la página)His parents were Sahei, a medical doctor, and Kaku Taniyama. Yutaka was born into a large family having two older brothers and three older sisters as well as a younger brother and a younger sister. Yutaka was a sickly child and suffered from tuberculosis which caused him to miss two years of high school. After graduating from the high school, he entered the University of Tokyo to study mathematics. During his undergraduate years he read Claude Chevalley's Theory of Lie groups and André Weil's Foundations of algebraic geometry as well as two other books by Weil on algebraic curves and abelian varieties. He attended algebra lectures by Masao Sugawara and these encouraged his interest in number theory. He graduated in March 1953.

He remained at the University of Tokyo as a 'special research student' in the Department of Mathematics, although he had no thesis advisor. Shimura writes in [1] about the apartment where Taniyama lived in Tokyo:

... he lived in a one-room apartment which consisted of 81 square feet of living space, a sink, and a tiny unfloored part behind the door. Running water, gas and electricity were provided separately in each room, but there was only one toilet on each floor of the two-storey building, shared by all the occupants of the dozen or so rooms of the floor. I remember that his was No. 20 on the second floor, close to the last. Thus it was more like a dormitory than an apartment, but it was more or less typical of the time. To take a bath, he had to go to a public bathhouse, a few minutes' walk from his apartment. The building, a rather shabby wooden structure, was named poetically 'Villa Tranquil Mountains'...

Taniyama's interests were in algebraic number theory and his fame is mainly due to two problems posed by him at the symposium on Algebraic Number Theory held in Tokyo and Nikko in 1955. His meeting with André Weil at this symposium was to have a major influence on Taniyama's work. These problems form the basis of a conjecture: every elliptic curve defined over the rational field is a factor of the Jacobian of a modular function field. This conjecture proved to be a major factor in the proof of Fermat's Last Theorem by Andrew Wiles. In the Proceeding of the conference he published the paper Jacobian varieties and number fields, then in the following year the paper L-functions of number fields and zeta functions of abelian varieties.

Other than these two papers the only other paper Taniyama published was Distribution of positive 0-cycles in absolute classes of an algebraic variety with finite constant field (1958). However, in addition to these papers, he wrote the book Modern number theory (1957) in Japanese, jointly with G Shimura.

With seemingly a great future in front of him, both in mathematics and his life (he was planning marriage to Misako Suzuki) he took his own life. In a long suicide note he left, he took great care to describe exactly where he had reached in the calculus and linear algebra courses he was teaching and to apologise to his colleagues for the trouble his death would cause them. As to the reason for taking his life he says:

Until yesterday I have had no definite intention of killing myself. But more than a few must have noticed I have been tired both physically and mentally. As to the cause of my suicide, I don't quite understand it myself, but it is not the result of a particular incident, nor of a specific matter. Merely may I say, I am in the frame of mind that I lost confidence in my future. There may be some to whom my suicide will be troubling or a blow to a certain degree. I sincerely hope that this incident will cast no dark shadow over the future of that person. At any rate I cannot deny that this is a kind of betrayal, but please excuse it as my last act in my own way, as I have been doing all my life.

About a month later his fiancé Misako Suzuki also committed suicide. She left a note which included the sentences:

We promised each other that no matter where we went, we would never be separated. Now that he is gone, I must go too in order to join him.

Shimura writes [1]:

...he was the moral support of many of those who came into mathematical contact with him, including of course myself. Probably he was never conscious of this role he was playing. But I feel his noble generosity in this respect even more strongly now than when he was alive. And yet nobody was able to give him any support when he desperately needed it. Reflecting on this, I am overwhelmed by the bitterest grief.

One might reasonably ask what Taniyama's interests were other than mathematics. He enjoyed listening to music, especially Beethoven's Eighth Symphony, and going to movies, his favourite film being 'The King and I'. His only hobby was writing articles which he never intended to publish, but he seemed to find writing them helped to organise his thoughts. Examples of the topics he wrote articles on included: reviews of books, ideas on how researchers should be trained, how to organise a new institute for mathematical sciences, and reviews of articles by others.

Full article by: J J O'Connor and E F Robertson
April 2009

Fórmula e identidad de Euler

Bueno hoy he trabajado con la exponencial compleja y me han presentado la fórmula considerada la más bella de las matemáticas, a saber
$$e^{i\pi}+1=0$$ …oportunamente mi profesor de variable compleja ha dicho que si esta fórmula no nos conmueve, simplemente debemos dedicarnos a otra cosa; toda persona con conocimientos matemáticos vulgares y silvestres nota las constantes matemáticas, las operaciones y los números naturales fundamentales en esa cortísima fórmula.

Bien, veamos qué dice esta ecuación; porqué es verdadera. Para ello primero introduciré la fórmula de Euler.

Sabemos que podemos desarrollar $\displaystyle{e^{x},x\in\mathbb{R}}$ en serie de potencias como
$$e^{x}=\sum_{k=0}^{\infty}\frac{x^{n}}{n!}$$ esto es,
$$e^{x}=1+x+\frac{x^{2}}{2!}+\frac{x^{3}}{3!}+\cdots$$ entonces, se sigue que
$$e^{ix}=1+ix+\frac{(ix)^{2}}{2!}+\frac{(ix)^{3}}{3!}+\cdots$$ recordemos ahora que
$$i=\sqrt{-1}\,\Longleftrightarrow\,{i^{2}=-1}$$ por lo que podemos simplificar nuestra ecuación como
$$e^{ix}=1+ix-\frac{x^{2}}{2!}-\frac{ix^{3}}{3!}+\frac{x^{4}}{4!}+\cdots$$ Recordemos que un número complejo z tiene la forma ${z=u+iv}$, con $\displaystyle{u,v\in\mathbb{R}}$, entonces convenientemente escribimos que
$$e^{ix}=\left(1-\frac{x^{2}}{2!}+\frac{x^{4}}{4!}-\cdots\right)+i\;\left(x-\frac{x^{3}}{3!}+\frac{x^{5}}{5!}-\cdots\right)$$ …hmmmm, esas series infinitas en lo que llamamos u,v parecen conocidas… tal vez sean las primeras series que nos enseñaron para ejemplificar las series de Taylor… Claro!
$$e^{ix}=\cos{x}+i\;{\sin{x}}$$ Y voilà! ésta es la llamada fórmula de Euler. Inmediatamente se sigue que ${e^{z}=e^{u}\left(\cos{v}+i\;\sin{v}\right)}$ , que es la exponencial compleja. En el plano complejo vemos esto así
Imagen (recarga la página)
La identidad de Euler es un caso particular de la fórmula de Euler, para la que según lo que hemos dicho, ${x=\pi}$, y entonces se tiene
$$e^{i\pi}=\cos(\pi)+i\;\sin(\pi)=-1$$ …es decir
$$e^{i\pi}+1=0$$ Esta fórmula es sumamente celebrada por la comunidad matemática y científica, y no podía ser para menos. El análisis complejo resulta una tentación enorme para quien disfruta de las matemáticas, el nombre complejo simplemente ha sido acarreado, mas sin embargo es una parte del análisis que nos revela una infinidad de hechos que generalmente ignoramos y que el análisis de variable real calla.

Ya habrá notado el lector que la exponencial compleja es una función periódica, mientras que para una exponencial real siempre hemos pensado en una función monótona, ésta y otras propiedades hacen en sumo interesante aventarse un clavado al análisis complejo.

Las curvas de nivel de la parte imaginaria de la exponencial compleja se ven algo así (la periodicidad es evidente).
Imagen (recarga la página)

Diversos factoriales!

Desde los cursos de matemáticas del bachillerato se habla del factorial, el cual sabemos tiene una definición muy sencilla, ${\forall n\in\mathbb{N}}$:
$$n!=n\cdot{(n-1)}\cdot\ldots\cdot{2}\cdot{1}$$ Generalizando un poco más, y para establecer 0!=1, podemos definirlo con la siguiente relación de recurrencia $\displaystyle{\forall{n}\in\mathbb{Z},n\geq{0}}$:

Sea ${a_n=n!}$, de modo que
$${a_0=1}\\
{a_{n+1}=(n+1)a_n}$$ En cualquier carrera de ciencias o ingeniería es muy común hacer uso del factorial, aparece en muchas formulaciones de utilidad y a veces cuando se quiere modelar o hacer un sistema efectivo, se necesita emplear de forma inteligente, por ejemplo:
  • La fórmula de Stirling para aproximación de factoriales grandes: ${n!\approx\sqrt{2\pi n}\left(\frac{n}{e}\right)^n}$.
  • Combinatoria, coeficientes binomiales: ${\left(\begin{array}{c}n\\k\end{array}\right)=\frac{n!}{k!(n-k)!}}$
  • La función Gamma, una generalización más amplia del factorial: ${\Gamma (z+1)=z\Gamma (z)}$
Pero bueno, aunque no conozcamos alguna de las anteriores, todos conocemos bien el factorial, pero, ¿es acaso éste el único tipo de factorial?, ¿ya conoces los multifactoriales o el subfactorial? Bien, comencemos con los multifactoriales introduciendo el doblefactorial. El doblefactorial de n, como era de esperarse, se denota n!!. La definición es muy sencilla, el doblefactorial de n es el producto de todos los enteros positivos menores o iguales a n que tienen la misma paridad que n. Es decir:
$$n!!=\begin{cases}n\cdot(n-2)\cdot\ldots\cdot4\cdot2,&\text{si n es par}\\n\cdot (n-2)\cdot\ldots\cdot 3\cdot1,&\text{si n es impar}\end{cases}$$ También por definición tenemos que 0!!=1 y en especial para multifactoriales (-1)!!=1. El doble factorial cumple varias relaciones que he encontrado en las referencias [1], [2] y [3] que resultan bastante interesantes, y son a considerar:
$$\int_0^{\pi/2}(\sin\theta)^{2n+1}\mathrm{d}\theta =\frac{(2n)!!}{(2n+1)!!}$$ y
$$\sum_{i=0}^{\infty}\frac{x^{2i}}{(2i)!!}=e^{x^{2}/2}$$ De la definición del doble-factorial, se sigue casi inmediatamente la del multifactorial. Podemos decir que el k-factorial de n, que denotaré ${n!^{(k)}}$, es el producto de enteros positivos menores o iguales a n congruentes con ${n\mod{k}}$. Es decir:
$$n!^{(k)}=n\cdot (n-k)\cdot (n-2k)\ldots$$ Donde los factores siempre son enteros mayores a cero. Por ejemplo podemos decir que ${8!!!=8\times{5}\times{2}}$ o que ${11!^{(4)}=11\times{7}\times{3}}$. Ésta por supuesto no es una definición muy rigurosa, pero puedes ahondar en el tema si deseas.

Para finalizar hablaré del subfactorial. En los últimos problemas que he resuelto de Project Euler me he topado con varios problemas que requieren programar combinatorias (o que una buena forma para resolverlos es usando combinatorias) y el factorial ha sido de gran ayuda. Los factoriales nos hablan de todas las formas en que un conjunto de n elementos distintos puede combinarse, que es precisamente n! formas. Ahora supón que no quieres que un elemento ocupe una posición que ya ha ocupado en alguna combinación previa, es decir no quieres que se quede en una misma posición. Pues entonces usas el subfactorial de n.

El subfactorial de n, denotado !n o de preferencia n¡, es el número de permutaciones de n objetos para las que ningún objeto aparece en su posición “natural”. Por ejemplo el conjunto {1,2,3} puede combinarse del modo descrito con las permutaciones {2,1,3}, {3,2,1}, por lo que decimos 3¡=2. Los subfactoriales satisfacen la siguiente relación de recurrencia:
$${n\text{<}=n\cdot(n-1)\text{<}+(-1)^{n}}$$ Por definición, también el subfactorial 0¡=1. Formas para calcular el subfactorial sencillamente, son las siguientes:
$$n\text{<}=\sum_{i=0}^{n} i\text{<}(-1)^{n-i}\left(\begin{array}{c}n\\i\end{array}\right)$$ para ${n\geq{1}}$, donde ${\left(\begin{array}{c}n\\i\end{array}\right)}$ es el coeficiente binomial. Y otra un poco más sencilla es:
$$n\text{<}=\left\lfloor{\frac{n\text{<}+1}{e}}\right\rfloor$$ para ${n\geq{1}}$, donde ${\lfloor{x}\rfloor}$ denota el mayor entero menor que x.

Los subfactoriales se implementan en Mathematica con la instrucción Subfactorial[n]. El comportamiento del subfatorial se puede observar mediante el siguiente gráfico de n¡, con la parte real marcada en color azul y la parte imaginaria en anaranjado.
Para más información sobre estos temas consulta las siguientes ligas:

[1] Weisstein, Eric W. "Double Factorial." From MathWorld--A Wolfram Web Resource.
[2] The Endeavour
[3] Gaussianos
[4]Weisstein, Eric W. "Subfactorial." From MathWorld--A Wolfram Web Resource.
[5]Weisstein, Eric W. "Multifactorial." From MathWorld--A Wolfram Web Resource.