[DFTB-Plus-User] Increasing the distance between two atoms

Timothy Kennedy timothyrkennedy27 at gmail.com
Sat Mar 30 03:33:35 CET 2024


Hi all,

I'm currently trying to replicate the results of a study in which they
simulated mechanical tensile force by moving two atoms apart in a molecule.
I'm not entirely sure what the best way to achieve this is using DFTB+. I'm
currently using this function in an attempt to increase the size of the
simulation box (the code was originally from my professor):

# Setup NPT dynamics
myNPT = NPT(
    atoms=atoms,
    timestep=0.1 * units.fs,
    externalstress=0,
    ttime=100 * 0.1 * units.fs,
    pfactor=1,
    temperature_K=300,
    mask=(0,0,0),
    logfile='mylog',
    trajectory='mytraj',
    append_trajectory=False,
    loginterval=1
)

# Define the strain rate
strain_rate = np.array([[.0001 / units.fs, 0, 0], [0, 0, 0], [0, 0, 0]])  #
strain x for Strauss unit cell
myNPT.set_strain_rate(strain_rate)

I'm not getting very useful results and I'm unsure what the problem is. I
can share my entire program if anybody would like to see it. Any help is
appreciated!
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.zfn.uni-bremen.de/pipermail/dftb-plus-user/attachments/20240329/ae2d8b15/attachment.htm>


More information about the DFTB-Plus-User mailing list