Calculation of normal contact forces between silica nanospheres

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.

The effect of collisions on the chemomechanics of ice-covered silica slabs: a molecular dynamics study

Using molecular dynamics simulation and the REAX potential, we study the collision of two planar silica surfaces covered by water ice. Without the ice cover, the two surfaces stick at all velocities investigated (160-1800 m s-1), due to the formation of chemical bonds between the colliding surfaces. A narrow ice cover - here of thickness 2 nm - prevents the sticking above a characteristic velocity, the bouncing velocity νb. During the collision, reactions occur at the silica-water interface; in particular, water molecules are dissociated and silanols are formed at the surface of the silica slabs. Passivation of the silica surface by H atoms is of little consequence to the magnitude of vb but reduces the number of surface reactions occurring.

Shock-Induced Microstructural Evolution, Phase Transformation, Sintering of Al-Ni Dissimilar Nanoparticles: A Molecular Dynamics Study

Molecular dynamic simulations have been performed to explore contact behavior, microstructure evolution and sintering mechanism of Al-Ni dissimilar nanoparticles under high-velocity impact. We confirmed that the simulated contact stress, contact radius, and contact force under low-velocity impact are in good agreement with the predicted results of the Hertz model. However, with increasing the impact velocity, the simulated results gradually deviate from the predicted results of the Hertz model due to the elastic-plastic transition and atomic discrete structure. The normalized contact radius versus strain exhibits a weak dependence on nanosphere diameter. Below a critical velocity, there are very few HCP atoms in the nanospheres after thermal equilibrium. There are two different sintering mechanisms: under low-velocity impact, the sintering process relies mainly on the dislocation slip of Al nanospheres, while the dislocation slip of Ni nanospheres and the atomic diffusion of Al nanospheres predominate under high-velocity impact.

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.

Dynamic penetration behaviors of single/multi-layer graphene using nanoprojectile under hypervelocity impact

Single/multilayer graphene holds great promise in withstanding impact/penetration as ideal protective material. In this work, dynamic penetration behaviors of graphene has been explored using molecular dynamics simulations. The crashworthiness performance of graphene is contingent upon the number of layers and impact velocity. The variables including residual velocity and kinetic energy loss under different layers or different impact velocities have been monitored during the hypervelocity impact. Results show that there exists deviation from the continuum Recht–Ipson and Rosenberg–Dekel models, but these models tend to hold to reasonably predict the ballistic limit velocity of graphene with increasing layers. Besides, fractal theory has been introduced here and proven valid to quantitatively describe the fracture morphology. Furthermore, Forrestal–Warren rigid body model II still can well estimate the depth of penetration of multilayer graphene under a certain range of velocity impact. Finally, one modified model has been proposed to correlate the specific penetration energy with the number of layer and impact velocity.

Numerical Simulation of the Rheological Behavior of Nanoparticulate Suspensions

Nanoparticles significantly alter the rheological properties of a polymer or monomeric resin with major effect on the further processing of the materials. In this matter, especially the influence of particle material and disperse properties on the viscosity is not yet understood fully, but can only be modelled to some extent empirically after extensive experimental effort. In this paper, a numerical study on an uncured monomeric epoxy resin, which is filled with boehmite nanoparticles, is presented to elucidate the working principles, which govern the rheological behavior of nanoparticulate suspensions and to simulate the suspension viscosity based on assessable material and system properties. To account for the effect of particle surface forces and hydrodynamic interactions on the rheological behavior, a resolved CFD is coupled with DEM. It can be shown that the particle interactions caused by surface forces induce velocity differences between the particles and their surrounding fluid, which result in increased drag forces and cause the additional energy dissipation during shearing. The paper points out the limits of the used simulation method and presents a correction technique with respect to the Péclet number, which broadens the range of applicability. Valuable information is gained for a future mechanistic modelling of nanoparticulate suspension viscosity by elucidating the interdependency between surface forces, shear rate and resulting drag forces on the particles.

A New Interaction Force Model of Gold Nanorods Derived by Molecular Dynamics Simulation

Interactions between nanoparticles is one of the key factors governing their assembly for ordered structures. Understanding such interactions between non-spherical nanoparticles and developing a quantitative force model are critical to achieving the ordered structures for various applications. In the present study, the non-contact interactions of two identical gold nanorods (AuNRs) with different aspect ratios have been studied by molecular dynamics simulation. A new interaction potential and force model for two nanorods approaching side-by-side has been proposed as a function of particle surface separation and their relative orientation. In addition, the interaction potentials of two nanorods approaching in other typical orientation configurations (i.e., crossed, head-to-head and head-to-side) have also been investigated.

Effect of nanoparticle size on the mechanical properties of nanoparticle assemblies

Nanoparticle assemblies (NPAs) have attracted tremendous interests of various research communities. The particle-size-effect on mechanical properties of NPAs is systematically studied. With decreasing the particle size d from 300 nm to 10 nm, the SiO2 NPAs become drastically harder (∼39×), stiffer (∼15×), and tougher (>3.5×). The results are consistent with the data scattered in the literature for various nanoparticle (NP) systems, indicating a fundamentally universal d-effect for all NPAs. A model is developed to correlate the hardness and the NP junction (NPJ) strength f. Here, f is mainly due to van der Waals and capillary interactions, roughly a constant (140 nN) for d = 100-300 nm, and then f decreases with decreasing d from ∼100 nm. The deformation mechanism of NPAs (for indentation depth ≫d) is shear plasticity involving shear breaking of NPJs. The fundamental mechanism for the d-effect is that, with decreasing d, the NPJ's planar density increases much faster than the decrease of f. Moreover, three deformation mechanisms of NPAs, (1) nanoparticle dislodging, (2) shear-band formation, and (3) cracking are naturally d-dependent. These new findings can provide important insights into the fundamental understanding of the inter-NP interaction, the mechanical behavior of the NPAs, and the design of robust NP-based devices.

The bouncing threshold in silica nanograin collisions

Using molecular dynamics simulations, we study collisions between amorphous silica nanoparticles. Our silica model contains uncontaminated surfaces, that is, the effect of surface hydroxylation or of adsorbed water layers is excluded. For central collisions, we characterize the boundary between sticking and bouncing collisions as a function of impact velocity and particle size and quantify the coefficient of restitution. We show that the traditional Johnson-Kendall-Roberts (JKR) model provides a valid description of the ingoing trajectory of two grains up to the moment of maximum compression. The distance of closest approach is slightly underestimated by the JKR model, due to the appearance of plasticity in the grains, which shows up in the form of localized shear transformation zones. The JKR model strongly underestimates the contact radius and the collision duration during the outgoing trajectory, evidencing that the breaking of covalent bonds during grain separation is not well described by this model. The adhesive neck formed between the two grains finally collapses while creating narrow filaments joining the grains, which eventually tear.

Enhanced charge transport and photovoltaic performance induced by incorporating rare-earth phosphor into organic-inorganic hybrid solar cells

In this work, dysprosium ion decorated yttrium oxide (Dy(3+):Y2O3) nanocrystal phosphors were incorporated into TiO2 acceptor thin film in a bid to enhance the light harvest, charge separation and transfer in the hybrid solar cells. The results show that the energy level offset between the donor (P3HT) and the acceptor (Dy(3+):Y2O3-TiO2) has been narrowed down, thus leading to the enhanced electron and hole transports, and also photovoltaic performances as compared to pure TiO2 without incorporating Dy(3+):Y2O3. By applying femtosecond transient optical spectroscopy, after the incorporation of dopant Dy(3+):Y2O3 into TiO2 at 6 wt%, both the hot electron and hole transfer lifetimes have been shortened, that is, from 30.2 ps and 6.94 ns to 25.1 ps and 1.26 ns, respectively, and an enhanced efficiency approaching 3% was achieved as compared to 2.0% without doping, indicating that the energetic charges are captured more efficiently benefitting a higher power conversion efficiency. Moreover, these results reveal that both the conduction band (CB) and valence band (VB) edges of the acceptor were elevated by 0.57 and 0.32 eV, respectively, after incorporating 6 wt% Dy(3+):Y2O3. This work demonstrates that distinct energy level alignment engineered by Dy(3+):Y2O3 phosphor has an important role in pursuing efficient future solar cells and underscores the promising potential of rare-earth phosphor in solar applications.

Energy level control: toward an efficient hot electron transport

Highly efficient hot electron transport represents one of the most important properties required for applications in photovoltaic devices. Whereas the fabrication of efficient hot electron capture and lost-cost devices remains a technological challenge, regulating the energy level of acceptor-donor system through the incorporation of foreign ions using the solution-processed technique is one of the most promising strategies to overcome this obstacle. Here we present a versatile acceptor-donor system by incorporating MoO₃:Eu nanophosphors, which reduces both the ‘excess' energy offset between the conduction band of acceptor and the lowest unoccupied molecular orbital of donor, and that between the valence band and highest occupied molecular orbital. Strikingly, the hot electron transfer time has been shortened. This work demonstrates that suitable energy level alignment can be tuned to gain the higher hot electron/hole transport efficiency in a simple approach without the need for complicated architectures. This work builds up the foundation of engineering building blocks for third-generation solar cells.

The dynamic effect on mechanical contacts between nanoparticles

The rich behaviors of high-speed mechanical contacts at the nanoscale have been studied. The seldom observed elastic-plastic transition governed by Hertz and Thornton models has been clearly unveiled, the origins of the hardening effect and the deformation mechanism of nanoscale plasticity have been discussed in terms of structural changes after compression and a series of physical quantities are measured including contact forces, contact radius, contact stress, coefficient of restitution and total impact time. Our simulation results closely resemble experiments and/or theoretical predictions: (i) when impact speed v is higher than Y/ρc0, the elastic-plastic deformation transition occurs, (ii) the yielded apparent elastic modulus and hardness are larger than those of the bulk, (iii) the initiating yield stress Y and hardness P0 still satisfy P0 ≈ 1.6Y, (iv) particle's volume decreases during compression, (v) contact radius a follows a [proportionality] v(2/5), (vi) at v ≥ 2000 m s(-1), the coefficient of restitution follows e [proportionality] v(-1/4) and (vii) the total time of impact follows Tc [proportionality] v(-1/5). However, there also exist many quantitative differences. The contact radius and final contact radius are underestimated by the continuum predictions while the total impact time is overestimated, but all of them reasonably agree with theoretical predictions with an increase of contact area and impact speed. The theoretical equation is adapted to predict the final contact radius during normal impact, in which the contact radius at zero load is also formulated.