Chain structure of liquid selenium investigated by a tight-binding Monte Carlo simulation

Crucial Role of S8-Rings in Structural, Relaxation, Vibrational and Electronic Properties of LiquidSulfur close to the λ Transition

Liquid sulfur has been studied by density-functional based molecular-dynamics simulations at different temperatures ranging from 400 K up to 700 K across the well-documented λ-transition. Structure models containing either a majority of Sn chains or S8 rings are considered and compared to experimental data from X-ray scattering. The comparison suggests a liquid structure of a majority of 2-fold sulfur at low temperature, dominated by S8 rings that open progressively upon temperature increase. Typical features associated with such rings are analyzed and indicate that they contribute to a specific third correlating distance in the pair correlation function and to a contribution at low wavevector k in reciprocal space. The vibrational properties of liquid sulfur are also considered and indicate a contribution at 60 meV that is associated with both chains and rings, albeit the latter lead to a more intense peak at this wavenumber. The underlying network structure also impacts the dynamic properties of the melts which display enhanced dynamic heterogeneities when S8 rings are present. The analysis of the electronic Kohn-Sham energies shows insulating character with a gap of about ≃2.0 eV, albeit the presence of localized mid-gap states is acknowledged that can be associated, in part, with the presence of S6rings.

Dynamic heterogeneity, cooperative motion, and Johari-Goldstein [Formula: see text]-relaxation in a metallic glass-forming material exhibiting a fragile-to-strong transition

We investigate the Johari-Goldstein (JG) [Formula: see text]-relaxation process in a model metallic glass-forming (GF) material ([Formula: see text]), previously studied extensively by both frequency-dependent mechanical measurements and simulation studies devoted to equilibrium properties, by molecular dynamics simulations based on validated and optimized interatomic potentials with the primary aim of better understanding the nature of this universal relaxation process from a dynamic heterogeneity (DH) perspective. The present relatively low temperature and long-time simulations reveal a direct correspondence between the JG [Formula: see text]-relaxation time [Formula: see text] and the lifetime of the mobile particle clusters [Formula: see text], defined as in previous DH studies, a relationship dual to the corresponding previously observed relationship between the [Formula: see text]-relaxation time [Formula: see text] and the lifetime of immobile particle clusters [Formula: see text]. Moreover, we find that the average diffusion coefficient D nearly coincides with [Formula: see text] of the smaller atomic species (Al) and that the 'hopping time' associated with D coincides with [Formula: see text] to within numerical uncertainty, both trends being in accord with experimental studies. This indicates that the JG [Formula: see text]-relaxation is dominated by the smaller atomic species and the observation of a direct relation between this relaxation process and rate of molecular diffusion in GF materials at low temperatures where the JG [Formula: see text]-relaxation becomes the prevalent mode of structural relaxation. As an unanticipated aspect of our study, we find that [Formula: see text] exhibits fragile-to-strong (FS) glass formation, as found in many other metallic GF liquids, but this fact does not greatly alter the geometrical nature of DH in this material and the relation of DH to dynamical properties. On the other hand, the temperature dependence of the DH and dynamical properties, such as the structural relaxation time, can be significantly altered from 'ordinary' GF liquids.

The Use of Scattering Data in the Study of the Molecular Organisation of Polymers in the Non-Crystalline State

Scattering data for polymers in the non-crystalline state, i.e., the glassy state or the molten state, may appear to contain little information. In this work, we review recent developments in the use of scattering data to evaluate in a quantitative manner the molecular organization of such polymer systems. The focus is on the local structure of chain segments, on the details of the chain conformation and on the imprint the inherent chemical connectivity has on this structure. We show the value of tightly coupling the scattering data to atomistic-level computer models. We show how quantitative information about the details of the chain conformation can be obtained directly using a model built from definitions of relatively few parameters. We show how scattering data may be supplemented with data from specific deuteration sites and used to obtain information hidden in the data. Finally, we show how we can exploit the reverse Monte Carlo approach to use the data to drive the convergence of the scattering calculated from a 3d atomistic-level model with the experimental data. We highlight the importance of the quality of the scattering data and the value in using broad Q scattering data obtained using neutrons. We illustrate these various methods with results drawn from a diverse range of polymers.

Self-healing effects in a semi-ordered liquid for stable electronic conversion of high-energy radiation

Radiation damage in solid-state semiconductors has, until now, placed strict limitations on the acceptable decay energies of radioisotopes in radiovoltaic cells. Relegation to low-energy beta-emitting isotopes has minimized the power output from these devices and limited the technology's ability to deliver greater energy densities and longer lifetimes than conventional batteries. We demonstrate the self-healing abilities of a liquid-phase semiconducting alloy which can withstand high-energy alpha radiation. Neutron diffraction of liquid selenium-sulfur shows the liquid phase repairing damage sustained in the irradiation of the solid phase. This self-healing behavior results in long-lived power output in a liquid selenium-sulfur alphavoltaic cell. To the best of our knowledge, this marks the only successful demonstration of resistance to high-energy radiation (>500 keV) in a semiconducting material. This new robustness can potentially allow increases to the available energy density in radiovoltaic cells near 1000 times the current state of the art.

Molecular dynamics simulation of structural and dynamic properties of selenium structures with different degrees of amorphization

We investigated Se structures of different degrees of disorder ranging from a 5% up to a 95% degree of amorphization. Starting from a trigonal crystalline structure we applied different strategies to introduce disorder into the Se configurations by irradiating atoms from their crystalline equilibrium positions. According to the symmetry of the trigonal phase, we introduced three types of disorder, i.e. the first type where only atoms forming layers of complete helical chains are shifted from their original positions (the thickness of these layers is chosen to represent the chosen degree of amorphicity), the second type where only atoms in planes-of respective thicknesses-lying perpendicular to the chains are displaced and the third type where only randomly chosen atoms are shifted from their crystalline equilibrium positions. After a thermal treatment of these disordered starting configurations, we calculated structural and dynamic properties (i.e. pair-correlation function and vibrational spectrum) and compared the results to both the original crystalline data and results obtained from corresponding glass structures.

Effect of through-space electron transfer on infrared spectrum of amorphous selenium

In this paper we present theoretical analyses on an infrared (IR) spectrum of amorphous selenium. The system is described by a 216-atom-chain model, and a set of molecular-dynamics simulations is performed to generate vitreous structures and vibrational modes. To describe an electronic structure of the system we employ a complete neglect of differential overlap model parametrized by ab initio cluster calculations. An IR intensity is evaluated with the Berry-phase formula for an electronic polarization. The effect of the through-space electron transfer on the IR spectrum is studied by artificially changing the magnitude of matrix elements associated with the electron transfer between nonbonded atoms in the chain. We find that the through-space electron transfer leads to (i) the enhancement of the bending IR peak at 135 cm(-1) and (ii) the appearance of a new low-frequency peak around 50 cm(-1), thus resulting in a good agreement with the experiment. The mechanism is discussed by a simple dipole model.

Raman scattering study on structural and dynamical features of noncrystalline selenium

We report on a detailed, temperature-dependent, off-resonant Raman scattering study of glassy and supercooled selenium. Raman spectra in the frequency regime of the first-order scattering (5-450 cm(-1)) have been recorded over a wide temperature range, i.e., 143-353 K. To facilitate the analysis, the spectra have intuitively been divided in three spectral regions. The analysis of the high frequency region (bond-stretching vibrational modes) yielded information on the rings-chains equilibrium. In particular, the polymer content was found to amount to more than 85% around the glass transition temperature, exhibiting a weak temperature dependence, which extrapolates nicely to the high-temperature dissolution data. The intermediate frequency range (representative of the medium-range structural order) was treated together with the low frequency regime (where low-energy excitations, i.e., the quasielastic line and the Boson peak are the dominant contributions) owing to their strong overlap. The study of the bond-bending regime revealed information which made it possible to clarify the role of ringlike and chainlike fragments incorporated in polymeric molecules. The temperature evolution of the Boson peak and the frequency dependence of the Raman coupling coefficient Comega were also determined. An attempt to decompose the partial contribution of the pure Boson peak to Comega revealed valuable information concerning the limiting (omega-->0) behavior of the coupling coefficient.