We use ab initio molecular dynamics simulations to investigate the properties of the dry surface of pure silicaand sodium silicate glasses. The surface layers are defined based on the atomic distributions along the direction(z direction) perpendicular to the surfaces. We show that these surfaces have a higher concentration of danglingbonds as well as two-membered (2M) rings than the bulk samples. Increasing the concentration of Na$_2$O reducesthe proportion of structural defects. From the vibrational density of states, one concludes that 2M rings have aunique vibrational signature at a frequency ≈850 cm$^{−1}$ , compatible with experimental findings. We also findthat, due to the presence of surfaces, the atomic vibration in the z direction is softer than for the two otherdirections. The electronic density of states shows clearly the differences between the surface and interior and wecan attribute these to specific structural units. Finally, the analysis of the electron localization function allows toget insight on the influence of local structure and the presence of Na on the nature of chemical bonding in the glasses.

We propose a novel model for a glass-forming liquid which allows to switch in a continuous manner froma standard three-dimensional liquid to a fully connected mean-field model. This is achieved by introducingk additional particle-particle interactions which thus augments the effective number of neighbors of eachparticle. Our computer simulations of this system show that the structure of the liquid does not changewith the introduction of these pseudo neighbours and by means of analytical calculations, we determine thestructural properties related to these additional neighbors. We show that the relaxation dynamics of thesystem slows down very quickly with increasing k and that the onset and the mode-coupling temperaturesincrease. The systems with high values of k follow the MCT power law behaviour for a larger temperaturerange compared to the ones with lower values of k. The dynamic susceptibility indicates that the dynamicheterogeneity decreases with increasing k whereas the non-Gaussian parameter is independent of it. Thus weconclude that with the increase in the number of pseudo neighbours the system becomes more mean-field like.By comparing our results with previous studies on mean-field like system we come to the conclusion that thedetails of how the mean-field limit is approached are important since they can lead to different dynamicalbehavior in this limit.

Using atomistic computer simulations we determine the roughness and topographical features of melt-formed (MS) and fracture surfaces (FS) of oxide glasses. We find that the topography of the MS is described well by the frozen capillary wave theory. The FS are significant rougher than the MS and depend strongly on glass composition. The height-height correlation function for the FS shows an unexpected logarithmic dependence on distance, in contrast to the power law found in experiments. We unravel the crucial role of spatial resolution on surface measurements and conclude that on length scales less than 10 nm FS are not self-affine fractals.

Using cyclic shear to drive a two dimensional granular system, we determine the structural char-acteristics for different inter-particle friction coefficients. These characteristics are the result of acompetition between mechanical stability and entropy, with the latter’s effect increasing with fric-tion. We show that a parameter-free maximum-entropy argument alone predicts an exponential cell order distribution, with excellent agreement with the experimental observation. We show thatfriction only tunes the mean cell order and, consequently, the exponential decay rate and the pack-ing fraction. We further show that cells, which can be very large in such systems, are short-lived,implying that our systems are liquid-like rather than glassy..

As a glass-forming liquid is cooled, the dynamics of its constituent particles changes from being liquid-like to more solid-like. The solidity of the resulting glassy material is believed to be due to a cage-formation process, whereby the motion of individual particles is increasingly constrained by neighbouring particles. This process begins at the temperature (or particle density) at which the glass-forming liquid first shows signs of glassy dynamics; however, the details of how the cages form remain unclear. Here we study cage formation at the particle level in a two-dimensional colloidal suspension (a glass-forming liquid). We use focused lasers to perturb the suspension at the particle level and monitor the nonlinear dynamic response of the system using video microscopy. All observables that we consider respond non-monotonically as a function of the particle density, peaking at the density at which glassy dynamics is first observed. We identify this maximum response as being due to cage formation, quantified by the appearance of domains in which particles move in a cooperative manner. As the particle density increases further, these local domains become increasingly rigid and dominate the macroscale particle dynamics. This microscale rheological deformation approach demonstrates that cage formation in glass-forming liquids is directly related to the merging of such domains, and reveals the first step in the transformation of liquids to glassy materials.

Using molecular dynamics simulations we investigate the dependence of the structural and vibrational properties of the surfaces of sodo-silicate glasses on the sodium content as well as the nature of the surface. Two types of glass surfaces are considered: A melt-formed surface (MS) in which a liquid with a free surface has been cooled down into the glass phase and a fracture surface (FS) obtained by tensile loading of a glass sample. We find that the MS is more abundant in Na and non-bridging oxygen atoms than the FS and the bulk glass, whereas the FS has higher concentration of structural defects such as two-membered rings and under-coordinated Si than the MS. We associate these structural differences to the production histories of the glasses and the mobility of the Na ions. It is also found that for Na-poor systems the fluctuations in composition and local atomic charge density decay with a power-law as a function of distance from the surface while Na-rich systems show an exponential decay with a typical decay length of ≈ 2.3 Å. The vibrational density of states shows that the presence of the surfaces leads to a decrease of the characteristic frequencies in the system. The two-membered rings give rise to a pronounce band at ≈ 880 cm−1 which is in good agreement experimental observations.

Disordered systems like liquids, gels, glasses, or granular materials are not only ubiquitous in daily life and in industrial applications, but they are also crucial for the mechanical stability of cells or the transport of chemical and biological agents in living organisms. Despite the importance of these systems, their microscopic structure is understood only on a rudimentary level, thus in stark contrast to the case of gases and crystals. Since scattering experiments and analytical calculations usually give only structural information that is spherically averaged, the three-dimensional (3D) structure of disordered systems is basically unknown. Here, we introduce a simple method that allows probing of the 3D structure of such systems. Using computer simulations, we find that hard sphere-like liquids have on intermediate and large scales a simple structural order given by alternating layers with icosahedral and dodecahedral symmetries, while open network liquids like silica have a structural order with tetrahedral symmetry. These results show that liquids have a highly nontrivial 3D structure and that this structural information is encoded in nonstandard correlation functions.