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

Redefined functions relating to collective diffusion coefficient and reduced ionic condictivity. #3

Open
wants to merge 2 commits into
base: master
Choose a base branch
from

Conversation

kitmccoll
Copy link

Left the notebook the same, since it will now be correct with the way the functions are defined.

Left the notebook the same, since it will now be correct with the way the functions are defined.
@bjmorgan
Copy link
Owner

bjmorgan commented Jun 25, 2020

  1. You have some rogue non-ascii characters in one of the docstrings.
  2. I've been comparing these to the definitions in the LLZO paper. In that paper, the "reduced ionic conductivity" is defined as D_J * the site occupancy x of the lattice. What is currently calculated is proportional to this, but is not exactly the same (although units for these are somewhat arbitrary). I am wondering whether we want to add a scaling factor of N_sites, which would make the updated code consistent with Eqn. 3.3, but would rescale the results presented in the paper.
  3. It might be nice to expand the docstrings so that these direct the reader to the relevant equations in the RSOS paper?

Fixed error in docstrings, updated example jupyter notebook and defined a reduced_ionic_conductivity function.
"""
Returns the reduced ionic conductivity, \u03C3\' ( D_J * n_natoms).
Returns the collective or "jump" diffusion coefficient multiplied by the number of atoms,
Copy link
Owner

Choose a reason for hiding this comment

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

This is potentially misleading. This method returns the collective diffusion coefficient, which is equal to the "jump" diffusion coefficient multiplied by the number of atoms, D_J * N.

Comment on lines +283 to +284
(Float): the collective or "jump" diffusion coefficient multiplied by the number of
atoms, ( D_J * number_of_atoms).
Copy link
Owner

Choose a reason for hiding this comment

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

Also needs changing here.

"""
if self.has_run:
return self.atoms.collective_dr_squared() / ( 6.0 * self.lattice.time )
else:
return None

@property
def collective_diffusion_coefficient( self ):
def collective_diffusion_coefficient_per_atom( self ):
"""
Returns the collective or "jump" diffusion coefficient, D_J.
Copy link
Owner

Choose a reason for hiding this comment

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

This is incorrect. Returns the collective diffusion coefficient per atom, which is equal to the "jump" diffusion coefficient.

@@ -298,7 +300,25 @@ def collective_diffusion_coefficient( self ):
(Float): The collective diffusion coefficient, D_J.
Copy link
Owner

Choose a reason for hiding this comment

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

Inconsistent. The collective diffusion coefficient is D_J * N_atoms.

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