Anais do Simpósio Brasileiro de Química Teórica
Proceedings of the XXI Brazilian Symposium on Theoretical Chemistry - Vol. 1 2021

Parabéns pelo trabalho, muito interessante seus sistemas. Como você obteve seus resultados ? Quais são as teorias e os modelos por trás dela ? Você não falou nada de sua metodologia. Obrigado, sucesso !

Olá Júlio, parabéns pelo trabalho, você foi bastante cuidadoso! Pretende expandir seus resultados para outras propriedades? Já pensou em realizar os cálculos de propriedades óticas não lineares para sólidos?

Oi Júlio, boa noite.

Quero parabenizá-lo pelo excelente trabalho. Adorei.

Oi Mateus, boa noite.

Quero parabenizá-lo pelo excelente trabalho. Adorei.

Oi Graziani, boa noite.

Quero parabenizá-lo pelo excelente trabalho. Adorei.

Boa tarde, Guilherme. Parabéns pelo trabalho!

Essas condições de contorno criadas têm alguma interpretação no seu problema físico?

Obrigado!

Parabéns pelo trabalho Incrível!!!!

Teria como adicionar efeito solvent explícito? Imagino que esse tipo de efeito mudaria a altura da barreira.

Eu tento trabalhar com efeito solvente. Pensei em colocar com densidade de cargas pela metodologia ASEC, caso o Sr. tenha interesse em colaborações. 

Otimo trabalho Guilherme! Fiquei pensando se voce consegue relacionar o overlap a I'm Delta learning na energia.

Olá Guilherme, parabéns pelo trabalho, muito interessante esta abordagem de aprendizado de máquina para estimar os descritores.

Gostaria de saber qual método de estrutura eletrônica foi utilizado para estimar o \Psi_minD.

Desde já agradeço e desejo muito sucesso na sua carreira!

Obrigado!

I am not 100% sure, but I feel like the problem of basis non-orthogonality can have to do with the fact of missing correlation in HF - i.e. it might be compensated if correlation is used.

If that happened to be true - your method might do better with any correlated calculation, e.g. with DFT (including hybrids) or CISD or CCSD.

In case you try it - I'd really like to know if it works.

Also, the DIIS joke gave me a good laugh (given I got it right).

Thanks.

Prezado Caio,

Achei bastante interessante o seu trabalho.

Gostaria de fazer alguns comentários e perguntas. Primeiramente, achei estranho você dizer que não conseguiu convergir o "SCF puro" para a água. Eu já fiz esse cálculo várias vezes e nunca tive problema de convergência. 

Já usei métodos de otimização global para átomos, mas nunca para moléculas. Existe alguma relação entre o seu método e métodos de otimização global? Posso lhe mandar um artigo se lhe interessar.

A função determinante de Slater, embora a forma mais usual de função de onda teste, por sua utilização em métodos Hartree-Fock e DFT, não tem a forma correta para tratamento de sistemas polieletrônicos. As funções DS não respeitam a simetria de permutação do hamiltoniano polieletrônico, ou seja, não são base para o grupo simétrico, como deveriam ser. Posso lhe mandar um artigo se interessar.

Teria como usar métodos de geometria diferencial para otimizar funções de onda que são bases para o grupo simétrico?

Olá Graziâni. Parabéns pelo trabalho! Tenho algumas dúvidas. 1) Pretende realizar TDDFT com campo elétrico para ver o efeito das energias de transição já que o gap homo-lumo não descreve adequadamente a energia de transição? 2) Depois da construção da célula com as moléculas, você rodou dinâmica para "termalizar" o sistema molecular?

Olá Lucelma, excelente trabalho. Os resultados que você mostrou foram para amostras de ZrO2:Yb3+, Er3+? Não ficou claro para mim. Parabéns pelos ótimos resultados.

Olá Felipe, gostei bastante do seu trabalho, principalmente por conta do tamanho do seu conjunto, parabéns. Existem alguns trabalhos que estudam a dopagem de C60 com ácidos carboxilicos para aplicações específicas. Minha pergunta é: Seria possível calcular o pKa de sistemas desse tipo C60-ácido?

Parabéns, Felipe. Belo trabalho. Vocês publicaram esse trabalho completo em algum lugar?

Abs,

Olá Guilherme, seu trabalho é muito interessante. Minha dúvida é, como ficam as propriedades partindo da sua função de onda (momento de dipolo, polarizabilidade...)?

Excelente trabalho e apresentação. Gostei muito do tema que vc está estudando. Algumas vezes o cálculo SCF falha, mesmo em sistemas que possuem a geometria bem definida. As alternativas ao DIIS disponíveis na literatura, como QC, XQC, aumentar o número de iterações ou mesmo reduzir a precisão também podem falhar causando uma situação decepcionante.

Eu gostaria de saber sobre o comportamento do determinante de Slater durante a otimização da função de onda. Será que alguns blocos de valores do determinante de Slater permanecem imutáveis durante o processo de otimização (principalmente os relacionados aos orbitais do caroço)? É possível que seu algoritmo evite modificar tais "blocos imutáveis" a fim de reduzir o custo computacional?

abraços

Olá Cayo, achei o seu trabalho super curioso e interessante! parabéns.   Eu não sou da matemática, então corro o risco de estar sugerindo algo a toa, mas será que daria para trocar o metodo de Newton pelo L-BFGS ou algum outro metodo mais rapido? Talvez com mais avanços incrementais da sua parte isso se torne algo ainda mais competitivo e possa ser implementado em outros codigos (vi na discussao abaixo que tem no Q-Chem). att

Olá, Caio. Seu trabalho é muito bonito! Parabéns pela pesquisa. Tenho uma pergunta: esse tipo de abordagem poderia ser aplicada para problemas de autovalor de matrizes esparças, como os que são tratados com o algoritmo de Davidson? Obrigado!

Caro Caio, muito interessante o trabalho

No HF usual, otimizamos os coeficientes de combinação linear dos orbitais atômicos para obter os orbitais moleculares do determinante de Slater

1) Quando você calcula gradientes e Hessianas o faz em reação aos mesmos coeficientes? Ou utiliza outras variáveis na otimização?

2) O HF usual procura o mínimo tipicamente seguindo o gradiente, enquanto você utiliza a geodésica, certo? Você poderia explicar em termos simples se haveria vantagem em seguir a geodésica?

Obrigado

Marcio

Ola Caio, parabens pela originalidade do seu trabalho. Precisamos de mais trabalhos de base desse tipo. Gostaria de entender melhor a questao do gradiente nulo. Voce aborda na apresentacao que a condicao de dE = 0 nos leva a tres possiveis condicoes (minimo, maximo e ponto de sela). Dai como podemos ter certeza se o calculo do gradiente nos levou ao minino?  

Boa noite, Caio.

Excelente apresentação e, além disso, muito didática. Conseguiu tornar o assunto altamente complexo em um pouco mais inteligível.

Eu gostaria de saber apenas quas foram as suas motivações para implementar um algoritmo para essa pesquisa.

Abraços.

Excelente trabalho, parabéns. Você saberia informar se sua metodologia poderia ser aplicada para análise do tunelamento em junções josephson? São sistemas onde dois supercondutores são separados por uma fina camada isolante e os elétrons tunelam através desta camada entre os supercondutores.

Gabriel, ótimo trabalh. Parabéns! Me tire duas dúvidas. 1) Quando você propõe as duas curvaturas, você imaginava que fossem aparecer espirais ou foi algo não planejado? 2)  Você pretende trabalhar com outras espirais ou propor novas curvaturas de diferentes curvas?

Felipe, muito bom. Quando você fala de conjuntos de treinamento e previsão é pra dizer que você utilizou um algoritmo de aprendizado de máquina?

Boa tarde, Guilherme, gostaria de parabenizá-lo pelo trabalho realizado, achei ele bem interessante. Tenho algumas dúvidas: 1) Estas funções de onda que você otimiza tem significado físico? 2) Vi que no tópico anterior você tinha dito que são gerados pontos aleatórios para a sua construção da função. Com as várias simulações que fez, conseguiu ver se esta função se ajusta ou parece com uma determinada forma funcional? Att

Olá Lucas, belo trabalho! Tenho duas dúvidas. 1) Você usou a rotina CREST do xTB ou pegou as estruturas das trajetórias das dinâmicas e aplicou a "cluster analysis"? 2) Por que você não usou os métodos compostos com funcionais "3c" do Grimme? Costuma ser rápido com resultados razoáveis para moléculas orgânicas e biológicas.

Olá Mariana, bastante interessante o seu trabalho. Tenho algumas questões. 1) Como era de se esperar, a presença de hidrogênio leva a um desvio maior do modelo harmônico devido a alta anarmonicidade de estiramentos C-H e similares. Você não gerou "scaling factor" em função do tipo de vibração ou da região da frequência? Isso seria interessante. 2) Por que você não usou funcionais híbridos duplos como o B2PLYP? Esse funcional apresenta em geral uma precisão maior para geometrias e frequências harmônicas. 3) Como você lidou com as frequências com menos de 100 cm-1 para o cálculo da entalpia? Essas frequências em geral possuem um grande erro numérico. 4) Poderia me dar mais detalhes do cálculo anarmônico? Como você fez para utilizá-lo na correção do harmônico? Pergunto isso porque existem diversos protocolos para isso como VSCF e VPT2.  Além disso, o cálculo anarmônico gera sempre mais bandas, em particular na região acima de 2500 cm-1 (bandas de combinação/overtones).

iOlá Eduardo, parabéns pelo interessante trabalho, sucesso ! E sua abordagem para supercondutores a altas temperaturas - ela pode ser adaptada ?

Parabéns pelo trabalho Mari! Vocês fizeram algum critério na escolha dos 10 funcionais ou os funcionais utilizados foram selecionados aleatoriamente?

Hi Carlos, I have been looking for some time to study something about ADC, of which I knew nothing, and this was a good way start, thanks for the talk!         I have some curiosities about it:         If I understood correctly, what you obtain from ADC are excitation energies directly; You don't have access to total electronic energies, is this right?         You say that it can be improved systematically, and make a parallel with perturbation theory. Does it converge to FCI excitation energies? Or it might diverge, like perturbation theory?          You answered in a question, now for the MR case, that MCSCF is used only to get a ground state, and on top of which ADC comes and give excitation energies. I don't know if this has any practical relevance, but what if you start with a excited MCSCF wave function? I guess that some of these excitation energies (relative to this chosen wave function) will be relative to new states (not obtained from ADC on top of ground state), some other to the same.... Do they compare well? Does it make sense trying ADC starting from an excited state?      Is there any restriction on which type of excited state can be reached, e.g , only dominated by single excitations?        Do you have a control on which state you are targeting? Or they just come in energy order?        Ops, I think that were too many questions... Perhaps the most important: which reference would you suggest to start learning this method?

Parabéns pelo trabalho Pedro! O Prof. Rogério nunca gostou mesmo da ideia de ter que usar o HLC. Você está ao lado de um cara genial, aproveite!

Parabéns Caio, muito interessante. Qual é a demanda computacional do seu método em comparação com o SCF/DIIS? A vantagem do seu método é combiná-lo com o SCF/DIIS para convergência mais rápido ?

Boa tarde, Caio. Parabéns pelo trabalho!

Você calcula as energias nos pontos de mínimo e de máximo usando algum método de estrutura eletrônica e depois faz a regressão para o seu potencial, certo? As energias utilizadas levam em consideração correção de ponto zero ou elas de certa forma se compensam? Isso influencia tanto quanto a "qualidade" das energias eletrônicas obtidas em cada ponto?

Obrigado!

Could the covalency issue be solved via adding dipoles for the bond overlap electron density, in addition to the functional group dipoles?
Large covalency should correspond to significant interatomic (bond) density.

Olá Letícia, tudo bem? Fiquei curioso para saber qual base você usou nos seu cálculos e qual o método de cálculo de integrais moleculares. Esse tempo que você compara, é apenas das iterações SCF, ou também do cálculo de integrais?

Olá, Letícia, tudo bem? Parabéns pela apresentação, foi realmente interessante. Estou para iniciar um trabalho semelhante ao seu e tenho algumas curiosidades.

1) Você implementou o método Hartree-Fock, incluindo a otimização das estruturas das moléculas, em Python e Julia? Se sim, que baita trabalho braçal, difícil, está de parabéns. Se não, o que exatamente fizeste?

2) Você não discute a precisão (ou exatidão) de seus algoritmos, apenas o tempo de computação. Como que ficaram seus resultados, em relação à energia e geometria das moléculas?

Novamente, parabéns pelo trabalho. Meus parabéns.

Oi Caio. Existe alguém que fez algo parecido na literatura?

Olá. Muito legal seu trabalho, parabéns! Não entendo muito desses cálculos de camada interna, mas você comenta que uma vantagem do método ADC é que você consegue calcular essas transições de camada interna sem ter que incluir esses orbitais no espaço ativo. Mas não ficou muito claro para mim o porquê disso ser uma vantagem em cima de incluir esses orbitais no espaço ativo e trabalhar com eles. 

Olá Carlos, ótimo trabalho! Tenho duas dúvidas. 1) Na construção do MR-ADC a energia é invariante às rotações dos orbitais? 2) Se você partir de uma referência CASSCF otimizada para um estado de camada interna, é possível realizar o procedimento com o MR-ADC para tentar simular um espectro Auger (pelo menos as energias)?

Oi Guilherme, tudo bem? Fiquei curioso para saber como é feita a minimização da energia com relação às funções K e W. Pelo que entendi essas funções estão definidas em um grid ao longo de xi/eta, é isso? Você usa sempre o mesmo número de pontos em xi e eta? É um método numérico, tipo gradiente descendente ou Newton? Como é calculada (se é que isso é necessário no seu método) a derivada da energia com relação a essas funções?

Olá João Gabrial, tudo bem? Muito interessante ver a resolução da eq. de Shcrödinger em espirais! Você considerou comparar o seu cálculo com outros métodos de integração numérica, como RK2 e Numerov (esse bastante usado para a resolução da Eq. de Schrödinger) para analisar a sua acurácia?

Oi Mariana, tudo bem? Esses fatores ótimos dependem muito da base escolhida (tanto os estabelecidos da literatura quanto os que você fez)? Qual base você usou nesse caso?

E cálculos com efeito de solvente, sabe se esses fatores mudam muito?

E para cálculos com métodos de função de onda? Há dados na literatura e/ou você planeja fazer esse tipo de estudo para métodos como MP2 ou CCSD?

Nice work, Carlos. I would like to know if the core electron excitation is taken into account in the MCSCF step. If so, how do you treat the problem of variational collapse?