7 from __future__
import print_function
16 minimum_particle_radius = 10
18 maximum_spring_constant = 1
20 maximum_diffusion_coefficient =\
25 expected_rotational_delta =\
27 time_step) * minimum_particle_radius
29 expected_spring_diffusion_length =\
31 .5 * maximum_spring_constant *
32 4 * expected_delta ** 2, time_step)
34 print(
"with a time step of", time_step,
"fs")
35 print(
"an object of radius", minimum_particle_radius,
"A will move",
36 expected_delta,
"A and a point on its surface will move",
37 expected_rotational_delta,
"A more")
38 print(
"the motion from fluctuations in the spring compression will be",
39 expected_spring_diffusion_length,
"A")
40 print(
"and a compression of 10% of the radius will induce a motion of",
42 .5 * maximum_spring_constant *
43 (.1 * minimum_particle_radius) ** 2,
Strings setup_from_argv(const Strings &argv, std::string description, std::string positional_description, int num_positional)
double get_diffusion_length(double D, double force, double t, double temp=-1)
double get_diffusion_angle(double D, double dtfs)
Functionality for loading, creating, manipulating and scoring atomic structures.
double get_einstein_diffusion_coefficient(double r, double temp=-1)