4.8.2. Linear Density — MDAnalysis.analysis.lineardensity
¶
A tool to compute mass and charge density profiles along the three cartesian axes [xyz] of the simulation cell. Works only for orthorombic, fixed volume cells (thus for simulations in canonical NVT ensemble).
-
class
MDAnalysis.analysis.lineardensity.
LinearDensity
(select, grouping='atoms', binsize=0.25, **kwargs)[source]¶ Linear density profile
Parameters: - select (AtomGroup) – any atomgroup
- grouping (str {'atoms', 'residues', 'segments', 'fragments'}) – Density profiles will be computed on the center of geometry of a selected group of atoms
- binsize (float) – Bin width in Angstrom used to build linear density histograms. Defines the resolution of the resulting density profile (smaller –> higher resolution)
- verbose (bool, optional) – Show detailed progress of the calculation if set to
True
-
results.x.
pos
¶ mass density in [xyz] direction
Type: numpy.ndarray
-
results.x.
pos_std
¶ standard deviation of the mass density in [xyz] direction
Type: numpy.ndarray
-
results.x.
char
¶ charge density in [xyz] direction
Type: numpy.ndarray
-
results.x.
char_std
¶ standard deviation of the charge density in [xyz] direction
Type: numpy.ndarray
Example
First create a
LinearDensity
object by supplying a selection, then use therun()
method. Finally access the results stored in results, i.e. the mass density in the x direction.ldens = LinearDensity(selection) ldens.run() print(ldens.results.x.pos)
New in version 0.14.0.
Changed in version 1.0.0: Support for the
start
,stop
, andstep
keywords has been removed. These should instead be passed toLinearDensity.run()
. Thesave()
method was also removed, you can usenp.savetxt()
ornp.save()
on theLinearDensity.results
dictionary contents instead.Changed in version 1.0.0: Changed selection keyword to select
Changed in version 2.0.0: Results are now instances of
Results
allowing access via key and attribute.