To cite this paper use one of the standards below:
Comparing Transition Potential and Linear Response TD-DFT Methods to Simulate X-ray Absorption Spectra for CO Desorbing from Ru(0001)
Gabriel Libanio Silva Rodrigues
Department of Physics / Stockholm University
Now you could share with me your questions, observations and congratulations
Create a topicTransition Potential DFT (TP-DFT) is often used for periodic XAS simulation because it does not need the calculations of so many excited-states, decreasing the computational cost.
TP-DFT has shown some drawbacks for simulating XAS for XFEL-related (pump and probe) experiments.
The TP-DFT has not been so investigated for this type of systems when you have desorbing adsorbates on surfaces.
Due to the local nature of the chemical bond finite-size clusters with linear response methods might be a viable alternative.
We showed that the cheap cluster model + linear response TD-DFT approach can be valid to simulate XAS on CO desorbing from Ru(0001).
This can open a new path to study these systems with more accurate linear-response based methods like ADC or EOM-CC in conjunction with cluster models.
Ricardo Oliveira
Hi Gabriel, very nice work! Probably TDDFT do not reproduce the correct vertical energy. Did you use any energy shift in your theoretical spectra? If so, do you use the same shift along the CO-surface distance?
BRENO RODRIGUES LAMAGHERE GALVAO
Hi Gabriel! I enjoyed a lot watching your video and learned a few things about Xray spectroscopic! You are going in an insightful research direction, instead of just crunching out numbers. I appreciate that, and frankly I wish I had the opportunity of doing that more often in my group. Congrats! The clusters are rigid at their experimental bulk equilibrium distances?
Gabriel Libanio Silva Rodrigues
Hi Breno! Thank you very much for the comments and I am often interested in giving (or trying to) this more insightful talks that could evolve into more technical aspects if necessary. Answering your question, the clusters are indeed frozen, but at their computationally optimized geometries calculated with periodic boundary conditions at RPBE level. Additionally, we have optimized the C-O bond length in the clusters at each calculated step of the desorption, as this is one of the most relevant parameters for the O1s-XAS, as it was shown in the presentation.
With nearly 200,000 papers published, Galoá empowers scholars to share and discover cutting-edge research through our streamlined and accessible academic publishing platform.
Learn more about our products:
This proceedings is identified by a DOI , for use in citations or bibliographic references. Attention: this is not a DOI for the paper and as such cannot be used in Lattes to identify a particular work.
Check the link "How to cite" in the paper's page, to see how to properly cite the paper
Gabriel Libanio Silva Rodrigues
Hi Ricardo! Thank you for the interest. Indeed, TDDFT always have an error in the absolute transition energies. Fortunately, this error is systematic and dependant on the chosen functional so we can always shift the spectra to compare with experimental data like you said. In this case it was not necessary because the observation of the peak's chemical shifts and the shape of the spectra was sufficient comparison. Therefore, all displayed TDDFT spectra is shown as it was computed with no shifts and the absolute values of the transition energies, which are underestimated if compared to the experiment. However, if we do an average of the peak's energies weighted by their intensities, we can see that the chemisorbed- to gas- phase chemical shift is quite constant, although a bit larger than experiment. We have all this data in the paper that it is already submitted, but let me know if you need anything else.
Ricardo Oliveira
Thank you! I have another question, how Tamm-Dancoff approximation can affect the intensities of CO+surface system?
Gabriel Libanio Silva Rodrigues
I have made calculations with both RPA and TDA and I have not seen any difference at all between the calculated spectra. At least not for the main 1s->pi* transition.