Using molecular dynamics simulations, we investigate the effect of uniaxial elastostatic compression on the potential energy, structural relaxation, and mechanical properties of binary glasses. We consider the three-dimensional Kob-Andersen binary mixture, which was initially cooled from the liquid state to the glass phase with a slow rate at zero pressure. The glass was then loaded with a static stress at the annealing temperature during extended time intervals. It is found that the static stress below the yielding point induces large-scale plastic deformation and significant rejuvenation when the annealing temperature is smaller than a fraction of the glass transition temperature. By contrast, aging effects become dominant at sufficiently small values of the static stress and higher annealing temperatures. The mechanical tests after the elastostatic compression have shown that both the elastic modulus and the yield stress decrease in rejuvenated samples, while the opposite trend was observed for relaxed glasses. These results might be useful for the thermomechanical processing of metallic glasses with optimized mechanical properties.
We present the results of the study of changes in liquid properties during ultrafast cooling. The molecular dynamics (MD) method is used, with aluminum melt as an example. Based on the changes in shear stress autocorrelation functions (SACF) with temperature in an ensemble of MD trajectories, we develop three approaches to the study of melt changes. In the first one, we investigate the asymptotic behavior of SACFs and a sharp increase in the melt viscosity, which is a conventional criterion. In the second approach, we present direct evidence of the transition of a metastable melt to a non-equilibrium state. In the third one, we show the appearance of transverse oscillations in a film of the melt. The most important observation is that all three phenomena occur in the same temperature range. On the basis of the current and the previous work, we conclude that there is a gap between the temperature of the splitting of the second peak of the pair correlation function and the temperature of the transition to a solid-like amorphous state. The splitting of the second peak occurs at a significantly higher temperature and this phenomenon is discussed.
The interaction and tunneling conductance between oppositely located ends of coaxial carbon nanotubes are studied by the example of two (11, 11) nanotubes with open ends terminated by hydrogen atoms. The Green function formalism is applied to determine the tunneling current through the nanotube ends as a function of the distance between the ends, relative orientation of the nanotubes and voltage applied. The energy favorable configuration of the coaxial nanotubes is obtained by the analysis of their interaction energy at different distances between the nanotube ends and angles of their relative rotation. Using these calculations, a general scheme of the force sensor based on the interaction between ends of coaxial nanotubes is proposed and the relation between the tunneling conductance and measured force is established for the considered nanotubes. The operational characteristics of this device as a magnetic field sensor based on measurements of the magnetic force acting on the coaxial nanotubes filled with magnetic endofullerenes are estimated.
The effect of cyclic loading on relaxation dynamics and mechanical properties of metallic glasses is studied using molecular dynamics simulations. We consider the Kob-Andersen three-dimensional binary mixture rapidly cooled across the glass transition and subjected to thousands of tension-compression cycles in the elastic range. It was found that during cyclic loading at constant pressure, the system is relocated to progressively lower levels of the potential energy, thus promoting greater densification and higher strength. Furthermore, with increasing stress amplitude, the average glass density increases and the minimum of the potential energy becomes deeper, while the elastic modulus is reduced. The typical size of clusters of atoms with large nonaffine displacements becomes smaller over consecutive cycles, which correlates with the gradual decrease in the potential energy. These results are important for thermomechanical processing of metallic glasses with improved mechanical properties.
The evolution of porous structure and mechanical properties of binary glasses under tensile loading were examined using molecular dynamics simulations. We consider vitreous systems obtained in the process of phase separation taking place after a rapid isochoric quench of a glass-forming liquid to temperatures below the glass transition point. The porous structure in undeformed samples varies from a connected porous network to a random distribution of isolated pores with increasing average density. We find that the elastic modulus follows a power-law dependence on the average glass density and the pore size distribution at small strain remains nearly the same as in quiescent samples. Upon further loading, the pores become significantly deformed and coalesce into larger voids that leads to formation of system-spanning empty regions associated with failure of the material.
The time evolution of the pore size distributions and mechanical properties of amorphous solids at constant pressure is studied using molecular dynamics simulations. The porous glasses were initially prepared at constant volume conditions via a rapid thermal quench from the liquid state to the glassy region and allowing for simultaneous phase separation and material solidification. We found that at constant pressure and low temperature, the porous network becomes more compact and the glassy systems relocate to progressively lower levels of the potential energy. Although the elastic modulus and the average glass density both increase with the waiting time, their dependence is described by the power-law function with the same exponent. Moreover, the results of numerical simulations demonstrated that under tensile loading at constant pressure, low-density porous samples become significantly deformed and break up into separate domains at high strain, while dense glasses form a nearly homogeneous solid material.
The influence of thermal processing on the potential energy, atomic structure, and mechanical properties of metallic glasses is examined using molecular dynamics simulations. We study the three-dimensional binary mixture, which was first relaxed near the glass transition temperature, and then rapidly cooled deep into the glass phase. It was found that glasses prepared at higher annealing temperatures are relocated to higher energy states and their average glass structure remains more disordered, as reflected in the shape of the pair correlation function. The results of mechanical testing demonstrate that both the shear modulus and yielding peak increase significantly when the annealing temperature approaches $T_g$ from above. Moreover, the shear modulus becomes a strong function of strain rate only for samples equilibrated at sufficiently high temperatures. Based on the spatial distribution of nonaffine displacements, we show that the deformation mode changes from brittle to ductile upon increasing annealing temperature. These results can be useful for the design and optimization of the fabrication processes of bulk glassy alloys with improved plasticity.
We study the effect of periodic, spatially uniform temperature variation on mechanical properties and structural relaxation of amorphous alloys using molecular dynamics simulations. The disordered material is modeled via a non-additive binary mixture, which is annealed from the liquid to the glassy state with various cooling rates and then either aged at constant temperature or subjected to thermal treatment. We found that in comparison to aged samples, thermal cycling with respect to a reference temperature of approximately half the glass transition temperature leads to more relaxed states with lower levels of potential energy. The largest energy decrease was observed for rapidly quenched glasses cycled with the thermal amplitude slightly smaller than the reference temperature. Following the thermal treatment, the mechanical properties were probed via uniaxial tensile strain at the reference temperature and constant pressure. The numerical results indicate an inverse correlation between the levels of potential energy and values of the elastic modulus and yield stress as a function of the thermal amplitude.
Non-equilibrium molecular dynamics simulations are performed to investigate the dynamic behavior of three-dimensional binary glasses prepared via an instantaneous quench across the glass transition. We found that with increasing strain amplitude up to a critical value, the potential energy approaches lower minima in steady state, whereas the amplitude of shear stress oscillations becomes larger. Below the yielding transition, the storage modulus dominates the mechanical response, and the gradual decay of the potential energy over consecutive cycles is accompanied by reduction in size of transient clusters of atoms with large nonaffine displacements. In contrast, above the yield strain, the loss modulus increases and the system settles at a higher level of potential energy due to formation of a system-spanning shear band after a number of transient cycles.
The telescopic contact between graphene layers with a dielectric spacer is considered as a new type of graphene-based nanoelectronic devices. The tunneling current through the contacts with and without an argon spacer is calculated as a function of the overlap length, stacking of the graphene layers and voltage applied using non-equilibrium Green function formalism. A negative differential resistance (similar to semiconductor tunnel diode) is found with the peak to valley ratio up to 10 and up to 2 for the contacts without any spacer and with the argon spacer, respectively. The capacitance of the contacts between the graphene layers with the argon spacer is calculated as a function of temperature taking into account the quantum contribution. The related RC time constant is estimated to be about 3 ps, which allows elaboration of fast-response nanoelectronic devices. The possibility of application of the contacts as memory cells is discussed.
The influence of the external pressure and surface energy on the wetting transition at nanotextured interfaces is studied using molecular dynamics and continuum simulations. The surface roughness of the composite interface is introduced via an array of spherical nanoparticles with controlled wettability. We find that in the absence of an external pressure, the liquid interface is flat and its location relative to the solid substrate is determined by the particle size and the local contact angle. With increasing pressure on the liquid film, the interface becomes more curved and the three-phase contact line is displaced along the spherical surface but remains stable due to re-entrant geometry. It is demonstrated that the results of molecular dynamics simulations for the critical pressure of the Cassie-Baxter wetting state agree well with the estimate of the critical pressure obtained by numerical minimization of the interfacial energy.