Para citar este trabalho use um dos padrões abaixo:
Simulating X-ray Photoelectron Spectra using Multireference Algebraic Diagrammatic Construction Theory
Carlos Eduardo Vieira de Moura
The Ohio State University
Agora você poderia compartilhar comigo suas dúvidas, observações e parabenizações
Crie um tópicoCore excited states in multireference systems require accurate ab initio methods. We developed a multireference algebraic diagrammatic construction (MR-ADC) for core ionization. Our implementation access excitations using the Core-Valence Separation (CVS) approximation. CVS-IP-MR-ADC presents good agreement with experiments for several multireference problems. Simulation of X-ray photoelectron spectra, including satellite states, can be accurately achieved.
Alexandre Rocha
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?
Amanda Dornela Torres
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.
Carlos Eduardo Vieira de Moura
Olá, Amanda. Obrigado pela pergunta! Esta vantagem é baseada na comparação com a aplicação de métodos como CASSCF especificamente para estados de camada interna. Para estes métodos, uma abordagem utilizada é incluir os orbitais de camada interna no seu espaço ativo e restringir a ocupação do mesmo. Restringir a ocupação com um elétron a menos, por exemplo, pode ser feito para descrever uma ionização de caroço. Este tipo de construção, por sua vez, demanda estratégias de otimização específicas (e aproximadas) para a função CASSCF. Como os estados de camada interna são estados excitados de altas energias, métodos variacionais mais diretos para minimização destas funções de onda precisariam do cálculo de muitas raízes antes de chegar nos estados de camada interna. Então, são necessárias restrições nas rotações dos orbitais durante o cálculo CASSCF, por exemplo. Desta maneira, você consegue recuperar energia de correlação eletrônica referente a esses elétrons, adicionando efeitos de relaxação, que são bastante importante para estes estados e significativos nos valores de potencial de ionização obtidos, por exemplo.
Ao aplicar métodos de propagadores, como o MR-ADC, os efeitos de correlação eletrônica são recuperados independente da inclusão destes orbitais no espaço ativo, o que permite a obtenção destes estados com boa qualidade e custo computacional menor do que métodos CC ou MRCI. A função CASSCF de referência pode ser escolhida, então, de forma a melhor reproduzir o caráter multireferência do seu sistema, como nos exemplos que mostrei na apresentação (O3 e benzinos).
É uma pergunta muito boa e a resposta ficou um pouco longa. Ainda tem outros aspectos que podem ser levados em consideração nesta comparação. Se ainda ficar alguma dúvida, fico feliz em tentar responder novamente!
Amanda Dornela Torres
Entendi! Muito obrigada pela resposta e parabéns de novo pelo trabalho.
Ricardo Oliveira
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)?
Carlos Eduardo Vieira de Moura
Olá, Ricardo! Obrigado pelas perguntas!
1) Os estados excitados obtidos devem sim ser invariantes às rotações dos orbitais do espaço ativo completo.
2) Essa é uma ideia bastante interessante. A implementação atual do método suporta apenas o CASSCF/CASCI como referência e uma generalização seria necessária para adicionar outras construções, como RASSCF. Isto também permitiria algumas outras aplicações. Uma forma de, de fato, simular os espectros Auger, seria através da implementação de um método como Feshbach-Fano, utilizando, por exemplo, ondas planas para descrever o contínuo. É um desenvolvimento que pretendo trabalhar, em algum futuro próximo!
Yuri Alexandre Aoto
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?
Carlos Eduardo Vieira de Moura
Hello, Yuri! Thank you for all your questions. I will try to answer them:
- You're right about ADC: it obtains directly the excitation energies and doesn't have access to the total energy. For the single-reference version of ADC, for example, the ground-state energy can be associated with MP2 energy, but there's some discussion about how accurate this assumption really is when you step back to the complete form of the propagator.
- About the perturbation theory: ADC converges to the FCI with increasing perturbation order. I think that the reason ADC is different from MP, in this case, is that ADC doesn't have the effect of intruder states during the perturbative expansion (but I will try to formulate it better)
- About the reference wavefunction in the excited state: The polarization propagator is defined in terms of the ground-state wavefunction. So, I am not sure if this approach is correct. On the other hand, it could maybe be an approach for processes like two-photon absorptions or Auger decay.
- About restrictions of excited states: there's no limitation related to single excitations, for example. We are interested in implementing and investigating satellite states, for example, which are essentially double-excitations.
- About the control of the target state: It's a multiple state method, so it doesn't a state-specific method. They're obtained according to the number of roots you ask. One actual limitation that we are improving in the MR-ADC is the guess dependence for the Davidson procedure, specially for high excitations.
- About the ADC Theory: I think a good start is this paper for Andreas Dreuw (doi:10.1002/wcms.1206) that can give you a general picture of the method. Then, the papers from Jochen Schirmer can give you more mathematical details. For MR-ADC, the first paper is very detailed about the method implementation and equations (doi: 10.1063/1.5055380).
Best regards!
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
Carlos Eduardo Vieira de Moura
Hello, Prof. Alexandre. Thanks for your question. In fact, the MCSCF reference calculation aims at the description of the ground state and, then, the excitation processes are generated by the creation and annihilation operators which acts overall orbitals, not only in the active space. Then, the inclusion of core orbitals to the active space is not required and there's no variational collapse, as usually found in self-consistent procedures for high energy excited states. It's possible to add the core orbitals to the MCSCF reference active space and then apply the CVS approach, but it probably doesn't give considerable improvements, since the MCSCF should not recover much correlation energy for them in the ground state. An advantage of the MR-ADC method is that the electronic correlation is heavily recovered, including orbital relaxation effects, which are crucial to core states, without specific self-consistent procedures in the MCSCF step.
Best regards!