Molecular dynamics of rolling and twisting motion of amorphous nanoparticles

Molecular dynamics simulations of head-on low-velocity collisions between particles

The particle contact model is important for powder simulations. Although several contact models have been proposed, their validity has not yet been well established. Therefore, we perform molecular dynamics (MD) simulations to clarify the particle interaction. We simulate head-on collisions of two particles with impact velocities less than a few percent of the sound velocity to investigate the dependence of the interparticle force and the coefficient of restitution on the impact velocity and particle radius. In this study, we treat particles with a radius of 10-100 nm and perform simulations with up to 0.2 billion atoms. We find that the interparticle force exhibits hysteresis between the loading and unloading phases. Larger impact velocities result in strong hysteresis and plastic deformation. For all impact velocities and particle radii, the coefficient of restitution is smaller than that given by the Johnson-Kendall-Robert theory, which is a contact model that gives the force between elastic spherical particles. A contact model of inelastic particles cannot reproduce our MD simulations. In particular, the coefficient of restitution is significantly reduced when the impact velocity exceeds a certain value. This significant energy dissipation cannot be explained even by the contact models including plastic deformation. We also find that the coefficient of restitution increases with increasing particle radius. We also find that the previous contact models including plastic deformation cannot explain the strong energy dissipation obtained in our MD simulations, although they agree with the MD results for very low impact velocities. Accordingly, we have constructed a new dissipative contact model in which the dissipative force increases with the stress generated by collisions. The new stress-dependent model successfully reproduces our MD results over a wider range of impact velocities than the conventional models do. In addition, we proposed another, simpler, dissipative contact model that can also reproduce the MD results.

An atomistic study of sticking, bouncing, and aggregate destruction in collisions of grains with small aggregates

Molecular dynamics simulations are used to study central collisions between spherical grains and between grains and small grain aggregates (up to 5 grains). For a model material (Lennard-Jones), grain-grain collisions are sticking when the relative velocity v is smaller than the so-called bouncing velocity and bouncing for higher velocities. We find a similar behavior for grain-aggregate collisions. The value of the bouncing velocity depends only negligibly on the aggregate size. However, it is by 35% larger than the separation velocity needed to break a contact; this is explained by energy dissipation processes during the collision. The separation velocity follows the predictions of the macroscopic Johnson-Kendall-Roberts theory of contacts. At even higher collision velocities, the aggregate is destroyed, first by the loss of a monomer grain and then by total disruption. In contrast to theoretical considerations, we do not find a proportionality of the collision energy needed for destruction and the number of bonds to be broken. Our study thus sheds novel light on the foundations of granular mechanics, namely the energy needed to separate two grains, the difference between grain-grain and grain-aggregate collisions, and the energy needed for aggregate destruction.

Collisions between CO, CO[Formula: see text], H[Formula: see text]O and Ar ice nanoparticles compared by molecular dynamics simulation

Molecular dynamics simulations are used to study collisions between amorphous ice nanoparticles consisting of CO, CO[Formula: see text], Ar and H[Formula: see text]O. The collisions are always sticking for the nanoparticle size (radius of 20 nm) considered. At higher collision velocities, the merged clusters show strong plastic deformation and material mixing in the collision zone. Collision-induced heating influences the collision outcome. Partial melting of the merged cluster in the collision zone contributes to energy dissipation and deformation. Considerable differences exist-even at comparable collision conditions-between the ices studied here. The number of ejecta emitted during the collision follows the trend in triple-point temperatures and increases exponentially with the NP temperature.

Bouncing and spinning of amorphous Lennard-Jones nanoparticles under oblique collisions

Collisions of Lennard-Jones nanoparticles (NPs) may be used to study the generic collision behavior of NPs. We study the collision dynamics of amorphous NPs for oblique collisions using molecular dynamics simulation as a function of collision velocity and impact parameter. In order to allow for NP bouncing, the attraction between atoms originating from differing NPs is reduced. For near-central collisions, a finite region of velocities – a ‘bouncing window’ – exists where the 2 NPs bounce from each other. At smaller velocities, energy dissipation and – at larger velocities – also NP deformation do not allow the NPs to surpass the attractive forces such that they stick to each other. Oblique collisions of non-rotating NPs convert angular momentum into NP spin. For low velocities, the NP spin is well described by assuming the NPs to come momentarily to a complete stop at the contact point (‘grip’), such that orbital and spin angular momentum share the pre-collision angular momentum in a ratio of 5:2. The normal coefficient of restitution increases with impact parameter for small velocities, but changes sign for larger velocities where the 2 NPs do not repel but their motion direction persists. The tangential coefficient of restitution is fixed in the ‘grip’ regime to a value of 5/7, but increases towards 1 for high-velocity collisions at not too small impact parameters, where the 2 NPs slide along each other.