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
LinearDensityobject 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, andstepkeywords has been removed. These should instead be passed toLinearDensity.run(). Thesave()method was also removed, you can usenp.savetxt()ornp.save()on theLinearDensity.resultsdictionary contents instead.Changed in version 1.0.0: Changed selection keyword to select
Changed in version 2.0.0: Results are now instances of
Resultsallowing access via key and attribute.