Computer Simulations of Amorphous Semiconductors

The subject of this thesis is the class of materials known as disordered solids. These solids are also referred to as non-crystalline solids, or amorphous solids or glasses. What these materials have in common is that their structure is not crystalline on any significant scale. The aim of this thesis is to provide insight into the atomic structure of disordered solids via computer simulation. Simply put, we try to find out where the atoms are located with respect to each other. For these materials, experimental techniques alone are not sufficient to determine the structure. As a result, computer simulations are often used. We describe a number of simulation techniques that allow for the investigation of large systems, containing over 100,000 particles. Moreover, we present a new technique to determine the configurational entropy

Phase separation in fluids exposed to spatially periodic external fields

When a fluid is confined within a spatially periodic external field, the liquid-vapor transition is replaced by a different transition called laser-induced condensation (LIC) [ Götze et al. Mol. Phys. 101 1651 (2003)]. In d=3 dimensions, the periodic field induces an additional phase, characterized by large density modulations along the field direction. At the triple point, all three phases (modulated, vapor, and liquid) coexist. At temperatures slightly above the triple point and for low (high) values of the chemical potential, two-phase coexistence between the modulated phase and the vapor (liquid) is observed; by increasing the temperature further, both coexistence regions terminate in critical points. In this paper, we reconsider LIC using the Ising model to resolve a number of open issues. To be specific, we (1) determine the universality class of the LIC critical points and elucidate the nature of the correlations along the field direction, (2) present a mean-field analysis to show how the LIC phase diagram changes as a function of the field wavelength and amplitude, (3) develop a simulation method by which the extremely low tension of the interface between modulated and vapor or liquid phase can be measured, (4) present a finite-size scaling analysis to accurately extract the LIC triple point from finite-size simulation data, and (5) consider the fate of LIC in d=2 dimensions

Large well-relaxed models of vitreous silica, coordination numbers and entropy

A Monte Carlo method is presented for the simulation of vitreous silica. Well-relaxed networks of vitreous silica are generated containing up to 300,000 atoms. The resulting networks, quenched under the BKS potential, display smaller bond-angle variations and lower defect concentrations, as compared to networks generated with molecular dynamics. The total correlation functions T(r) of our networks are in excellent agreement with neutron scattering data, provided that thermal effects and the maximum inverse wavelength used in the experiment are included in the comparison. A procedure commonly used in experiments to obtain coordination numbers from scattering data is to fit peaks in rT(r) with a gaussian. We show that this procedure can easily produce incorrect results. Finally, we estimate the configurational entropy of vitreous silica.Comment: 7 pages, 4 figures (two column version to save paper

Towards device-size atomistic models of amorphous silicon

The atomic structure of amorphous materials is believed to be well described by the continuous random network model. We present an algorithm for the generation of large, high-quality continuous random networks. The algorithm is a variation of the "sillium" approach introduced by Wooten, Winer, and Weaire. By employing local relaxation techniques, local atomic rearrangements can be tried that scale almost independently of system size. This scaling property of the algorithm paves the way for the generation of realistic device-size atomic networks.Comment: 7 pages, 3 figure

Energy landscape of relaxed amorphous silicon

We analyze the structure of the energy landscape of a well-relaxed 1000-atom model of amorphous silicon using the activation-relaxation technique (ART nouveau). Generating more than 40,000 events starting from a single minimum, we find that activated mechanisms are local in nature, that they are distributed uniformly throughout the model and that the activation energy is limited by the cost of breaking one bond, independently of the complexity of the mechanism. The overall shape of the activation-energy-barrier distribution is also insensitive to the exact details of the configuration, indicating that well-relaxed configurations see essentially the same environment. These results underscore the localized nature of relaxation in this material.Comment: 8 pages, 12 figure

Wang-Landau study of the 3D Ising model with bond disorder

We implement a two-stage approach of the Wang-Landau algorithm to investigate the critical properties of the 3D Ising model with quenched bond randomness. In particular, we consider the case where disorder couples to the nearest-neighbor ferromagnetic interaction, in terms of a bimodal distribution of strong versus weak bonds. Our simulations are carried out for large ensembles of disorder realizations and lattices with linear sizes $L$ in the range $L=8-64$. We apply well-established finite-size scaling techniques and concepts from the scaling theory of disordered systems to describe the nature of the phase transition of the disordered model, departing gradually from the fixed point of the pure system. Our analysis (based on the determination of the critical exponents) shows that the 3D random-bond Ising model belongs to the same universality class with the site- and bond-dilution models, providing a single universality class for the 3D Ising model with these three types of quenched uncorrelated disorder.Comment: 7 pages, 7 figures, to be published in Eur. Phys. J.

Structural Information in Two-Dimensional Patterns: Entropy Convergence and Excess Entropy

We develop information-theoretic measures of spatial structure and pattern in more than one dimension. As is well known, the entropy density of a two-dimensional configuration can be efficiently and accurately estimated via a converging sequence of conditional entropies. We show that the manner in which these conditional entropies converge to their asymptotic value serves as a measure of global correlation and structure for spatial systems in any dimension. We compare and contrast entropy-convergence with mutual-information and structure-factor techniques for quantifying and detecting spatial structure.Comment: 11 pages, 5 figures, http://www.santafe.edu/projects/CompMech/papers/2dnnn.htm

Monte Carlo Methods for Estimating Interfacial Free Energies and Line Tensions

Excess contributions to the free energy due to interfaces occur for many problems encountered in the statistical physics of condensed matter when coexistence between different phases is possible (e.g. wetting phenomena, nucleation, crystal growth, etc.). This article reviews two methods to estimate both interfacial free energies and line tensions by Monte Carlo simulations of simple models, (e.g. the Ising model, a symmetrical binary Lennard-Jones fluid exhibiting a miscibility gap, and a simple Lennard-Jones fluid). One method is based on thermodynamic integration. This method is useful to study flat and inclined interfaces for Ising lattices, allowing also the estimation of line tensions of three-phase contact lines, when the interfaces meet walls (where "surface fields" may act). A generalization to off-lattice systems is described as well. The second method is based on the sampling of the order parameter distribution of the system throughout the two-phase coexistence region of the model. Both the interface free energies of flat interfaces and of (spherical or cylindrical) droplets (or bubbles) can be estimated, including also systems with walls, where sphere-cap shaped wall-attached droplets occur. The curvature-dependence of the interfacial free energy is discussed, and estimates for the line tensions are compared to results from the thermodynamic integration method. Basic limitations of all these methods are critically discussed, and an outlook on other approaches is given

Critical aspects of the random-field Ising model

We investigate the critical behavior of the three-dimensional random-field Ising model (RFIM) with a Gaussian field distribution at zero temperature. By implementing a computational approach that maps the ground-state of the RFIM to the maximum-flow optimization problem of a network, we simulate large ensembles of disorder realizations of the model for a broad range of values of the disorder strength h and system sizes  = L3, with L ≤ 156. Our averaging procedure outcomes previous studies of the model, increasing the sampling of ground states by a factor of 103. Using well-established finite-size scaling schemes, the fourth-order’s Binder cumulant, and the sample-to-sample fluctuations of various thermodynamic quantities, we provide high-accuracy estimates for the critical field hc, as well as the critical exponents ν, β/ν, and γ̅/ν of the correlation length, order parameter, and disconnected susceptibility, respectively. Moreover, using properly defined noise to signal ratios, we depict the variation of the self-averaging property of the model, by crossing the phase boundary into the ordered phase. Finally, we discuss the controversial issue of the specific heat based on a scaling analysis of the bond energy, providing evidence that its critical exponent α ≈ 0−

Revisiting the scaling of the specific heat of the three-dimensional random-field Ising model

We revisit the scaling behavior of the specific heat of the three-dimensional random-field Ising model with a Gaussian distribution of the disorder. Exact ground states of the model are obtained using graph-theoretical algorithms for different strengths = 268 3 spins. By numerically differentiating the bond energy with respect to h, a specific-heat-like quantity is obtained whose maximum is found to converge to a constant in the thermodynamic limit. Compared to a previous study following the same approach, we have studied here much larger system sizes with an increased statistical accuracy. We discuss the relevance of our results under the prism of a modified Rushbrooke inequality for the case of a saturating specific heat. Finally, as a byproduct of our analysis, we provide high-accuracy estimates of the critical field hc = 2.279(7) and the critical exponent of the correlation exponent ν = 1.37(1), in excellent agreement to the most recent computations in the literature