Comments (7)
Ewald force should include k-space (reciprocal space) contribution
http://www.pages.drexel.edu/~cfa22/msim/node49.html
from mcmd.
name md_NVE_example
mode md ! mc or md
checkpoints_option off
input_atoms input.pdb ! PDB ONLY
potential_form lj
com_option off ! center of mass
carbasis 25 25 25 90 90 90
ensemble nve
! MOLECULAR DYNAMICS ONLY ===============================
md_pbc on
md_dt 0.2 ! fs
md_ft 2 ns
md_corrtime 1000 ! steps
!md_init_vel 0 ! A / fs^2
md_mode molecular
md_rotations off
from mcmd.
ATOM 1 HT h2o M 1 -1.752770 7.859250 -4.727420 default 0.4170
ATOM 2 OT h2o M 1 -1.871140 8.575620 -4.113410 default -0.8340
ATOM 3 HT h2o M 1 -1.746980 9.388710 -4.590600 default 0.4170
ATOM 4 HT h2o M 2 -2.090070 6.957620 -4.921120 default 0.4170
ATOM 5 OT h2o M 2 -2.854070 7.014120 -5.484420 default -0.8340
ATOM 6 HT h2o M 2 -3.625470 6.980620 -4.929320 default 0.4170
ATOM 7 HT h2o M 3 -2.665930 6.111700 -3.476070 default 0.4170
ATOM 8 OT h2o M 3 -1.816200 5.786300 -3.199900 default -0.8340
ATOM 9 HT h2o M 3 -1.744820 4.883110 -3.488810 default 0.4170
ATOM 10 HT h2o M 4 -0.068850 6.740530 -2.796270 default 0.4170
ATOM 11 OT h2o M 4 -0.817170 7.262010 -2.527460 default -0.8340
ATOM 12 HT h2o M 4 -1.348260 7.424680 -3.299340 default 0.4170
from mcmd.
quaternion dynamics tests
lj, no rotations, conserves in nve.
from mcmd.
LJES issue fixed via ae0614c
Actually, 0800139
k-space force needed to be doubled to compensate for half Ewald sphere
Now, need to add LJES routine to GPU CUDA
...done; matches cpu force.
NOW rotations...
from mcmd.
theta should maybe define angular position w/respect to ORIGINAL position; not previous.
so theta is recursive by integration function but the actual rotation should be performed on original relative coordinates r-r_com
from mcmd.
or,
U_rotation = torque dot theta , add this to normal potential?
from mcmd.
Related Issues (20)
- add surface dynamics (turn on pbc on one side)
- add Si, Cd UFF4MOF params
- MSD/diffusion calc HOT 2
- add auto-parameters for optimization/flexible sim by the given structure HOT 1
- do ewald ES energy for periodic optimization, too
- OpenMP for bonded energy/forces HOT 1
- make element determination function more selective
- box dynamics
- thermostat for flexible/atomic/molecular MD modes HOT 1
- long range polarization force
- make output for user-selected pair energies
- annealing for MD HOT 1
- make option for individual function datafiles
- make VACF calculator HOT 1
- lammps write bond/angle information HOT 1
- add impropers to default bond parameters detection HOT 2
- compilation error with gpu HOT 2
- Error while running "bash compile.sh" in cygwin HOT 2
- Feature to control max bondlength for specific atom names
- add common cations that bind to MOFs
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from mcmd.