DPDSimulator

class pylimer_tools_cpp.DPDSimulator(self: pylimer_tools_cpp.DPDSimulator, universe: pylimer_tools_cpp.Universe, crosslinker_type: int = 2, slipspring_bond_type: int = 9, is_2d: bool = False, seed: str = '')

Bases: pybind11_object

A quick-and-dirty implementation of the dissipative particle dynamics (DPD) simulation with slip-springs as presented by Langeloth et al. [LMBohmMullerPlathe13] and Schneider et al. [SFKarimiVarzanehMullerPlathe21].

Get an instance of this class

Methods Summary

assume_box_large_enough(self)	Configure whether to run PBC on the bonds or not.
config_a(self[, A])	Configure the force-field (pair-style) parameter A.
config_allow_relocation_in_network(self[, ...])	Configure whether a relocation step may happen when a slip-spring has ended at a crosslink.
config_auto_correlator_output(self, values)	Set which values to compute multiple-tau autocorrelation for.
config_average_output(self, values)	Set which values to compute averages for.
config_bond_formation(self, ...[, ...])	Configure how to do bond formation during the run.
config_box_deformation(self, target_box)	Configure where to (incrementally) deform the box to during the next simulation run.
config_lambda(self[, l])	Configure the modified velocity verlet integration parameter lambda.
config_num_steps_dpd(self[, num_steps])	Configure the number of steps to do in one DPD sequence.
config_num_steps_mc(self[, num_steps])	Configure the number of steps to do in one MC sequence.
config_restart_output(self, file[, output_every])	Set when to output a restart file.
config_shift_one_at_a_time(self[, ...])	Configure whether to shift atoms one at a time.
config_shift_possibility_empty(self[, ...])	Configure the possibility of shifting to empty positions.
config_sigma(self[, sigma])	Configure the force-field (pair-style) parameter sigma.
config_slipspring_high_cutoff(self[, cutoff])	Configure the lower cut-off of how far a pair may be distanced for a slip-spring to be created.
config_slipspring_low_cutoff(self[, cutoff])	Configure the higher cut-off of how far a pair may be distanced for a slip-spring to be created.
config_spring_constant(self[, k])	Configure the force-field (bond-style) parameter k, the spring constant.
config_step_output(self, values)	Set which values to log.
create_slip_springs(self, num[, bond_type])	Randomly add the specified number of slip-springs to neighbours within the specified cut-offs.
get_bond_lengths(self)
get_coordinates(self)
get_current_timestep(self)	Get the current timestep number.
get_nr_of_atoms(self)
get_nr_of_bonds(self)
get_nr_of_bonds_to_form(self)	Get the number of bonds that are configured to have to be formed.
get_nr_of_extra_atoms(self)
get_nr_of_extra_bonds(self)
get_nr_of_slip_springs(self)
get_nr_of_steps_dpd(self)
get_nr_of_steps_mc(self)
get_shift_one_at_a_time(self)
get_shift_possibility_empty(self)
get_slip_spring_bond_type(self)
get_spring_constant(self)
get_stress_tensor(self)
get_temperature(self)
get_timestep(self)	Get the timestep used in the simulation.
get_universe(self[, with_slipsprings])	Get a universe instance from the current coordinates (and connectivity).
get_volume(self)
read_restart_file(file)	Read a restart file in order to continue a simulation.
refresh_current_state(self)	After re-configuring the force-field parameters, this method should be called to update the current stress tensor etc.
run_simulation(self, n_steps[, dt, with_MC])
start_measuring_msd_for_atoms(self, atom_ids)	Set a new origin for measuing the mean square displacement for a specified set of atoms
validate_neighbour_list(self, arg0)	Validate the neighbor list consistency.
validate_state(self)	Validate the current simulation state.
write_restart_file(*args, **kwargs)	Overloaded function.

Methods Documentation

assume_box_large_enough(self: pylimer_tools_cpp.DPDSimulator) → None

Configure whether to run PBC on the bonds or not.

If your bonds could get larger than half the box length, this must be kept false (default). Otherwise, you can set it to true and therewith get some securities.

config_a(self: pylimer_tools_cpp.DPDSimulator, A: float = 25.0) → None

Configure the force-field (pair-style) parameter A.

config_allow_relocation_in_network(self: pylimer_tools_cpp.DPDSimulator, allow_relocation_in_network: bool = False) → None

Configure whether a relocation step may happen when a slip-spring has ended at a crosslink.

Side-effect: if true, the relocations may also happen to a slip-spring next to a crosslink.

Parameters:

(bool) (allow_relocation_in_network) – Whether to allow relocation in the network or not.

config_auto_correlator_output(self: pylimer_tools_cpp.DPDSimulator, values: list[pylimer_tools_cpp.OutputConfiguration], num_corr_in: int = 32, p: int = 16, m: int = 2) → None

Set which values to compute multiple-tau autocorrelation for. If you use this, you should cite Ramírez et al. [RamirezSVL10].

Parameters:

• values – a list of OutputConfiguration structs

• num_corr_in – Number of correlations in

• p – Parameter p for the autocorrelator

• m – Parameter m for the autocorrelator

config_average_output(self: pylimer_tools_cpp.DPDSimulator, values: list[pylimer_tools_cpp.OutputConfiguration]) → None

Set which values to compute averages for.

Parameters:

values – A list of OutputConfiguration structs specifying what to average

config_bond_formation(self: pylimer_tools_cpp.DPDSimulator, num_bonds_to_form: int, max_bonds_per_atom_type: dict[int, int], bond_formation_dist: float = 1.0, attempt_bond_formation_every: int = 50, atom_type_form_from: int = 2, atom_type_form_to: int = 1) → None

Configure how to do bond formation during the run.

Parameters:

• (int) (atom_type_form_to) – The nr of bonds to form in total. Use 0 to stop bond formation.

• (dict) (num_bonds_per_atom_type) – The nr of bonds each atom type may have at most (e.g., 2 for strand atoms, 4 for a tertiary crosslinkers)

• (float) (bond_formation_dist) – The maximum distance allowed to form bonds

• (int) – attempt to form bonds every this many steps during the simulation run

• (int) – The atom type to start forming bonds from.

• (int) – The atom type to start forming bonds to.

config_box_deformation(self: pylimer_tools_cpp.DPDSimulator, target_box: pylimer_tools_cpp.Box) → None

Configure where to (incrementally) deform the box to during the next simulation run.

config_lambda(self: pylimer_tools_cpp.DPDSimulator, l: float = 0.65) → None

Configure the modified velocity verlet integration parameter lambda.

config_num_steps_dpd(self: pylimer_tools_cpp.DPDSimulator, num_steps: int = 500) → None

Configure the number of steps to do in one DPD sequence.

Parameters:

num_steps – Number of DPD steps per sequence

config_num_steps_mc(self: pylimer_tools_cpp.DPDSimulator, num_steps: int = 500) → None

Configure the number of steps to do in one MC sequence.

Parameters:

num_steps – Number of Monte Carlo steps per sequence

config_restart_output(self: pylimer_tools_cpp.DPDSimulator, file: str, output_every: int = 50000) → None

Set when to output a restart file.

Note

The filename determines the type of serialization: .json, .xml are supported; other file endings will lead to binary serialization (fastest!).

Caution

This method may not be backwards- nor forward-compatible. Use the same version of pylimer-tools if you want to be sure that things work.

Parameters:

• file – The file path to the restart file to write

• output_every – How often to write the restart file (default: 50000)

config_shift_one_at_a_time(self: pylimer_tools_cpp.DPDSimulator, shift_one_at_a_time: bool = False) → None

Configure whether to shift atoms one at a time.

This setting affects Monte Carlo move behavior in the simulation.

config_shift_possibility_empty(self: pylimer_tools_cpp.DPDSimulator, shift_possibility_empty: bool = True) → None

Configure the possibility of shifting to empty positions.

This setting affects Monte Carlo moves in the simulation.

config_sigma(self: pylimer_tools_cpp.DPDSimulator, sigma: float = 3.0) → None

Configure the force-field (pair-style) parameter sigma.

config_slipspring_high_cutoff(self: pylimer_tools_cpp.DPDSimulator, cutoff: float = 2.0) → None

Configure the lower cut-off of how far a pair may be distanced for a slip-spring to be created.

config_slipspring_low_cutoff(self: pylimer_tools_cpp.DPDSimulator, cutoff: float = 0.5) → None

Configure the higher cut-off of how far a pair may be distanced for a slip-spring to be created.

config_spring_constant(self: pylimer_tools_cpp.DPDSimulator, k: float = 2.0) → None

Configure the force-field (bond-style) parameter k, the spring constant.

config_step_output(self: pylimer_tools_cpp.DPDSimulator, values: list[pylimer_tools_cpp.OutputConfiguration]) → None

Set which values to log.

Parameters:

values – a list of OutputConfiguration structs

create_slip_springs(self: pylimer_tools_cpp.DPDSimulator, num: int, bond_type: int = 9) → int

Randomly add the specified number of slip-springs to neighbours within the specified cut-offs.

get_bond_lengths(self: pylimer_tools_cpp.DPDSimulator) → numpy.ndarray[numpy.float64[m, 1]]

get_coordinates(self: pylimer_tools_cpp.DPDSimulator) → numpy.ndarray[numpy.float64[m, 1]]

get_current_timestep(self: pylimer_tools_cpp.DPDSimulator) → int

Get the current timestep number.

Returns:

The current timestep index

get_nr_of_atoms(self: pylimer_tools_cpp.DPDSimulator) → int

get_nr_of_bonds(self: pylimer_tools_cpp.DPDSimulator) → int

get_nr_of_bonds_to_form(self: pylimer_tools_cpp.DPDSimulator) → int

Get the number of bonds that are configured to have to be formed.

get_nr_of_extra_atoms(self: pylimer_tools_cpp.DPDSimulator) → int

get_nr_of_extra_bonds(self: pylimer_tools_cpp.DPDSimulator) → int

get_nr_of_slip_springs(self: pylimer_tools_cpp.DPDSimulator) → int

get_nr_of_steps_dpd(self: pylimer_tools_cpp.DPDSimulator) → int

get_nr_of_steps_mc(self: pylimer_tools_cpp.DPDSimulator) → int

get_shift_one_at_a_time(self: pylimer_tools_cpp.DPDSimulator) → bool

get_shift_possibility_empty(self: pylimer_tools_cpp.DPDSimulator) → bool

get_slip_spring_bond_type(self: pylimer_tools_cpp.DPDSimulator) → int

get_spring_constant(self: pylimer_tools_cpp.DPDSimulator) → float

get_stress_tensor(self: pylimer_tools_cpp.DPDSimulator) → numpy.ndarray[numpy.float64[3, 3]]

get_temperature(self: pylimer_tools_cpp.DPDSimulator) → float

get_timestep(self: pylimer_tools_cpp.DPDSimulator) → float

Get the timestep used in the simulation.

Returns:

The simulation timestep value

get_universe(self: pylimer_tools_cpp.DPDSimulator, with_slipsprings: bool = True) → pylimer_tools_cpp.Universe

Get a universe instance from the current coordinates (and connectivity).

Parameters:

with_slipsprings – Whether to include slip-springs in the returned universe (default: True)

get_volume(self: pylimer_tools_cpp.DPDSimulator) → float

static read_restart_file(file: str) → pylimer_tools_cpp.DPDSimulator

Read a restart file in order to continue a simulation.

Parameters:

file – The file path to the restart file to read

refresh_current_state(self: pylimer_tools_cpp.DPDSimulator) → None

After re-configuring the force-field parameters, this method should be called to update the current stress tensor etc.

run_simulation(self: pylimer_tools_cpp.DPDSimulator, n_steps: int, dt: float = 0.06, with_MC: bool = False) → None

start_measuring_msd_for_atoms(self: pylimer_tools_cpp.DPDSimulator, atom_ids: list[int]) → None

Set a new origin for measuing the mean square displacement for a specified set of atoms

validate_neighbour_list(self: pylimer_tools_cpp.DPDSimulator, arg0: float) → None

Validate the neighbor list consistency.

Returns:

True if the neighbor list is valid

validate_state(self: pylimer_tools_cpp.DPDSimulator) → None

Validate the current simulation state.

Returns:

True if the simulation state is valid

write_restart_file(*args, **kwargs)

Overloaded function.

1. write_restart_file(self: pylimer_tools_cpp.DPDSimulator, file: str) -> None

   Explicitly force the writing of a restart file, now!

   param file:

   The file path and name of the restart file to be written. Can end in .xml, .json or anything else (-> binary)

2. write_restart_file(self: pylimer_tools_cpp.DPDSimulator, file: str) -> None

   Explicitly force the writing of a restart file, now!

   param file:

   The file path and name of the restart file to be written. Can end in .xml, .json or anything else (-> binary)