Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Feature: Ground state projection tool for RT-TDDFT #5477

Open
wants to merge 17 commits into
base: develop
Choose a base branch
from

Conversation

ESROAMER
Copy link
Collaborator

Background

For Rt-TDDFT calculations, the ground-state projection is an important feature. By projecting the excited-state wavefunction at time t onto the reference state at time 0, real-time excitation infos can be obtained. Below, we use N2 as an example to demonstrate this feature.
By applying a Gaussian pulse around 2 fs, the electrons in N2 are excited and begin to oscillate. This oscillation can be observed by monitoring the electric dipole.
image
Meanwhile, the number of electrons in the ground state begins to decrease. By summing the number of ground-state electron
image
By analyzing the real-time projection matrix, we can clearly observe three main excitations:
HOMO-1 -> LUMO+1
HOMO-2 -> LUMO
HOMO -> LUMO+2
result
These three excitations can be clearly seen in the absorption spectrum obtained via Fourier transformation.
image

What's changed?

  1. Created new folder 'rt-tddft-tools' The previous absorption tools have been moved there.
  2. Uploaded the ground-state projection tool.
  3. Uploaded the tool for calculating overlap matrices. In cases where atoms are moved, this script must be used to regenerate the overlap matrix before performing the ground-state projection correctly.
  4. Uploaded an example for Si to demonstrate the use of the ground-state projection tool.

Linked Issue

#4925

@ESROAMER ESROAMER changed the title Feature: Ground state projection Feature: Ground state projection tool for RT-TDDFT Nov 13, 2024
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

why do we need this file, and the file is huge.

Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This image is used to display the results of the example. I can convert it to JPEG format, or maybe I should just delete it?

Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

the file is large, too

Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This one will keep only 10 bands and 100 steps of data.

Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

file too large

tools/rt-tddft-tools/overlap_generator.py Outdated Show resolved Hide resolved
symbols = atoms0.get_chemical_symbols()
element_counts = Counter(symbols)
unique_elements = set(symbols)
assert orb_file_num == len(unique_elements), f"轨道数目{orb_file_num}与元素数目{len(unique_elements)}不一致!"
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

do not use Chinese

tools/rt-tddft-tools/overlap_generator.py Outdated Show resolved Hide resolved
tools/rt-tddft-tools/overlap_generator.py Outdated Show resolved Hide resolved
@mohanchen mohanchen added the Features Needed The features are indeed needed, and developers should have sophisticated knowledge label Nov 14, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Features Needed The features are indeed needed, and developers should have sophisticated knowledge
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants