E nossos filhos nos chamarão de bárbaros

Rookie

election

É difícil explicar a cadeia de sentimentos que me atingiu com a eleição de Donald Trump para a presidência dos Estados Unidos. Em primeiro lugar, feriu meu brio estatístico, eu estava disposto a apostar dinheiro alto em Hillary pelas projeções e pesquisas feitas exaustivamente por matemáticos e cientistas americanos, eles adoram essas coisas; e foi outro Brexit. Em segundo lugar, pelo desespero trazido pela realização que, Europa ou EUA, a esperança não é mais capaz de vencer o medo.

Sou, e sempre fui, um otimista. Sempre acho que meu ano presente é o melhor de minha vida, que a sociedade caminha em uma direção mais aberta, pluralista, receptiva e tolerante, que diversos direitos conquistados por minorias são agora inalienáveis e que a questão do racismo, da misoginia e da intolerância será naturalmente resolvida pela passagem inexorável das horas, pela morte, de velhice e pacífica, dessa geração criada com outros modos e outros costumes. Depois de hoje, não sei se sou capaz de levantar essa bandeira com tanta convicção.

É difícil generalizar esses eleitores, mas vou tentar assim mesmo. Sinto como se eles quisessem surrar com violência o establishment, Washington no caso dos americanos, e Trump era o punho. Sinto também, na Inglaterra e nos EUA, como se esse voto fosse o último grito de uma maioria branca que percebe no horizonte o fim de sua hegemonia e grita, esperneia, porque “cultura do país” e “sua cultura” não serão sinônimos óbvios. “Cultura branca” não é algo que existe, chamamos isso nesses países, infelizmente, de apenas cultura.

Porque ainda que tanto o Brexit quanto a eleição de Trump estejam fortemente correlacionados a raça, escolaridade e idade dos eleitores (brancos, sem ensino superior e acima de 50 anos votaram em bloco para o sucesso de ambos os cataclismos),  sei que o silêncio de um dos lados é tão culpado quanto o barulho do outro.

Eu me considero um liberal, em um sentido bem vago da palavra. Posto bastante sobre política, então evito falar de minhas opiniões para não perder credibilidade, mas Donald Trump não é exatamente uma questão de opinião. Essa vitória chocante é um evento nauseante na história da democracia americana, um homem carregado ao poder por forças de nativismo, nacionalismo, misoginia e racismo. Trocam o primeiro presidente negro por um charlatão que sugou apoio de áreas xenofóbicas e de supremacia racial branca. É impossível reagir a esse evento com algo menor que repulsa. E ainda, apesar de todos os escândalos, apesar de todos os absurdos, os eleitores odiaram mais o que ele dizia combater do que seus hábitos horripilantes, e pagaram para ver o que estava “atrás da porta número 3”, como se em um programa de auditório, ao invés de uma mulher competente com história de vida pública, talvez história até demais.

Li uma vez um artigo no Le Monde cujo título me marcou bastante, ele dizia: “Et nos enfants nous appelleront: ‘barbares'”, ou, na minha tradução fraca, o título desse post. Ele me faz pensar em qual é o papel de alguém que se diz progressista, alguém que deseja estar do lado certo da história. Não é apenas, como alguns pensam, defendem direitos estabelecidos e claros, mas buscar novos. Nossa tarefa, nessas derrotas tenebrosas, nesses tempos sombrios, é erguer a cabeça e continuar lutando contra os males de hoje e de amanhã. Meu objetivo é não me tornar como minha avó, uma pessoa maravilhosa e mulher progressista em seu tempo, mas que não era capaz de admitir que Rui Barbosa era pardo, tampouco Machado de Assis.

Pensamos na escravidão, na ditadura militar, na repressão da era Vargas, na segregação americana e no Apartheid africano, no papel das mulheres durante toda a antiguidade até mais ou menos 1970, em homossexuais sendo castrados quimicamente, condenados à morte, em transexuais atravessando uma vida de repressão e culpa, e pensamos: “Bárbaros!”. Nossa tarefa não é apenas combater os males de hoje, não é apenas ter um amigo negro, ou um chefe gay, e achar que racismo e homofobia não se aplicam a você. Nossa tarefa é pensar: “por que meus filhos me chamarão de bárbaro?”, como eu chamo minha avó, e ela chamava, provavelmente, sua avó. Respondendo essa pergunta, você descobre o mal a combater.

Eu queria poder postar algo científico, estatístico, sobre a vitória de Trump, analisando detalhes da sondagem e apontando os pontos fracos, mas o fato é: ninguém sabe o que aconteceu. Explicações sobram, mas quando já se sabe a resposta fica bem mais fácil. Talvez as pessoas que mais votariam no Trump também seriam as que mais desconfiariam da mídia golpista que ligaria para fazer uma sondagem, o que criou uma sub-representação de trumpistas em três ou quatro estados chave, o suficiente para o sistema americano colapsar. Prometo um post sobre a forma de governo americana, e sobre como ela é particularmente ruim, para aparecer semana que vem.

É bom estar de volta, queria que as circunstâncias fossem outras.

Ondas gravitacionais

Rookie

Estou atrasado para comentar algo sobre ondas gravitacionais, mas esse atraso tem dois motivos. O primeiro é a saturação da mídia, é difícil achar algo para escrever que já não tenha sido escrito, e melhor, sobre o tema; e o segundo é eu já ter comentado isso no alarme falso de ondas gravitacionais que tivemos em 2014.

Dessa outra vez, como todo o resto da mídia, empolguei-me antecipadamente, caindo no grande erro de cantar vitória antes do peer-review. Esse episódio serve de lição no que é o centro da ciência: a ideia da autocorreção, e de não ter medo de dizer que está errado. Comento isso em outro post, nesse vamos ver o experimento LIGO, o que fez e o que mediu.

O artigo você encontra aqui. Vou tentar contar a história que o artigo apresenta sem a linguagem técnica, que nem eu entendo tão bem, não sou astrofísico, mas meu curso de relatividade geral no mestrado vai me servir de algo. Em setembro de 2015, o observatório LIGO (Laser Interferometer Gravitational-Wave Observatory), encarregado de medir alterações no campo gravitacional em sua volta, detectou um forte sinal. Ele possui dois detectores, um em Washington e um em Louisiana, e é importante ter dois detectores para diferenciar o que é ruído e o que é sinal. Todo detector, como antena de televisão antiga, possui ruído, e ondas gravitacionais são tão difíceis de medir que apenas eventos que produzam alguma alteração considerável na gravidade podem vencer o ruído e ser medidos. O sinal medido foi esse:

LIGO_data_1O eixo x, cortado do gráfico, representa intervalos de 0.05s nos ticks maiores. Na esquerda você vê o sinal medido em Washington e na direita, Washington e Louisiana sobrepostos. Pelo gráfico da direita você consegue ver que o começo do sinal ainda está cheio de ruído, há áreas em que eles não concordam nem um pouco, mas a partir do meio do gráfico há forte coincidência nos sinais, uma oscilação acelerando, crescente, e uma queda brusca na intensidade. Esses 0.1s de sinal são o que físicos esperavam faz 100 anos, desde que Einstein formulou a relatividade geral e propôs que a gravidade, como o eletromagnetismo, deve produzir ondas.

Essas ondas não são algo fácil de explicar ou entender sem alguma matemática, mas vou tentar. Gravidade, como descrita pela relatividade geral, é uma deformação do espaço-tempo. Eu sei que sem nenhuma equação isso parece mais uma frase de terapia de autoajuda quântica que ciência séria, então tente imaginar que gravidade acontece porque a massa dos corpos afeta o espaço em torno do corpo e atrai outros corpos para si. Quanto maior a massa, maior é esse efeito. Se um corpo massivo acelera, ele afeta o espaço em torno de si de um jeito um pouco diferente, ele perturba esse espaço com sua gravidade e essa perturbação pode se propagar.

Em uma analogia, uma grande explosão de dinamite libera energia que atinge as moléculas no ar, elas se propagam atingindo outras moléculas no ar, que atingem outras, e outras, até atingirem o frágil osso em nosso tímpano e registramos isso como um barulho alto. A gravidade se comporta de um jeito parecido, mas, diferentemente do som, ela é extremamente fraca, o que pode ser entendido na analogia como se fôssemos praticamente surdos. Apesar de nossa audição comprometida, Einstein previu faz cem anos que o som existe, e estamos tentando medi-lo ou detectá-lo de vários jeitos, tentando criar um aparelho auditivo cada vez melhor. Um dia, uma explosão acontece, e temos um excelente aparelho auditivo (o LIGO); finalmente, e pela primeira vez, conseguimos ouvir algum som do universo. Esse som você vê na figura acima.

E isso é diferente de simplesmente “detectar gravidade”, porque isso detectamos o tempo todo. Se você pular do quinto andar, certamente detectará mais gravidade do que gostaria, essa medida é de uma natureza diferente. Ainda na nossa analogia com som, é como se soprassem forte em seu rosto. Você consegue sentir vento e sopro sem problemas, detectar ar e movimento de ar não é a questão. Mas som é outra coisa. Som é a propagação de vibrações no ar causadas por objetos à distância, e ele carrega informação sobre esses objetos, é bem diferente de vento e brisa. Eu sinto por qualquer astrofísico que esteja lendo essa analogia entre gravidade e vento, mas foi o melhor que pude fazer.

Vale a pena discutir um pouco como esse tal de LIGO funciona. Ele tem dois braços (cada um de 4 km cada), um perpendicular ao outro, e cada braço possui um laser que vai e volta esses 4 km. Quando uma onda gravitacional passa, ela deforma o braço, porque altera a gravidade sutilmente e essa alteração faz a interferência entre os lasers mudar ligeiramente, mas o suficiente para medir uma alteração muito pequena. A sensibilidade de LIGO é excepcional, ele é capaz de detectar uma diferença no tamanho dos braços de uma distância menor que um núcleo atômico.

Depois de medir esse sinal, os cientistas usaram uma boa análise comparativa para diminuir o ruído e dar aquela alisada nos dados, todos nós físicos somos versados na arte de alisamento, chapinha e escova de dados experimentais. O resultado é:

LIGO_2Agora sim isso está bonito. Tendo o sinal bonito e alisado, que são as linhas cinza desses gráficos, resta fazer o trabalho de detetive. Pela primeira vez escutamos um som com nossos ouvidos, mas que barulho é esse? O que produziu esse som? Parece ambicioso querer saber a natureza do barulho logo na primeira vez que escutamos qualquer coisa, mas poucas coisas no universo são capazes de produzir um ruído desse calibre, então podemos listar as mais prováveis e ver qual melhor se encaixa no que escutamos.

Ondas gravitacionais são produzidas por objetos massivos acelerados, isso sabemos. No espaço, a única fonte de aceleração possível para um corpo massivo desses (tirando a expansão do universo e tal) é a presença de outro corpo massivo que, pela sua gravidade, puxa aquele primeiro. E o primeiro puxa também o segundo. Assim como o Sol puxa a Terra, a Terra também puxa o Sol, mas é mais ou menos um cabo de guerra entre uma criança de cinco anos e um levantador olímpico de peso, o Sol nem se mexe. Quando dois corpos bem massivos se encontram, de massas comparáveis, nenhum ganha no cabo de guerra e, ao invés de um orbitar o outro, ambos orbitam um ponto comum em uma bela dança conhecida como sistema binário, algo assim:

E isso causa alterações na gravidade em forma de ondas periódicas, exatamente o que o LIGO mediu, então é uma boa sugestão que esse sinal tenha vindo de um sistema binário. Mas a frequência e a intensidade do sinal aumentam e colapsam, terminando em uma linha praticamente reta como o sinal cardíaco de alguém em seu último suspiro. O que poderia ter causado esse fim ao sistema binário? A solução mais simples, os corpos do sistema binário se encontraram e viraram um só.

Um sistema binário costuma ser bem estável, assim como a Terra girando em torno do Sol é um sistema, felizmente, bem estável. Mas ao longo das centenas de milhões de anos esse sistema vai perdendo energia, muito em parte por emitir as tais ondas gravitacionais, e essa perda de energia faz os corpos irem se aproximando nessa dança. Chega uma hora que os parceiros desse baile se encontram, e a explosão monumental que ocorre foi o que o LIGO mediu. Pelos dados que temos, e pela teoria que sabemos, podemos medir a massa e o tamanho desses corpos massivos que, no final da dança, se fizeram ouvir pelos quatro cantos do universo: são, sem dúvida, dois buracos negros em um sistema binário que colidiram e se tornaram um.

Não há analogia possível para esse evento. Posso apresentar números, mas nada que qualquer cérebro humano seja capaz de conceber, de tentar entender. Um evento nessa proporção foge do que nosso cérebro foi condicionado para imaginar. Posso dizer quantas centenas de bilhões de bombas nucleares isso representa, mas essa analogia tanto é sem sentido quanto ainda está longe de representar a grandeza em termos de energia dessa colisão. Temos sorte de não estar nas redondezas, de não ver de mais perto o evento, essa explosão monumental, essa mostra de poder do universo, esse Krakatoa do cosmos.

Com os grandes progressos da análise computacional e de modelos calculáveis de relatividade geral, podemos simular em nossos computadores os valores desse explosão, programar em C++ o nosso próprio Krakatoa, e ver que tipo de onda gravitacional sairia disso. É com bastante alegria que digo que isto são as linhas coloridas do gráfico anterior, alinhando perfeitamente com o que foi medido. A história desse evento, que aconteceu todo em menos de um segundo, com buracos negros de massa dezenas de vezes maiores que o Sol se movendo a metade da velocidade da luz, é dada por esse gráfico:

LIGO_3E foi isso que foi medido entre Washington e Louisiana em setembro de 2015.

Eu gosto muito de física, e é difícil segurar o fascínio, e um pouco da emoção (sou clichê assim mesmo), quando contemplo algo parecido. A medição traz algo de pessoal a esse fato. Eu certamente conseguiria conceber uma colisão de estrelas, de buracos negros ou de galáxias teoricamente, eu já fiz exercícios na faculdade de coisas parecidas, mas isso tudo é impessoal, teórico demais; é muito diferente de de fato medir algo assim, saber que aconteceu, onde e quando.

Poderia parar o post por aqui, mas tenho recebido várias perguntas de “para que isso serve?”. Tenho minha resposta favorita a essa pergunta, mas resisto porque sei que nem todo mundo tem esse gosto pela coisa. Mantenho a analogia com a audição: ganhamos um sentido novo. Podemos medir coisas que antes não eram possíveis. Há diversos objetos que não interagem de forma eletromagnética, em especial matéria escura e energia escura (são “escuras” com um motivo), mas que interagem de forma gravitacional. Podemos agora ter esperança de medir eventos e objetos que apenas quem tem ouvidos pode escutar. Claro, ainda estamos na infância disso, nosso aparelho de audição pode melhorar bastante, ele ouviu algo porque a explosão foi bem grande, mas é um começo. O universo está rico em sons gravitacionais de várias fontes, e de várias formas, e, pela primeira vez, estamos escutando.

5 sextas, 5 sábados, 5 domingos

Rookie

Chega.

Já chega.

Por favor, chega de enviar correntes de e-mails me dizendo que Janeiro terá cinco sextas, cinco sábados e cinco domingos. Já não consigo mais receber essas mensagens sem espumar de raiva, por diversos motivos, e tomo meu tempo para escrever esse post, um pouco fora do tema desse blog, para dissertar sobre um de meus combates épicos: o spam mal feito.

Tenho vários combates, todos perfeitamente justificados. Coxinha comida pelo lado errado, e-mail escrito em azul (ou qualquer outra cor além de preto), e-mails @hotmail, gente que vai “guardar mesa” enquanto outra pessoa está na fila, spoilers na internet, gente que conta o final da piada quando eu estou perto de terminar de contar, mensagens-imagens de “boa quinta-feira” no grupo de WhatsApp da família, astrologia, Windows Vista, café solúvel, francês que deseja “bonne année!” em fevereiro, teclados com switch de borracha, notebooks da HP, tomada de três pinos (sou clichê assim), texto “científico” escrito em Word, gente que critica arte moderna dizendo que “fazia pintura mais bonita quando tinha cinco anos”, emissora de televisão que acha que perguntar para gente aleatória na rua o que eles acharam daquilo configura “jornalismo”, qualquer uso de caps lock que não seja para siglas; todos esses eventos são bandeiras que carrego com mais ou menos intensidade. Poucos deles, no entanto, provocam em mim meus instintos mais primitivos quanto spam mal feito, mal pesquisado, mal compartilhado.

E o e-mail das cinco sextas, sábados e domingos é um exemplo primoroso dessa ojeriza que me é evocada. Sim, janeiro 2016 é um mês com essas características, mas isso não é surpreendente ou impressionante porque não só é um evento comum, mas acontece pela simples passagem do tempo que, por definição, é a coisa mais previsível que existe. Podemos fazer uma conta rápida para estimar a frequência com que isso acontece. Para que um mês possua cinco desses dias, basta que comece em uma sexta-feira e tenha 31 dias. As chances de um mês aleatório começar em uma sexta, pela isotropia da semana, são \frac{1}{7}. São sete meses no ano com 31 dias (usando os ossos da mão para contar fica fácil), então a chance de um mês ter essas duas características é \frac{1}{7}\times\frac{7}{12}=\frac{1}{12}.

Esse evento acontece, em média, uma vez por ano, não uma vez a cada 823 anos, de onde tiraram esse número?! Ele nem funciona! Janeiro de 2839 não terá cinco sexta, cinco sábados e cinco domingos! O primeiro dia de janeiro de 2839 é um sábado, não uma sexta. Se você quer repetição, um número desses seria formado com o produto da duração de todos os ciclos envolvidos nessa repetição. Por exemplo, a semana tem sete dias e ano bissexto acontece uma vez a cada quatro anos, então tenho certeza de que o mês de janeiro de daqui a 4\times 7=28 anos terá também cinco sextas, cinco sábados e cinco domingos. Isso não funciona sempre porque ano bissexto tem problemas quando o ano é múltiplo de 100 e 400, o que me forçaria a multiplicar por coisa demais, mas isso é um jeito de garantir a periodicidade e é conhecido como teorema de Lagrange. Se você começar uma atividade no domingo e repeti-la de dois em dois dias, depois de sete repetições você voltará a fazê-la no domingo. Se repetir de três em três dias, de quatro em quatro dias, de cinco em cinco dias, não importa; depois de sete repetições você estará sempre no domingo, com certeza absoluta.

O número correto de anos seria um número formado do produto de diversos outros números, cada um correspondente à duração dos ciclos envolvidos nessa periodicidade (semana, ano bissexto). Mas 823 não é apenas o número errado, ele é primo! E não é tão fácil chutar um primo tão alto assim, então parabéns aos criadores do spam, eles conseguiram, por azar, escolher um dos números mais errados possíveis, se existe uma gradação de errado para esse problema.

E o compartilhamento desse spam me causa tanto desgosto porque uma simples volta ao calendário poderia resolver o problema. Esse mesmo evento aconteceu em maio do ano passado, meu aniversário foi em um desses “sábados sagrados dos chineses”. E o que chineses têm a ver com isso? Eles nem usavam o calendário gregoriano!

Por fim, me acalmo. Sei que o combate é vão, que spam são meus gigantes em moinhos de vento, bem como e-mail em azul e a tomada de três pinos, não vencerei. Fica esse desabafo, e esse comentário sobre essa interessante propriedade de repetição de teoria dos grupos. Sei que esse post é inútil, a intersecção de leitores desse blog com quem compartilha esse tipo de abominação é perto de zero, meu alcance é baixo, meu poder evangelístico de apresentar as boas novas de pegar o calendário deste ano e verificar que a mesma coisa acontecerá em julho é pequeno.

A biblioteca de Babel

Rookie

Hoje não tem modelo ou coisa muito complicada, vamos falar de conhecimento, e de um de meus contos favoritos: a biblioteca de Babel. Publicado por Jorge Luis Borges em 1941, esse conto apresente um universo distópico de excesso de informação fascinante. Passando por ele, comento um pouco do propósito desse blog, e do valor da educação convencional e análise crítica e científica em uma era dominada pela consulta imediata a todo conhecimento humano na ponta dos dedos em um celular.

Imagine uma grande biblioteca, grande mesmo. A descrição de Borges é intrincada e interessante, a estrutura da biblioteca é feita de hexágonos intrincados com prateleiras repletas de livros em ordem aparentemente aleatória. A biblioteca contém apenas livros de exatas 410 páginas, com cada página tendo 40 linhas e cada linha contendo 80 símbolos. Para simplificar, esses símbolos são apenas 22 letras, espaço, vírgula e ponto; estamos usando um alfabeto simplificado. Essa biblioteca é única por uma característica curiosa, ela contém todos os livros possíveis de 410 páginas.

Uma representação da biblioteca.

É um projeto ambicioso, mas seguramente podemos dizer que a biblioteca de Babel é finita. Uma conta rápida nos mostra que ela contém 25^{410\times 40\times 80}\approx 10^{1.834.097} livros. Em comparação, a maior biblioteca do mundo, a British Library, tem 1.5.10^8 livros. Em termos de bits, a biblioteca de babel contém aproximadamente 2^{6.092.738} bits, ou bytes, não lembro a diferença, mas chega nessa escala de números essas definições são pálidas para descrever a grandeza da biblioteca de Babel.

O conto é muito interessante e descreve o universo dos cidadãos atormentados por esse sonho sadista de um deus cruel criador do maior pesadelo de um bibliotecário. Existem puristas que desejam queimar todos os livros que não fazem sentido para preservar apenas os puros, os sagrados, que contam algo com nexo; ainda que haja intensos argumentos e debates, pois um livro aparentemente sem sentido pode ter nexo, coerência e coesão em alguma língua desconhecida. Talvez seja o trabalho deles descobrir essa língua, mas até que ponto isso não seria, ao mesmo tempo, definir o conteúdo do livro?

Mais do que isso. Nessa biblioteca certamente há um livro relatando exatamente o seu futuro, e o meu, e todos os nossos futuros possíveis que podem ser descritos em 410 páginas. Nessa biblioteca há uma lista de todos os seus romances, suas aventuras, seus desejos, e esse post do blog. Ainda, é uma coleção finita de livros, bastaria você procurar com carinho e você acharia, após algum tempo, todas as respostas que busca, e todas as perguntas que teve.

Não quero lançar spoilers, paro por aqui. A ideia da biblioteca me fascina por motivos diferentes, eu vejo partes de Babel em meu dia a dia como divulgador e como quem se interessa por muitas coisas que não conhece. O maior problema da biblioteca, e talvez sua principal característica, é sua inutilidade completa. Se você quer alugar um livro, digamos “História do Brasil“, de Boris Fausto, dizer título e autor não basta. Não apenas esse livro está na biblioteca, mas todas as 410\times 40\times 80 variantes com um erro de digitação, ou as variantes com dois erros de digitação, três, quatro, cinco e assim por diante. Para obter o livro “História o Brasil“, você precisaria especificar letra por letra o livro que deseja, ferindo o propósito do aluguel. Em outras palavras, para especificar um livro, você precisa fornecer a mesma informação contida no livro, tornando a obtenção do livro inútil.

E observo Babel na vida real diversas vezes na Internet. Como disse um amigo, quando eu cheguei na Internet era tudo mato. Quinze anos depois, essa soma do conhecimento humano está longe de atingir os padrões de Babel em termos de bits ou bytes, mas está crescendo, tornando-se gigantesca e por vezes lembra os hexágonos da biblioteca do pesadelo. Imagine que você busca algo sobre evolução, vacinas, dietas, curas, física quântica ou sobre a importância das estrelas em nosso dia a dia. Esses temas possuem tanto ruído, tantos “livros com erros de digitação”, que você precisa fornecer mais informação para obter o conhecimento que deseja. Você precisa já saber e conhecer alguns desses temas para navegar nas sugestões de sua consulta, correndo o risco de cair em um dos livros da categoria que os puristas adorariam queimar.

Muita coisa na área científica dessa grande Internet está mal explicada, distorcida, deturpada, por uma mistura de charlatães e leigos bem-intencionados com pouca paciência para aprender ciência como se deve. Meu trabalho, e o seu, é iluminar os bons livros, permitindo que sejam encontrados com menos informação do que eles contêm. É uma cruzada contra o ruído, contra a pseudociência, abrindo caminho nessa internet que ainda tem muito mato para se cortar.

 

Os aluguéis de Berlim

Geek Hardcore

Minha namorada está buscando um apartamento para alugar em Paris começando no mês de novembro, e a busca tem sido uma tragédia em diversos atos. Cubículos pequenos, muita competição, preços altos e até imobiliárias que mais parecem esquema de pirâmide que um estabelecimento sério têm acompanhado a jornada dessa física teórica nas últimas duas semanas. Quanto aos preços, ao menos há um fator positivo, uma lei recente baixada inspirada na solução que Berlim tenta implementar para resolver seu problema de aluguéis: um teto de preço baseado no valor médio das redondezas.

A lei diz que o valor de um aluguel não pode ser superior à média do bairro de uma diferença maior que 10%. Ou seja, se a média do bairro é 22 euros o metro quadrado, um aluguel não pode ser mais alto que 24,2 euros o metro quadrado. Simples assim. Essa é uma ideia bem divertida, que me fez considerar a seguinte questão: o problema da Berlim gananciosa.

Nessa nova Berlim, a lei é um pouco mais rígida. Uma pessoa não pode impor um aluguel maior que a média de seus vizinhos cuja diferença entre a média e o aluguel seja maior que 10%. A pergunta é: é possível que exista uma cidade com essa lei em que todas as pessoas cobram o maior aluguel possível pela lei? Ou seja, obrigando as pessoas a cobrarem exatamente 110% da média de seus vizinhos, é possível haver uma cidade em que todas as pessoas conseguem cobrar essa quantia?

Para bem definir o problema, vamos estabelecer hipóteses fundamentadas na realidade. A primeira: minha cidade é unidimensional, parecida com Cuernavaca, porque quero uma noção de vizinhança bem definida e porque D=1 é meu número favorito de dimensões. Segunda: casas na periferia não têm esperança de ganhar dinheiro no aluguel em relação aos vizinhos e possuem preço fixo. Isso vem do fato, inspirado em minha experiência parisiense, de que morar longe do centro só é bom se for na direção de Versailles, e olhe lá. Isso fixa as condições de contorno do problema. Dado um preço fixo nas casas da periferia, é possível que o resto da cidade cobre um aluguel superior em \alpha porcento à média de seus vizinhos?

A resposta, e já mando com spoilers, é: depende do \alpha. Há valores em que isso é possível, mas existe um limite para a extorsão de seus inquilinos se você precisa respeitar a regra dos \alpha porcento. Para entender o que acontece, vamos começar com o caso de uma cidade pequena, a cidade de Berlim de três habitantes, chamada B_3. Vejamos a cara dessa cidade.

Cidade de Berlim 3A regra geral das casas de Berlim gananciosa é:

 x_n = \frac{(1+\alpha)}{2}(x_{n+1}+x_{n-1}).

Seja (1+\alpha)/2=\lambda, para simplificar nossa vida. A única exigência da primeira cidade de Berlim que vamos estudar, B_3, é que x_2=\lambda (1+1)=2\lambda, então o valor de \lambda pode ser qualquer coisa, não há limites para a ganância de B_3

O caso de B_4, no entanto, é diferente. Vejamos como é essa cidade:

Cidade de Berlim 4Neste caso, os valores de x_2 e x_3 serão fixos pelo sistema

 \begin{cases} x_2=\lambda(1+x_3) \\ x_3 = \lambda(1+x_2) \end{cases} .

Vocês sabem que eu gosto de matrizes, então vamos colocar isso em notação de gente grande. Os aluguéis de B_4 serão definidos pelo sistema

 \left(\begin{matrix} 1 & -\lambda \\ -\lambda & 1 \end{matrix}\right)\left(\begin{matrix} x_2 \\ x_3 \end{matrix}\right) = \left(\begin{matrix} \lambda \\ \lambda \end{matrix}\right).

Esse é um sistema linear bem bonito, mas vocês devem se lembrar das aulas do ensino que diziam, com razão, que a vida fica mais cinza, o mundo perde a cor e você perde a vontade de cantar uma bela canção quando o determinante da matriz que define o sistema linear se anula. Quando isso acontece, não há mais solução possível para o sistema. E note que o determinante que nossa querida matriz é

 \det \left(\begin{matrix} 1 & -\lambda \\ -\lambda & 1 \end{matrix}\right) = 1-\lambda^2.

Ou seja, algo muito errado acontece quando \lambda=1, que acontece no mesmo momento em que \alpha=1. Não apenas o sistema não tem solução, não há mais valor de alguel possível para satisfazer tanta ganância, mas se você tentar impor um \lambda mais alto que um, não encontrará a solução que deseja. Se \lambda > 1, a solução do sistema apresentará valores negativos de x_2 e x_3, e acho que é seguro acrescentar como hipótese baseada na realidade que preços de aluguel devem ser positivos.

Assim, percebemos que com quatro casas o valor máximo que a lei pode impor para a extorsão dos inquilinos é de \lambda=1, o que equivale a \alpha=1, ou o máximo que podemos impor é que uma casa cobre o dobro da média dos vizinhos. Vejamos o que acontece quando a cidade é maior, visitemos B_5.

Cidade de Berlim 5Nesse caso o sistema linear que define os valores de x_2, x_3 e x_4 será

 \left(\begin{matrix} 1 & -\lambda & 0 \\ -\lambda & 1 & -\lambda \\ 0 & -\lambda & 1 \end{matrix}\right)\left(\begin{matrix} x_2 \\ x_3 \\ x_4 \end{matrix}\right) = \left(\begin{matrix} \lambda \\ 0 \\ \lambda \end{matrix}\right).

Vamos chamar a matriz que define esse sistema 3 x 3 de A_3. A partir de agora, a matriz que define o sistema n \times n será A_n. O determinante dela pode ser calculado com aquela técnica que você aprendeu no segundo colegial e errava toda vez que tentava fazer de cabeça sem usar aquela técnica das duas colunas fantasmas que você copia ao lado da última.

 \det A_3 =1-2\lambda^2

E o determinante se anula no valor \lambda = \sqrt{2}/2, o que equivale ao valor de \alpha=\sqrt{2}-1\approx 0.71, que é menor que um. É fácil perceber a tendência desse problema: quanto maior a cidade, menor é a taxa de ganho extra que podemos cobrar e ainda exigir que todos cobrem mais aluguel que seus vizinhos. O valor de \alpha diminui com o número de casas N da cidade! Mas… diminui como?

E aqui vou descrever como foi meu processo para resolver esse problema. Eu queria achar a relação entre o \alpha máximo e N, o tamanho da cidade. Eu sabia que \alpha seria decrescente em N, mas não sabia como. Seria \alpha proporcional a 1/N? Ou a 1/N^2? Ou talvez uma exponencial negativa de N? A exponencial eu podia ignorar, vi matrizes vezes demais no olhos e sei que esse tipo de problema terá uma relação polinomial, mas não vou excluir a possibilidade. Para resolver esse problema, a primeira coisa que fiz foi pensar como físico, ou seja, trigger warning aos que gostam de um raciocínio matemático rigoroso. Vamos pegar uma cidade muito grande, grande o suficiente para ignorar as bordas, e vamos voltar à equação principal do modelo:

 x_n = \frac{(1+\alpha)}{2}(x_{n+1}+x_{n-1}).

E podemos jogar um pedaço para lá e outro para cá para obter uma outra forma dessa equação de recorrência:

 \alpha (x_{n+1}+x_{n-1}) = x_{n+1}+x_{n-1}-2x_n.

O termo à direita é um velho amigo nosso, é uma segunda derivada discreta do aluguel. Nós brincamos com algo parecido em um dos primeiros posts desse blog, sobre a interpretação física do Laplaciano, e chegou a hora de usá-la. Eu tenho que dividir esse termo da direita por um diferencial ao quadrado para obter uma derivada, e esse diferencial é o tamanho de meu grid na derivada discreta. Em termos do modelo, eu defino que a posição de cada casa é dada pela variável l e a distância entre duas casas é dada por \Delta l^2. Dividindo ambos os lados da equação por \Delta l eu obtenho:

 \frac{\alpha}{\Delta l^2}(x_{n+1}+x_{n-1}) = \frac{x_{n+1}+x_{n-1}-2x_n}{\Delta l^2}.

O lado direito dessa equação, quando a cidade é bem grande e as casas são bem próximas, converge para \frac{d^2x}{dl^2}, enquanto o termo da esquerda pode ser simplificado lembrando que a soma dos vizinhos é \frac{2x_n}{(1+\alpha)}. Isso tudo resulta na seguinte equação diferencial:

 \sigma.x(l) = \frac{d^2x}{dl^2}.

Sendo

 \sigma = \frac{2\alpha}{(1+\alpha)dl^2}.

Eu sei que esse diferencial no denominador é revoltante a seu espírito matemático, mas não é mais minha responsabilidade zelar por seu coração sensível. Esse dl é a distância entre duas casas. Suponhamos que a cidade toda tenha tamanho 1, em unidades de tamanho de cidade, então a distância média entre duas casas é dl=1/N. Para que essa equação diferencial faça algum sentido, eu preciso que \sigma não seja nem zero, nem infinito, então, como exigência, eu preciso que \alpha diminua conforme dl diminui, ou seja, preciso que:

 \sigma=\frac{2\alpha}{(1+\alpha)dl^2}\sim 1\Rightarrow\alpha\sim\frac{1}{N^2}

Sendo esse \sim o símbolo oficial de “tem o mesmo comportamento assimptótico de”, mas costumamos ler como: “se comporta como”, “é mais ou menos parecido em algum limite que eu espero que você tenha entendido” ou “não exatamente, mas tipo quase lá”.

Então o valor máximo da extorsão deve decrescer com o quadrado do número de casas. O problema está razoavelmente resolvido, mas o algebrista dentro de mim grita por uma solução melhor que essa. Vamos ver como seria aquela matriz A_n para uma cidade bem grande. Teremos algo da forma:

 A_n = \left(\begin{matrix} 1 & -\lambda & 0 & 0 & \cdots & 0 \\ -\lambda & 1 & -\lambda & 0 & \cdots & 0 \\ 0 & -\lambda & 1 & -\lambda & \cdots & 0\\ \vdots & \ddots & \ddots & \ddots & \ddots & \vdots \\ 0 & 0 & 0 & \cdots & -\lambda & 1 \end{matrix}\right).

A_n é uma matriz tridiagonal com diagonal repleta de 1’s e co-diagonais preenchidas por -\lambda. Eu não apenas quero o determinante dessa criatura, eu quero encontrar as raízes do determinante e, mais especificamente, a primeira raiz que seja maior que \frac{1}{2}. Eu sei que \lambda=\frac{1}{2} equivale ao caso \alpha=0, onde todos aluguéis são iguais à média dos vizinhos e, pelas condições nas bordas, isso representaria o caso em que todos os aluguéis são iguais a 1. Valores de \lambda menores que meio representariam alguém cobrar um aluguel menor que a média dos vizinhos e, por mais que esse seja o sonho de minha namorada, não serão considerados. A primeira raiz do determinante que é maior que meio equivale ao primeiro valor de \lambda que torna o sistema impossível, valores acima disso tornam alguns aluguéis negativos e não me interessam. Eu quero saber como a primeira raiz do determinante acima de meio se aproxima de meio quando o tamanho da matriz aumenta, e isso não é um problema muito fácil.

O determinante de uma matriz tridiagonal tem fórmula, e nesse caso essa fórmula é:

 \det A_n=\det A_{n-1}-\lambda^2\det A_{n-2} .

Isso não me ajuda muito, aliás, quase me desespera. Isso significa que o determinante será um polinômio de grau cada vez mais alto em \lambda e eu conheço minhas limitações, resolver um polinômio de grau maior que quatro exatamente quase nunca é possível e eu começo a perder a esperança quando escrevo os primeiros cinco casos desse polinômio, eles não parecem ter nenhum padrão reconhecível e minha busca parece vã. Eu vou atrás das raízes desses polinômios, e lá as coisas começam a ficar interessantes.

No caso de A_2, encontramos \lambda=\pm 1 como raízes. Para A_3, temos as raízes \lambda=\pm\frac{\sqrt{2}}{2}. O determinante A_4 ainda tem raízes, é um polinômio de grau 4 que pode ser reduzido a um polinômio de grau 2 (a origem dessa redução é a simetria do problema, a segunda casa, por simetria, deve ter o mesmo valor da penúltima, por exemplo). As raízes são \lambda=\pm \frac{1\pm\sqrt{5}}{2}, que, aliás, são a razão áurea e sua inversa.

Espera… 1, \frac{\sqrt{2}}{2} e razão áurea… eu conheço esses valores. Eles são todos cossenos de alguém (ou quase)! Isso não pode ser coincidência. O determinante de A_n é um polinômio cujas raízes possuem uma ligação com cossenos. E quando eu percebi isso, fechei os punhos e soprei por entre os dentes cerrados a palavra Chebyshev

Eu e você nos lembramos da aula de Física-Matemática II, quando Domingos Marchetti mandou aquela lista de exercícios que pedia para provar propriedades dos polinômios de Chebyshev. Entre elas, o fato de que polinômios de Chebyshev possuem \cos(\frac{k\pi}{n+1}) como raízes. Se eu quero um polinômio com essas raízes, preciso torturar Chebyshev. E eu devia ter desconfiado, veja a relação de recorrência que gera esses polinômios. Seja T_n(x) o polinômio de Chebyshev (de segundo tipo) de grau n, teremos:

 T_{n}(x) = 2xT_{n-1}(x)-T_{n-2}(x)

E me diz se isso não é muito parecido com minha relação de recorrência dos determinantes de A_n! O determinante de A_n não é exatamente T_n(x), mas eu consigo ligar os dois com algumas manobras. O post já está longo e vou escrever direto a relação, é fácil chegar nela operando a relação de recorrência. Teremos

 \det A_n=(2\lambda)^nT_n\left(\frac{1}{2\lambda}\right).

Ou seja, se queremos as raízes do determinante de A_n, basta tomar o inverso das raízes do polinômio de Chebyshev! Sabemos que essas raízes são x_k=\cos(\frac{k\pi}{n+1}). Assim, as raízes de nosso determinante tridiagonal são dadas por:

 \lambda_k=\frac{1}{2\cos\left(\frac{k\pi}{n+1}\right)}

Sem surpresas, esses valores coincidem exatamente com 1, \frac{\sqrt{2}}{2} e razão áurea. Queremos a primeira raiz acima de meio, ela é dada quando calculamos o valor k=1. Dessa forma, o maior valor possível da Berlim gananciosa será

 \lambda_{\text{max}}=\frac{1}{2\cos\left(\frac{\pi}{n+1}\right)}\sim\frac{1}{2}+\frac{\pi^2}{4N^2},

onde expandi para grandes valores de N para ver o comportamento assimptótico de uma cidade grande. Como \lambda=\frac{1}{2}+\frac{\alpha}{2}, confirmamos que, para que esse sistema tenha soluções positivas, existe um limite máximo no lucro dos proprietários que diminui conforme o tamanho da cidade aumenta, e é dado por:

 \alpha<\frac{\pi^2}{2N^2}

Confirmando o que nossa intuição física já sabia naquele cálculo sem rigor acima, a dependência ocorre com o inverso do quadrado do número de habitantes na cidade unidimensional.

E no caso da cidade bidimensional? Não consegui resolver a matriz dessa, não houve Chebyshev que salvasse, mas o cálculo não rigoroso, que deu certo uma vez, parece indicar uma dependência no inverso de N, considerando uma cidade quadrada cujas bordas contém \sqrt{N} casas enfileiradas.

Eu imagino que esse problema tenha sido um pouco mais difícil que o problema médio que apresento aqui, mas é fascinante. Nunca antes eu tinha encontrado polinômios de Chebyshev in the wild, tampouco com uma motivação tão interessante. Vou guardar esse na manga, esperando algum dia algum aluno preguiçoso, como eu era, perguntar-me de que servem esses tais polinômios cujo nome eu tenho que dar ctrl-c ctrl-v para não errar.

E uma nota final: minha namorada encontrou um apartamento finalmente. Fica em uma cidade simpática ao sul de Paris, com acesso ao metrô, mas não tão perto, com alguns metros quadrados, mas não tantos. Foi difícil achar, então comemoramos; mas confesso que, se perguntasse um pouco pelos corredores do prédio, provavelmente ficaria estarrecido em descobrir que o \alpha desse lugar não respeita nenhum majorante.