Para citar este trabalho use um dos padrões abaixo:
A geometric interpretation for coupled-cluster theory being more accurate than corresponding configuration interaction
Yuri Alexandre Aoto
Universidade Federal do ABC
Agora você poderia compartilhar comigo suas dúvidas, observações e parabenizações
Crie um tópicoWatch this next:
The geometry associated to the wave function methods of electronic structure is being studied.
The manifold of CC wave functions tends to be closer to the exact wave function than the manifold of CI wave functions.
Newton method is used to find the coupled-cluster wave function closest to the full configuration interaction wave function.
Wagner Richter
Nice presentation indeed! I was wondering if switching from CISD to QCISD should improve the results. I believe the curve for QCISD (as it has the quadratic term) would be very close to the CCSD curve, or at least closer to CCSD than to CISD. Did you (or will you) test this?
Marcio Varella
Caro Yuri,
achei o trabalho muito interessante, não conhecia a abordagem
Eu teria duas perguntas/curiosidades:
1) Não entendi o argumento de porque a excitação quádrupla, existente apenas no CC, pode ser descartada como explicação para seu melhor desempenho em relação ao CI
2) O quão difícil seria generalizar a abordagem para problemas de maior dimensionalidade? (maiores espaços ativos)
Obrigado
Marcio
Yuri Alexandre Aoto
Oi Marcio, obrigado pelas perguntas!
1) o argumento isolado de o CC possuir um determinante quadruplamente excitado não pode ser usado como argumento pois esse determinante a mais na função de onda não possui um coeficiente que pode ser variado livremente, de maneira independente dos outros. Então não existe um grau de liberdade a mais, em relação ao CI, na definição da função de onda. Sendo assim que garantia existe que esse determinante a mais contribuirá para melhorar (e não piorar) a descrição da função de onda? É isso que o trabalho visa elucidar, se essa contribuição a mais de fato ajuda, fazendo que as funções de onda CC (isto é, aquela superfície associada ao CC) esteja mais próxima da função de onda FCI.
2) Bom, a abordagem está generalizada para, a princípio, qualquer dimensão, as imagens e a discussão no início são apenas ilustrativas. Os exemplos Li2 e H4 são dois casos em que temos uma dimensão muito maior. A principal dificuldade para fazer isso é que obter a expressão para os coeficientes para os determinantes da função de onda CC em termos das amplitudes é trabalhoso e custoso computacionalmente. O fato de ser trabalhoso, para minha sorte, foi superado usando um trabalho que mostra como fazer essa transformação (a referência 3 do resumo, achei que não valia a pena entrar nesse detalho no vídeo). O fato de ser custoso computacionalmente faz com que não dê para ir para sistemas ou bases muito grandes (8 elétrons e base 6-31G foi o máximo que consegui). Como isso tudo depende de calcular antes a função de onda FCI, que não é viável para sistemas muito grandes de qualquer jeito, não é um problema tão drástico.
Yuri
Antonio Gustavo Sampaio de Oliveira Filho
Oi Yuri,
Muito bom o trabalho. Parabéns.
Você chegou a avaliar a geometria da solução por teoria de perturbação MP?
Esse estudo da forma do subespaço da solução aproximada poderia levar a proposta novos ansätze que sejam mais próximos da solução exata?
Abraços,
Antonio
Yuri Alexandre Aoto
Oi Antonio! Uma abordagem para estudar teoria MP vai na direção da pergunta do Wagner, porque aí a maneira como as amplitudes são obtidas está muito entranhada na proposta da teoria. Por exemplo, na função de onda MP1 apenas os determinantes duplamente excitados entram na função de onda, assim como no CCD. Então é necessário também realizar a interpretação geométrica para as equações para calcular a função de onda em si. Isso certamente é um dos alvos, e estamos estudando CC nessa direção.
Para funções de onda em ordens mais altas de perturbação, onde começa a haver relação entre amplitudes de diferentes ordens, isso começa a ficar mais interessante, mas ainda muito relacionado às equações para determinar a função de onda.
Essa sua pergunta seria a aplicação mais atrativa desse estudo, isto é, desenvolver novos métodos justamente focado nas características desses conjuntos, por exemplo estarem mais próximos da solução FCI. Ainda não chegamos lá, mas esperamos que nossos estudos sugiram rotas para isso.
Abraços!
Annika Bande
Hello Yuri and greetings to Brazil,
In the slide on 2 x H_2 (maybe also for Li_2) you show the comparison in energies of CI and CC with FCI. I understand this as the lack of correlation energy captured in the respective truncation scheme. Is that so? Then you also show the comparison of the distances that you found. For CI the two curves came out very similarly, but for the CC case at intermediate distances the curvatures differed. Do you know the reason for this?
Cheers from Berlin, Annika
Yuri Alexandre Aoto
Hi Annika, thank you for the question! "the lack of correlation energy captured in the respective truncation scheme. Is that so?": Yes, that is right. Concerning the behaviour at intermediate distance for 2 H2, it is indeed different, and this increases with basis set (for even smaller basis sets this is less prominent). The easy answer is that these curves are associated with different quantities, energy and distance of wave functions, and thus a very fine comparison of details must be made with care, for instance there are other metrics that could be used for the distance, and the curvature would change. With this in mind, the top graph, with the distance, shows that the FCI wave function approaches the CCSD manifold rather slowly, but the CCSD energy is corrected more quickly. This CCSD energy is regular, traditional CCSD. The point of the CCSD manifold that minimises the distance to FCI (whose distance to FCI is what is plotted) is another point, let's call it CCSD_minD. I actually calculated the distance from CCSD_minD to the regular CCSD wave function, and this is relatively large in this region, what indicates that, although CCSD_minD is the point closest to the FCI wave function (under the metric we used), is not necessarily the "best" CCSD for the energy, and the amplitudes equation can take this into account. I can't (yet) calculate the energy for the CCSD_minD with the same integrals as done for the CCSD and FCI energies, but my guess is that, if done, the correlation energy recovered by CCSD_minD would behave much more like the distance.
Com ~200 mil publicações revisadas por pesquisadores do mundo todo, o Galoá impulsiona cientistas na descoberta de pesquisas de ponta por meio de nossa plataforma indexada.
Confira nossos produtos e como podemos ajudá-lo a dar mais alcance para sua pesquisa:
Esse proceedings é identificado por um DOI , para usar em citações ou referências bibliográficas. Atenção: este não é um DOI para o jornal e, como tal, não pode ser usado em Lattes para identificar um trabalho específico.
Verifique o link "Como citar" na página do trabalho, para ver como citar corretamente o artigo
Yuri Alexandre Aoto
Hi Wagner, thank you for the great question. The problem with QCI is that its equations are not obtained as in CI or CC, where one starts with a well defined wave function, and thus with a corresponding manifold. In QCI one starts with the equations for CI (or CC), and add some terms (or remove some terms) in the equations to determine the amplitudes. This means that a geometric analysis must include a geometric interpretation for the amplitudes equations as well. It is possible, but more involved. I have been doing this for different variants of CC only, but definitely considering QCI is also a great idea.