The Debye-Waller factor of liquid silica: Theory and simulation

We show that the prediction of mode-coupling theory for a model of a network-forming strong glass-former correctly describes the wave-vector dependence of the Debye-Waller factor. To obtain a good description it is important to take into account the triplet correlation function c_3, which we evaluate from a computer simulation. Our results support the possibility that this theory is able to accurately describe the non-ergodicity parameters of simple as well as of network-forming liquids.Comment: 5 pages of Latex, 3 figure

Computer investigation of the energy landscape of amorphous silica

The multidimensional topography of the collective potential energy function of a so-called strong glass former (silica) is analyzed by means of classical molecular dynamics calculations. Features qualitatively similar to those of fragile glasses are recovered at high temperatures : in particular an intrinsic characteristic temperature $T_c\simeq 3500$K is evidenced above which the system starts to investigate non-harmonic potential energy basins. It is shown that the anharmonicities are essentially characterized by a roughness appearing in the potential energy valleys explored by the system for temperatures above $T_c$.Comment: 5 pages; accepted for publication in PR

Matrix controlled channel diffusion of sodium in amorphous silica

To find the origin of the diffusion channels observed in sodium-silicate glasses, we have performed classical molecular dynamics simulations of Na$_2$O--4SiO$_2$ during which the mass of the Si and O atoms has been multiplied by a tuning coefficient. We observe that the channels disappear and that the diffusive motion of the sodium atoms vanishes if this coefficient is larger than a threshold value. Above this threshold the vibrational states of the matrix are not compatible with those of the sodium ions. We interpret hence the decrease of the diffusion by the absence of resonance conditions.Comment: 5 pages, 4 figure

Static and Dynamic Properties of a Viscous Silica Melt Molecular Dynamics Computer Simulations

We present the results of a large scale molecular dynamics computer simulation in which we investigated the static and dynamic properties of a silica melt in the temperature range in which the viscosity of the system changes from O(10^-2) Poise to O(10^2) Poise. We show that even at temperatures as high as 4000 K the structure of this system is very similar to the random tetrahedral network found in silica at lower temperatures. The temperature dependence of the concentration of the defects in this network shows an Arrhenius law. From the partial structure factors we calculate the neutron scattering function and find that it agrees very well with experimental neutron scattering data. At low temperatures the temperature dependence of the diffusion constants $D$ shows an Arrhenius law with activation energies which are in very good agreement with the experimental values. With increasing temperature we find that this dependence shows a cross-over to one which can be described well by a power-law, D\propto (T-T_c)^gamma. The critical temperature T_c is 3330 K and the exponent gamma is close to 2.1. Since we find a similar cross-over in the viscosity we have evidence that the relaxation dynamics of the system changes from a flow-like motion of the particles, as described by the ideal version of mode-coupling theory, to a hopping like motion. We show that such a change of the transport mechanism is also observed in the product of the diffusion constant and the life time of a Si-O bond, or the space and time dependence of the van Hove correlation functions.Comment: 30 pages of Latex, 14 figure

Relationship between Structure, Entropy and Diffusivity in Water and Water-like Liquids

Anomalous behaviour of the excess entropy ($S_e$) and the associated scaling relationship with diffusivity are compared in liquids with very different underlying interactions but similar water-like anomalies: water (SPC/E and TIP3P models), tetrahedral ionic melts (SiO$_2$ and BeF$_2$) and a fluid with core-softened, two-scale ramp (2SRP) interactions. We demonstrate the presence of an excess entropy anomaly in the two water models. Using length and energy scales appropriate for onset of anomalous behaviour, the density range of the excess entropy anomaly is shown to be much narrower in water than in ionic melts or the 2SRP fluid. While the reduced diffusivities ($D^*$) conform to the excess entropy scaling relation, $D^* =A\exp (\alpha S_e)$ for all the systems (Y. Rosenfeld, Phys. Rev. A {\bf 1977}, {\it 15}, 2545), the exponential scaling parameter, $\alpha$, shows a small isochore-dependence in the case of water. Replacing $S_e$ by pair correlation-based approximants accentuates the isochore-dependence of the diffusivity scaling. Isochores with similar diffusivity scaling parameters are shown to have the temperature dependence of the corresponding entropic contribution. The relationship between diffusivity, excess entropy and pair correlation approximants to the excess entropy are very similar in all the tetrahedral liquids.Comment: 24 pages, 4 figures, to be published in Journal of Physical Chemistry

The electronic structure of amorphous silica: A numerical study

We present a computational study of the electronic properties of amorphous SiO2. The ionic configurations used are the ones generated by an earlier molecular dynamics simulations in which the system was cooled with different cooling rates from the liquid state to a glass, thus giving access to glass-like configurations with different degrees of disorder [Phys. Rev. B 54, 15808 (1996)]. The electronic structure is described by a tight-binding Hamiltonian. We study the influence of the degree of disorder on the density of states, the localization properties, the optical absorption, the nature of defects within the mobility gap, and on the fluctuations of the Madelung potential, where the disorder manifests itself most prominently. The experimentally observed mismatch between a photoconductivity threshold of 9 eV and the onset of the optical absorption around 7 eV is interpreted by the picture of eigenstates localized by potential energy fluctuations in a mobility gap of approximately 9 eV and a density of states that exhibits valence and conduction band tails which are, even in the absence of defects, deeply located within the former band gap.Comment: 21 pages of Latex, 5 eps figure

Frequency dependent specific heat of viscous silica

We apply the Mori-Zwanzig projection operator formalism to obtain an expression for the frequency dependent specific heat c(z) of a liquid. By using an exact transformation formula due to Lebowitz et al., we derive a relation between c(z) and K(t), the autocorrelation function of temperature fluctuations in the microcanonical ensemble. This connection thus allows to determine c(z) from computer simulations in equilibrium, i.e. without an external perturbation. By considering the generalization of K(t) to finite wave-vectors, we derive an expression to determine the thermal conductivity \lambda from such simulations. We present the results of extensive computer simulations in which we use the derived relations to determine c(z) over eight decades in frequency, as well as \lambda. The system investigated is a simple but realistic model for amorphous silica. We find that at high frequencies the real part of c(z) has the value of an ideal gas. c'(\omega) increases quickly at those frequencies which correspond to the vibrational excitations of the system. At low temperatures c'(\omega) shows a second step. The frequency at which this step is observed is comparable to the one at which the \alpha-relaxation peak is observed in the intermediate scattering function. Also the temperature dependence of the location of this second step is the same as the one of the $\alpha-$peak, thus showing that these quantities are intimately connected to each other. From c'(\omega) we estimate the temperature dependence of the vibrational and configurational part of the specific heat. We find that the static value of c(z) as well as \lambda are in good agreement with experimental data.Comment: 27 pages of Latex, 8 figure

Molecular structural order and anomalies in liquid silica

The present investigation examines the relationship between structural order, diffusivity anomalies, and density anomalies in liquid silica by means of molecular dynamics simulations. We use previously defined orientational and translational order parameters to quantify local structural order in atomic configurations. Extensive simulations are performed at different state points to measure structural order, diffusivity, and thermodynamic properties. It is found that silica shares many trends recently reported for water [J. R. Errington and P. G. Debenedetti, Nature 409, 318 (2001)]. At intermediate densities, the distribution of local orientational order is bimodal. At fixed temperature, order parameter extrema occur upon compression: a maximum in orientational order followed by a minimum in translational order. Unlike water, however, silica's translational order parameter minimum is broad, and there is no range of thermodynamic conditions where both parameters are strictly coupled. Furthermore, the temperature-density regime where both structural order parameters decrease upon isothermal compression (the structurally anomalous regime) does not encompass the region of diffusivity anomalies, as was the case for water.Comment: 30 pages, 8 figure

Path integral for half-binding potentials as quantum mechanical analog for black hole partition functions

The semi-classical approximation to black hole partition functions is not well-defined, because the classical action is unbounded and the first variation of the uncorrected action does not vanish for all variations preserving the boundary conditions. Both problems can be solved by adding a Hamilton-Jacobi counterterm. I show that the same problem and solution arises in quantum mechanics for half-binding potentials.Comment: 6 pages, proceedings contribution to "Path integrals - New Trends and Perspectives", Dresden, September 200