Onde as leis da mecânica quântica se quebram em simulações?

29

Como alguém que é bacharel em física, fiquei um pouco escandalizado quando comecei a trabalhar com simulações moleculares. Foi um pouco chocante descobrir que mesmo as simulações mais detalhadas e computacionalmente caras não podem reproduzir quantitativamente todo o comportamento da água a partir dos primeiros princípios.

Anteriormente, eu tinha a impressão de que as leis básicas da mecânica quântica eram um problema resolvido (além da gravidade, que geralmente se supõe irrelevante em escala molecular). No entanto, parece que uma vez que você tenta escalar essas leis e aplicá-las a algo maior ou mais complexo que um átomo de hidrogênio, seu poder preditivo começa a quebrar.

Do ponto de vista matemático, entendo que as funções de onda rapidamente se tornam muito complicadas para resolver e que aproximações (como Born-Oppenheimer) são necessárias para tornar as funções de onda mais tratáveis. Entendo também que essas aproximações introduzem erros que se propagam cada vez mais à medida que o tempo e as escalas espaciais do sistema em estudo aumentam.

Qual é a natureza do maior e mais significativo desses erros de aproximação? Como posso obter uma compreensão intuitiva desses erros? Mais importante, como podemos avançar para um método ab-initio que nos permita simular com precisão moléculas inteiras e populações de moléculas? Quais são os maiores problemas não resolvidos que estão impedindo as pessoas de desenvolver esses tipos de simulações?

tel
fonte
6
Er ... o que todos fizeram você pensar que "as leis básicas da mecânica quântica eram um problema resolvido" equivalia a ser capaz de "reproduzir o comportamento completo da água a partir dos primeiros princípios [em simulação]" ? É um problema de treze corpos.
dmckee
@ DMCkee ver, isso é exatamente o que eu estou confuso. 13 problema corporal não significa solução analítica, com certeza, mas o que está nos impedindo de encontrar uma solução numérica de precisão arbitrária? Simplesmente você atinge a parede do que é computacionalmente possível? Você já está no ponto em que uma computação exige que a vida de um sol seja concluída? Em caso afirmativo, que tipos de aproximações você pode fazer para simplificar o problema? Você consegue entender essas aproximações em um nível intuitivo? Existem maneiras de melhorar as aproximações, reduzir o nível de erro que elas introduzem? Quebrá-lo para mim
tel
@dmckee pelo que me fez pensar que a água deveria ser simples em primeiro lugar ... Eu culpo os simuladores de proteínas. Eles me fizeram sonhar com o que foi possível :)
tel

Respostas:

16

Tanto quanto sei, os métodos mais precisos para cálculos estáticos são a interação de configuração completa com um Dirac Hamiltonian de quatro componentes totalmente relativístico e um conjunto de bases "suficientemente completo". Eu não sou especialista nessa área em particular, mas pelo que sei do método, resolvê-lo usando um método variacional (em vez de um método baseado em Monte-Carlo) aumenta muito, pois acho que o número de determinantes de Slater que você tem para incluir em sua matriz dimensiona algo como . (Há um artigo sobre o custo computacional aquiO(norbsCne).) Os métodos de Monte-Carlo relacionados e os métodos baseados neles usando "andadores" e redes de determinantes podem fornecer resultados mais rapidamente, mas, como está implícito acima, não são variacionais. E ainda são terrivelmente caros.

As aproximações atualmente em uso prático apenas para energias de mais de dois átomos incluem:

  • Nascido em Oppenheimer, como você diz: isso quase nunca é um problema, a menos que seu sistema envolva o tunelamento de átomos de hidrogênio ou a menos que você esteja muito perto de uma travessia estadual / travessia evitada. (Veja, por exemplo, interseções cônicas.) Conceitualmente, existem métodos não adiabáticos para a função de onda / densidade, incluindo a CPMD, e também o Path-Integral MD, que pode explicar os efeitos do tunelamento nuclear.
  • Cálculos não-relativísticos e aproximações de dois componentes da equação de Dirac: você pode obter uma formulação exata de dois componentes da equação de Dirac, mas mais praticamente a aproximação regular de ordem de Zeroth (ver Lenthe et al, JChemPhys, 1993) ou Douglas- O Hamiltonian de Kroll-Hess (veja Reiher, ComputMolSci, 2012) é comumente usado e frequentemente (provavelmente geralmente) negligenciando o acoplamento de rotação-órbita.
  • Conjuntos de bases e LCAO: os conjuntos de bases não são perfeitos, mas você sempre pode torná-los mais completos.
  • Funcionais de DFT, que tendem a tentar fornecer uma tentativa suficientemente boa de troca e correlação sem o custo computacional dos métodos mais avançados abaixo. (E que vêm em alguns níveis diferentes de aproximação. A LDA é a de nível básico, GGA, metaGGA e a troca exata inclui mais do que isso, e incluir o RPA ainda é uma técnica bastante cara e nova, tanto quanto eu. Também estou ciente de que existem funcionais que usam técnicas diferentes em função da separação e alguns que usam a vorticidade que acho que têm aplicação em estudos magnéticos ou de aromaticidade.) (B3LYP, o funcional que algumas pessoas amam e outras adoram odiar, é um GGA que inclui uma porcentagem da troca exata.)
  • Truncamentos de interação de configuração: CIS, CISD, CISDT, CISD (T), CASSCF, RASSCF, etc. Essas são todas aproximações ao IC que assumem que os determinantes excitados mais importantes são os menos excitados.
  • Interação de configuração de várias referências (truncamentos): Idem, mas com alguns estados de referência inicial diferentes.
  • Método Coupled-Cluster: não pretendo entender corretamente como isso funciona, mas obtém resultados semelhantes aos truncamentos da Interação de configuração com o benefício da consistência de tamanho (por exemplo, ) (em grande separação)).E(H2)×2=E((H2)2

Para dinâmica, muitas das aproximações referem-se a coisas como o tamanho limitado de um sistema tratável e a opção prática de escalonamento de tempo - são coisas bastante padrão no campo de simulação numérica de tempo. Também há manutenção de temperatura (consulte termostatos Nose-Hoover ou Langevin). No entanto, isso é principalmente um conjunto de problemas de mecânica estatística, como eu o entendo.

De qualquer forma, se você tem uma mente física, pode ter uma idéia muito boa do que é negligenciado observando as formulações e os papéis sobre esses métodos: os métodos mais comumente usados ​​terão pelo menos um ou dois papéis que não são a especificação original explicando sua formulação e o que ela inclui. Ou você pode apenas conversar com pessoas que os usam. (As pessoas que estudam sistemas periódicos com DFT estão sempre resmungando sobre o que os diferentes funcionais fazem e não incluem e explicam.) Muito poucos dos métodos têm omissões surpreendentes específicas ou modos de falha. O problema mais difícil parece ser o tratamento adequado da correlação eletrônica, e qualquer coisa acima do método Hartree-Fock, que não explica nada disso, é uma tentativa de incluí-lo.

Pelo que entendi, obter a precisão do IC relativístico completo com conjuntos completos de bases nunca será barato sem reinventar dramaticamente (ou jogar fora) os algoritmos que usamos atualmente. (E para as pessoas que dizem que DFT é a solução para tudo, estou aguardando suas formulações de densidade pura sem orbitais.)

Há também o problema de que quanto mais preciso você faz sua simulação, incluindo mais contribuições e formulações mais complexas, mais difícil é realmente fazer alguma coisa. Por exemplo, o acoplamento de órbita de rotação às vezes é evitado apenas porque torna tudo mais complicado de analisar (mas às vezes também porque tem efeito desprezível), e os orbitais canônicos Hartree-Fock ou Kohn-Sham podem ser bastante úteis para entender os recursos qualitativos de um sistema sem estratificar a saída adicional de métodos mais avançados.

(Espero que algo disso faça sentido, provavelmente seja um pouco irregular. E provavelmente perdi a aproximação ou a coisinha favorita de alguém.)

Aesin
fonte
11

O(Ne3.7)NeNe=104N=39

A questão principal será que, além do aumento da potência computacional, você precisará criar algoritmos melhores que possam derrubar o expoente 3.7 para algo mais gerenciável.

aeismail
fonte
6
O(Ne3.7)
11
Eu realmente gosto e aprecio essa discussão!
Abrir o caminho
Meu entendimento é que a mecânica quântica (ou pelo menos a teoria da estrutura eletrônica) seria considerada um problema resolvido se os métodos mais precisos fossem escalados como O (N ^ 3). O problema é que são essencialmente apenas os piores métodos, aproximações médias de campo, que se aproximam dessa escala e algo como o IC completo é escalonado exponencialmente com o número de elétrons (ou mais tipicamente as funções básicas).
Tyberius
6

O problema é amplamente equivalente à diferença entre computadores clássicos e computadores quânticos. Computadores clássicos trabalham com valores únicos de uma só vez, pois apenas um futuro / histórico é possível para uma entrada determinística. No entanto, um computador quântico pode operar em todas as entradas possíveis simultaneamente, porque pode ser colocado em uma superposição de todos os estados possíveis.

Da mesma forma, um computador clássico precisa calcular todas as propriedades individualmente, mas o sistema quântico que está simulando possui todas as leis do universo para calcular todas as propriedades simultaneamente.

O problema é exacerbado pela maneira como temos que passar dados quase em série através de uma CPU, ou no máximo alguns milhares de CPUs. Por outro lado, o universo tem um conjunto quase ilimitado de cálculos simultâneos acontecendo ao mesmo tempo.

Considere como exemplo 3 elétrons em uma caixa. Um computador precisa escolher um passo no tempo (primeira aproximação) e continuar recalcular as interações de cada elétron entre si, por meio de um número limitado de CPUs. Na realidade, os elétrons têm um número incognoscível de partículas de troca reais e virtuais em trânsito, sendo absorvidos e emitidos, como um processo contínuo. Cada partícula e ponto no espaço tem alguma interação acontecendo, o que precisaria de um computador para simular.

Simulação é realmente a arte de escolher suas aproximações e seus algoritmos para modelar o assunto da melhor maneira possível com os recursos disponíveis. Se você quer perfeição, receio que seja a matemática de galinhas esféricas nos aspiradores; só podemos simular perfeitamente o muito simples.

Phil H
fonte
3
muito bom "Simulação é realmente a arte de escolher suas aproximações e seus algoritmos para modelar o assunto, bem como possível com os recursos que você tem disponível"
Abrir o caminho
É verdade que apenas fetichistas esféricos de frango se preocupam com a perfeição. A verdadeira questão é o que nos impede de chegar ao "bom o suficiente"? Para muitos problemas de interesse biológico (ou seja, todos os problemas de ligação a medicamentos de todos os tempos), precisos o suficiente seriam calcular as energias para aproximadamente 1 kT ou mais. Isso às vezes é chamado de "precisão química".
tel
11
@tel: Depende da área. Para algumas coisas, temos mais precisão nos modelos do que podemos alcançar na prática, por exemplo, modelando orbitais de elétrons de hidrogênio. Para outros, geralmente sistemas não lineares de muitos corpos, onde múltiplos efeitos entram em jogo, lutamos para combinar o experimento; química quântica para coisas como energias de ligação (consulte Teoria Funcional da Densidade), dobra de proteínas, esses são lugares onde ainda não podemos reproduzir de maneira confiável o experimento com os recursos disponíveis. Computadores quânticos de tamanho razoável fariam o trabalho.
Phil H
4

Não sei se o seguinte ajuda, mas para mim foi muito esclarecedor visualizar o comportamento de escala dos sistemas quânticos:

O principal problema vem do fato de que o espaço de Hilbert dos estados quânticos cresce exponencialmente com o número de partículas. Isso pode ser visto com muita facilidade em sistemas discretos. Pense em um par de poços em potencial que estão conectados entre si, podem apenas dois: poço 1 e poço 2. Agora adicione bósons (por exemplo, Rubidium 87, apenas como exemplo), primeiro apenas um. Quantos vetores de base possíveis existem?

  • vetor base 1: bóson no poço 1
  • vetor base 2: bóson no poço 2

|1 1,0 0|0 0,1 1

Agora, suponha que o bóson possa saltar (ou escavar um túnel) de um poço para o outro. O Hamiltoniano que descreve o sistema pode então ser escrito como notação matricial, como

H^=(ϵ1 1ttϵ2)

Onde ϵ1 1,2|1 1,0 0|0 0,1 1

Esse problema é tão simples que pode ser resolvido manualmente.

Agora, suponha que tenhamos mais poços em potencial e mais bósons, por exemplo, no caso de quatro poços com dois bósons, há 10 possibilidades diferentes para distribuir os bósons entre os poços. Então o hamiltoniano teria 10x10 = 100 elementos e 10 eigenstates.

número de eigenstates=(número de poços+número de bósons-1 1número de bósons)

92,3782

2.7105310107 elementos, ocupando tanto espaço que precisaríamos de todas as partículas de 10 milhões de universos como o nosso apenas para codificar essa informação.

Robert
fonte
2

n3n

pontos-1 12

A teoria funcional da densidade é outra maneira de lidar com esse problema, mas é uma aproximação. É uma aproximação muito boa em alguns casos, mas em outros casos pode ser surpreendentemente ruim.

Dan
fonte
-1

Esse problema é resolvido pela teoria funcional de densidade. A essência está substituindo muitos graus de liberdade do corpo por vários campos, um deles sendo a densidade de elétrons. Para uma grande exposição, consulte a palestra nobel de um dos fundadores da DFT: http://www.nobelprize.org/nobel_prizes/chemistry/laureates/1998/kohn-lecture.pdf

Artan
fonte
Você poderia dar algum contexto ao link que está fornecendo? Nós desencorajamos respostas que apenas fornecem um link sem qualquer tipo de explicação, e esse tipo de resposta é excluído, a menos que seja editado.
precisa saber é o seguinte
10
e, a propósito, você realmente deve tomar cuidado com "Este problema é resolvido por ....". Uma vez que existem limites para DFT que alguém deve mencionar
Abrir o caminho
3
O DFT fornece uma aproximação muito útil, mas não 'resolve' nada! Não é exato sem funcionais exatos para a troca e correlação e, mesmo assim, não produz as funções de onda, mas a densidade de elétrons.
Phil H
Muitos QM do corpo não se decompõem como uma teoria, é apenas NP difícil. DFT é uma teoria com complexidade polinomial que resolve com a mesma precisão dos princípios básicos QM a estrutura eletrônica de todos os elementos químicos. É por isso que ganhou o Prêmio Nobel de química. Ele forneceu excelentes resultados para grandes sistemas quando comparado a experimentos.
25412 Artan
2
Você está errado. A DFT não resolve "o problema" com a mesma precisão. Ele "resolve" um caso específico (estado fundamental), introduzindo uma funcionalidade de correlação de troca completamente desconhecida.
M13