help出力 md_packets_FuncAnimation.py --help
usage: md_packets_FuncAnimation.py [-h] [--n N] [--a A] [--m0 M0] [--k2 K2]
[--k4 K4] [--method {euler,verlet}]
[--dt DT] [--nstep NSTEP]
[--interval INTERVAL]
[--save_interval SAVE_INTERVAL]
[--boundary {fixed}] [--A A]
[--sigma SIGMA] [--k0 K0]
[--initial_mode {left_only,two_packets}]
[--x0_left X0_LEFT] [--x0_right X0_RIGHT]
[--scatterer_mode {none,heavy_mass,light_mass,soft_bond,hard_bond,nonlinear_bond,barrier}]
[--scatterer_center SCATTERER_CENTER]
[--mass_width MASS_WIDTH]
[--heavy_mass_factor HEAVY_MASS_FACTOR]
[--light_mass_factor LIGHT_MASS_FACTOR]
[--bond_width BOND_WIDTH]
[--soft_bond_k2_factor SOFT_BOND_K2_FACTOR]
[--soft_bond_k4_factor SOFT_BOND_K4_FACTOR]
[--hard_bond_k2_factor HARD_BOND_K2_FACTOR]
[--hard_bond_k4_factor HARD_BOND_K4_FACTOR]
[--nonlinear_bond_k2_factor NONLINEAR_BOND_K2_FACTOR]
[--nonlinear_bond_k4_factor NONLINEAR_BOND_K4_FACTOR]
[--barrier_half_width BARRIER_HALF_WIDTH]
[--barrier_mass BARRIER_MASS]
[--barrier_k2 BARRIER_K2]
[--barrier_k4 BARRIER_K4] [--save {0,1}]
[--gif_fps GIF_FPS] [--gif_dpi GIF_DPI]
1D coupled oscillator wave-packet scattering simulation with linear + quartic
springs.
options:
-h, --help show this help message and exit
--n N number of masses (default: 200)
--a A lattice spacing (default: 1.0)
--m0 M0 background mass (default: 1.0)
--k2 K2 background linear spring constant (default: 5.0)
--k4 K4 background nonlinear spring constant (default: 500.0)
--method {euler,verlet}
time-integration method (default: verlet)
--dt DT time step (default: 0.05)
--nstep NSTEP number of MD/animation steps (default: 1000)
--interval INTERVAL animation interval [ms] (default: 1)
--save_interval SAVE_INTERVAL
interval for FFT and energy snapshots (default: 10)
--boundary {fixed} boundary condition (default: fixed)
--A A packet amplitude (default: 0.2)
--sigma SIGMA packet width (default: 2.0)
--k0 K0 central wave number (default: 0.1)
--initial_mode {left_only,two_packets}
initial packet configuration (default: left_only)
--x0_left X0_LEFT left packet center. If omitted, n*a*0.25 is used.
(default: None)
--x0_right X0_RIGHT right packet center. If omitted, n*a*0.75 is used.
(default: None)
--scatterer_mode {none,heavy_mass,light_mass,soft_bond,hard_bond,nonlinear_bond,barrier}
scatterer type (default: heavy_mass)
--scatterer_center SCATTERER_CENTER
scatterer center index. If omitted, n//2 is used.
(default: None)
--mass_width MASS_WIDTH
mass-defect width (default: 3)
--heavy_mass_factor HEAVY_MASS_FACTOR
heavy mass factor (default: 10.0)
--light_mass_factor LIGHT_MASS_FACTOR
light mass factor (default: 0.2)
--bond_width BOND_WIDTH
bond-defect width (default: 4)
--soft_bond_k2_factor SOFT_BOND_K2_FACTOR
soft-bond k2 factor (default: 0.2)
--soft_bond_k4_factor SOFT_BOND_K4_FACTOR
soft-bond k4 factor (default: 1.0)
--hard_bond_k2_factor HARD_BOND_K2_FACTOR
hard-bond k2 factor (default: 5.0)
--hard_bond_k4_factor HARD_BOND_K4_FACTOR
hard-bond k4 factor (default: 1.0)
--nonlinear_bond_k2_factor NONLINEAR_BOND_K2_FACTOR
nonlinear-bond k2 factor (default: 1.0)
--nonlinear_bond_k4_factor NONLINEAR_BOND_K4_FACTOR
nonlinear-bond k4 factor (default: 10.0)
--barrier_half_width BARRIER_HALF_WIDTH
barrier half width (default: 3)
--barrier_mass BARRIER_MASS
barrier mass value (default: 8.0)
--barrier_k2 BARRIER_K2
barrier k2 value (default: 20.0)
--barrier_k4 BARRIER_K4
barrier k4 value. If omitted, background k4 is used.
(default: None)
--save {0,1} save GIF with ani.save after the animation window is
closed (default: 0)
--gif_fps GIF_FPS frames per second for GIF save (default: 30)
--gif_dpi GIF_DPI DPI for GIF save (default: 100)
