Servicios Personalizados
Revista
Articulo
Indicadores
- Citado por SciELO
- Accesos
Links relacionados
- Similares en SciELO
Compartir
Revista Boliviana de Física
versión On-line ISSN 1562-3823
Revista Boliviana de Física v.17 n.17 La Paz 2010
OBSERVACIÓN DE NEUTRONES SOLARES Y EFECTOS SOBRE EL CAMPO MAGNÉTICO TERRRESTRE ASOCIADOS CON LA FULGURACIÓN SOLAR DEL 2 DE NOVIEMBRE DE 2003
OBSERVATION OF SOLAR NEUTRONS AND THE EFFECTS ON THE EARTHS MAGNETIC FIELD ASSOCIATEDWITH THE SOLAR FLARE OF THE 2ndNOVEMBER 2003
Chamani W.1, Ticona R.2 y Ricaldi E.1
1 Carrera de Física
2 Instituto de Investigaciones Físicas
Universidad Mayor de San Andrés
La Paz-Bolivia
RESUMEN
El 2 de noviembre de 2003 a las 17:03:00 (GMT) se produjo una fulguración solar que fue verificada por el satélite GOES y se registró por el Monitor de Neutrones (NM-64) a las 17:17:24 (GMT) del laboratorio de Física Cósmica en el monte Chacaltaya. En particular, se observó que la fulguración produjo la variación del campo geomagnético local (F); se produjeron variaciones de la intensidad de la componente horizontal (H), de la declinación (D) y de la componente vertical (Z). Se observaron los registros de los incrementos con respecto a los valores promedio (registrados diariamente) de dichas componentes estimados en 71 nT, 10 grados y \9nT respectivamente en el Observatorio Geomagnético de Patacamaya, así como el incremento del campo neto F estimado en 66 nT registrado por el Observatorio Geomagnético de Villa Remedios. Las duraciones de la señal de neutrones y de la señal magnética fueron aproximadamente de 40 y 54 minutos respectivamente, lo que muestra que la fulguración fue intensa y la emisión de neutrones fue de larga duración.
Descriptores: fulguración solar neutrones solares crochets magnéticos
ABSTRACT
On the 2ndNovember 2003 at 17:03:00 (GMT) there was a solar fiare recorded by the GOES satellite. The event was also registered at 17:17:24 (GMT) by the Neutrón Monitor (NM-64) of the Cosmic Physics Laboratory located at Mount Chacaltaya in Bolivia. The fiare produced changes in the intensity of the local geomagnetic field (F) affecting all three vector components horizontal component (H), declination (D) and the vertical component (Z). Records of the increase in relation to average valúes of the aforementioned components were observed at 71 nT, 10 degrees and \9nT at the Geomagnetic Observatory at Patacamaya (Bolivian Altiplano) and the net increase of the F field estimated at 66nT registered by the Geomagnetic Observatory at Villa Remedios. The duration of the neutrón and magnetic signáis were between 40 and 54 minutes approximately. This demónstrates that the fiare was intense and that the emission of neutrons took place over a long period of time.
Subject headings: solar fiares solar neutrons magnetic crochet
1. INTRODUCCIÓN
El estudio y la observación de las fulguraciones solares y los efectos que produce sobre el campo geomagnético han merecido la atención de muchos científicos sobre todo a partir de las observaciones de Carrington y Hodgson en 1859, cuando encontraron una interrelación entre un crochet magnético y una fulguración solar. Las fulguraciones son los eventos más violentos y energéticos en el sistema solar, los cuales están ñsicamente relacionados con las manchas solares; pueden acelerar iones a altas energías. Estos eventos se producen cuando la energía magnética acumulada en la atmósfera solar, cercana a una mancha, es repentinamente liberada. Seguidamente, generan ondas electromagnéticas a lo ancho de todo el espectro, desde ondas de radio hasta frecuencias de rayos gamma, que son producidas por emisión sincrotónica en el campo magnético solar de la corona. Mientras los rayos X (RX) son generados en la atmósfera solar mediante el fenómeno del Bremsstrahlung, los neutrones son producidos por la interacción de iones acelerados con la atmósfera solar. Este tipo de eventos pueden ser observados por distintos métodos como: 1) la observación de la emisión con espectro-Heliógrafos, 2) la observación del flujo de rayos X con satélites y también pueden ser reconocidos y observados indirectamente sobre los efectos que producen en el campo geomagnético, ya que ocasionan variacio-
FlG. 1. Incremento del flujo de rayos X detectado por el satélite GOES entre el 1 y el 3 de Noviembre. El tiempo de inicio de la fulguración X8.3 fué a las 17:03 UT. La línea roja expresa el flujo de rayos X en el rango de longitud de onda de 1.0 8.QÁ y la azul para el rango de 0.5 41 (color online).
nes y perturbaciones en sus tres componentes: la horizontal (H), la vertical (Z) y la de declinación (D), debido a que un frente de partículas cargadas provenientes de la fulguración chocan con la magnetósfera de la Tierra, produciendo su deformación y compresión, seguidas de tormentas magnéticas en la atmósfera terrestre con un retraso de hasta algunos días. Otro de los efectos que producen son los crochets magnéticos que son variaciones magnéticas rápidas debidas a una desviación del sistema de corriente ionosférica, donde la variación geomagnética en el hemisferio iluminado por el sol se intensifica por un corto tiempo; el efecto es más grande cerca del punto subsolar (zenit). Por otro lado, los neutrones solares producidos son inmunes a los campos magnéticos, lo que hace posible que puedan alcanzar la superficie de la Tierra y por lo tanto nos permitan estudiar los mecanismos por los cuales han sido acelerados. Debido al interés por estudiar y observar este tipo de fenómenos que produce el sol, en este trabajo presentamos resultados preliminares del evento de una fulguración solar producida el 2 de noviembre de 2003 y reportamos los efectos que ha tenido sobre el campo magnético de la Tierra. Se observa una correlación del evento entre los datos del experimento del Monitor de Neutrones 12NM-64 de Chacaltayay dos observatorios geo-magnéticos, ubicados en Patacamaya y en Villa Remedios, de la Universidad Mayor de San Andrés (UMSA).
2. LA PROPAGACIÓN DE LOS NEUTRONES SOLARES
Los neutrones producidos en una fulguración solar no son afectados por el campo magnético del sol ni por el campo magnético terrestre; pero dependiendo de su energía logran alcanzar la superficie de la Tierra antes de que decaigan durante su viaje. Por ello la probabilidad de sobrevivencia, P(E), de un neutrón es:
Donde t es su tiempo de vuelo entre el Sol y la Tierra, yes el factor de Lorentz y t = 8865, el tiempo de vida del neutrón.
FlG. 2. El contaje registrado cada 5 minutos por el Monitor de Neutrones (12-NM64) localizado en el monte Chacaltaya, Bolivia el 2 de Noviembre de 2003. Los datos después de las 18:00 UT no fueron registrados debido a un corte de electricidad. El tiempo (inicio de la fulguración) registrado fue a horas 17:17:24 UT el cual correlaciona con el tiempo registrado por el satélite GOES. Se observa que la emisión neutrónica solar intensa máxima fue a horas 17:30:01 UT.
3. MÉTODO DEL TIEMPO DE VUELO
Los neutrones solares no pueden viajar a la velocidad de la luz debido a que tienen masa. El tiempo de su llegada a la Tierra es más largo para neutrones de bajas energías y si existe un amplio rango de energía de los neutrones producidos al mismo tiempo en el Sol, entonces existe también un amplio rango de tiempos de llegada. El tiempo de retraso Át, está dado por:
Donde R es la distancia entre el Sol y la Tierra (\UA = 1.496 x \0um), c es la velocidad de la luz y v la velocidad del neutrón. Combinando la ecuación (2) con la (3) se obtiene:
Entonces la energía cinética del neutrón Tn es:
Donde M es la energía de la masa en reposo del neutrón (940 MeV).
4. ARREGLO EXPERIMENTAL
4.1. Monitor de Neutrones NM-64
El monitor de neutrones es un instrumento para monitorear la intensidad de los rayos cósmicos galácticos y solares. Consiste principalmente de tres partes: el reflector, el productor y el contador. El monitor MN-64 usa polietileno como reflector. El rol del reflector es el de rechazar los neutrones de bajas energías producidos en la atmósfera y en las sustancias cercanas al monitor de neutrones. El plomo que rodea al reflector es usado como productor. Un contador de BF3
<
Fig. 3. Incremento de la componente horizontal H del campo magnético de la Tierra, el tiempo de inicio fue a horas 17:43 UT. La intensidad de campo magnético del crochet fue aproximadamente 71 nT. La tormenta magnética se identifica como un cambio abrupto de la componente H al tercer día después de la fulguración.
Fig. 4. Decremento de la componente vertical Z del campo magnético de la Tierra que se inicia a horas 17:35 UT. La intensidad de dicha componente file de 19nT.
está localizado dentro del productor. El moderador interno es usado para desacelerar neutrones hasta energías térmicas, el cual está localizado entre el contador y el productor.
4.2. Observatorio Geomagnético de Patacamaya
En este laboratorio se cuenta con un arreglo óptico-mecánico (Variógrafos La Cour), un teodolito magnético y un magnetómetro de protones. El variógrafo toma los datos de las variaciones diarias de los elementos H, D, Z del campo magnético terrestre, siendo la más sensible la componente H en latitudes muy cercanas al ecuador y menos sensible la componente Z.
4.3. Observatorio Geomagnético de Villa Remedios
Este laboratorio cuenta con un arreglo de bobinas que trabajan a modo Flux Gate, modelo FGE, de tres componentes X, Y, ZyF, proveniente de un magnetómetro fijo Overhauser. Los datos son registrados por un Data Logger. Los registros de este observatorio son rígidos y están sujetos a un mínimo
Fig. 5. Incremento de la componente de declinación D del campo magnético de la Tierra registrada a horas 17:51. La intensidad de dicha componente fue de 10 grados.
Fig. 6. Variación del campo magnético neto F que se inicia a horas 17:16 UT. La intensidad del campo neto F fue de aproximadamente 66nT.
de fuentes de error, siendo la variación de la temperatura del recinto del registro la más significativa.
5. REGISTROS Y OBSERVACIONES
Una fulguración solar intensa ocurrió el 2 de noviembre de 2003. La emisión del flujo de rayos X blandos observados por el satélite GOES (Fig. 1), comenzó a incrementarse a las 17:03 UT, alcanzó su máximo a las 17:25 UT y decayó a la mitad de su valor máximo a las 17:39 UT. El pico del flujo de RX fue clasificado como X8.3. La mancha solar corresponde a la región 10486, ubicada en las coordenadas S\4°, W56°. En la misma región se produjo el evento solar de neutrones el 28 de Octubre de 2003. La Fig. 2 muestra el contaje registrado cada 5 minutos por el Monitor de Neutrones NM-64. La emisión neutrónica comenzó a incrementarse a las 17:17:24 UT, la señal intensa alcanzó su máximo a las 17:30:01 UT, el que correlaciona con el tiempo registrado por el GOES. La duración del evento fue de aproximadamente de 40 minutos. Las Figs. 4, 5 y 6 muestran el efecto que produce la fulguración solar sobre el campo magnético de la Tierra. Se producen variaciones respecto de los valores,
registrados diariamente, de las componentes horizontal (H), vertical (Z) y de declinación (D) del campo geomagnético por el Observatorio Geomagnético de Patacamaya. La Fig. 7 muestra la variación del campo geomagnético neto (F) registrado por el Observatorio Geomagnético de Villa Remedios. La Fig. 4 muestra una anomalía notable (crochet) en la componente H del campo, que alcanzó su máximo a las 18:04 UT con una intensidad de 71 nT. A partir de los datos que se leen de este registro, la duración de la señal magnética fue de 54 minutos, por otro lado, la nube de plasma solar llegó aproximadamente un día y medio después de la fulguración. La tormenta magnética se identifica como un cambio abrupto de la componente H al tercer día después de la fulguración, en este mismo lapso de tiempo ocurre la fase de compresión y deformación de la magnetósfera terrestre. La Fig. 5 muestra el decremento de la intensidad de la componente vertical Z del campo estimada en 19 nT que se inició a horas 17:35 UT. Por el contrario, en este caso se observa un incremento de la componente de declinación D del campo estimado en 10 grados, que es el ángulo formado entre el polo norte geográfico y la componente H del campo.
6. CONCLUSIONES
Los resultados obtenidos en este estudio permiten concluir que la fulguración solar fue intensa, de larga duración y de emisión de neutrones parecida al evento registrado el 7 de septiembre de 2005 (Ricaldi 2007). Además se ha visto que el efecto que ha producido ha sido el cambio drástico de la intensidad del campo geomagnético ocasionando los llamados crochets magnéticos o anomalías magnéticas que producen la variación diaria solar (Sq), el cual se manifestó significativamente en el incremento de la componente H, lo cual es de esperarse para latitudes bajas como la de Chacaltaya, cercanas al ecuador magnético. En el instante que ocurrió la fulguración solar, el Sol estaba situado casi en el cénit sobre el observatorio de Chacaltaya (292.0°E, \6.2°S, 5250m.s.n.m.), por lo tanto se encontró en una muy buena posición para la observación de los neutrones solares. El rango de energías de los neutrones solares obtenido mediante el método del tiempo de vuelo estuvo entre 50 y \94MeV, con velocidades relativistas entre 0.31cy 0.56 c. Por otro lado, se han estimado los incrementos y decrementos del campo neto F, las componentes D y Z del mismo en 66nT, 10 grados y \9nT, respectivamente. Se ha visto que la tormenta geomagnética se produce al tercer día después de la fulguración, la cual desprendió una nube de plasma que tardó 1 d 12h y 46 m inutos en llegar a la Tierra, lo hizo a una velocidad de casi 1133 km/s y tuvo una extensión de 13.7x 106fcm. Este evento produjo la disminución en la intensidad del campo magnético neto F de la Tierra por un intervalo de tiempo de 3 horas 21 minutos en el observatorio geomagnético de Villa Remedios. Por último, los resultados presentados son preliminares ya que falta observar el registro del evento por el Detector de Neutrones Solares y así poder confirmar o rechazar alguna correlación con los datos del Monitor de Neutrones.
REFERENCIAS
1.- Space Environment Center-N.O.A.A. 2010, http ://www. sec .noaa.gov/today.html .xrav. [ Links ]
2.- Rastogi, R. G. et al. 1999, PhysEarth PlamísSpace, 51, 947. [ Links ]
3.- Ricaldi, E. 2007, Observación simultánea de neutrones solares en asociación con una fulguración solar del 7 de septiembre de 2005 (HF-UMSA). [ Links ]
4.-Ticona, R 2004, El Monitor de Neutrones (IIF-UMSA). [ Links ]
5.- Watanabe, K. 2005, Solar Neutrón Events Associated with Solar Fiares (PhD Thesis). [ Links ]
6.- Zirin, H. 1965, Solar Fiares: Observations, Chapter 6 of Stellar and Solar Magnetic Fields (IAU Symposium No 22) [ Links ]^rND^sRastogi^nR. G.^rND^nPedro^sMiranda Loza^rND^nFabricio G.^sÁvila Elías^rND^nPedro^sMiranda Loza^rND^nFabricio G.^sÁvila Elías^rND^nPedro^sMiranda Loza^rND^nFabricio G^sÁvila Elías
DISEÑO DE UN OSCILOSCOPIO MULTICANAL CON FPGA (PROYECTO LAGO)
THE DESIGN AND CONSTRUCTION OF A FPGA MULTI CHANNEL OSCILOSCOPE (LAGO PROJECT)
Pedro Miranda Loza1 & Fabricio G. Ávila Elías2
1 Instituto de Investigaciones Físicas-UMSA
2 Carrera de Electrónica, Facultad de Ingeniería-UMSA
RESUMEN
Se diseñaron dos instrumentos basados en FPGA (Field Programmable Gate Array) para el proyecto LAGO (Large Aperture GRB Observatory) que facilitan la observación de señales eléctricas producidas por variaciones de rayos gamma de alta energía. Estos instrumentos son: osciloscopio de almacenamiento digital básico de alta velocidad con el que se puede observar las variaciones de las señales eléctricas en el tiempo, analizador multicanal que aporta datos cuantitativos acerca de la frecuencia con que las señales eléctricas alcanzan determinados niveles de voltaje. Ambos instrumentos se implementaron en hardware digital (re)programable de alta velocidad, utilizando Lenguaje de Descripción de Hardware (HDL). El despliegue y almacenamiento de los datos se realiza por un ordenador que facilita el acceso a funciones de control de los instrumentos. Para esto se desarrolló un software propio en visual C++. Se ha verificado el funcionamiento del diseño final conectado a un detector real en el Laboratorio de Física Cósmica en el monte Chacaltaya.
Descriptores: sistemas computacionales instrumentos para rayos gamma
ABSTRACT
Two instruments based on FPGA (Field Programmable Gate Array) were designed and constructed for the LAGO (Large Aperture GRB Observatory) project. The instruments allow the observation of electric signáis producedby variations in high energy gamma rays. The instruments developed were: a high velocity basic digital storage oscilloscope which enables the observation of variations in electrical signáis over time, a multi channel analyzer that supports quantative data and picks up on the frequency at which electrical signáis reach certain voltages. Both instruments opérate with Language Description Hardware, a high velocity digital and reprogrammable hardware. The collection and storage of data is performed using a computer that facilítates access to the instruments' control functions. For this purpose original software was developed in C++. The design and function of the instruments were tested by connecting them to a real detector at the Cosmic Physics Laboratory at Chacaltaya, Bolivia.
Subjectheadings: computer systems gamma ray instruments
1. INTRODUCCIÓN
En los experimentos de ñsica de partículas, las herramientas de adquisición y tratamiento de datos están en constante evolución. Son varios los factores que favorecen ésto, desde el aumento de la frecuencia de trabajo de los conversores análogo-digitales hasta la disponibilidad de realizar prototipos basados en FPGA (Field Programmable Gate Array), haciendo posible realizar un circuito integrado a medida, sin los riesgos económicos asociados a las otras opciones tecnológicas. Aunque existen módulos especializados, disponibles comercialmente, éstos son de un costo mucho mayor al desarrollo de sistemas específicos.
El diseño de instrumentos de nueva generación con FPGA, está destinado a su utilización en el proyecto LAGO (Large Aperture GRB Observatory)1, el cual está implementado para la observación de las variaciones de rayos gamma, Gamma Ray Burst (GRB), de alta energía. (Laboratorio de Física Cósmica-Chacaltaya)2.
Las señales analógicas producidas por los detectores de rayos cósmicos de alta energía, tienen características especiales tales como: la duración del pulso, el valor máximo de éste, el tiempo promedio entre pulsos; siendo posible estudiar éstas señales mediante un Osciloscopio de alta velocidad. También es de interés el número de pulsos que se han producido en un
1 Proyecto LAGO,
http://partículas.enea.gov.ar/experiments/lago.
2 Proyecto LAGO-Bolivia,
http://www.f iumsa.edu.bo/iif/lago/inf ormes.
Pedro Miranda Loza & Fabricio G. Ávila Elías
Fig. 1. Diagrama general de conexiones.
Fig. 2. Diagrama funcional del diseño general.
intervalo de tiempo; esta tarea es realizada por un Analizador Multicanal, el cual tiene el propósito de contar dichos pulsos, además de clasificarlos por amplitudes.
Los detectores de Rayos Cósmicos (RC) del proyecto LAGO entregan una señal eléctrica de corta duración y amplitud proporcional a la energía de la partícula. Para este efecto el detector está formado por las siguientes partes: Detector Cherenkov con agua (Water Cherenkov Detector-WCD), fotomultiplicador (PMT) y alimentación de alto voltaje (HV).
Una de las principales características de la tecnología FPGA es la (re)programación del hardware digital por el usuario. Por lo que luego es posible diseñar nuevos instrumentos con una inversión mínima. Esto permite que un solo circuito digital en PCB (Tarjeta de Circuito Impreso) pueda desempeñar diferentes funciones según lo requiera la ocasión. Por ejemplo, en la implementación de los instrumentos descritos más arriba, podría programarse el hardware como sólo Osciloscopio, o como Osciloscopio y Analizador Multicanal.
La velocidad en la captación y el proceso de datos, hacen del FPGA un dispositivo ideal para el experimento realizado; en vez de realizar la codificación de algoritmos en un microprocesador, éstos se mapean directamente en hardware. Esta última característica permite ejecutar instrucciones en paralelo; logrando así que dos instrumentos, Osciloscopio y Analizador Multicanal, procesen datos al mismo tiempo.
El costo económico estimado en la implementación del prototipo, comparado con el costo de una solución comercial es cinco veces menor. Por ello el Instituto de Investigaciones Físicas de la U.M.S.A. creyó que este proyecto se justificaría por su menor costo y aporte a la investigación en física cósmica.
Se elaboraron las GUFs (Graphic User Interface) en el ordenador para la visualización de las señales detectadas por los instrumentos. Además de proveer las funciones gráficas, el ordenador provee la mayoría de las funciones de control.
2. HARDWARE
En la Fig. 1 se muestra el diagrama general de conexiones entre el Detector, el Inversor de Alta Velocidad, el Dispositivo de Adquisición de Datos y el PC. El Inversor de Alta Velocidad es utilizado para invertir el pulso negativo generado por el Detector, está construido con elementos discretos y tiene un ancho de banda de 300MHz con una impedancia de entrada de 1MCX
La señal (negativa) proveniente del Detector (fotomultiplicador), es conducida a través de cable coaxial hasta el Inversor de Alta Velocidad, luego la señal (positiva) de salida del inversor es introducida al Dispositivo de Adquisición de Datos, éste digitaliza las señales analógicas y transmite los datos al PC, a través de una interfaz serial RS-232 con conectares DB9. Para la interfaz serial, es utilizado un cable de nueve hilos no cruzado (tipo módem).
El Dispositivo de Adquisición de Datos está conformado por una tarjeta de digitalización que contiene un conver-sor análogo-digital de 100MSPS (ADC08100-NationalInstruments), y una tarjeta de desarrollo con el FPGA ACEX EP1K10TC100 (Altera).
Los bloques de la Fig. 1 a implementarse en el FPGA son: la memoria RAM donde se guardarán las muestras, la lógica de control y los módulos de recepción y transmisión de datos del puerto serial.
FlG. 3. Ventanas de la interfaz gráfica.
3. SOFTWARE
3.1. Diseño general
Tomando en cuenta la capacidad de los FPGAs de ejecutar instrucciones en paralelo, se han implementado el Oscilos-copio de Almacenamiento Digital y el Analizador Multicanal en un mismo diseño escrito en lenguaje de descripción de hardware (HDL). Ya que los dos instrumentos ejecutan conjuntos de procesos por separado, y estos procesos acceden a recursos compartidos, se ha visto de gran utilidad describir el diseño mediante diagramas de acceso, éstos son diagramas de interconexión funcionales para diseños concurrentes. En la Fig. 2 se tiene la representación funcional del diseño general, donde los componentes activos o procesos son presentados como círculos y los componentes pasivos como rectángulos, las flechas representan derechos de acceso o conexiones.
Los módulos OSC y PMCA corresponden a los dos instrumentos, Osciloscopio y Analizador Multicanal de Pulsos, respectivamente. Ambos módulos agrupan procesos y realizan el control de bancos de memoria RAM propios. También se puede observar que los procesos representados por dos círculos concéntricos son procesos que responden requerimientos de interrupción (IRQ), que deben entenderse desde el punto de vista del dispositivo: éste se encuentra realizando un proceso y a partir de un evento se pone en marcha otro proceso.
Ambos módulos utilizan recursos comunes, el mismo canal de comunicación con el ordenador(puerto serial) para enviar datos y recibir comandos, y los datos entregados por el conversor análogo digital; éste se encuentra configurado en modo libre, entregando el resultado de una conversión cada 10 ns.
El proceso FADC realiza un seguimiento continuo de la señal digital de entrada, una vez que esta señal cumple ciertas condiciones, éste proceso avisa al módulo OSC que comience la grabación de los datos digitales en posiciones consecutivas de la memoria RAM. Una vez que se han grabado todas las posiciones de la memoria disponibles, el módulo OSC accede al proceso Rx/Tx, y envía el contenido total de la memoria RAM al ordenador a través de la interfaz serial.
De manera similar ocurre el funcionamiento del módulo PMCA, el proceso FADC identifica el máximo (digital) de un pulso, con este dato accede al módulo PMCA, éste a su vez actualiza la cuenta del respectivo canal (en la memoria del espectro) y devuelve el control al proceso FADC. Si el módulo PMCA recibe el comando enviar todos los datos de la memoria de espectro al ordenador, se detiene el proceso de actualización del PMCA, se accede al proceso Rx/Tx y se envía el contenido total de la memoria de espectro al ordenador a través de la interfaz serial.
Las opciones implementadas para el instrumento Osciloscopio son las siguientes: velocidad de muestreo, nivel de dis-
FlG. 4. Curva de platean para PMT 9353KB.
paro, pendiente de disparo, predisparo y autodisparo. Las opciones implementadas para el instrumento Analizador Mul-ticanal son: nivel de discriminación o umbral, control de inicio/fin, borrado de toda la memoria de espectro, envió de datos al ordenador.
3.2. ínter faz gráfica
El software escrito en visual C++ realiza las siguientes funciones:
Acceso al puerto de comunicación serial del ordenador para adquirir los datos.
Proceso y despliegue gráfico de los datos capturados por el instrumento.
Operación de las funciones del instrumento a través de una interfaz gráfica de usuario.
El acceso al puerto serial se realizó a través de funciones de API de Windows, que facilitó el manejo de buffers de datos de envío y recepción y la modificación de los parámetros de transmisión. Teniendo en cuenta la gran utilidad de visualizar los datos capturados por el instrumento, el despliegue de gráficos se realizó utilizando librerías OpenGL3, que facilitan la representación de datos y su animación en dos dimensiones, logrando de esta manera una representación adecuada para los instrumentos Osciloscopio y Analizador Multicanal.
Un ejemplo del despliegue de datos en la interfaz gráfica se puede observar en la Fig. 3. En la parte superior está la ventana correspondiente al instrumento Osciloscopio con un pulso generado por el detector de Rayos Cósmicos, en la parte inferior está la ventana del Analizador Multicanal de pulsos con un espectro característico de Rayos Cósmicos y en la parte derecha está el menú de opciones de los instrumentos.
4. RESULTADOS
Una vez efectuadas varias pruebas con un generador de funciones, se verificó que la interfaz gráfica en la PC muestra correctamente los datos captados por el instrumento Osciloscopio. Los controles proveídos en el diseño hacen posible que se puedan observar señales en el rango de frecuencias de IKHz hasta 50MHz; el límite superior es dado por el
3 OpeGL UserInterfaceLibrary, http://glui.sourceforge.net.
teorema de Nyquist ya que el ADC utilizado trabaja a una frecuencia de 100MHz. Aunque en la práctica, para obtener una forma de onda más detallada es necesario tener por lo menos nueve muestras, en vez de dos muestras por periodo. Esto hace que la máxima frecuencia, presente en la señal a analizar debe ser de jq la frecuencia de muestreo. Los pulsos producidos por un PMT tienen una frecuencia aproximada de 40MHz, estos pulsos son detectados por el instrumento Osciloscopio (dos o tres muestras por pulso) y registrados en un archivo para posterior análisis.
Se aplicó el instrumento Analizador Multicanal de pulsos a la determinación de la curva de platean para un fotomulti-plicador9353KB deElectron Tubes4. Elfotomultiplicadorse encuentra instalado conjuntamente a un centelleador de agua, denominado WCD3, en el Laboratorio de Física Cósmica, en el monte Chacaltaya.
Se polarizó con alto voltaje (HV) el fotomultiplicadory se realizó el conteo de pulsos durante un tiempo fijo. Se polarizó el PMT de 900 a 1700 V, en intervalos de 100 V.
Los resultados de la aplicación se observan en la Fig. 4 en escala lineal y logarítmica. En esta ultima se puede apreciar mejor la región en la que un aumento del alto voltaje causa una menor diferencia entre las cuentas. Entre 1100 V y 1400 V está un segmento con menor pendiente; por lo tanto en ese segmento, a diferencia de otros, las cuentas no aumentarán significativamente conforme vaya aumentando el alto voltaje. Según la Fig. 4, el valor sugerido para polarizar el PMT es de 1200 V, y éste es el valor5 elegido para la operación de WCD3 en Chacaltaya.
Además de las cuentas totales de pulsos, también se dispone del espectro de picos, correspondiente a cada alto voltaje aplicado al PMT de WCD3. Se han comparado estos espectros con las mediciones independientes realizadas por el proyecto LAGO para WCD3, y se ha verificado que los espectros son similares para cada alto voltaje.
Los instrumentos descritos han sido diseñados para su uso específico en la observación de señales eléctricas rápidas, producidas por rayos cósmicos, en general de alta energía,
4 Electrón Tubes Enterprises Ltd., //www. electrontubes. com.
5 Este valor fue elegido después de varias pruebas con equipos especiali
zados instalados en Chacaltaya.
del proyecto LAGO, pero éstos también tienen un uso general. En el caso del Osciloscopio, puede utilizarse como un instrumento de laboratorio de nivel universitario, para la observación de señales periódicas y no periódicas, dada su alta tasa de conversión (100 MSPS). En el caso del Analizador Multicanal, éste es muy útil en análisis de espectroscopia nuclear, en el cual la distribución de probabilidad de la radiación detectada está en función de su energía. El Analizador Multicanal también puede utilizarse en aplicaciones en las que se requiera el conteo de pulsos.
REFERENCIAS
1.- Rivera, B. H. 2007, Implementación y Calibración de un Prototipo de Detector Cherenkov de Chubascos Atmosféricos (Tesis de Grado, Universidad Mayor de San Andrés) [ Links ]^rND^sRivera^nB. H.^rND^sJosé Antonio^nZelaya^rND^sJosé Antonio^nZelaya^rND^sJosé Antonio^nZelaya
DOS PASIONES: LA FÍSICA CÓSMICA Y LAS NIEVES ANDINAS
TWO PASSIONS: COSMIC PHYSICS AND ANDEAN SNOWY PEAKS
Gral. Ing. José Antonio Zelaya Academia Nacional de Ciencias de Bolivia
RESUMEN
De la historia del Laboratorio de Física Cósmica de Chacaltaya. El Ing. José Antonio Zelaya rememora su experiencia como investigador en el laboratorio en la época inicial cuando el Dr. Ismael Escobar, fundador del Laboratorio, era su Director.
Descriptores: historia de la ciencia observatorios rayos cósmicos
ABSTRACT
A history of the Cosmic Physics Laboratory at Chacaltaya. José Antonio Zelaya reminisces about his experiences as a researcher at the Laboratory during its beginnings when Dr. Ismael Escobar, founder of the Laboratory, was Director.
Subjectheadings: history of science observatories cosmic rays
"La vida atrae la vida, nada ocurre por azar" Javier Sierra.
La atracción de las nieves andinas y por el esquí, inspiraron al Profesor universitario Dr. Ismael Escobar, para promover las investigaciones de Física Cósmica en el nevado de Chacaltaya. Su entusiasmo consiguió despertar el interés del rectorado de la UMSA, y también el interés de otras universidades del exterior. Su plan inicial fue instalar una estación meteorológica y poco a poco las instalaciones. A partir de las primeras habitaciones fueron agregándose otras para recibir las expediciones extranjeras que solicitaban llegar a este laboratorio que poco a poco fue adquiriendo fama internacional.
Cuando en 1964, el Ing. Ramón Schulczewski iniciaba otro experimento internacional sobre estudios de la ionosfera, llegó la primera antena parabólica a Bolivia, fui invitado para orientar astronómicamente mediante un azimut geodésico. En dos noches, utilizando la estrella polar del Sur, la Sigma Octantis, logramos cumplir con el pedido. Nos pidieron la cuenta por el trabajo realizado y mi respuesta fue que se trataba de una colaboración a la UMSA, donde ya dictaba la cátedra de astronomía geodésica. Cuando pasaron algunos años recibí un llamado telefónico para visitar al Director del "Laboratorio de Física Cósmica" Dr. Escobar. En muy breves palabras me invitó para integrar la planta de investigadores del Laboratorio. Por mi grado militar de Mayor, le dije que tendría que solicitarse siguiendo el conducto regular mediante una nota al Comandante del Ejército. Así lo hizo y a la semana me volvió a llamar y me dio la nota de respuesta del Comando del Ejército, negando mi participación. Recuerdo sus palabras: "No entiendo, cómo una institución como el Ejército se niega conseguir el prestigio de participar en la investigación científica en el ya famoso Laboratorio de Física Cósmica de Chacaltaya, otra institución militar se sentiría halagada por esta invitación - Deje que yo aplique mis métodos." Yo levanté los hombros, dando así mi consentimiento para que ponga en práctica su idea. Llamó a su secretaria y pidió hablar con el Dr. Paz Estenssoro, Presidente de la República, a los pocos minutos escuché la conversación y me dijo "esperaremos la respuesta", unos minutos más y evidentemente recibió la respuesta. Y, volviéndose me dijo, "bienvenido Mayor como Investigador asociado, en comisión del Ejército cumpliendo las leyes de oficiales en comisión." Esta anécdota, explica la capacidad emprendedora del Dr. Escobar. Cuando llegué por primera vez a Chacaltaya se me asignó una habitación en el segundo piso donde ya estaban instalados varios equipos de investigación; dormí en el piso porque no habían catres. Al comedor llegábamos en tres turnos, donde se reunían americanos, japoneses e investigadores bolivianos. Durante un año realicé en el nevado, a 5.420 metros de altura, las investigaciones sobre las "Nubes de libración en las constelaciones opuestas a la trayectoria de la Luna para evitar interferencias en las placas de emulsión en la cámara con telescopio incluido. El Dr. Stevenson, del Servicio Geológico de los Estados Unidos, me dio las instrucciones para las operaciones en la obtención de las placas con exposiciones hasta de 10 minutos escogiendo una estrella de seguimiento, en los mapas estelares en los lugares aproximados a los puntos de libración de la Luna. En dos noches de trabajo obtuvimos 8 placas, las que procesamos en la cámara oscura del laboratorio. Satisfecho con mi trabajo, el Dr Stevenson, me dejó que realizara solo todo el proceso, sin importarme la fatigosa ascensión al nevado, ni el frío (23° bajo cero) y los fuertes vientos que había que soportar. La cúpula plateada que adornaba la cumbre de la pista de esquí, presentó un inconveniente, reflejaba el brillo de algunas estrellas en la sensible emulsión cuando el obturador estaba abierto. Recibí una carta para que detectara el inconveniente y, al comprobar que se trataba de la cúpula, avisé al Dr. Escobar y cuando me preguntó la solución le dije que había que construir otro mojón de soporte para la cámara a la intemperie algo alejada de la cúpula. Muy rápidamente se realizó el traslado y las placas obtenidas ya no presentaban interferencias, pero la falta de protección aumentaba la incomodidad del frío y el viento. Cuando la motivación del trabajo llegó a despertar mi curiosidad personal, aumenté las exposiciones hasta 20 minutos y luego hasta 30 minutos. Las placas así obtenidas me mostraron la inmensidad del misterioso cosmos. El motivo principal de la investigación era descubrir las rocas en suspensión o equilibrio gravitacional llamadas las rocas de Kordalewsky un astrónomo polaco que escribió su teoría de las nubes de libración en las que se encontraban rocas de hasta 200m de longitud, provenientes de las colisiones de meteoros o meteoritos en la superficie lunar, las mismas que dejaron las huellas como cráteres volcánicos. La teoría agregaba que de las altas montañas estas rocas incluso se las podría ver a simple vista. Cuando cumplí la fecha de mi compromiso fui preguntado por carta si había visto las rocas a simple vista, mi respuesta fue clara, siguiendo el dictado que me enseñaron en la EMI. La ciencia es la búsqueda de la verdad. Y dije, "no, no he visto". Al poco tiempo el hombre llegó a la Luna. Recibí una nota de agradecimiento por mi trabajo, indicándome que en una reunión de científicos en el Servicio Geológico, había ocupado el segundo lugar entre los trabajos científicos más sacrificados que se hacían en el mundo, el primer lugar lo ocupó un investigador de los peces ciegos en Laponia donde se cavaron pozos en el hielo y los investigadores zambullían para agarrar los peces con las manos y subirlos a superficie.
Durante mi estadía en Chacaltaya logré la donación de 20 catres del Ministerio de Defensa, fui nuevamente invitado para colaborar en la investigación del ozono en superficie, en el manejo de los telescopios para la detección de mesones pesados en interior mina y en el procesamiento de datos meteorológicos y paralelamente, mientras trabajaba en las noches, durante el día colaboré como guía para atender la gran cantidad de turistas que visitaban las instalaciones incomodando muchas veces la atención de los investigadores. El Dr. Escobar agradeció mis servicios y se despidió cuando fue distinguido para ocupar un cargo en el exterior. Guardo para él los más gratos recuerdos del mismo modo para los colegas extranjeros y bolivianos con los que tuve la oportunidad de iniciar una permanente amistad.
Los trabajos realizados en Chacaltaya fueron tomados en cuenta para ocupar la silla de Astronomía en la Academia de Ciencias, que había dejado vacante el padre De Cotte del observatorio del Colegio San Calixto. Con el tiempo, ocupé la dirección del Instituto de Investigaciones Físicas de la UMSA y logré que la sección de Construcciones de la Universidad aprobara los planos arquitectónicos para un nuevo laboratorio en Chacaltaya. Estos planos deben estar en los archivos correspondientes.
XXI REUNIÓN NACIONAL DE LA SOCIEDAD BOLIVIANA DE FÍSICA
DEL 26 AL 31 DE OCTUBRE DE 2009
COCHABAMBA-BOLIVIA
Sociedad Boliviana de Física
RESUMEN
Se presenta la relación de ponencias de la XXI Reunión Nacional de Física realizada en Cochabamba entre el 26 y el 31 de Octubre de 2009.
Subject headings: physics workshops
Título:
EVALUACIÓN DEL MODELO CLIMÁTICO REGIONAL PRECIS PARA EL ÁREA DE BO-LIVIA: COMPARACIÓN CON ESTACIONES DE SUPERFICIE
Autor:
Marcos Froilán Andrade Flores
Resumen:
Los posibles efectos del cambio climático sobre regiones como Bolivia son de mucho interés tanto a nivel científico como político. Una de las herramientas más importantes para estudiar estos efectos son los modelos numéricos del clima. A fin de establecer la confiabilidad de las proyecciones de la evolución del clima futuro producidas por dichos modelos es necesario, sin embargo, validar las salidas de los mismos. En este trabajo se estudia el comportamiento del modelo climático regional PRECIS (Providing REgional Climates for Impacts Studies) del Hadley Centre de la Oficina Meteorológica Inglesa, sobre la región de Boli-via. La comparación de resultados del modelo con datos observados en superficie a nivel mensual de temperatura y precipitación en el periodo 1961-1990 muestran que el modelo tiene un buen desempeño en zonas bajas (menos de 500 m snm) pero que sobreestima precipitación tanto en zonas de altura intermedia (entre 500 y 3500 m snm) como en la zona andina (situada a una altura mayor a 3500 m snm). A su vez la comparación muestra que el modelo subestima la temperatura (media, máxima y mínima) en la mismas dos regiones mientras que realiza un trabajo aceptable en las zonas bajas. Aunque la parametrización de nubes y precipitación en las zonas de los valles y andina parece ser el elemento que más podría contribuir a las diferencias observadas es necesario un estudio más profundo para entender los mecanismos involucrados.
Título:
HERRAMIENTAS TECNOLÓGICAS Y EL PROCESO DE ENSEÑANZA APRENDIZAJE EN AULA PARA LA MATERIA DE FÍSICA
Autor:
Milka Mónica Torrico Troche
Resumen:
La Física como ciencia es indispensable en la formación de los futuros profesionales de Ciencias e Ingenierías, considerando como población los estudiantes que cursan los primeros semestres de las carreras de Ciencias e Ingenierías de las Universidades publica y privadas de la provincia Cercado del Departamento de Cochabamba, las estadísticas nos muestran muchas deficiencias en la respuesta de los estudiantes hacia la materia de Física Básica, promoviendo un análisis de la situación que se centra en el proceso de enseñanza en aula, identificando a través de encuestas a docentes y estudiantes involucrados en el proceso, tres problemas, que son: (i) deficiencia de material bibliográfico adecuado; (ii) falta de visualización del fenómeno físico en los estudiantes; (iii) deficiencia en la comunicación docente - estudiante. De esta manera el presente trabajo tiene por objetivo: diseñar una guía metodológica que involucren el uso de nuevas herramientas tecnológicas, para mejorar el proceso de enseñanza aprendizaje en aula en la asignatura de Física Básica de las carreras de Ciencias e Ingenierías, tomando en cuenta los diferentes métodos didáctico y recursos didácticos, donde se incorpora el uso de herramientas computacionales en función a su aplicabi-lidad con el campo de la Física. Una vez definidos los métodos y recursos se planifica el desarrollo de las clases, estructurando las actividades, material, contenido, herramientas y recursos en un PLAN DE CLASES, que se constituye en una guía para el docente. El Plan de Clases lo elabora cada docente definiendo el rumbo que tomara la clase, de esta manera la introducción de las nuevas herramientas tecnológicas promueve interés en el estudiante que esta habituado a la nueva tecnología y la gama de alternativas computacionales permite que el docente se convierta en investigador en el proceso de enseñanza, recuérdese que enseñar es un arte donde el artista es el docente. Finalmente se muestra como ejemplo parte de la aplicación de la guía metodológica propuesta a la materia de Física II, tema Estática de Fluidos. generales y sobre la que no existe aún un concepto universal o una ciencia unificadora que valga para todos los ámbitos donde ella aparece. En el presente trabajo, el autor concibe una base teórica firme sobre la cual construye una teoría que explica el fenómeno de la complejidad de una manera coherente como un planteamiento científico inédito.
Título:
MAPA DE RADIACIÓN SOLAR PARA SOLIVIA
Autor:
Marcelo Javier Lucano
Resumen:
En este trabajo se estimaron los valores de la radiación solar global a través de modelos de Sistemas de Información Geográfica e imágenes satelitales obtenidos a través de la metodología SSE de la NASA para un periodo de registros de 22 años. Con estos valores de radiación, se obtuvo la distribución temporal y espacial de la radiación solar para Bolivia mediante un proceso de interpolación utilizando el paquete de análisis geoestadístico ILWIS, considerando datos de 182 grillas de Io x Io que abarcan la zona de estudio y utilizando el método de superficie móvil con distancia inversa, exponente de peso 2, distancia esférica y segundo grado parabólico. Se observó que los niveles de radiación más bajos corresponden a los meses de mayo a julio que pertenecen a la estación de invierno y los niveles de radiación más altos corresponden a los meses de octubre a diciembre que pertenecen a las estaciones de primavera y verano. Asimismo, se observó que la parte noreste de Bolivia, a partir de la Cordillera Oriental presenta los valores más bajos de media anual de radiación global, entre 3.9 - 5.1 kW-h/m2-día, en tanto que, en la zona suroeste la radiación solar aumenta con valores entre 5.1 - 7.2 kW-h/m 2-día.
Título:
LA TEORÍA DE LA ORGANIZATIVIDAD: 'EL FINAL DE LA TERMODINÁMICA CLÁSICA Y EL COMIENZO DE LA ORGANIZATIVIDAD' (UNA TEORÍA QUE EXPLICA LA COMPLEJIDAD)
Autor:
Andrés Trepp del Carpió
Resumen:
La complejidad es una cualidad abstracta de ciertos sistemas de la naturaleza y del mundo humano cultural, de la que hoy sólo se conoce sus características
Título:
BÚSQUEDA DE ANTIMATERIA EN EL UNIVERSO
Autor:
Kenny Choque Talavera
Resumen:
El modelo estándar de la física de partículas establece que tanto partículas como antipartículas fueron creadas simultáneamente en cantidades iguales. Según el modelo estándar del Big Bang, el Universo fue creado en un tiempo primordial. Es natural pensar que en el inicio del Universo, materia y antimateria fueron creadas en cantidades iguales. Sin embargo, las observaciones astronómicas muestran que el Universo observable está constituido por materia. La aparente asimetría entre materia y antimateria es un enigma fundamental de la ñsica de partículas y también de la astronomía. Si la antimateria existe en algún lugar del Universo, la búsqueda de señales de rayos gamma provenientes de su aniquilación, o aun, de rayos cósmicos de antipartículas, serían indicadores de la existencia de dominios de antimateria.
Título:
¿ES LA HORCA DEL INCA UN OBSERVATORIO ASTRONÓMICO?
Autor:
Gonzalo Pereira Quiroga
Resumen:
A propósito de la conmemoración del Año Internacional de la Astronomía, algunos estamentos medios del Gobierno de Bolivia, se han enfrascado en tramitar la declaratoria de "Patrimonio Astronómico de la Humanidad" ante la UNESCO, al monumento arqueológico de la Horca del Inca, situado en la localidad de Copa-cabana. Uno de sus principales argumentos para obtener la declaratoria de patrimonio astronómico es una serie de trabajos arqueoastronómicos y de datación atribuidos al Laboratorio de Física Cósmica de la Universidad Mayor de San Andrés. La ponencia que proponemos presentar es el resultado de un análisis de los trabajos, tanto arqueológicos como astronómicos, que se emplean para apoyar la propuesta: La Horca del Inca es un Observatorio Astronómico.
Título:
PRUEBAS DE ALINEALTDAD POR DATOS SUBROGADOS
Autor:
Y Miguel Peñafiel N.
Resumen:
Se analiza un registro de las tres componentes del campogeomagnético (Patacamaya) mediante datos subrogados de Fourier y subrogados restringidos para decidir sobre sus características de linealidad. Se compara los resultados obtenidos con los correspondientes coeficientes maximales de Liapunov.
Título:
MÉTODO DE DEPROYECCIÓN Y OBTENCIÓN DE LA ESTRUCTURA DE BRAZOS DE GALAXIAS ESPIRALES
Autor:
Deterlino Urzagasti
Resumen:
Se desarrolla un método para la deproyección de imágenes de discos de galaxias espirales sin barra y la obtención de su estructura de brazos mediante el modelo de espirales logarítmicas. El método se basa en la obtención de los puntos de máximos y mínimos locales de las fluctuaciones de la intensidad en las imágenes galácticas a través de un proceso de suavizado con el filtro de Savitzky-Golay. Estos puntos luego sirven para determinar los ángulos de inclinación y de posición del disco galáctico así como el pitch angle de sus brazos espirales. El método se aplicado a varias galaxias espirales, encontrándose un buen acuerdo con otros resultados de la literatura.
Título:
ANÁLISIS DE PATRONES EN FRACTALES
Autor:
Flavio Ghezzi
Resumen:
Se discute la dimensión frac tal de sistemas deter-minísticos y no-determinísticos haciendo énfasis en la relación entre la dimensión fractal y diferentes parámetros ñsicos. Finalmente, se discute la invarian-cia de escala de sistemas coloidales.
Título:
AUTOSIMILARLDADES EN EL ESPACIO DE PARÁMETROS PARA EL CmCUITO DE CHUA
Autores:
Marcelo Ramírez Ávila, Jason Gallas
Resumen:
Se muestran autosimilaridades en el espacio de parámetros del circuito de Chua, considerando tanto una función lineal por partes, así como una función cúbica. También se muestra el isomorfismo entre los planos de fase para ambas no-linealidades.
Título:
ASPECTOS TEÓRICOS DE LA TEORÍA ESPE CIAL DE LA RELATrVTDAD
Autor:
Luis Panfilo Yapu Quispe
Resumen:
Se describe la estructura del grupo de Lorentz que permite relacionar la teoría especial de la relatividad con otras áreas de la física moderna. Para ilustrar esta descripción se deducirán algunas fórmulas "clásicas".
Título:
OSCILACIONES DE BLOCH EN UN CIRCUITO CUÁNTICO INDUCTIVO ACOPLADO A LAS
VIBRACIONES DE UNA RED
Autor:
Marcelo Calcina Nogales
Resumen:
Se acopla un circuito puramente inductivo con carga discreta a los modos vibracionales de una red unidimensionales (fonones). Los modos vibracionales deforman el anillo y por lo tanto se modifica su inductan-cia (sistema electromecánico-cuántico). Cuando este anillo se conecta a una fuerza electromotriz de valor £o, el periodo usual de la corriente (qe £o/V) se ve afectado por un armónico con periodo igual al doble de su valor original (2qe £o/V). Este fenómeno es equivalente al de las oscilaciones de Bloch en un cristal en el cual el electrón experimenta una frecuencia de oscilación aFIW. En nuestro caso, la carga eléctrica fundamental qe juega el rol de la constante de red.
Título:
LOS COMETAS METEOROLÓGICOS
Autores:
Grupo Astronáutico Boliviano
Resumen:
Se estudiaron dos tipos de cometas de papel (delta y rombo) y su desempeño bajos las condiciones atmosféricas de La Paz, con el fin de diseñar cometas capaces de elevar un equipo de medición atmosférica. La gran preferencia de usar cometas en lugar de globos aerostáticos se basa en la facilidad de recuperar los instrumentos de medición así como en la economía para construirlos. Comprobamos que el método de unir varios cometas es adecuado poder elevar una carga de prueba y alcanzar una altura mayor sin perder estabilidad. Lo primero que se hace es elevar un cometa piloto que cuando está estable jala a un segundo cometa que lleva la carga de prueba, y así se siguen uniendo varios cometas hasta alcanzar alturas ideales de unos 250m.
Título:
PROYECTO "LAGO"
Autor:
Alfonso Ve larde Chave z
Resumen:
Se describe la instalación y puesta en marcha de la estación LAGO-Bolivia.
Título:
"SISTEMATIZANDO" LA INVESTIGACIÓN CIENTÍFICA
Autor:
Marco Antonio Córdova Olivares
Resumen:
Por mi experiencia a lo largo de mi formación académica, he visto "desvanecerse" trabajos de investigación tanto de docentes como de estudiantes, ya que caen en "saco roto". Ello se debe a que no existe una política institucional ni gubernamental que pueda "agarrar" estos trabajos de investigación. En este sentido, propongo un método que facilita la digitalización y el control de todos los trabajos de investigación científica, tanto a nivel personal como institucional, enfocando mi propia experiencia sobre la construcción de mi trabajo de tesis.
Título:
DOSIMETRÍA DE NEUTRONES - MONITOR DE NEUTRONES
Autor:
Rolando Ticona Peralta
Título:
ESTUDIO DE LA SINCRONIZACIÓN DE PÉNDULOS ACOPLADOS MEDIANTE SIMULACIONES COMPUTACIONALES
Autores:
Armando Reynaldo Ticona Bustillos, Marcelo Gonzalo Ramírez Ávila
Resumen:
Usando un modelo discreto de péndulo amortiguado, estudiamos las condiciones de sincronización para dos péndulos acoplados. En primer lugar analizamos las cuencas de atracción estudiando la razón entre los periodos de los dos péndulos (Ti/T{) para diferentes condiciones iniciales, es decir, diferentes ángulos desde los cuales los péndulos inician su movimiento, considerando para ésto diferentes relaciones de largo, masa, acoplamiento entre péndulos y disipación. A partir de estos resultados podemos caracterizar las lenguas de Arnold, analizando nuevamente el comportamiento de la relación entre periodos en función a las condiciones iniciales de uno de los péndulos y los diferentes valores de masa, largo, acoplamiento y disipación. En estos gráficos podemos identificar sincronización 1 : 1 así como sincronizaciones de otros órdenes (m : n) y también regiones cuasi periódicas, las cuales podrían conducirnos también a regiones caóticas en trabajos futuros; además se tienen todas las condiciones para ampliar este trabajo a más de dos péndulos.
Título:
CAMAS SOLARES
Autores:
Emiliana Juaniquina, Fernando Ve larde, Flavio Ghezzi, Armando Ticona
Resumen:
Se realiza la caracterización del espectro de emisión electromagnética de una lámpara fluorescente usada en camas solares, para esto usamos un espectro fotómetro, que nos permite analizar el espectro desde aproximadamente 180nm, hasta aproximadamente 880nm, dividiendo este intervalo en 3500 canales. Debido a que el espectro fotómetro adquiere los datos mediante una fibra óptica, las limitaciones de ésta no nos permiten analizar todo el espectro, dejándonos sólo información entre las bandas UVB y el visible. Las lámparas usadas en estas camas solares deberían tener una emisión rica en la banda del UVA, lo que permitiría el bronceado; y el resto de emisión debería estar en el visible. Nuestros resultados muestran un contenido pequeño de UVB en la emisión de dichas lámparas, lo cual representa un peligro debido a lo nocivo de este tipo de radiación, por lo cual se recomienda un uso muy cuidadoso y supervisado de dichos instrumentos cosmetológicos.
^rND^nA. R.^sTlCONA BUSTILLOS^rND^nG. M.^sRAMÍREZ ÁVILA^rND^nA. R.^sTlCONA BUSTILLOS^rND^nG. M.^sRAMÍREZ ÁVILA^rND^nA. R^sTlCONA BUSTILLOS^rND^nG. M^sRAMÍREZ ÁVILADETERMINACIÓN DEL COEFICIENTE DE DISIPACIÓN DE UN PÉNDULO
DETERMINATION OF A PENDULUM'S DISSIPATION COEFFICIENT
A. R. TlCONA BUSTILLOS & G. M. RAMÍREZ ÁVILA
Instituto de Investigaciones Físics
Universidad Mayor de San Andrés La Paz - Bolivia
RESUMEN
Usando un sensor de posición se toman datos de las oscilaciones de un péndulo amortiguado. Analizando las características del péndulo se determina su coeficiente de disipación.
Descriptores: procedimientos de laboratorio — dinámica de la partícula
ABSTRACT
Using a position gauge we register data from the oscillations of a damped pendulum. We analize the pendulums characteristics and determine its dissipation coefficient.
Subjectheadings: laboratory procedures — dynamics of the particle
El movimiento armónico simple (MAS) encuentra aplicaciones en muchas ramas de la física; sin embargo, la primera dificultad con la que se tropieza al comparar sus resultados con oscilaciones reales, es la disipación de energía mecánica debida a la interacción de estos sistemas con el medio viscoso que los rodea. Esto hace que los sistemas reales sólo obedezcan por un tiempo relativamente corto los resultados de las soluciones del MAS. El problema de introducir el término de disipación (que puede ser función de la velocidad, de la masa y forma del objeto, así como de la viscosidad del medio) en las ecuaciones está en que éstas ya no son sencillas pues pierden su carácter lineal y por tanto, la solución de las mismas no es inmediata y se deben recurrir a aproximaciones o directamente a la resolución numérica de las ecuaciones. En este trabajo, se presenta un ejemplo de cómo medir directamente este coeficiente mediante las oscilaciones de un péndulo cuya lenteja puede ser considerada como un objeto puntual, utilizando un sensor de movimiento que nos da la posici´on de la lenteja del péndulo en función del tiempo. Este trabajo surge como parte de la investigación realizada en el estudio de péndulos acoplados (Ticona Bustillos & Ramirez Avila 2008).
2. PÉNDULO CON DISIPACIÓN
Se considera un péndulo con masa puntual, como se muestra en la Fig. 1. Las fuerzas que actúan sobre este objeto son el peso y la fuerza de disipación por parte del aire. Aplicando la segunda ley de Newton en la dirección del movimiento, se tiene:Pero la aceleración se puede expresar en función del desplazamiento angular como:
Además, considerando desplazamientos angulares pequeños para los cuales se cumple senf9 « 9, la ecuación se reduce a:
que tiene una solución conocida y detallada en libros de física general (Halliday & Krane 1998; Serway 1997; Tipler 1995) dada por:
angular dadapor:
FIG. 1.— Diagrama de fuerzas y aceleración para el péndulo con disipación.
FIG. 4.— Ángulo del péndulo en función del tiempo.
La frecuencia angular es menor a la del péndulo sin amortiguamiento, debido a la interacción con el medio viscoso.
3. MEDIDA DEL COEFICIENTE DE DISIPACIÓN
Utilizamos un sensor de posición PASCO, el cual toma datos de la distancia de un objeto a partir de un nivel de referencia, como se muestra en la Fig. 2. Estos datos son enviados a una computadora. Se utilizó un péndulo de 53.0cm de largo y 150.Og de masa. La aceleración debida a la gravedad en la ciudad de La Paz se considera igual a 9.775m/s2. Los datos obtenidos se muestran en la Fig. 3. En estos datos, se puede observar que las distancias muy cercanas al detector no son obtenidas correctamente debido a la propia sensibilidad del detector; además, existen algunos casos en los que el péndulo no es detectado cuando está en la máxima distancia, dándonos el dato de la pared que está a mayor distancia.
Como se conoce el largo del péndulo, se pueden transformar estos datos a ángulos con lo que se obtiene el gráfico de la Fig. 4, donde sólo se consideran los datos a partir del movimiento del péndulo; en estos datos también se realizó un desplazamiento desde el nivel de referencia mostrado en la Fig. 3.
En la ecuación (4) se pueden considerar solamente los máximos de las oscilaciones, es decir, los puntos en lo cuales el coseno es igual a la unidad. En ese caso sólo nos quedamos con la parte exponencial del comportamiento.
De igual manera, esto puede ser realizado con los datos obtenidos con el detector, tanto en las distancias mínimas como en las máximas. Como mencionamos antes, las distancias menores tienen un problema debido al detector, por lo cual sólo consideramos las distancias máximas en cada oscilación. De esta manera, se puede obtener el ajuste exponencial mostrado en la Fig. 5.
Para obtener la Fig. 5, algunos de los puntos han sido eliminados, ya que quedaban fuera de la línea de ajuste, debido a las circunstancias ya explicadas. A partir de este ajuste podemos determinar que el valor del coeficiente de disipación es igual a 0,013 kg/s.
4. CONCLUSIONES
Este método nos brinda una forma muy sencilla de calcular el coeficiente de disipación de un objeto
relativamente pequeño, sin tomar en cuenta medidas de los factores geométricos ni la viscosidad del fluido.
FIG. 5.— Ajuste exponencial del ángulo en función del tiempo.
REFERENCIAS
1.-Halliday, R. & Krane. 1998, Física, vol. 1, cuarta ed. (CECSA).
[ Links ]2.- Serway, R. 1997, Física, vol. 1, cuarta ed. (McGraw-Hill).
[ Links ]3.- Ticona Bustillos, A. R. & Ramírez Avila, G. M. 2008, RBF, 14, 121.
[ Links ]4.-Tipler. 1995, Física (Reverte)
[ Links ]^rND^sHalliday^nR.^rND^sTicona Bustillos^nA. R.^rND^sRamírez Avila^nG. M.^rND^sChoque Saire^nM. P.^rND^sChoque Saire^nM. P.^rND^sChoque Saire^nM. PLA OCARINA DE ZANAHORIA A CARROT OCARINA
Choque Saire M. P.
Laboratorio de Física 121 Universidad Mayor de San Andrés
RESUMEN
En los instrumentos prehispánicos aerófonos el sonido es producido por la vibración de una columna de aire como es el caso de la quena, la zampona y la tarkha. Se suele olvidar, sin embargo, a otros instrumentos pertenecientes a esta misma familia, como es el caso de la ocarina, cuya afinación y fabricación es dificultosa. Este instrumento musical es de tamaño reducido, su sonido es agudo (523 - 1046 kHz) y se la encuentra en múltiples diseños. Su construcción es ardua debido a que habitualmente las ocarinas están hechas de cerámica o madera, con un centro hueco y orificios con áreas iguales o desiguales que son cubiertas para producir diferentes notas; esta dificultad se puede remediar usando un material dócil. Este material tiene que ser de forma alargada, cilindrica y de consistencia maciza para evitar una deformación posterior; además, tal material debe permitir labrar el instrumento cómodamente. Entre las posibles opciones disponibles se eligió a una zanahoria (aunque parezca inusual) pues reúne todos los requisitos para la construcción de la ocarina; el único inconveniente es su sensibilidad al paso del tiempo ya que su volumen se puede reducir por la deshidratación.
Descriptores: instrumentos musicales propiedades acústicas de sólidos
ABSTRACT
In pre Colonial wind instruments sound is produced through the vibration of an air column such as in the Quena, Zampona and Tarka. A less well known wind instrument of the same family is the Ocarina whose fabrication and tuning are particularly difficult. The Ocarina is a small wind instrument with a high pitch (523 - 1046 kHz) and is found in many forms and designs. The instrument is arduous to make given that it is usually crafted from materials such as ceramic or wood and is composed of a hollow centre and finger holes of varying or equal diameter This difficulty can be remedied by using a more pliable and yet durable material that can be formed into an elongated, and cylindrical form. Looking into the possible available alternative materials the carrot was chosen. The carrot met all the requirements although the material shrinks over time due to dehydration.
Subject headings: musical instruments acoustical properties of solids
1. INTRODUCCIÓN
La ocarina se remonta a la América prehispánica, a culturas como la aymara y la quechua que utilizaban este instrumento musical acompañados de quenas, zamponas y otros instrumentos. Se usa en países como Perú, Bolivia, Venezuela, norte de Chile, norte de Argentina y toda Me-soamérica.
Cabe recalcar que aunque existe teoría que se refiere al origen de este instrumento en Europa, se trata de otro instrumento de similares características ya que la Ocarina Sud-Americana tiene orígenes de cientos de años atrás.
La podemos encontrar en modelos sencillos, dobles o triples y con diferente número de orificios, como muestra la Fig. 1. Su tono depende de la relación entre el volumen del aire y el área del agujero destapado. Por ejemplo, en una ocarina de cuatro agujeros de áreas iguales, cubriendo todos los
FIG. 1. Los tipos de ocarinas dependen del tamaño y del agujeros que posee.
FIG. 2.-temático.
Ocarina artesanal hecha de arcilla usada para el modelado maorificios previamente para luego pasar a quitar el dedo de uno de ellos, se obtiene la nota Sol, independientemente de qué agujero se destape. Si el área descubierta es igual, la nota producida es la misma.
2. ESCALA PENTATÓNICA
En nuestro sistema musical occidental, es conveniente utilizar sólo unas frecuencias concretas, a las cuales se las llama notas.
Las frecuencias de las notas se dividen en porciones llamadas "octavas", y cada octava se divide en 12 porciones llamadas notas. Cada nota de una octava tiene exactamente la mitad de frecuencia que la misma nota en la octava superior.
Con el oído humano solamente se pueden captar notas que sean superiores a 18/fz y por debajo de los 20kHz (muy aproximadamente). Es así que sólo podemos oír unas diez octavas como mucho, con doce notas cada una. Por debajo de 18/fz se llaman infrasonidos y por encima, ultrasonidos. El margen auditivo de las personas varía según la edad y otros factores. Los animales tienen un margen auditivo diferente; así, es muy conocido el hecho que los perros pueden sentir frecuencias mucho más altas, dentro del margen de los ultrasonidos.
Para que nuestra ocarina produzca sonidos agradables al oído humano es necesario ajustaría a una escala musical pentatónica porque es la escala más simple e intuitiva. La pentafonía sistema musical en que se usa la escala pentatónica se usa en la generalidad de los sistemas musicales tradicionales por ser usados por los músicos antiguos.
Nuestro oído humano tiene una "construcción" tal, que los sonidos cuyas frecuencias están en la proporción simple (2/1, 3/2, 4/3, etc.), suenan juntos de una manera agradable. Por ejemplo, la nota Do en su quinta octava tiene aproximadamente una frecuencia de 1046.5Hz y junto a un Sol en la misma octava de frecuencia 1567.9 Hz suena agradable porque están a una proporción de 3:2, respectivamente. Esta proporción es la combinación de sonidos más pura y se la suele llamar quinta justa.
En un sistema musical pentafónico se necesita una nota base que es la nota más grave que puede producir nuestro instrumento y una configuración de frecuencias de las notas;
FIG. 3. Relación lineal que existe entre el área descubierta y la frecuencia.
esta configuración es:
Tono, 9/8 Tono, (3/2)(8/9) = 4/3 Tono, 3/2 Tono, (3/2) (9/8) = 27/16 Tono, (3/2)2(8/9) = 2 Tono
La ocarina propuesta consta de una escala pentatónica, de cinco notas, con una nota base en Do sostenido (Do#). Y obedeciendo las reglas de esta escala su sistema musical tiene que ser como indica la Tabla 1.
3. MODELO MATEMÁTICO PARA EL NÚMERO DE ORIFICIOS DE LA OCARINA
Los orificios de la ocarina son la parte más importante en el diseño de la ocarina ya que aportan la variedad tonal propuesta en la Tabla 1.
Así que tuvimos que diseñar un modelo que pueda relacionar el área del orifico con la frecuencia que deseamos obtener, recordando que la nota depende de la suma de orificios destapados.
Primeramente se tomó una ocarina artesanal (Fig. 2) hecha de arcilla que constaba de seis orificios de \mm de diámetro cada uno y se aumentó nueve orificios más; teniendo un total de quince orificios. La posición de los orificios fue relevante porque no perturbó en gran magnitud a las frecuencias obtenidas.
Luego se procedió a tomar la frecuencia del total de orificios descubiertos con ayuda del software Analysis Center 2010. Su diseño nos da la facilidad de poder obtener la frecuencia con el micrófono en tiempo real.
Las frecuencias que se obtuvieron en relación con el área destapada (Tabla 2) nos dan una idea de que existe una relación lineal. Esta relación se puede observar en la Fig. 3.
TABLA 2
FlG. 4. Pasos realizados para la construcción de la ocarina. Las primeras dos fotos muestran la parte resonadora y la boquilla. Las dos últimas fotos muestran la parte que permite el flujo del aire.
Así que procedemos hacer el ajuste de curvas (1): Usando la regresión lineal observamos que>> es el área y x es la frecuencia. Con todo esto podemos concluir que nuestro modelo está hecho y es:
donde A es el área del orificio y i7 es la frecuencia deseada. O, en función del diámetro D del orificio, la cual emplearemos más adelante para la construcción de la ocarina de zanahoria:
4. CONSTRUCCIÓN DE LA OCARINA DE ZANAHORIA
En la construcción se necesitaron dos zanahorias lo más uniformes posible y de mayor volumen para poder evitar posibles fisuras al afollar. Y los pasos fueron los siguientes:
Tomamos la primera zanahoria y la cortamos por la parte delgada, después con un taladro ahuecamos la zanahoria por la porción central. Esta sirve para la parte resonante de la ocarina (Fig. 4a).
Hacemos un pequeño orificio cuadrangular en forma descendente en la parte central de su superficie con ayuda de un cincel y un estilete pequeño (Fig. 4b).
Agarramos la segunda zanahoria y la cortamos de tal
TABLA 3
manera que pueda ingresar en la primera zanahoria (Fig. 4c).
Cortamos en la parte que se inserta en la primera zanahoria un canal cuadrangular que permite el flujo de aire (Fig. 4d).
Cuando pasamos al último paso, el tallado de los orificios, tenemos que hacer uso del modelo (2) y la Tabla 1.
Reemplazamos la frecuencia de la segunda nota (porque la primera la obtenemos sin necesidad de orificios) en (2) y obtenemos el diámetro del primer orificio, es decir, con F = 775.2 se obtiene D = 6.7mm.
Para el segundo orificio debemos modificar (2) porque si lo usáramos directamente obtendríamos el área necesaria más el área del primer orificio obteniendo datos que no nos sirven de manera inmediata. Es por ello que usaremos:
Donde F es la frecuencia requerida y F0 es la frecuencia de la nota previa. Con este proceso obtenemos los diámetros requeridos para las frecuencias de las notas deseadas (Tabla 3). Se puede deducir que esta escala tiene una tendencia a orificio de 6mm. Lamentablemente estas medidas no se pueden conseguir a cabalidad porque se necesitaría algún instrumento que pueda perforar con una exactitud de micrómetros
5. PRUEBA DE SONIDO Y ANÁLISIS ESPECTRAL
Ya construida la ocarina ahora pasamos a analizar su afinación, es decir, con un análisis espectral observamos los picos más altos de frecuencia frente a su amplitud. Estos picos obtenidos gracias a las herramientas del software MATLAB 7.4.0 como la transformada rápida de Fourier (FFT) nos dicen su acercamiento a la nota musical que esperamos obtener.
Inicialmente de manera un tanto obvia podemos asegurar que las notas más graves, de manera específica las tres primeras notas (Fa, Sol y La# ) suenan mejor afinadas que las ultimas tres (Do, Re y Fa), posiblemente por la variación en la octava a que corresponden. Pero más específicamente se puede ver por las gráficas obtenidas que se muestran en la Figs. 5 a 10, que las dos primeras notas tienen una mayor precisión respecto al tono.
La frecuencia de la primera nota Fa de la cuarta octava (nuestra nota base) tuvo un rango de frecuencia entre 660 Hz y 680 Hz, se acercó bastante al resultado esperado, exactamente unos 9AHz de diferencia (Fig. 5).
En la segunda nota Sol (Fig. 6) se esperó que diera 115.2Hz y obtuvimos un rango entre 160Hz y 780Hz. El resultado requerido se encuentra dentro de nuestro rango. Sin embargo, se observa que el rango va en aumento.
.
FlG. 9. Esta gráfica corresponde a la frecuencia de la quinta nota Re.
Los resultados de las notas La sostenido y Do tienen una similitud bastante interesante: Para la nota La sostenido calculamos una frecuencia de 918.8/fz, y en la Fig. 7 vemos
FlG. 10. Esta gráfica corresponde a la frecuencia de la sexta nota Fa.
que la mayor amplitud se encuentra entre 900Hz y 950Hz; en la cuarta nota Do se calculó una frecuencia de 1033.6Hzy como se muestra en la Fig. 8 se observa una mayor amplitud entre 1000Hz y 1050Hz. Las Figs. 7 y 8 muestran un intervalo de mayor amplitud, de 50Hz y la frecuencia calculada está dentro de los rangos.
Cuando las notas de nuestra ocarina subieron a una quinta octava el sonido de las notas Re y Fa simplemente ya no se escuchó y esto se puede ver en la gran dispersión de las Figs. 9 y 10.
6. CONCLUSIONES
La ocarina, y en especial esta ocarina hecha de un vegetal muy común, tiene una cierta ventaja en la construcción pero una seria desventaja de duración puesto que mientras el tiempo se prolongue su sonido irá decayendo en claridad ya que la boquilla se deshidrata cada vez más perdiendo el corte en el flujo que hace el sonido. Otro problema con este instrumento es que, como explicamos en su construcción, la zanahoria que dirige el flujo del aire tiene que caber exactamente sin dejar que escape el aire del resonador, y por el mismo problema de la deshidratación este va perdiendo vo-lumeny cada vez deja escapar más aire, cambiando las notas y el tono del instrumento.
Pero al estar éste fresco, ésto es, aproximadamente unos dos días después de la inmediata construcción, puede darnos sonidos muy claros, variados y afinados de la manera que se requiera. Si usamos notas que pertenezcan a la cuarta o tercera octava obtendremos sonidos mucho más claros porque el instrumento hecho de zanahoria no alcanza frecuencias tan altas demandadas para la quinta octava en adelante.
Y se puede decir que este es el instrumento más sencillo, de bajo costo y sorprendente, porque al no depender de la resonancia, es decir, siendo su forma relevante a la hora de la construcción, se puede formar de materiales orgánicosaunque parezca humorístico como pepinos manzanas o hasta muy difícilmente de huevo, porque estos tienen un gran volumen y una gran versatilidad al momento de su construcción.
REFERENCIAS
1.- Halliday, R. & Krane. 1998, Física, vol. 1, cuarta ed. (CECSA) [ Links ]^rND^sHalliday^nR.^rND^nA. R.^sTicona Bustillos^rND^nG. M.^sRamírez Ávila^rND^nA. R.^sTicona Bustillos^rND^nG. M.^sRamírez Ávila^rND^nA. R^sTicona Bustillos^rND^nG. M^sRamírez Ávila
CONDICIONES DE SINCRONIZACIÓN EN DOS PÉNDULOS ACOPLADOS
Sinchronization Conditions of Two Coupled Pendula
A. R. Ticona Bustillos1, G. M. Ramírez Ávila1,2
1Instituto de Investigaciones Físicas, Universidad Mayor de San Andrés
Casilla 8635, La Paz, Bolivia
2AG Nichtlineare Dynamik (S) / Kardiovasculäre Physik
Institut für Physik
Humboldt-Universität zu Berlin, Robert-Koch-Platz 4, 10115 Berlin, Alemania
(Recibido 5 de febrero de 2010; aceptado 7 de marzo de 2010)
Abstract
Based on a damped pendulum discrete model, we studied the synchronization conditions for two coupled pendula, varying both the pendula's features and coupling conditions. We found the basis for attraction in several situations in which the control parameters were fixed. Varying the control parameters (length, mass and damping coefficient), we found phase diagrams related to the initial conditions of one of the pendula; in these diagrams we identified synchronization regions. We emphasize the synchronization with a winding number ρ ≈ 1 (synchronization 1:1); nevertheless, other synchronization orders are possible (ρ ≠ 1).
Subject headings: dynamical systems (non-linear) - synchronization - coupled oscillators Código(s) PACS: 05.45._a, 05.45.Xt
Resumen
Con base en un modelo discreto de péndulo amortiguado, se estudian las condiciones de sincronización para dos péndulos acoplados, variando las características propias de los péndulos, así como las condiciones de acoplamiento. Se encuentran las cuencas de atracción para diferentes situaciones en las que se fijan los parámetros de control. Variando los valores de los parámetros de control (longitud, masa y coeficiente de disipación), se encuentran diagramas de fase relacionados con las condiciones iniciales de uno de los péndulos, mediante los cuales es posible identificar regiones de sincronización. Se hace énfasis en la sincronización 1:1 aunque sincronizaciones de otros órdenes son también posibles.
Descriptores: sistemas dinámicos no-lineales - sincronización - osciladores acoplados
1 Introducción
El fenómeno de sincronización es muy común en la naturaleza y muchos sistemas de diversa índole exhiben este comportamiento [212003Strogatz,172003Rosenblum & Pikovsky], destacándose entre ellos los sistemas biológicos [42001Glass], las reacciones químicas [192003Shabunin et al.Shabunin, Astakhov, Demidov, Provata, Baras, Nicolis, & Anishchenko,32005Fukuda et al.Fukuda, Morimura, & Kai], los circuitos electrónicos [91998Kittel et al.Kittel, Parisi, & Pyragas,132003 Ramírez Ávila et al. Ramírez Ávila , Guisset, & Deneubourg], los láseres [181994Roy & Thornburg,112009López-Gutiérrez et al.López-Gutiérrez, Posadas-Castillo, López-Mancilla, & Cruz-Hernández] y por supuesto, los péndulos [71986Huygens,202003Smith et al.Smith, Blackburn, & Baker] que desde el punto de vista histórico, constituye el primer sistema en el cual se observó sincronización [61673Huygens]; los péndulos, a pesar de ser sistemas en apariencia simples siguen concitando la atención y son sujetos de investigación tanto desde el punto de vista teórico como experimental [22009Baker & Blackburn]. En este trabajo, nos abocamos a encontrar las condiciones de sincronización para dos péndulos acoplados. Se sabe que para sistemas no lineales disipativos, es posible la existencia de más de un atractor, por lo que diferentes condiciones iniciales pueden evolucionar hacia cualquiera de los atractores coexistentes; así, el conjunto de condiciones iniciales que se aproxima a un atractor, es llamado la cuenca de atracción de este atractor [81999Kapitaniak & Bishop]. En §2 se describe el modelo utilizado para estudiar el sistema de dos péndulos acoplados con disipación, en la primera parte de §3 se muestran diferentes cuencas de atracción del sistema que nos dan una idea de cómo obtener condiciones de sincronización; posteriormente, se muestran las regiones de sincronización, que pueden asociarse a las llamadas lenguas de Arnold (en el sentido de que representan regiones de sincronización de una manera similar a la utilizada en [162001Pikovsky et al.Pikovsky, Rosenblum, & Kurths,122004 Ramírez Ávila ]1) para diferentes situaciones en las que los parámetros varían. Se darán finalmente en §4 las conclusiones y las perspectivas de este trabajo.
2 Modelo
Al igual que en [142008 Ticona Bustillos & Ramírez Ávila ], utilizamos el modelo discretizado de péndulo y lo adaptamos para la situación en la cual se tienen dos péndulos acoplados con disipación, cuyas ecuaciones están dadas por:
|
donde los superíndices 1 y 2 identifican a cada uno de los péndulos acoplados2 y w representa las condiciones de acoplamiento, la cual contiene a la distancia de separación entre los péndulos y a la rigidez del material usado para acoplar los mismos; además:
| (2) |
Figure 1: Cuenca de atracción para identificar las regiones de sincronización 1:1 (regiones oscuras) cuando la relación de longitud de los péndulos acoplados es l2 = 2l1 con ρ = 1.0000 ± 0.0020. Se observa una simetría en la región de sincronización que nos lleva a interpretar que la sincronización 1:1 es favorecida cuando las condiciones iniciales son bastante diferentes.
Figure 2: Cuenca de atracción delimitando las situaciones en las que ρ = 1.0000 ± 0.0001 (región oscura), cuando se considera que las masas de las lentejas de los péndulos son de 0.050 kg y 1.000 kg respectivamente. La sincronización 1:1 es favorecida cuando las condiciones iniciales de ambos péndulos tienen valores similares y con el mismo signo (extremos inferior izquierdo y superior derecho), tienen valores similares con diferente signo (región casi rectilínea con pendiente aproximada de -1) o tienen valores diferentes pero con el mismo signo y con la particularidad que uno de los péndulos debe tener una condición inicial no muy lejana a cero (el resto de las regiones de sincronización).
Este modelo fue comparado con valores experimentales en [142008 Ticona Bustillos & Ramírez Ávila ], dando muy buenos resultados, con lo que pudimos estudiar las principales características de transmisión de movimiento en función de las características del sistema. Además, se pudo verificar la fuerte dependencia de este sistema con las condiciones iniciales y algunas características de sincronización, las cuales estudiamos con más detalle a continuación.
3 Resultados
Primeramente, se determinaron las cuencas de atracción para posibles situaciones en las cuales se pueden controlar ciertos parámetros, como la masa de la lenteja y la longitud del péndulo, la disipación y el factor de acoplamiento. Para el trabajo numérico, se utilizó el valor de la aceleración debida a la gravedad en La Paz: g = 9.775 m/s2; este valor ya fue utilizado en [142008 Ticona Bustillos & Ramírez Ávila ], para comparar los resultados con medidas experimentales, las cuales fueron realizadas en esta ciudad. Analizamos el comportamiento de la relación ρ = T2/T1 = 1 de los períodos de los péndulos, variando las condiciones iniciales de los mismos. Esta consideración, nos permite identificar las condiciones iniciales para las cuales es posible la sincronización 1:1 entre los péndulos. Tomamos condiciones iniciales desde -0.7 hasta 0.7 rad. Nuestro primer estudio (Figs. 1-2) implica el establecimiento de las regiones de sincronización 1:1 (regiones oscuras) de las que no siguen esta relación (regiones claras), utilizando los siguientes valores para el factor resultante de disipación y el factor de acoplamiento: b(1) = b(2) = 0.9999 y w = 1×10−6 respectivamente. Cada punto mostrado en estas figuras corresponde a un promedio de ρ obtenido sobre las 50 últimas oscilaciones consideradas en las simulaciones para cada péndulo. En la Fig. 1 mostramos un caso en el cual ambos péndulos tienen una masa de 0.0500 kg y los largos son de 1.00 m y 2.00 m, respectivamente. Los puntos nos muestran los valores obtenidos para ρ a partir de las diferentes condiciones iniciales, donde hacemos la distinción entre las regiones en las cuales se presenta la sincronización 1:1 (oscura) de las que no presentan esta característica (clara). Podemos observar que la región en forma de cruz corresponde a las condiciones de sincronización 1:1 ya que se obtienen valores de ρ muy cercanos a la unidad. Es interesante observar la simetría en la Fig. 1, lo que nos indica por una parte que el conjunto de condiciones iniciales que nos llevan a ρ ≈ 1 es bastante reducido y por otra parte que para obtener esta situación, las condiciones iniciales de ambos péndulos deben ser muy próximas a cero o en su defecto, uno de los péndulos debe tener una condición inicial cercana a cero y el otro un ángulo inicial considerable, situación que da lugar a que la región de sincronización 1:1 se ensanche. Esto último implica que es más fácil alcanzar una sincronización 1:1 entre los péndulos cuando uno de ello empieza en el reposo relativo (θ0 = 0) y el otro con una amplitud grande.
Figure 3: Cuenca de atracción cuando l1 = l2 = 1.00 m y m1 = m2 = 0.050 kg, fijando b(1) = b(2) = 0.999 con ρ = 1.0000 ± 0.0004 en las regiones de sincronización, las cuales presentan cierta simetría en el sentido que los péndulos deben tener condiciones iniciales de signo contrario cuando los ángulos iniciales son pequeños. Es notable también el ensanchamiento de las regiones de sincronización a medida que se consideran condiciones iniciales mayores. Se observa finalmente la aparición de una nueva región de sincronización no simétrica para valores pequeños del ángulo inicial del péndulo 2.
Figure 4: Cuenca de atracción para longitudes y masas iguales para los péndulos, con b(1) = b(2) = 0.9999, cuando se considera un factor de acoplamiento w = 5×10−6. Las regiones de sincronización oscuras corresponden al valor ρ = 1.0000 ± 0.0021. Existe una simetría entre estas regiones y se nota que la densidad de las mismas es mayor que en los anteriores casos.
Para los resultados de la Fig. 2 el largo de ambos péndulos es de 1.00 m y las masas son de 0.050 kg y 1.000 kg respectivamente. En este caso, se observa que también existe una simetría en las regiones de sincronización 1:1, donde se pueden resaltar tres aspectos interesantes, a saber: (i) Los extremos inferior izquierdo y superior derecho presentan una simetría que implica que las condiciones iniciales de ambos péndulos deben tener valores similares y del mismo signo. (ii) La región rectilínea de pendiente negativa muestra que los péndulos deben comenzar cada uno en extremos opuestos para favorecer la sincronización 1:1. (iii) Las otras regiones de sincronización muestran que los valores de las condiciones iniciales pueden ser bastante diferentes pero deben tener el mismo signo y además uno de los péndulos debe comenzar con un ángulo pequeño no muy alejado de cero. En las Figs. 3 y 4 utilizamos valores iguales para los largos y las masas de ambos péndulos (l1 = l2 y m1 = m2), fijando el valor de b(1) = b(2) = 0.999. En la Fig. 3, la densidad de regiones de sincronización 1:1 es mayor, resaltando las regiones en las cuales los ángulos iniciales deben ser pequeños para ambos péndulos y de diferente signo (parte central), notándose luego un ensanchamiento considerable a medida que los ángulos iniciales son mayores; se observa también la aparición de nuevas regiones de sincronización no simétricas como la que aparece en la parte central y derecha de la Fig. 3, en la cual, ambos ángulos iniciales tienen valores positivos pero el del péndulo 2 puede permanecer cercano a cero, en tanto que el del péndulo 1 puede crecer linealmente. Finalmente, en la Fig. 4 fijamos el factor de acoplamiento en w = 5×10−6. La distribución de las regiones de sincronización es parecida a la de la Fig. 3 pero estas regiones se ensanchan, lo que nos lleva a pensar que gracias al mayor acoplamiento entre los péndulos, la sincronización 1:1 se ve favorecida. Nuevamente, la lógica de condiciones iniciales de la misma magnitud pero de signos contrarios son las que llevan a la sincronización 1:1; sin embargo, como se señaló anteriormente, existe un aumento notorio de estas regiones con una tendencia a ocupar gran parte del plano θ02 vs. θ01. También se pueden observar manchas difusas representando regiones de sincronización 1:1 lo que muestra una vez más la tendencia que tienen las regiones de sincronización a hacerse más extensas, máxime si el factor de acoplamiento aumenta su valor. A continuación, se estudian la sincronización cuando se hacen variar un parámetro de control y las condiciones iniciales del segundo péndulo, fijando los valores del primer péndulo en: m1 = 0.050 kg y l1 = 1.00 m, siendo la condición inicial del primer péndulo en todos los casos θ01 = 0.5 rad. De esta manera, se obtienen las regiones que guardan similitud con las lenguas de Arnold como se mencionó anteriormente, por lo que se utilizará en lo sucesivo esta terminología.
Figure 5: (Color online) Plano de fases l2 vs. θ02 para el cual se muestra la lengua de Arnold correspondiente a ρ = 1.0000 ± 0.0289 (región oscura identificada con 1), cuando m1 = m2 = 0.050 kg. Las otras regiones que se identifican con 0, 2-7 corresponden a zonas en las que se pueden encontrar otros órdenes de sincronización con ρ = (0,1), ρ = (1,2), ... , ρ = (6,7). Es claro que en estos intervalos, ρ ⊂ P+ por lo que no sólo existen regiones de sincronización sino también cuasi-periodicidades y caos.
Figure 6: (Color online) Plano de fases l2 vs. θ02 para el cual se muestra la lengua de Arnold correspondiente a ρ = 1.0000 ± 0.0201 (región oscura identificada con 1), cuando m1 = 5m2 = 0.050 kg. Las regiones identificadas con los números 0, 2, 3, 4 y 5 corresponden a regiones en las cuales ρ = (0,1), ρ = (1,2), ρ = (2,3), ρ = (3,4) y ρ = (4,5) respectivamente. Se nota el desplazamiento del eje de simetría hacia la izquierda (θ02 < 0), lo que repercute en el hecho que la región 5 sólo se presente en la parte derecha del plano. Las regiones blancas en forma de "dedos" al interior de la lengua de Arnold es posible que sean el resultado de los errores del cálculo numérico pues si se considera un intervalo ligeramente mayor para el valor de ρ de la lengua de Arnold, estas desaparecen.
Figure 7: Plano de fases m2 vs. θ02 para el cual se muestra la lengua de Arnold correspondiente a. ρ = 1.0000 ± 0.0005 (región oscura identificada con 1), cuando l1 = l2 = 1.00 m. Se evidencia la presencia de regiones para ρ = (0,1) y ρ = (1,2), (regiones identificadas con 0 y 2) aunque muchas de las regiones, sobre todo las identificadas con 0, en realidad pueden ser consideradas como parte de la lengua de Arnold ρ ≈ 1 puesto que la precisión numérica es mayor a la que se puede obtener experimentalmente.
En primer lugar tomamos en cuenta la variación del largo del segundo péndulo, consideramos valores desde 0.10 m hasta 50.00 m, la masa de la lenteja de este péndulo es mantenida en 0.050 kg. En las Figs. 5-11, se muestran las regiones equivalentes a las lenguas de Arnold, donde nos concentramos en la región con ρ ≈ 1, aunque, otras regiones con otros órdenes de sincronización pueden estar también presentes. En la Fig. 5, se muestra el diagrama de fase l2 vs. θ02, en el cual aparece claramente la lengua de Arnold con ρ ≈ 1 (región oscura identificada con el número 1). Además, se identifican regiones en las que ρ toma valores comprendidos en los intervalos (0,1), (1,2), ... , (6,7) que se identifican con los números 0, 2, ... . , 7 respectivamente. Es interesante notar que estas regiones, en general, están bien definidas y son simétricas respecto de la situación en la cual θ02 = 0. Modificando el valor del coeficiente de disipación, no se notan cambios apreciables, pero cambiando el valor de la masa del segundo péndulo, por ejemplo a 0.010 kg, constatamos que la lengua de Arnold para ρ ≈ 1 es más extensa que en el caso de masas iguales de las lentejas de ambos péndulos, como se puede observar en la Fig. 6. Además, se constata que la simetría existente en el caso anterior (Fig. 5) se rompe y se podría pensar que el eje de simetría ya no pasa por θ02 = 0, sino que está desplazado a la izquierda del plano (θ02 < 0). Como en el anterior caso, se identifican otras regiones en las cuales ρ = (0,1), ρ = (1,2), ... , ρ = (4,5), identificando las mismas con los números 0, 2, ... , 5 respectivamente, siendo ρ real y positivo lo que puede dar lugar a regiones en las cuales el comportamiento de los péndulos es periódico, cuasi-periódico o caótico. Un aspecto interesante es la aparición de regiones claras al interior de la lengua de Arnold; en principio, esto significaría que regiones con ρ < 1 rompen la simetría de la lengua de Arnold: Sin embargo, estas regiones en particular si bien poseen ρ < 1, este valor es muy cercano a 1 por lo que se podría pensar que los errores propios del cálculo numérico hacen que no se consideren estas regiones en el intervalo elegido para ρ. Otro aspecto a resaltar es la aparición de una pequeña región ρ = (4,5) sólo en la parte derecha del plano; esto se justifica fácilmente por el hecho de la ruptura de simetría referida líneas arriba. Ahora, fijando l1 = l2 = 1.00 m y variando el valor de m2, podemos observar en la Fig. 7 que la lengua de Arnold para ρ ≈ 1 ya no tiene una estructura sencilla como era el caso en las Figs. 5-6 y ocupa varias regiones del plano, el cual contiene también regiones con ρ = (0,1) y ρ = (1,2). Nuevamente, nuestros resultados son conservadores en cuanto a la región de sincronización 1:1 puesto que la precisión numérica que se considera supera la precisión que se podría tener en las mediciones experimentales, por lo que es posible que las regiones de sincronización 1:1 ocupen una extensión mayor en el plano de fases m2 vs. θ02. Por otra parte, es interesante notar que en la zona en la cual m2 ≈ m1, se tiene una tendencia a que la región de sincronización esté presente casi para todos los valores de θ02, lo que en cierta manera confirma los resultados obtenidos en la Fig. 3.
Figure 8: Plano de fases m2 vs. θ02 para el cual se muestra la lengua de Arnold correspondiente a. ρ = 1.0000 ± 0.0006 (región oscura identificada con 1), cuando l2 = 2l1 = 2.00 m. También existen regiones en las cuales, ρ = (0,1) y ρ = (1,2), (regiones identificadas con 0 y 2), pero al igual que en los anteriores casos, ρ = (0,1) puede ser absorbida por la región de sincronización 1:1.
En cambio si l2 = 2l1 = 2.00 m, la lengua de Arnold en el plano m2 vs. θ02 es muy pequeña, como se ve en la Fig. 8. Es claro que la mayor parte del plano corresponde a valores ρ = (1,2) por lo que en estas regiones bien podrían haber sincronizaciones de otros órdenes, cuasi-periodicidades o caos. Se observa también que cuando m2 → m1, las únicas condiciones iniciales θ02 que preservan la sincronización 1:1 son aquellas en las que θ02 → 0, en concordancia con los resultados de la Fig. 1. Los anteriores resultados indican que tanto la masa de la lenteja como la longitud del péndulo juegan papeles importantes en la sincronización de los dos péndulos acoplados. De igual manera podemos analizar la sincronización en función del factor de acoplamiento. En la Fig. 9 se nota que para valores pequeños del factor de acoplamiento, existe sincronización casi para todas las condiciones iniciales, en cambio para valores grandes de dicho factor, la sincronización parece desaparecer dando lugar a una "granularidad" en la que se tienen valores de ρ = (0,2). En la Fig. 9, se consideraron masas de las lentejas de los péndulos iguales pero no las longitudes de los mismos; sin embargo, el comportamiento cuando se tiene l1 = l2 no se modifica sustancialmente. Finalmente, si variamos el valor del coeficiente de disipación λ = λ(1) = λ(2) tenemos que para péndulos con las mismas características (l1 = l2 y m1 = m2), la región de sincronización 1:1, abarca casi todo el plano de fases, como se puede observar en la Fig. 10. Regiones con ρ = (0,1) y ρ = (1,2) están también presentes en el plano pero al igual que en los anteriores casos, las regiones con ρ = (0,1), en general pueden ser absorbidas por las regiones de sincronización 1:1, lo que aumentaría aún más la densidad de regiones con ρ ≈ 1. Otro aspecto interesante es el referido a los valores elevados de λ, para los cuales, la sincronización se pierde; esto está en relación con el hecho de que cuando la disipación es muy grande, los péndulos tienden a volver a su situación de reposo relativo y las medidas de comportamiento síncrono de los péndulos ya no es posible hacerlas con un transientes adecuado.
Figure 9: Plano de fases w vs. θ02 para el cual se muestra la lengua de Arnold correspondiente a. ρ = 1.0000 ± 0.0248 (región oscura identificada con 1), cuando l2 = 2l1 = 2.00 m y m1 = m2 = 0.050 kg. La región de sincronización 1:1 está bien definida para valores w ≈ 5x10−6; para valores mayores, ya no se encuentra una región bien definida para la sincronización 1:1 y se tiene una especie de "granularidad" con valores para ρ comprendidos en el intervalo (0,2).
En cambio si las longitudes de los péndulos son diferentes, la región de sincronización 1:1 se reduce drásticamente, como se puede observar en la Fig. 11. Trabajando con ρ = 1.0000 ± 0.0013 se observa que la lengua cerca de la parte central contiene regiones blancas y grises que en principio significan valores de ρ = (0,1) y ρ = (1,2) respectivamente; sin embargo, si se considera ρ = 1.0000 ± 0.0097, estas regiones desaparecen de la lengua, la cual queda bien definida y además, la región correspondiente a valores elevados del coeficiente de disipación considerados en el plano, pasa a ser considerada como región de sincronización 1:1. En este último caso, se puede pensar también que el hecho de trabajar con valores elevados del factor de disipación λ, no permite cuantificar adecuadamente las oscilaciones tal como se explicó para la Fig. 10. El hecho que la región de sincronización 1:1 sea relativamente pequeña concuerda con los anteriores resultados para situaciones en las que l1 ≠ l2 (Figs. 1 y 8), en las cuales se tiene la misma situación; sin embargo, no se debe olvidar los resultados mostrados en las Figs. 5 y 6 en las cuales, el intervalo de condiciones iniciales θ02 que conducen a la sincronización 1:1 crece con la longitud l2.
4 Conclusiones y Perspectivas
Se determinaron las cuencas de atracción para varios valores de parámetros, pudiéndose observar que las mismas presentan diferentes formas según el valor que toman los parámetros. Nos concentramos en las condiciones iniciales de los péndulos para las que puede existir sincronización 1:1.
Figure 10: Plano de fases λ vs. θ02 para el cual se muestra la lengua de Arnold correspondiente a. ρ = 1.0000 ± 0.0006 (región oscura identificada con 1), cuando l2 = l1 = 1.00 m y m1 = m2 = 0.050 kg. Otras regiones con ρ = (0,1) y ρ = (1,2) aparecen pero es notable la predominancia de la región de sincronización 1:1. En la parte superior del plano (valores elevados de λ, existe una región en la que ρ ≈ (0.96,1)
Figure 11: Plano de fases λ vs. θ02 para el cual se muestra la lengua de Arnold correspondiente a. ρ = 1.0000 ± 0.0013 (región oscura identificada con 1), cuando las longitudes de los péndulos son diferentes: l2 = 2l1 = 2.00 m y m1 = m2 = 0.050 kg. Las regiones blancas y grises que aparecen sobre la lengua, así como la región blanca en la parte superior, son regiones con ρ ≈ 1 pero con menor precisión de la establecida al principio; por lo que se espera que sean parte de la lengua con sincronización 1:1.
Se han identificado las regiones de sincronización para varias condiciones de los péndulos, en función a las condiciones iniciales del péndulo 2. Si bien nos abocamos a la determinación de la sincronización de orden 1:1, es claro que existen otros órdenes de sincronización. Además, como se tienen valores de ρ ⊂ P+, pueden existir también comportamientos cuasi periódicos y caóticos. Se puede observar que la longitud y la masa de los péndulos ejercen una gran influencia en la sincronización. Se ve claramente en todos los casos estudiados, la importancia de la precisión numérica con la que se trabaja y su relación con la exactitud que se puede lograr en mediciones experimentales. A partir de estos resultados, podemos analizar cuasi periodicidades, diagramas de bifurcación, rutas hacia el caos, en las regiones donde no existe sincronización. Evidentemente, este modelo es susceptible a ampliarse a más de dos péndulos lo que constituye la continuación natural de esta investigación. Una parte importante de nuestra investigación es trabajar tanto con mediciones experimentales como con modelos que justamente son validados por estas mediciones. Si bien, por el momento se tienen concordancias importantes entre los resultados experimentales con los obtenidos a partir del modelo, no se descarta para el futuro considerar el sistema de ecuaciones diferenciales para describir el sistema de péndulos acoplados. Finalmente, se pretende que este tipo de trabajos pueda a posteriori encontrar su aplicación también en aspectos académicos y didácticos, tal como ocurrió con la primera parte de nuestra investigación con péndulos [152010 Ticona Bustillos & Ramírez Ávila ].
Agradecimientos
Agradecemos los valiosos comentarios de los árbitros de la revista, gracias a los cuales pudimos expresar con mayor claridad las ideas de este trabajo. GMRA agradece al Deutscher Akademischer Austausch Dienst (DAAD) por la beca de investigación otorgada.
References
- [12002Bai et al.Bai, Lonngren, & Sprott]
-
1.- Bai, E.-W., Lonngren, K. E., & Sprott, J. C. 2002, Chaos, Solitons & Fractals, 13, 1515 [ Links ]
- [22009Baker & Blackburn]
-
2.- Baker, J. L. & Blackburn, J. L. 2009, The Pendulum: A Case Study in Physics (Oxford: Oxford University Press) [ Links ]
- [32005Fukuda et al.Fukuda, Morimura, & Kai]
-
3.- Fukuda, H., Morimura, H., & Kai, S. 2005, Physica D: Nonlinear Phenomena, 205, 80 [ Links ]
- [42001Glass]
-
4.- Glass, L. 2001, Nature, 410, 277 [ Links ]
- [52003Grosu]
-
5.- Grosu, I. 2003, Master-slave synchronization and control of chaotic Sprott circuit (11th International IEEE Workshop on Nonlinear Dynamics of Electronic Systems, Scuol/Schuls, Switzerland) [ Links ]
- [61673Huygens]
-
6.- Huygens, C. 1673, Horologium Oscillatorium (Paris: Apud F. Muguet) [ Links ]
- [71986Huygens]
-
7.- -. 1986, The pendulum clock or geometrical demonstrations concerning the motion of pendula as applied to clocks (Iowa: Ames, Iowa : Iowa State University Press) [ Links ]
- [81999Kapitaniak & Bishop]
-
8.- Kapitaniak, T. & Bishop, S. R. 1999, The Illustrated Dictionary of Nonlinear Dynamics and Chaos Chichester (John Wiley & Sons, Inc.) [ Links ]
- [91998Kittel et al.Kittel, Parisi, & Pyragas]
-
9.- Kittel, A., Parisi, J., & Pyragas, K. 1998, Physica D: Nonlinear Phenomena, 112, 459 [ Links ]
- [102004Li et al.Li, Chen, Cai, & Zhao]
-
10.- Li, Y.-N., Chen, L., Cai, Z.-S., & Zhao, X.-Z. 2004, Chaos, Solitons & Fractals, 22, 767 [ Links ]
- [112009López-Gutiérrez et al.López-Gutiérrez, Posadas-Castillo, López-Mancilla, & Cruz-Hernández]
-
11.- López-Gutiérrez, R. M., Posadas-Castillo, C., López-Mancilla, D., & Cruz-Hernández, C. 2009, Chaos, Solitons & Fractals, 42, 277 [ Links ]
- [122004 Ramírez Ávila ]
-
12.- Ramírez Ávila, G. M. 2004, Synchronization phenomena in light-controlled oscillators (PhD thesis) [ Links ]
- [132003 Ramírez Ávila et al. Ramírez Ávila , Guisset, & Deneubourg]
-
13.- Ramírez Ávila, G. M., Guisset, J. L., & Deneubourg, J. L. 2003, Physica D: Nonlinear Phenomena, 182, 254 [ Links ]
- [142008 Ticona Bustillos & Ramírez Ávila ]
-
14.- Ticona Bustillos, A. R. & Ramírez Ávila, G. M. 2008, Revista Boliviana de Física, 14, 121 [ Links ]
- [152010 Ticona Bustillos & Ramírez Ávila ]
-
15.- -. 2010, Revista Boliviana de Física, 16, 41 [ Links ]
- [162001Pikovsky et al.Pikovsky, Rosenblum, & Kurths]
-
16.- Pikovsky, A., Rosenblum, M., & Kurths, J. 2001, Synchronization: a universal concept in nonlinear sciences, vol 12 (New York: Cambridge University Press) [ Links ]
- [172003Rosenblum & Pikovsky]
-
17.- Rosenblum, M. & Pikovsky, A. 2003, Contemporary Physics, 44, 401 [ Links ]
- [181994Roy & Thornburg]
-
18.- Roy, R. & Thornburg, K. S. J. 1994, Physical Review Letters, 72, 2009 [ Links ]
- [192003Shabunin et al.Shabunin, Astakhov, Demidov, Provata, Baras, Nicolis, & Anishchenko]
-
19.- Shabunin, A., Astakhov, V., Demidov, V., Provata, A., Baras, F., Nicolis, G., & Anishchenko, V. 2003, Chaos, Solitons & Fractals, 15, 395 [ Links ]
- [202003Smith et al.Smith, Blackburn, & Baker]
-
20.- Smith, H. J. T., Blackburn, J. A., & Baker, G. L. 2003, International Journal of Bifurcation and Chaos, 13, 7 [ Links ]
- [212003Strogatz]
-
21.- Strogatz, S. H. 2003, Sync: The Emerging Science of Spontaneous Order (New York: Hyperion Press) [ Links ]
- [221994Sugawara et al.Sugawara, Tachikawa, Tsukamoto, & Shimizu]
-
22.- Sugawara, T., Tachikawa, M., Tsukamoto, T., & Shimizu, T. 1994, Physical Review Letters, 72, 3502 [ Links ]
Footnotes:
1Estas regiones de sincronización se diferencian de las lenguas de Arnold, en el sentido que en lugar de representarlas en un plano intensidad de acoplamiento vs. frecuencia, son representaciones en el plano parámetro de control vs. condiciones iniciales. 2Por razones de comodidad en la notación, se utilizará a lo largo del artículo mi = m(i) y li = l(i).
File translated from TEX by TTH, version 3.89.
^rND^sKittel^nA.^rND^sParisi,^nJ.^rND^sPyragas^nK.^rND^nLi, Y.-N.^sChen^rND^sCai^nZ.-S^rND^sZhao^nX.-Z.^rND^sLópez-Gutiérrez^nR. M.^rND^sPosadas-Castillo^nC.^rND^sLópez-Mancilla^nD.^rND^sCruz-Hernández^nC.^rND^sRamírez Ávila^nG. M.^rND^sGuisset^nJ. L.^rND^sDeneubourg^nJ. L.^rND^sTicona Bustillos^nA. R.^rND^sRamírez Ávila^nG. M.^rND^sPikovsky^nA.^rND^sRosenblum^nM.^rND^sKurths^nJ.^rND^sRosenblum^nM.^rND^sPikovsky^nA.^rND^sRoy^nR.^rND^sThornburg^nK. S. J.^rND^sShabunin^nA.^rND^sAstakhov^nV.^rND^sDemidov^nV.^rND^sProvata^nA.^rND^sBaras^nF.^rND^sNicolis^nG.^rND^sAnishchenko^nV.^rND^sSmith^nH. J. T.^rND^sBlackburn^nJ. A.^rND^sBaker^nG. L.^rND^sStrogatz^nS. H.^rND^sSugawara^nT.^rND^sTachikawa^nM.^rND^sTsukamoto^nT.^rND^sShimizu^nT.^rND^nFrancesco^sZaratti^rND^nRicardo^sForno^rND^nGonzalo^sGutiérrez^rND^nRoger^sApaza^rND^nFernando^sVelarde^rND^nFrancesco^sZaratti^rND^nRicardo^sForno^rND^nGonzalo^sGutiérrez^rND^nRoger^sApaza^rND^nFernando^sVelarde^rND^nFrancesco^sZaratti^rND^nRicardo^sForno^rND^nGonzalo^sGutiérrez^rND^nRoger^sApaza^rND^nFernando^sVelarde
On 8 Nov 2010, 09:37.\uppercase{Medidas de albedo en UV-B en el Salar de Uyuni}\ \vspace{1em}{ Albedo's Measurements in UV-B at Uyuni's Salt Lake (Bolivia)} MEDIDAS DE ALBEDO EN UV-B EN EL SALAR DE UYUNI
Albedo's Measurements in UV-B at Uyuni's Salt Lake (Bolivia)Francesco Zarattia, Ricardo Forno, Gonzalo Gutiérrez,
Roger Apaza & Fernando VelardeLaboratorio de Física de la Atmósfera
Instituto de Investigaciones Físicas, FCPN, UMSA
aEmail:fzaratti@fiumsa.edu.bo
(Recibido el 28 de julio de 2010; aceptado el 11 de agosto de 2010)
Abstract
We report the final results of a biannual project concerning two measurement campaigns of ultraviolet radiation (RUV) at Uyuniâs Salt Lake: in May and November 2008. The gathered data and their respective analysis resulted in important findings relating to the projectâs principle objectives: the measurement of superficial albedo of the UV band and the quantification of the extent of albedo in the area surrounding the Uyuni salt lake., In the first case we found that the albedo depends not only on the time of year, such as, the rainy season between January and March but also on climate change effects in the region which could lead to a decrease in albedo, a sign of the salt lake's environmental contamination in particular from anthropogenic activities. In addition we measured the extent of the albedo effect in the salt lake and surrounding area however the analysis of these results has so far given inconclusive findings. Subject headings: UV radiation - radiative processes - solar radiation (albedo) Código(s) PACS: 61.80.Ba, 92.60.Vb Resumen
Se reportan los resultados finales del proyecto bianual, centrado en dos campañas de medición de la Radiación Ultravioleta (RUV) en el Salar de Uyuni: la primera en mayo de 2008 y la otra en noviembre de 2008. Los datos recogidos en ambas campañas y el análisis respectivo han permitido llegar a importantes resultados en los dos objetivos principales del proyecto: medida del albedo superficial en la banda UV y cuantificación del alcance del mismo en la región circundante al Salar. En el primer objetivo, se han encontrado indicios sólidos de que el albedo depende no sólo de la estación del año, como se podría esperar debido al régimen de lluvias localizado en los meses de enero a marzo, sino también de los cambios climáticos que se están dando en la región, lo que podría conllevar una disminución del albedo, señal del "ensuciamiento" del Salar. Adicionalmente se ha medido el alcance de los efectos de albedo en el Salar, pero el análisis no muestra aún resultados concluyentes. Descriptores: radiación UV - procesos radiativos - radiación solar (albedo)
1 Antecedentes
Figure 1: Arreglo experimental utilizado para medir el albedo. El sistema permite la variación de altura entre los domos y el suelo y, además, permite calibrar la horizontalidad.
El año 2004 se realizó una campaña, junto a colegas de la Universidad e Munich en el Salar de Uyuni para medir el albedo "plano" en UV y realizar mediciones de la irradiancia solar en la banda UV. Si bien esa campaña resultó exitosa en cuanto se pudo medir el albedo y la irradiancia UV en diferentes sitios del Salar, quedaron algunas preguntas abiertas que el presente trabajo intenta responder. En especial nos interesaba conocer el alcance que tiene el albedo UV del Salar en su entorno, o sea hasta dónde afecta a la irradiancia total el cambio de albedo UV del suelo, con la mirada puesta en la problemática del cambio de uso de suelo. El desarrollo, dificultades y resultados de la I campaña en el Salar (mayo de 2008), han sido reportadas en el artículo "Resultados preliminares del proyecto: Medidas de albedo en el Salar de Uyuni" presentado en noviembre de 2009. A partir de los resultados de esa campaña, en el presente artículo analizaremos los datos colectados en la II campaña (noviembre de 2008).
Figure 2: La imagen muestra la dependencia del valor del offset del radiómetro YES_139 con el valor del voltaje proporcionado por la batería de alimentación en un experimento realizado en laboratorio.
2 Objetivos
El proyecto fue planificado para un período de dos años, mediante dos campañas, una exploratoria y otra final. La campaña exploratoria, realizada en mayo de 2008, tenía el objetivo de medir el albedo del Salar en diferentes lugares, diferentes por su distancia a la orilla o por las características de la superficie [42008Zaratti et al.Zaratti, Forno, Gutiérrez, Apaza, & Velarde]. La hipótesis para verificar era que la superficie del Salar no es homogénea y le corresponde más de un valor de albedo. Asimismo, se asume que el albedo depende de la época del año, razón por la cual se planeó repetir la campaña en noviembre, aplicando las lecciones aprendidas. El otro objetivo del proyecto era medir la irradiancia en función de la distancia a la orilla, suponiendo que, a medida que nos adentramos en el Salar la irradiancia aumenta y a medida que nos alejamos hacia el altiplano la irradiancia disminuye. La pregunta interesante es: ?`a qué distancia de la orilla se deja de percibir el efecto de borde?
3 Planificación de la campaña
En la campaña de noviembre de 2008 participaron dos investigadores, un tesista y un técnico y en ella se realizaron varias innovaciones, como consecuencia de los resultados de la campaña de mayo. La lección principal aprendida en la campaña de mayo fue la dificultad de realizar mediciones alimentando los equipos con baterías. En un primer momento se optó por usar baterías de motos, por su fácil manejo y recarga, pero los resultados no fueron buenos, debido al elevado consumo de carga de los equipos durante las más de seis horas de trabajo dentro del Salar, sin posibilidad de recarga. Las consecuencias principales se manifestaban en variaciones de los valores del offset de los radiómetros, en función del voltaje de las baterías.
Figure 3: Variación del Offset y de la temperatura del domo del instrumento YES 139 con el valor del voltaje de la batería.
Figure 4: Datos originales (en milivoltios) tomados por el instrumento YES 139 conectado a un datalogger Campbell que tenía un divisor de voltaje para evitar saturaciones. La curva superior es la irradiancia medida con el domo mirando al cenit, la inferior mirando al nadir.
El equipo instrumental consistía de dos radiómetros de radiación ultravioleta YES UVB-1. Se trata de dos radiómetros de banda ancha (280-320 nm) fabricados por Yankee Environmental System Inc. que miden la radiación solar global UV-B. Usan filtros de vidrio colorado y fósforo UV-B que convierte la radiación UV-B en luz verde, la cual es medida por fotodetectores de estado sólido. Sus principales características de operación son [21997 Yankee Environmental System (YES) ]: respuesta espectral característica (Diffey simulada), respuesta de coseno ±5% (para 0° <SZA,<60° ), tiempo de respuesta ±0.1 s y temperatura de operación entre -40° C y +40° C. El primero de los radiómetros funciona en el LFA desde el año 1998 (YES 138) y otro ha sido recientemente adquirido (llamado YES 139). Asimismo se contó con dos dataloggers Campbell CR-200 que fueron adecuados para funcionar con los radiómetros y dos GPS para registrar la posición de los sitios de medición. Un aspecto preocupante era el desgaste que presentaba el YES 138 después de más de 10 años de operación casi continua y en ambientes extremos, como Chacaltaya. Por esa razón se planificó realizar medidas "en paralelo" de los dos radiómetros, en el sentido que cada instrumento tomara datos independientemente del otro, como se explicará más adelante.
Figure 5: Curva de irradiancia UVB, representada en mV, del YES 139 con el domo mirando al cenit (día 17/11/08) interpolada con un polinomio de cuarto grado.
Figure 6: El equivalente de la Fig. 5 con el domo mirando hacia abajo (día 17/11/08).
4 Realización de la campaña
Las medidas de la segunda campaña se realizaron durante cinco días, del 15 al 19 de noviembre, con el apoyo logístico de una agencia de turismo local. Los cinco días fueron despejados, registrándose la insolación máxima anual, debido al paso del sol por el cenit de esa localidad. La rutina diaria fue ajustada a los resultados del día y a la emergencia de algunos problemas que se presentaron con las medidas. En efecto, los primeros dos días se realizaron medidas de albedo en el mismo sitio de las medidas de albedo de la anterior campaña (20° 27' 41" S, 67° 15' 24" W, 3670 msnm), esto es, casi en el centro del Salar, pero con los dos radiómetros en posición fija (uno mirando al zenit y el otro al nadir), utilizando el mismo arreglo experimental ya probadodo durante la campaña de mayo (ver Fig. 1). El tercer día, debido a algunas dificultades de intercalibración bajo condiciones térmicas y eléctricas diferentes, como se explicará más adelante, se repitieron las medidas alternando la posición de los radiómetros cada 20 minutos. El día martes 18 se realizó el cruce del Salar midiendo la irradiancia RUV mediante uno de los radiómetros en sitios distantes unos 15 km entre sí, a medida que nos adentrábamos en el Salar hasta el centro del mismo. Luego, a la vuelta, se repitieron las medidas hasta unos 10 km desde la orilla en el altiplano. El segundo radiómetro se quedó midiendo en un lugar fijo (el Hotel de Sal a unos 6 km de la orilla), como referencia. Finalmente, el día miércoles 19 se realizó una medida de absorción sobre superficies con diferente inclinación. Cabe añadir que, en cada sitio donde se realizaron medidas, se tomaron muestras de la superficie de sal para ser analizadas posteriormente.
Figure 7: Valor del albedo UVB, medido mediante comparación de las curvas de irradiancia con domo mirando hacia arriba (eje abscisas) y hacia abajo (eje ordenadas) para el día 17/11/08. Los valores, en mV, representan la mitad del valor real.
Figure 8: Curva del albedo UVB instantáneo del YES 139 (día 17/11/09).
5 Resultados de las medidas de albedo
Con base a las lecciones aprendidas de la primera campaña, se puso particular atención a los problemas del offset (y de corriente oscura) de cada instrumento. En efecto, el estudio realizado en el LFA ha mostrado la dependencia del offset del voltaje de la batería, siendo estable entre 12 y 12.5 V, condiciones que no siempre se puede mantener cuando se trabaja en el campo. Asimismo se puso en evidencia un transitorio en el tiempo de estabilización de la temperatura del sensor de aproximadamente 20 minutos, lo que invalidaba algunos datos que se tomaron durante ese tiempo, después de conectar los instrumentos a la batería (ver Figs. 2 y 3). En particular, se llegó a la conclusión de que era preferible utilizar medidas de cada instrumento por separado (como se hizo el día 17) que intentar corregir y recalcular los datos con radiómetros en posición fija (datos de los días anteriores). Por tanto, sobre la base de los datos originales, mostrados en la Fig. 4, se procedió a reconstruir, por interpolación las curvas con el domo arriba y con el domo abajo.
Figure 9: Variación temporal del albedo medido en la campaña de mayo de 2008. Es evidente el desgaste del instrumento YES 138.
Figure 10: Intercomparación de las señales de los domos YES 139 (en dirección al nadir) y YES 137 (en dirección al cenit) durante la campaña de mayo 2008. La pendiente de la recta de regresión nos muestra el albedo del lugar.
Para simplificar el cálculo se normalizaron los datos de ambas series, se hizo un ajuste de curva de cuarto grado, lo más apropiado en torno al mediodía, se interpolaron las curvas y se volvió a los datos en mV. Las figuras siguientes muestran el resultado final. Una vez obtenidas las curvas anteriores, sin necesidad de mayores ajustes, tratándose del mismo instrumento, se calculó el valor del albedo en UV con diferentes métodos. El primero consiste en graficar las dos curvas una contra la otra y deducir, de la pendiente, el valor del albedo A. De ese modo, se halla el valor: A= 0.387 ± 0.001 (Fig. 7). Si se obligara a la recta de ajuste a pasar por el origen, lo que parece lógico tratándose de un mismo instrumento, el valor del albedo resultaría: A'= 43 %. El otro método consiste en calcular el valor instantáneo del albedo, para cada dato, y luego calcular el valor medio. La Fig. 8 muestra los valores del albedo "instantáneo" medidos. En el intervalo temporal 10am - 4pm, las horas centrales del día, el albedo resulta A"= 41.0% ± 1.2%. En resumen, el albedo medido con el radiómetro YES 139 está entre 39 y 41 %. Cuando analizamos los datos del YES 138, tomados al mismo tiempo que los del YES 139 pero en sentido invertido, obtenemos los valores respectivos: A = 39% (R2 =1). A' = 47.6% (R2 = 0.92) entre 10am y 4pm. A" = (49 ± 3)% entre 10am y 4 pm. Las significativas diferencias en el albedo instantáneo del YES 138, se deben posiblemente a variaciones de valores del offset con la temperatura y el voltaje de la batería a lo largo del día, pero es importante recordar, además, que el YES 138 se ha demostrado muy poco confiable, debido a su largo tiempo de vida (10 años). Una prueba de esta afirmación es la curva temporal del albedo del YES 138, comparada con la equivalente del YES 139 (Fig. 9).
6 Discusión de las medidas de albedo
Como se mencionó en la sección 4, en mayo y en noviembre de 2008 las medidas de albedo se realizaron en el mismo sitio, localizado con los GPS y con señales dejadas en el lugar. La Fig. 10 muestra el resultado obtenido con dos radiómetros, el 137 y el 139, en esa campaña. Además se cuenta con las medidas de albedo realizadas, con otra clase de radiómetro (Scie-tech), en mayo del año 2005 [32007Reuter et al.Reuter, Ghezzi, Palenque, Torrez, Andrade, & Zaratti]. Por tanto, podemos intentar una comparación somera entre las tres medidas, de acuerdo a la Tabla 1:
Table 1: Valores del albedo medidos en las tres campañas realizadas por el LFA-UMSA.
Fecha Albedo medido mayo 2005 69% mayo 2008 56% noviembre 2008 39% Figure 11: Variación temporal del albedo medido en la campaña de mayo de 2005 en torno al medio día (cero en las abscisas). La asimetría cenital es despreciable.
La Tabla 1 nos muestra dos fenómenos. Por un lado no debería sorprender que el albedo en noviembre sea inferior al albedo en mayo, por la simple razón que en mayo el Salar resiente todavía de la época de lluvias, mientras en noviembre la superficie viene de soportar más de seis meses de sequía y vientos que suelen transportar polvo y tierra hacia el Salar desde el Altiplano. Desde luego, el exacto valor de esa disminución del albedo dependerá de las particularidades de la estación de lluvias primero y de vientos después. En efecto, cuando se analizan cortes de la costra del Salar es posible reconocer diferencias en la acumulación de tierra en el subsuelo del Salar, dependientes del clima. Más interesante es la disminución de un 20% del albedo en el mes de mayo a distancia de tres años. Si bien se trata de instrumentos diferente, y por tanto con una diferente respuesta coseno y espectral, es posible asumir que esas diferencias no son suficientemente significativas como para explicar las diferencias de los albedos, particularmente para pequeños ángulos cenitales, como a los que se hicieron las medidas en el mes de noviembre. Asimismo los lugares de observación fueron distintos; sin embargo, en vista de que el año 2005 se realizó la medida a menso de 10 km de la orilla del Salar (y, por tanto, dentro del rango de influencia de la frontera Salar-Altiplano), es más llamativo que esa medida resulte bastante mayor a la realizada en mayo de 2008. Por tanto, sin descartar fluctuaciones locales del clima, es posible avanzar la hipótesis, con cargo a comprobación futura, de un real "ensuciamiento", o obscurecimiento, del Salar, debido no tanto, como a veces se escucha, a los efectos del turismo masivo sino a un verdadero cambio climático de la región que se manifiesta, de acuerdo a algunos datos meteorológicos del Servicio Nacional de Meteorología e Hidrología (SENAMHI) en lluvias más concentradas en los primeros tres meses del años, con períodos de sequía más prolongados. De hecho los resultados de los modelos climáticos más acreditados apuntan a patrones climáticos como el que acabamos de describir (mayor duración de la época seca) en la región de los Andes Occidentales [12007Magrin et al.Magrin, Gay García , Cruz Choque , Giménez, Moreno, Nagy, Nobre, & Villamizar].
7 Medidas del alcance del albedo
El otro objetivo del proyecto era medir el alcance de la influencia de una región con determinado albedo (el altiplano) sobre la irradiancia medida a diferente distancia de aquella en una región de diferente albedo (Salar). Para ese fin se diseñó un experimento para medir, estáticamente, la irradiancia durante períodos de 20 minutos en diferentes puntos dentro del Salar (en lo posible sobre el mismo paralelo, para evitar correcciones de latitud), y así poner en evidencia el porcentaje de disminución de la irradiancia al acercarse a la orilla del Salar. Para fines de comparación, uno de los radiómetros (YES 138) se quedó midiendo en un punto fijo, cerca del Hotel de Sal, a 6 km de la orilla. El análisis de los resultados se ha revelado más complejo de lo esperado y requiere más tiempo y dedicación. De hecho, las diferencias que se quieren medir no son muy significativas, debido a que están tomadas cerca al mediodía y las correcciones temporales son bastante delicadas. Adicionalmente se tuvieron algunos percances experimentales que obligaron a un trabajo aún más complejo de reconstrucción de la curva con base a los valores del día anterior. Por todas estas razones estimamos que este tema merece un artículo a parte, el cual ya está bastante avanzado.
8 Conclusiones
Las dos campañas, llevadas a cabo en el marco del proyecto bianual, han permitido aprender algunas lecciones y han logrado importantes resultados. Las lecciones se refieren principalmente a la dificultad, no siempre considerada, de realizar medidas ópticas en condiciones extremas, con fuentes de potencia no confiables, con instrumentos que requieren particular cuidado, debido a su tiempo de uso y al tiempo limitado que se dispone. No obstante, se han logrado resultados que ponen en discusión anteriores conclusiones, por la dinámica climática de la región del Salar. El principal es el ensuciamiento progresivo del Salar, debido posiblemente a la acumulación de polvo por la reducción del período anual de lluvias, aspecto que deberá ser confirmado por el análisis de datos meteorológicos. Desafortunadamente la estación situada en la isla central del Salar (Incahuasi) ha estado parada y abandonada en los últimos años, pero se tienen datos de la estación de Uyuni. En todo caso, los resultados obtenidos, alientan al Laboratorio de Física de la Atmósfera a emprender otra campaña en un futuro próximo para confirmar y cuantificar las anomalías encontradas.
Agradecimientos
Nuestros sinceros agradecimientos al personal del Instituto de Investigaciones Físicas de la Universidad Mayor de San Andrés (UMSA) por el apoyo generoso en la ejecución de las campañas, empezando por el Director Dr. Wilfredo Tavera. Asimismo agradecemos el apoyo logístico proporcionado satisfactoriamente por las empresas Peru Bolivian Tours y Licancabur. Finalmente, agradecemos a la Empresa Ferroviaria Andina por el transporte gratuito de los equipos a Uyuni.
References
- [12007Magrin et al.Magrin, Gay García , Cruz Choque , Giménez, Moreno, Nagy, Nobre, & Villamizar]
-
1.- Magrin, G., Gay García, C., Cruz Choque, D., Giménez, J., Moreno, A., Nagy, G., Nobre, C., & Villamizar, A. 2007, Latin America. Climate Change 2007: Impacts, Adaptation and Vulnerability. Contribution of Working Group II to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change (M.L. Parry, O.F. Canziani, J.P. Palutikof, P.J. van der Linden and C.E. Hanson, Eds., Cambridge University Press, Cambridge, UK, 581-615) [ Links ]
- [21997 Yankee Environmental System (YES) ]
-
2.- Yankee Environmental System (YES). 1997, UVB-1 Ultraviolet Pyranometer (Installation and User Guide, Version 2.0) [ Links ]
- [32007Reuter et al.Reuter, Ghezzi, Palenque, Torrez, Andrade, & Zaratti]
-
3.- Reuter, J., Ghezzi, F., Palenque, E., Torrez, R., Andrade, M., & Zaratti, F. 2007, Journal of Photochemistry and Photobiology B: Biology, 87, 1 [ Links ]
- [42008Zaratti et al.Zaratti, Forno, Gutiérrez, Apaza, & Velarde]
-
4.- Zaratti, F., Forno, R., Gutiérrez, G., Apaza, R., & Velarde, F. 2008, Resultados preliminares del proyecto Medidas de albedo en el Salar de Uyuni (IIF-UMSA, Informe gestión 2008) [ Links ]
File translated from TEX by TTH, version 3.89.
^rND^sMagrin^nG.^rND^sGay García^nC.^rND^sCruz Choque^nD.^rND^sGiménez^nJ.^rND^sMoreno^nA.^rND^sNagy^nG.^rND^sNobre^nC.^rND^sVillamizar^nA.^rND^sReuter^nJ.^rND^sGhezzi^nF.^rND^sPalenque^nE.^rND^sTorrez^nR.^rND^sAndrade^nM.^rND^sZaratti^nF.^rND^sZaratti^nF.^rND^sForno^nR.^rND^sGutiérrez^nG.^rND^sApaza^nR.^rND^sVelarde^nF.^rND^sWinkelmann^nR.^rND^sWinkelmann^nR.^rND^sWinkelmann^nR
On 8 Nov 2010, 09:35.EFECTO DEL CAMPO EL\'{E}CTRICO ATMOSF\'ERICO SOBRE EL CONTEO DE PART\'ICULAS CARGADAS DETECTADAS POR UN CENTELLADOR EN EL LABORATORIO DE F\'ISICA C\'OSMICA DE CHACALTAYA\ \vspace{1em}{ The Effect of the Atmospheric Electric Field on the Counting of Charged Particles Detected by a Scintillator in the Laboratory of Cosmic Ray Physics at Mount Chacaltaya (Bolivia)} EFECTO DEL CAMPO ELÉCTRICO ATMOSFÉRICO SOBRE EL CONTEO DE PARTÍCULAS CARGADAS DETECTADAS POR UN CENTELLADOR EN EL LABORATORIO DE FÍSICA CÓSMICA DE CHACALTAYA
The Effect of the Atmospheric Electric Field on the Counting of Charged Particles Detected by a Scintillator in the Laboratory of Cosmic Ray Physics at Mount Chacaltaya (Bolivia)Winkelmann R.
Carrera de Física
Universidad Mayor de San Andrés
La Paz-Bolivia
(Recibido el 30 de julio de 2010; aceptado el 15 de agosto de 2010)
Abstract
Secondary particles produced by primary cosmic rays are affected by diverse atmospheric factors such as pressure and temperature. The Atmospheric Electric Field is one such factor. There is evidence that during electric storms the counting rate of particles can increase considerably. In this work we look at the relationship between the Atmospheric Electric Field and the secondary particles, at the Chacaltaya Cosmic Ray Physics Laboratory, under two distinct weather conditions: good and unsettled weather.
Subject headings: atmospheric electric field - cosmic rays Código(s) PACS: 92.60.Pw, 95.50.S-, 95.45._i
Resumen
Las partículas secundarias que son producidas por rayos cósmicos primarios se ven afectadas por diversos factores en la atmósfera, como presión, temperatura, y otros. El Campo Eléctrico Atmosférico (CEA) es uno de estos factores. Se observa que durante tormentas eléctricas la tasa de conteo de las partículas puede incrementarse considerablemente. En el presente trabajo se busca observar en el Laboratorio de Física Cósmica de Chacaltaya (5220 msnm) la relación existente entre el CEA y las partículas secundarias producidas por rayos cósmicos primarios, bajo dos distintas condiciones climáticas: buen tiempo meteorológico y tiempo meteorológico perturbado.Descriptores: campo eléctrico atmosférico - rayos cósmicos.
1 Introducción
El conjunto de partículas generadas por un rayo cósmico primario posee tres componentes: hadrónica, electromagnética y muónica; éste conjunto de partículas se ven inmersas en una serie de condiciones atmosféricas, como modulación de campo eléctrico [52000Feynman et al.Feynman, Leighton, & Sands], presión y temperatura [32008Alvarez] entre otros, presentes desde el momento de su formación, pero no todas las partículas se verán afectadas de la misma forma. La idea de una posible influencia del Campo Eléctrico Atmosférico (CEA) durante tormentas eléctricas sobre partículas cargadas, en principio electrones, fue planteada por primera vez en la década de los 20 [111925Wilson]; durante las siguientes décadas varios fueron los intentos por observar una aceleración en partículas cargadas antes y/o durante tormentas eléctricas, los resultados de dichas observaciones fueron contradictorios y nada claros. En 1985 Alexeenko et al. realiza la primera inspección minuciosa entre medidas del CEA y variaciones en el registro de rayos cósmicos encontrando correlaciones de corta duración (alrededor de 8 a 16 minutos) con un incremento en la tasa de conteo en el registro de rayos cósmicos de 0.2 a 0.3%; trabajos posteriores como el de [11999Aglietta et al.], [42000Brunetti et al.] y [102001Vernetto et al.], entre muchos otros, demostraron que el CEA afecta tanto a muones como electrones de carga positiva y negativa, además de mostrar un incremento en la radiación X y γ. Vernetto reporta un incremento en la tasa de conteo de un Chubasco Atmosférico Extenso (EAS por sus siglas en inglés) de 10 a 15% en un periodo de 10 a 20 minutos además, [21985Alexeenko et al.] reportan incrementos en la tasa de conteo de hasta un 20% antes de una descarga eléctrica.
Figure 1: Datos correspondientes al 06 de Octubre. (a) Registro al segundo de un día para el centellador. (b) Registro al segundo de un día para el CEA. (c) Registro del centellador promediado al minuto, la curva negra corresponde a los datos alisados. (d) Registro del CEA promediado al minuto, la curva negra corresponde a los datos alisados.
La generación del CEA se debe a la diferencia de potencial existente entre la Tierra y la Electrósfera, otras influencias atmosféricas como las nubes cargadas eléctricamente son capaces de producir intensos campos eléctricos. El modelo clásico de la distribución de cargas en las nubes establece que ésta puede llegar a tener hasta tres centros de carga: uno ubicado en la cima de la nube, con una carga de +40 coulombs, alrededor de los 10 a 12 Km desde la superficie terrestre, una segunda concentración de carga, con -40 coulombs, ubicada alrededor de los 5 a 7 Km, también medida desde la superficie, y una tercera carga no tan significativa que puede o no estar presente, con +3 coulombs, ubicada en la base de la nube alrededor de los 2 Km, éstos valores solo son referenciales, ya que dependerán del tipo de nube, su desarrollo, topología del terreno, etc. Las partículas que se ven inmersas en un CEA son aceleradas bajo el mecanismo Electrón Fugitivo (Runaway Breakdown) propuesto por [61999Gurevich et al.Gurevich, Zybin, & Roussel], el cual establece que las partículas que son aceleradas dentro del campo eléctrico atmosférico experimentan un incremento en su energía y pueden llegar a producir nuevas partículas mediante interacciones con núcleos de la atmósfera, a su vez éstas producen nuevas partículas por el mismo mecanismo hasta que las últimas partículas no tengan la energía suficiente para continuar con el proceso, de manera que el número de partículas aumenta de forma exponencial. El presente trabajo esta destinado a realizar un estudio de la influencia del CEA en el conteo de partículas generadas a partir de rayos cósmicos detectados en el Laboratorio de Física Cósmica de Chacaltaya (ubicado a 5230 m.s.n.m.) por un centellador, durante dos diferentes escenarios climáticos: días con buen tiempo meteorológico y días con tiempo meteorológico perturbado por nubes, en un período de tres meses.
2 Detectores de Partículas en Chacaltaya
Existen diversos arreglos experimentales instalados en el Laboratorio de Física Cósmica de Chacaltaya destinados a detectar partículas secundarias, las cuales, al igual que en cualquier parte del mundo, son afectadas por factores atmosféricos. Los registros de partículas cargadas son constantemente corregidos por presión y temperatura, ya que se cuenta con instrumentos para realizar éstas medidas. Variaciones en el conteo de partículas durante tormentas eléctricas son observadas constantemente, algunas de ellas son reportadas en [72002Huaygua et al.Huaygua, Velarde, & Saavedra], encontrando un número estadísticamente significativo de tales eventos mediante correlaciones con equipos improvisados para detectar descargas eléctricas atmosféricas; no es sino hasta Septiembre del año 2008 que el Laboratorio empieza a registrar medidas del CEA con un instrumento específico para éste propósito, es así que por primera vez, en éste trabajo, se realiza un estudio de las variaciones en el registro de partículas acompañadas por medidas del CEA.
Figure 2: Modulaciones promedio. (a) Campo eléctrico atmosférico. (b) Cuentas de partículas.
Figure 3: (a) Correlación promedio con los datos de las curvas alisadas. (b) Ajuste cuadrático en la región de descenso del número de partículas.
3 Obtención y Metodología de Datos
En el Laboratorio de Física Cósmica de Chacaltaya se encuentran instalados los equipos utilizados en éste trabajo. El primero, un sensor de campo eléctrico atmosférico de la marca Boltek, modelo EFM- 100 con un tiempo de respuesta de 0.1 segundos, un rango de resolución digital de 0.01 kV/m, y un rango de medición desde -20 kV/m hasta +20 kV/m, sí los valores medidos superan éstos límites el sensor se satura; para evitar saturaciones es posible disminuir la sensibilidad del EFM- 100 a través de resistencias con valores adecuados. Debido al intenso campo eléctrico registrado durante las pruebas preliminares a la instalación del EFM- 100, se opta por reducir la sensibilidad del mismo a un 25% del valor original. El segundo es un detector de partículas compuesto por un centellador de plástico de 0.25 m2 de área y 0.1 m de espesor. Como se sabe, los centelladores sólo detectas partículas cargadas, fundamentalmente electrones y muones, por lo tanto resulta muy adecuado para nuestro propósito.
3.1 Selección de datos
Los datos utilizados para éste trabajo son los registrados durante los meses de Septiembre a Noviembre de 2008; se seleccionan y separan primero los registros diarios del CEA que presenten una modulación suave y de baja intensidad, luego se escogen días en los cuales la intensidad del CEA muestre alteraciones bruscas y saturaciones en el equipo. El primer caso es denominado "días con buen tiempo meteorológico", ya que presentan un cielo con poco o nada de nubosidad, esto se corrobora con imágenes satelitales obtenidas en [82010NASA]. El segundo grupo es llamado "días con tiempo meteorológico perturbado" por tratarse de días con alteraciones eléctricas en la atmósfera debidas a la presencia de nubes cargadas eléctricamente, precipitaciones, rayos y otros, así, el primer grupo suma un total de 14 casos, en cambio, el segundo tiene 24 casos.
Figure 4: Variación porcentual en el registro del centellador vs. variación porcentual del CEA. Cada punto representa un día.
Figure 5: Incremento en el número de partículas registrado el 07 de Octubre de 2008. Las líneas verticales muestran la coincidencia temporal entre ambos registros: CEA (a) y N (b).
3.2 Metodología
3.2.1 Días con buen tiempo meteorológico
Los datos obtenidos por ambos detectores son registrados al segundo, lo que hace un total de 86400 valores al día para cada detector. La Fig. 1 (a) y (c) muestran dicho registro para el centellador y el CEA respectivamente (los datos corresponden al 06 de Octubre). En el caso del centellador (Fig. 1(a)) puede apreciarse tan sólo una franja, sin una modulación, sin embargo, si se toma el promedio de los valores sobre un minuto reducimos el número de registros a 1440 y posteriormente se realiza un alisado de la curva a través de una media móvil, aparece claramente una modulación. La Fig. 1(b) muestra el resultado de dicho proceso, donde la curva negra es el alisado del registro promediado. Un proceso idéntico se lleva acabo para el caso del CEA, la Fig. 1(d) muestra el resultado para este caso. Una vez que se obtienen las curvas alisadas de ambos registros para los 14 casos se realiza una gráfica promedio de ambas modulaciones (Fig. 2). Se ve que ambas modulaciones tienen inicio a la misma hora aproximadamente, entre las 7 y 8 de la mañana hora local.
Figure 6: Dos incrementos notables en el registro de partículas coincidentes con alteraciones en el CEA registrados el 13 de Octubre de 2008. Las líneas verticales muestran la coincidencia temporal entre ambos registros: CEA (a) y N (b).
Figure 7: Variaciones notables en el registro de partículas no coincidentes temporalmente con alteraciones en el CEA registrados el 10 de Noviembre de 2008. El desfasaje es de aproximadamente de siete horas. Las líneas verticales muestran el tiempo aproximado de duración de las perturbaciones en cada registro: CEA (a) y N (b).
A continuación se realiza una correlación entre ambas modulaciones (Fig. 2), la Fig. 3(a) muestra el resultado, donde es posible distinguir el ciclo debido a que ambas modulaciones son opuestas, mientras una se incrementa, la otra disminuye y viceversa. En la etapa en que los valores del registro del centellador disminuyen, es posible realizar un ajuste cuadrático (Fig. 3(b)) con un coeficiente de correlación de 0.99, esta correlación entre ambos registros se mantiene por el lapso de 7 horas aproximadamente.
Figure 8: Variación de corta duración registrada el 18 de Octubre de 2008. (a) Intensas variaciones en el registro del CEA. (b) Registro de partículas con dos variaciones de corta duración.
3.2.2 Días con tiempo meteorológico perturbado
Primero realizamos un análisis similar al caso anterior, reducimos el número de datos al minuto, luego se restan las modulaciones diarias obtenidas para el caso de días con buen tiempo meteorológico para así obtener ambos registros libres de las modulaciones. Se encuentran tres casos en los que destacan incrementos de larga duración en el registro de partículas acompañado de intensos y bruscos cambios en la intensidad del CEA, sin embargo, el procedimiento descrito anteriormente hace imperceptibles las variaciones de corta duración en ambos registros, sobre todo en el de partículas, por lo que para detectar éstas pequeñas variaciones se trabaja con los registros tomados al segundo. Para encontrar éstas variaciones de corta duración se analizan las variaciones porcentuales en el registro de partículas que coincidan temporalmente con saturaciones prolongadas o súbitos cambios en la intensidad del CEA. Analizar las variaciones porcentuales no basta para poder identificar éstos incrementos de corta duración, por lo que se realizó una media móvil al registro de partículas, esto ayuda a minimizar las fluctuaciones, poniendo en evidencia los incrementos de corta duración. Las gráficas correspondientes a éstos análisis se muestran en la siguiente sección.
4 Resultados
4.1 Días con buen tiempo meteorológico
Para verificar una variación en el número de cuentas de las partículas se restan las modulaciones diarias de ambos registros, una vez hecho esto es posible analizar las variaciones porcentuales máximas. La Fig. 4 muestra el resultado de dicho tratamiento, además, pone en evidencia la independencia de ambas modulaciones, es decir, no importa cuanto pueda crecer el CEA, como puede verse hasta casi un 200% respecto al mínimo diario, las partículas no sienten un incremento considerable, siempre están entre 2 y 5% aproximadamente. Lo que podría explicarse porque la intensidad del CEA en días con buen tiempo meteorológico no es lo suficientemente grande como para iniciar el proceso del Electrón Fugitivo.
4.2 Días con tiempo meteorológico perturbado
Como resultado del análisis de variaciones de larga duración con tiempo meteorológico perturbado se obtuvieron 3 casos. El primero de ellos (Fig. 5) muestra un salto desde valores negativos hasta valores positivos respecto a la media en el registro de partículas que coincide con un descenso en el valor del CEA desde valores mayores a los 80 kV/m hasta un mínimo de 40 kV/m aproximadamente, con una duración cerca de 30 minutos empezando alrededor de las 13 horas. El siguiente caso (Fig. 6) muestra dos notables variaciones del CEA que coinciden también con dos incrementos en el registro de partículas. La primera variación que alcanza un poco más de 20 kV/m con una duración de 3 horas aproximadamente que coincide con dos pequeños picos en el registro de partículas, pero que no se ajustan del todo, temporalmente hablando, a los picos del CEA. La segunda variación en el CEA alcanza los 33 kV/m, dura alrededor de 2 horas y también coincide con incrementos en el registro de partículas. Cabe destacar que el registro del CEA no sufre saturaciones, es decir, no llega al máximo valor del sensor que es de 80 kV/m, y aún así existe incrementos notables en el registro de partículas. El registro realizado el 11 de Noviembre de 2008 (Fig. 7) muestra un CEA sumamente perturbado, desde variaciones de corta duración, presumiblemente causadas por descargas eléctricas, hasta saturaciones prolongadas, siendo la más larga cerca de 2 horas. A diferencia de los dos casos anteriores, la variación en el registro de partículas no es coincidente en tiempo, sino que después de varias horas de haber ocurrido las alteraciones en el CEA las partículas empiezan a sentir este efecto, y les toma alrededor de una hora más regresar al número de partículas promedio, es decir que la perturbación del CEA dura cerca de 10 horas, en cambio el incremento en el número de partículas dura aproximadamente 11 horas.
Figure 9: Variación de corta duración registrada el 13 de Noviembre de 2008. (a) Intensas variaciones en el registro del CEA. (b) Registro de partículas con dos variaciones de corta duración.
A continuación se muestran los casos con variaciones de corta duración, en los cuales la media móvil, que es la línea obscura y con poca variación colocada sobre el registro de partículas, ayuda a identificar los incrementos reales en el registro de partículas. El primero de éstos casos es el registro tomado el día 18 de Octubre de 2008 (Fig. 8), en el cual podemos destacar las intensas variaciones y saturaciones prolongadas en el registro del CEA (Fig. 8(a)), al mismo tiempo encontramos dos incrementos en el registro de partículas (Fig. 8(b)), el primero con una duración de 5 minutos aproximadamente y un incremento del 28.9 %, y el segundo, con una duración de alrededor de 15 minutos y un incremento del 33.4 %. El segundo caso fue registrado el 13 de Noviembre de 2008 (Fig. 9). En el panel superior podemos ver saturaciones en el registro que cambian rápidamente de signo, claro ejemplo de descargas eléctricas, existen tres secciones que podemos identificar con claridad, la primera que ocurre alrededor de las 11 horas, una segunda que ocurre entre las 13 y 16 horas aproximadamente, y la tercera alrededor de las 18 horas. En cambio, en el panel inferior, que es el registro de partículas, podemos observar tan sólo dos alteraciones que corresponden temporalmente a la segunda variación del CEA entre las 13 y 16 horas, el primero de ellos logra una variación del 31.2 % y dura aproximadamente 10 minutos, el segundo, un poco más intenso y duradero, logra un incremento del 36.1 % por un tiempo no menor a los 20 minutos. Finalmente, el registro realizado el 19 de Noviembre de 2008 muestra un caso similar al anterior, en el registro del CEA tenemos dos regiones con cambios súbitos, el primero entre las 10 y 13 horas, y el segundo, menos intenso, entre las 16 y 17 horas, mientras que el registro de partículas (Fig. 10(b)) muestra solo una variación entre las 16 y 17 horas, en otras palabras, ésta variación coincide con la variación del CEA menos intensa. La variación en el registro de partículas posee un incremento del 29.4 % y dura al rededor de 10 minutos.
5 Conclusiones
Se logra realizar un primer estudio del CEA en el Laboratorio de Física Cósmica de Chacaltaya encontrando un perfil promedio para días con buen tiempo meteorológico (Fig. 2(a)), encontrando su valor máximo cerca a las 15 horas en tiempo local, que es un resultado coincidente con el expuesto en [52000Feynman et al.Feynman, Leighton, & Sands], el cual establece que sin importar la ubicación geográfica la máxima variación del potencial atmosférico en días con buen tiempo meteorológico ocurre a las 19 horas del meridiano de Greenwich, además, Chacaltaya es un lugar bastante particular en este sentido, ya que el mínimo valor para días con buen tiempo meteorológico es del orden de 3 kV/m, a diferencia de los 0.1 kV/m mostrados en [92007Ramachandran et al.] y [52000Feynman et al.Feynman, Leighton, & Sands]. Esta característica es debida principalmente a dos peculiaridades que posee el laboratorio: la primera se debe a la altura a la que se encuentra Chacaltaya (5230 m.s.n.m.), además de tratarse de una montaña, y la segunda se debe a la particular ubicación del sensor EFM -100, ya que éste se encuentra en una torre que a su vez se halla sobre depósitos de galena que intensifican el CEA, lo cual pudo comprobarse durante la calibración del sensor.
Figure 10: Variación de corta duración registrada el 19 de Noviembre de 2008. (a) Variaciones en el registro del CEA con dos regiones muy bien diferenciadas. (b) Registro de partículas con tan solo una variación de corta duración.
Por otra parte, el análisis realizado para encontrar alguna correlación entre el CEA y el registro de partículas durante días con buen tiempo meteorológico reveló que ambas modulaciones diarias se originan al mismo tiempo y mantienen una alta correlación por varias horas (en promedio un coeficiente de correlación cuadrático de 0.99 por un lapso de aproximadamente 7 horas). Sin embargo, el resultado final de éste análisis lo muestra la Fig. 4, en el cual se aprecia una completa independencia entre ambas modulaciones. Para el caso de días con tiempo meteorológico perturbado se encuentran dos tipos de resultados, variaciones de larga y corta duración, que precisamente concuerdan con los reportados en [11999Aglietta et al.]. De 24 casos estudiados, tan solo 4 muestran variaciones de larga duración, siendo éstos poco homogéneos en su comportamiento, es decir, los incrementos en el registro de partículas no siempre ocurren simultáneamente a las perturbaciones en el CEA, como es el caso del registro realizado el 10 de Noviembre de 2008 (Fig. 7), el cual muestra un retraso de varias horas en el incremento del registro de partículas comparado al tiempo de inicio en las perturbaciones en el CEA. Para el caso de variaciones de corta duración se encontraron registros de 3 días que concuerdan con éste perfil, los cuales revelan un incremento en la razón de conteo de entre 29 y 36% con una duración de 5 a 20 minutos aproximadamente. Al parecer las variaciones de corta duración siempre son coincidentes, temporalmente hablando, con bruscos incrementos y saturaciones en el registro del CEA. Sin embargo, no siempre se presentan cuando el CEA sufre intensas y prolongadas fluctuaciones, claro ejemplo de ello son los registros realizados los días 13 y 19 de Noviembre (Figs. 9 y 10), el registro del CEA muestra marcadas regiones con intensas perturbaciones, pero el incremento en la razón de conteo de las partículas solo es coincidente con algunas de éstas regiones. En conclusión podemos decir que los incrementos en las razones de conteos durante perturbaciones en el CEA, no son un efecto que siempre se presente, ni siquiera que tengan duraciones o intensidades iguales en todos los casos. Una razón para que éste fenómeno no sea constante podría radicar en el detector de centelleo, ya que al ser de una superficie pequeña la fluctuacíon sobre la media es muy grande (de más de 4σ), por lo que las variaciones debidas al CEA podrían estar inmersas dentro de las fluctuaciones, es por esto que un futuro análisis debe usar un arreglo de detectores para el registro de partículas.
References
- [11999Aglietta et al.]
-
1.- Aglietta, M. et al. 1999, Proc. 26-th ICRC, 7, 351 [ Links ]
- [21985Alexeenko et al.]
-
2.- Alexeenko, V. V. et al. 1985, Proc. 19-th ICRC, 5, 352 [ Links ]
- [32008Alvarez]
-
3.- Alvarez, J. 2008, Efecto de las Tormentas Eléctricas en los Rayos Cósmicos Detectados en la Superficie Terrestre (Tesis de Maestría, UNAM) [ Links ]
- [42000Brunetti et al.]
-
4.- Brunetti, M. et al. 2000, Geophys. Res. Letter, 27, 1599 [ Links ]
- [52000Feynman et al.Feynman, Leighton, & Sands]
-
5.- Feynman, Leighton, & Sands. 2000, Física Vol II (Addison-Wesley-Longman) [ Links ]
- [61999Gurevich et al.Gurevich, Zybin, & Roussel]
-
6.- Gurevich, A. V., Zybin, K. P., & Roussel, R. A. 1999, Dupre Phys.Lett. A, 254, 79 [ Links ]
- [72002Huaygua et al.Huaygua, Velarde, & Saavedra]
-
7.- Huaygua, N., Velarde, A., & Saavedra, O. 2002, Revista Boliviana de Física, 8, 23 [ Links ]
- [82010NASA]
-
8.- NASA. 2010, http://rapidfire.sci.gsfc.nasa.gov [ Links ]
- [92007Ramachandran et al.]
-
9.- Ramachandran, V. et al. 2007, The South Pacific Journal of Natural Science, 12, 70 [ Links ]
- [102001Vernetto et al.]
-
10.- Vernetto, S. et al. 2001, Proc. 27-th ICRC, 10, 4165 [ Links ]
- [111925Wilson]
-
11.- Wilson, C. T. 1925, Proc. Cambridge Phil. Soc., 22, 534 [ Links ]
File translated from TEX by TTH, version 3.89.
On 8 Nov 2010, 09:33.