O que acontece quando uma molécula não tem dipolo permanente
A grande maioria das moléculas orgânicas comuns — metano, dióxido de carbono, benzeno, Cl — são apolares na prática. Elas não têm carga positiva e negativa permanentemente separadas. Quando você aproxima um íon ou uma molécula polar dessas substâncias, o campo elétrico externo distorce a nuvem eletrônica. O centro de cargas negativas se desloca em relação ao centro de cargas positivas. Isso gera um dipolo induzido. Simples assim. Eu aprendi isso na prática num projeto de simulação de dispersão de gases nobres em zeólitas. A gente tentava modelar a interação do xenônio com a estrutura porosa e o software simplesmente não convergia. O problema era que a parametrização padrão do campo de forças tratava o Xe como átomo sem carga, sem dipolo permanente, e ignorava completamente a polarizabilidade. A energia de adsorção saía totalmente errada — tipo 40% abaixo do valor experimental. A solução foi ativar o termo de indução no modelo, usando o potencial de Drude oscilante ou, numa abordagem mais simples, adicionar o termo de London (dispersão) explicitamente com coeficientes C6 derivados da polarizabilidade experimental de cada espécie. Aí sim a simulação casou com os dados de isotermas reais.
Como calcular a energia de interação com dipolos induzidos
A energia de interação entre um íon de carga q e uma molécula com polarizabilidade , a uma distância r, segue a lei clássica: E_induzido = q² / (8²r)
Para a interação entre duas moléculas apolares, o termo de dispersão de London domina: E_disp = 3/2 · (II)/(I+I) · ()/((4)²r)
Onde I é a energia de ionização primeira de cada espécie. Em simulações de dinâmica molecular, isso costuma ser implementado via campos de forças polarizáveis (Drude, AMOEBA) ou corrigido pós-simulação com a abordagem de indução contínua. Num computador comum, calcular esses termos diretamente para cada par de partículas escala como O(N²) e pode.triplicar o tempo de simulação. A workaround prática que eu uso é o método de reação de campo (reaction field), que trata o meio exterior como um dielétrico contínuo com constante . Isso reduz a complexidade para O(N log N) usando listas de vizinhança e cortes de corte de 12 a 14 Å — suficiente para a maioria dos sistemas líquidos. Um detalhe que os livros não enfatizam bastante: a polarizabilidade não é constante. Ela depende do ambiente eletrônico local. Num solvente polar como a água, a polarizabilidade efetiva de uma molécula de O dissolved pode ser até 15% maior do que no vácuo, porque as moléculas de água ao redor comprimem a nuvem eletrônica e facilitam a deformação. Se você estiver fazendo simulações de alta precisão, use valores de polarizabilidade derivadas de cálculosQM level (CCSD(T) ou DFT com funcional adequade como B97X-D) e não valores tabelados do Handbook.
Aplicações práticas que importam de verdade
Dipolos induzidos são responsáveis por forças de van der Waals, que por sua vez governam pontos de ebulição de gases nobres, solubilidade de não-polares em solventes apolares, e a própria estabilidade de proteínas via empacotamento hidrofóbico. Sem eles, o metano seria gasoso em temperaturas extremamente mais baixas e a água não formaria gotas — literalmente, a tensão superficial desapareceria. Na indústria farmacêutica, negligenciar a indução leva a erros sistemáticos na previsão de afinidade de ligação. Um inibidor enzimático que interage com resíduo de fenilalanina via empacotamento - depende essencialmente de dipolos induzidos mútuos. Campos de forças fixos (como o antigo CHARMM22) subestimam essa interação em média 2-3 kcal/mol por contacto. Versões mais recentes com polarizabilidade embutida (CHARMM-Drude) corrigem isso, mas o custo computacional dobra ou triplica. A decisão entre um e outro depende do que você está disposto a perder em troca de acurácia.
Erros comuns que eu já vi acontecerem
O primeiro erro frequente é confundir dipolo induzido com dipolo permanente. CO é apolar porque os dois vetores C=O se cancelam geométricamente. Mas ele tem polarizabilidade anisotrópica — maior perpendicular ao eixo molecular do que ao longo dele. Se você modelar CO como uma esfera isotrópica com médio, estará errando interações com superfícies carregadas. O correto é usar um tensor de polarizabilidade de 3x3. Outro erro crasso: assumir que dipolos induzidos são sempre atrativos. Eles são, desde que o campo externo seja estático. Mas com campos oscilantes — luz, radiação EM — a resposta da nuvem eletrônica tem fase defasada, e surgem forças de dispersão retardada (forças de Casimir-Polder), que decaem como 1/r em vez de 1/r. Isso é relevante em nanoindentação e microscopia de força atômica, mas a maioria dos pacotes de simulação moleculares ignora esse regime completamente.
Se você está começando agora e quer só uma referência rápida para consultar valores de polarizabilidade, o banco de dados do NIST Chemistry WebBook tem uma seção dedicada com dados experimentais e calculados para mais de 3.000 espécies. É o ponto de partida mais confiável que eu conheço. O resto é exercício de aplicação.
