Skip to content

add the gamma-point hessian with phonopy - #860

Open
puzhichen wants to merge 1 commit into
pyscf:masterfrom
bytedance:feat/pbc-gamma-hessian
Open

add the gamma-point hessian with phonopy#860
puzhichen wants to merge 1 commit into
pyscf:masterfrom
bytedance:feat/pbc-gamma-hessian

Conversation

@puzhichen

Copy link
Copy Markdown
Collaborator

No description provided.

@puzhichen
puzhichen marked this pull request as draft August 17, 2026 08:29
@puzhichen
puzhichen marked this pull request as ready for review August 18, 2026 06:26
dm0 = mf.make_rdm1()
try:
for index, displaced in enumerate(displaced_cells, 1):
cell.set_geom_(displaced.positions, unit="Angstrom")

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

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

cell = cell.set_geom_(displaced.positions, unit="Angstrom", inplace=False)

)
return self.force_constants
finally:
cell.set_geom_(

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

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

By using inplace=False above, restoring original geometry is not needed

# -----------------------------------------------------------------
# VASP | Angstrom AMU eV/Angstrom eV/Angstrom^2
unitcell = PhonopyAtoms(
symbols=[cell.atom_symbol(i) for i in range(natm)],

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

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

symbols=cell.elements

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants