Uma abordagem diferente na estatística do experimento Millikan

A different approach to the statistics of the Millikan experiment

ARTIGOS GERAIS

Uma abordagem diferente na estatística do experimento Millikan

A different approach to the statistics of the Millikan experiment

N. VeissidI, ¹; L.A. Pereira^I; A.F.V. Peña^II

^ILaboratório Associado de Sensores e Materiais, Instituto Nacional de Pesquisas Espaciais, São José dos Campos, SP, Brasil

^IIDepartamento de Física, Química e Biologia, Faculdade de Ciências e Tecnologia, Universidade Estadual Paulista "Júlio de Mesquita Filho", Presidente Prudente, SP, Brasil

RESUMO

A experiência clássica de Millikan para a determinação da carga do elétron é muito interessante do ponto de vista de estudos estatísticos de medidas experimentais. Em uma experiência didática de laboratório de física, é sugerido dividir os valores da carga elétrica obtida pelo valor da carga do elétron e, desta forma, estabelecer experimentalmente que as gotas de óleo carregam números inteiros de elétrons. Este trabalho apresenta um algoritmo estatístico novo para tratamento das medidas e para o cálculo do melhor valor da carga do elétron em um experimento Millikan.

Palavras-chave: experimento Millikan, estatística, carga do elétron.

ABSTRACT

The classic Millikan experiment to determine the electron charge is very interesting from the point of view of statistics of experimental measurements. In a physics laboratory it is suggested to divide the obtained values of the electric charge by the value of an electron charge, and it is possible to experimentally establish that the oil droplets carry integer numbers of electrons. This work presents a new statistical algorithm for the treatment of measurements and for determining the best value of the electron charge in a Millikan experiment.

Keywords: Millikan experiment, statistics, electron charge.

1 Introdução

Em 1909, R.A. Millikan realizou uma série de experimentos que demonstraram que as cargas elétricas sempre acontecem em múltiplos inteiros de uma unidade elementar. Millikan também determinou o valor desta carga com uma precisão de 0,1%. Atualmente o experimento de Millikan é um sistema borrifador de gotículas de óleo ionizadas dentro de um campo elétrico orientado em sentidos alternados. A velocidade de ascensão de uma determinada gota de óleo dentro de um meio viscoso permite determinar a carga elétrica que a gota carrega. Ao mudar a polaridade do campo elétrico entre as placas, é possível inverter o sentido de movimento de uma gota, permitindo que a gota analisada fique na visada do observador por tempo suficiente para realizar várias medidas. Supondo apenas que a velocidade terminal da gota é proporcional à força agindo sobre ela é possível verificar que as cargas sempre ocorrem em múltiplos de uma unidade fundamental, cujo valor real é 1,602 x 10^-19 C [1].

O livro de Crease [2] relata várias curiosidades a respeito do experimento. Por exemplo, a defesa do físico Ehrenhaft a respeito da existência de subelétrons baseado nos dados do próprio Millikan. Outra curiosidade foi a indicação do eminente físico Ernest Rutherford que a maior dificuldade do experimento residia na taxa de evaporação das gotículas de água. Desta forma, para resolver o problema, Millikan ou Fletcher [3] desenvolveram um método de borrifação de gotículas de óleo. Gotas minúsculas observadas desciam sob o efeito da gravidade e com aplicação do campo elétrico elas invertiam o sentido. Ligando e desligando o campo, a mesma gota pode ser observada varias vezes "passeando" pela visada do microscópio. Esta gota equilibrada [2] permite muitas medidas e, consequentemente, a possibilidade de realizar um estudo estatístico nos dados do experimento.

Em laboratórios de física, a execução do experimento Millikan tem como motivação mostrar que as gotículas de óleo possuem números inteiros de elétrons. Os valores obtidos experimentalmente das cargas elétricas são divididos pelo valor da carga elementar do elétron, ou seja, 1,602 x 10^-19 C e, deste modo, prova-se que essa carga elétrica está quantizada. Ela é um produto entre um número natural n e uma quantidade bem definida de carga e (e, 2e, 3e,..,ne). A estatística se restringe em fazer histogramas e verificar que as colunas mais altas estão centradas em valores múltiplos de e.

Para todo conjunto de pontos experimentais é necessário uma análise, ou seja, fazer um ajuste da melhor função que define estes pontos. Para definir o melhor conjunto de valores dos parâmetros ajustáveis de uma função é necessário partir de um critério. A teoria estatística estabelece a máxima probabilidade, que é entendida como o valor mínimo de χ² (chi-quadrado) [4]. Este critério é usado no trabalho e, também, é conhecido como método dos mínimos quadrados. Para o ajuste de um histograma o valor do chi-quadrado é dado pela Eq. (1) [5].

onde N_i é a altura da coluna do histograma e n_i é o valor obtido no ajuste de uma função distribuição tipo gaussiana.

2 Fundamentação estatística do experimento Millikan

Os valores da Tabela 1 são de 100 medidas tomadas no dia primeiro de setembro de 2011 no laboratório de Estrutura da Matéria da Faculdade de Ciências e Tecnologia da Unesp de Presidente Prudente por um dos autores deste trabalho (L.A.P.) que na época era estudante universitário.

A Tabela 1 mostra o número da medida e o valor da respectiva medida. A coluna "num" é o número de cargas inteiras da gotícula e foi calculado pela razão arredondada para inteiro entre a medida e o valor de 1,602 x 10^-19. A coluna "carga e" é o valor de cada carga obtido pela divisão do "valor" pelo "num". Por exemplo, a medida 6 provavelmente corresponde a cinco elétrons e o valor da carga e é próximo do esperado 1,602 x 10^-19 C. As medidas 49 (2,603 x 10^-19 C) e 50 (2,445 x 10^-19 C), foram removidas porque a análise dos dados revela que provavelmente estas medidas estão acompanhadas de alguma incoerência, talvez algum tipo de erro, e atrapalha o método usado no trabalho. A seção "Critério para Rejeição de Medidas Experimentais" explica o motivo da eliminação destes dois pontos e, desta forma, a quantidade de medidas agora fica reduzida para 98 medidas.

A Fig. 1 mostra o histograma feito com as medidas da Tabela 1. O histograma mostra que a variação dos valores das cargas medidas foi entre 1 x 10^-19 C e 16 x 10^-19 C. Teoricamente, os picos deste histograma estariam centrados nos múltiplos da carga do elétron, 1,602 x 10^-19 C, mas devido a existência de erros sistemáticos e aleatórios estes picos estão mascarados e, provavelmente, centrados em múltiplos de algum outro valor.

Sabendo-se que a carga do elétron é de 1,602 x 10^-19 C e dividindo-se com arredondamento para um número inteiro os valores da Tabela 1, obtém-se o histograma dado pela Fig. 2. No entanto, este artifício não é didático no sentido do cálculo da carga de um elétron porque este valor foi usado a priori. Os valores medidos estão acompanhados de erros sistemáticos e aleatórios e, portanto, o melhor valor da carga do elétron, do ponto de vista da estatística da medida, deve ser calculado. A equação da gaussiana da Fig. 2 é definida como y em função de x, onde y é a frequência e x é a carga elétrica. Nesta equação xc é o ponto central da gaussiana, w é a largura a meia altura e A é a área sob a gaussiana. Os valores são mostrados no quadro da Fig. 2 e todos estes parâmetros foram obtidos pelo método dos mínimos quadrados no ajuste das colunas do histograma. Observa-se na Fig. 2 que apesar de termos usado o valor de 1,60 x 10 C, o ajuste calculou um melhor valor de q para 1,5837 x 10 C. Neste ponto do trabalho vale ressaltar que o símbolo q está sendo usado para a variável independente da carga do elétron para se calcular o valor de chi-quadrado e o símbolo e é usado para o resultado da carga do elétron.

O algoritmo usado neste trabalho é o cálculo do chi-quadrado da Eq. (1) em função da variação do valor da carga q obtido no ajuste. Onde o valor da carga q é a variável independente e foi escolhido variar entre 1,60 x 10^-19 C e 1,75 x 10^-19 C, com passo de 0,01 x 10^-19 C. Por exemplo, o histograma da Fig. 2 apresenta um chi-quadrado com valor de 57,6. Segundo a teoria de erros [4] o melhor valor para a carga do elétron é aquele que apresenta a maior probabilidade de acontecer, que corresponde ao menor valor de chi-quadrado.

Variando o valor da carga q temos que os valores das colunas "num" e "carga e" da Tabela 1 sofrem alterações. Por exemplo, a medida 89 com carga de 10,572 x 10^-19 C pode ter as seguintes cargas e:

• para 1,60 x 10^-19 C, temos "num"= 7 elétrons e, portanto, e=1,510 x 10 C;

• para 1,65 x 10^-19 C, temos "num"= 6 elétrons e, portanto, e=1,762 x 10 C;

• para 1,70 x 10^-19 C, temos "num"= 6 elétrons e, portanto, e=1,762 x 10 C e

• para 1,75 x 10^-19 C, temos "num"= 6 elétrons e, portanto, e=1,762 x 10 C.

Portanto, o histograma mostrado na Fig. 2 se modifica gradualmente e de modo sutil com a variação da carga q e alcança uma configuração gaussiana boa para um determinado valor de chi-quadrado, ou melhor, para o menor valor do chi-quadrado que é dado pela Eq. (1).

3 Critério para rejeição de pontos experimentais

Segundo o algoritmo deste trabalho espera-se uma curva contínua ou bem comportada do valor do chi-quadrado em função do valor da carga do elétron. Desta forma, o valor ótimo da carga do elétron é calculado no ponto mínimo desta curva. No entanto, devido ao algoritmo proposto existem valores de cargas que estão no limiar n elétrons ou n - 1 elétrons. Por exemplo, a medida 50 que é 2,445 x 10^-19 C quando dividido pela carga q de 1,63 x 10^-19 C fornece de maneira arredondado o número inteiro de 2 e, em seguida, para q de 1,64 x 10^-19 C fornece um número arredondado de uma unidade de elétron. Os valores de chi-quadrado para os dois casos são de 11,7 e 3273, respectivamente. A Fig. 3(B) mostra um histograma com uma coluna afastada no lado direito e a Fig. 3(A) não tem esta coluna. Este fato é causado pelo ponto 50 que sai da parte esquerda do histograma (A) e vai para a parte direita de (B), quando o valor de q muda de 1,63 x 10^-19 C para 1,64 x 10^-19 C. Consequentemente, em função deste fato, o valor do chi-quadrado sofre uma descontinuidade de 11,7 para 3273.

De maneira análoga, a Fig. 4 mostra o caso idêntico que o ponto 49 sofre quando o valor da carga q muda de 1,73 x 10^-19 C para 1,74 x 10^-19 C. O gráfico da Fig. 4(B) apresenta mais uma coluna afastada do histograma no lado direito. Em função destas explicações que estão ilustradas nas Figs. 3 e 4, foi eliminado da estatística os pontos de números 49 e 50 da Tabela 1.

O gráfico da variação do chi-quadrado em função de q, sem a remoção dos pontos 49 e 50, é mostrado na Fig. 5. Os picos da curva existentes nos valores de q em 1,65 x 10^-19 C e 1,78 x 10^-19 C são explicados por estes pontos, conforme mostrados nas Fig. 3 e 4.

4 Resultado

O critério numérico a ser adotado para se obter o melhor valor de q ou a melhor gaussiana é o método da máxima probabilidade [4], que corresponde o menor valor de chi-quadrado dado pela Eq. (1). Quanto menor for o valor de chi-quadrado significa que mais perfeita é a distribuição gaussiana. Vários valores para q foram experimentados e a Fig. 6 mostra o gráfico da variação do chi-quadrado em função do valor de q. Neste ponto, convêm lembrar que os valores de q escolhidos foi entre 1,60 x 10^-19 C e 1,86 x 10^-19 C, com passos regulares de 0,01 x 10^-19 C, mas o valor de q obtido no ajuste de cada gaussiana é próximo do valor escolhido, por isso os pontos da Fig. 6 não estão igualmente espaçados. Por exemplo, o histograma da Fig. 2 foi feito para q de 1,60 x 10^-19 C e o ajuste revelou um valor de q de 1,587 x 10^-19 C.

A carga do elétron é determinada no mínimo da parábola da Fig. 6, que é dado pela Eq. (2). O ajuste ou "fitting" da parábola da Fig. 6 foi feito tomando os pontos de q entre 1,68 x 10^-19 C e 1,86 x 10^-19 C e fornecem os valores dos parâmetros (A, B e C) juntamente com seus respectivos desvios padrões, conforme mostrado no quadro da Fig. 6. Estes desvios padrões foram obtidos pelo método dos mínimos quadrados [4].

A incerteza da carga do elétron é obtida por propagação de erros usando os valores dados no quadro da Fig. 6. Esta incerteza é de 0,073 x 10^-19 C. Portanto, o valor da carga do elétron, com seu respectivo desvio padrão, obtido neste trabalho através do ajuste dos valores da Tabela 1 é dado pela Eq. (3).

A Fig. 7 ilustra de maneira bem didática que o histograma para q próximo de 1,81 x 10^-19 C é mais parecido com gaussiana do que o equivalente para q próximo de 1,60 x 10^-19 C.

5 Conclusão

A Fig. 6 mostra que o melhor valor para a carga do elétron no experimento Millikan realizado no laboratório da FCT - UNESP é 1,808 x 10^-19 C, com desvio padrão de 0,073 x 10^-19 C. Evidentemente que ele não é o valor real da carga do elétron, este valor é o valor obtido neste trabalho, usando os dados da Tabela 1, e o algoritmo mostrado. Este resultado para a carga q revela que o equipamento talvez tenha um erro sistemático. A diferença maior do que 0,2 x 10^-19 C entre o valor esperado e o valor calculado pelas medidas é superior do que dois desvios padrão e, desta forma, é impossível afirmar com certeza que esta diferença é um erro sistemático. Somente novas medidas podem definir melhor se o equipamento e o método estão produzindo erro sistemático ou é uma oscilação aleatória.

A Fig. 7 mostra a essência do trabalho que é procurar um valor para q que forneça um histograma parecido com uma curva gaussiana. Desta forma, pode-se observar de maneira clara que o valor próximo de 1,81 x 10^-19 C é melhor do que o valor esperado de 1,60 x 10^-19 C. Esta nova proposta de cálculo, para o valor da carga do elétron em um experimento Millikan, adiciona informação estatística ao nível universitário sem a necessidade de acrescentar modificações em aparatos experimentais já existentes em laboratórios de física.

Recebido em 1/3/2013

Aceito em 21/9/2013

Publicado em 6/2/2014

• [1] P.A. Tipler e R.A. Llewellyn, Física Moderna (LTC, Rio de Janeiro, 2001).
• [2] R.P. Crease, Os Dez Mais Belos Experimentos Científicos (Ed. Jorge Zahar, Rio de Janeiro, 2006).
• [3] C.A. dos Santos, Revista Brasileira de Ensino de Física 17, 107 (1995).
• [4] J.H. Vuolo, Fundamentos da Teoria de Erros (Ed. Edgard Blücher Ltda, São Paulo, 1992).
• [5] W.H. Press, B.P. Flannery, S.A. Teukolsky and W.T. Vetterling, Numerical Recipes in Pascal (Cambridge University Press, Cambridge, 1992).

Datas de Publicação

• Publicação nesta coleção
  29 Abr 2014
• Data do Fascículo
  Mar 2014

Histórico

• Aceito
  21 Set 2013
• Recebido
  01 Mar 2013