Numerical study of the low-frequency atomic dynamics in a Lennard-Jones glass

Aging effects on thermal conductivity of glass-forming liquids

Thermal conductivity of a model glass-forming system in the liquid and glass states is studied using extensive numerical simulations. We show that near the glass transition temperature, where the structural relaxation time becomes very long, the measured thermal conductivity decreases with increasing age. Second, the thermal conductivity of the disordered solid obtained at low temperatures is found to depend on the cooling rate with which it was prepared. For the cooling rates accessible in simulations, lower cooling rates lead to lower thermal conductivity. Our analysis links this decrease of the thermal conductivity with increased exploration of lower-energy inherent structures of the underlying potential energy landscape. Further, we show that the lowering of conductivity for lower-energy inherent structures is related to the high-frequency harmonic modes associated with the inherent structure being less extended. Possible effects of considering relatively small systems and fast cooling rates in the simulations are discussed.

Boson peak and Ioffe-Regel criterion in amorphous siliconlike materials: The effect of bond directionality

The vibrational properties of model amorphous materials are studied by combining complete analysis of the vibration modes, dynamical structure factor, and energy diffusivity with exact diagonalization of the dynamical matrix and the kernel polynomial method, which allows a study of very large system sizes. Different materials are studied that differ only by the bending rigidity of the interactions in a Stillinger-Weber modelization used to describe amorphous silicon. The local bending rigidity can thus be used as a control parameter, to tune the sound velocity together with local bonds directionality. It is shown that for all the systems studied, the upper limit of the Boson peak corresponds to the Ioffe-Regel criterion for transverse waves, as well as to a minimum of the diffusivity. The Boson peak is followed by a diffusivity's increase supported by longitudinal phonons. The Ioffe-Regel criterion for transverse waves corresponds to a common characteristic mean-free path of 5-7 Å (which is slightly bigger for longitudinal phonons), while the fine structure of the vibrational density of states is shown to be sensitive to the local bending rigidity.

Anharmonicity, mechanical instability, and thermodynamic properties of the Cr-Re σ-phase

Using density-functional theory in combination with the direct force method and molecular dynamics we investigate the vibrational properties of a binary Cr-Re σ-phase. In the harmonic approximation, we have computed phonon dispersion curves and density of states, evidencing structural and chemical effects. We found that the σ-phase is mechanically unstable in some configurations, for example, when all crystallographic sites are occupied by Re atoms. By using a molecular-dynamics-based method, we have analysed the anharmonicity in the system and found negligible effects (~0.5 kJ/mol) on the Helmholtz energy of the binary Cr-Re σ-phase up to 2000 K (~0.8T(m)). Finally, we show that the vibrational contribution has significant consequences on the disordering of the σ-phase at high temperature.

Vibrational spectrum of topologically disordered systems

The topological nature of the disorder of glasses and supercooled liquids strongly affects their high-frequency dynamics. In order to understand its main features, we analytically studied a simple topologically disordered model, where the particles oscillate around randomly distributed centers, interacting through a generic pair potential. We present results of a resummation of the perturbative expansion in the inverse particle density for the dynamic structure factor and density of states. This gives accurate results for the range of densities found in real systems.

Experimental evidence of the acousticlike character of the high frequency excitations in glasses

The dynamic structure factor S(Q,E) of glassy glycerol has been measured by inelastic x-ray scattering as a function of momentum transfer Q and at constant energy transfer E. This allows one to establish independently from specific models of S(Q,E) the following: (i) Propagating collective excitations exist in glasses at high Q. (ii) Their dispersion up to E higher than E(BP) (the boson peak energy) confirms that E(BP) is not the onset of modes localization. (iii) The observation of an almost Q independent plateau on the high Q side of the Brillouin peak supports numerical simulations on glasses, describing the vibrational eigenvectors in terms of acousticlike and "random" components.

Dynamical structure factor in disordered systems

We study the spectral width as a function of the external momentum for the dynamical structure factor of a disordered harmonic solid, considered as a toy model for supercooled liquids and glasses. In the contexts of both the single-link coherent potential approximation and a single-defect approximation, two different regimes are clearly identified: if the density of states at zero energy is zero, the usual p(4) law is recovered for small momentum. On the contrary, if the disorder induces a nonvanishing density of states at zero energy, a linear behavior is obtained. The dynamical structure factor is numerically calculated in lattices as large as 96(3) and satisfactorily agrees with the analytical computations.

Relaxation processes in harmonic glasses?

A relaxation process, with the associated phenomenology of sound attenuation and sound velocity dispersion, is found in a simulated harmonic Lennard-Jones glass. We propose to identify this process with the so-called microscopic (or, instantaneous) relaxation process observed in real glasses and supercooled liquids. A model based on the memory function approach accounts for the observation and allows one to relate to each other (1) the characteristic time and strength of this process, (2) the low frequency limit of the dynamic structure factor of the glass, and (3) the high frequency sound attenuation coefficient, with its observed quadratic dependence on the momentum transfer.

Determination of the short-wavelength propagation threshold in the collective excitations of liquid ammonia

The dynamics structure factor S(Q,E) of liquid ammonia l-NH3 at T = 200 K and at its vapor pressure has been measured by inelastic x-ray scattering (IXS) in the 1-15 nm(-1) momentum transfer ( Q) range. Contrary to previous IXS studies on other associated liquids and glasses, in l-NH3 a large inelastic signal is observed up to Q = 15 nm(-1). This, enabling S(Q,E) measurements as a function of Q at constant E transfer, allows us to demonstrate experimentally the transition from a propagating dynamics regime, where the acoustic excitation energy linearly disperses with Q, to a high-Q regime, where it is no longer possible to observe a dominant excitation in the S(Q,E).

Evolution of vibrational excitations in glassy systems

The equations of the mode-coupling theory (MCT) for ideal liquid-glass transitions are used for a discussion of the evolution of the density-fluctuation spectra of glass-forming systems for frequencies within the dynamical window between the band of high-frequency motion and the band of low-frequency-structural-relaxation processes. It is shown that the strong interaction between density fluctuations with microscopic wavelength and the arrested glass structure causes an anomalous-oscillation peak, which exhibits the properties of the so-called boson peak. It produces an elastic modulus which governs the hybridization of density fluctuations of mesoscopic wavelength with the boson-peak oscillations. This leads to the existence of high-frequency sound with properties as found by x-ray-scattering spectroscopy of glasses and glassy liquids. The results of the theory are demonstrated for a model of the hard-sphere system. It is also derived that certain schematic MCT models, whose spectra for the stiff-glass states can be expressed by elementary formulas, provide reasonable approximations for the solutions of the general MCT equations.