skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Disorder by design: A data-driven approach to amorphous semiconductors without total-energy functionals
Abstract X-ray diffraction, Amorphous silicon, Multi-objective optimization, Monte Carlo methods. This paper addresses a difficult inverse problem that involves the reconstruction of a three-dimensional model of tetrahedral amorphous semiconductors via inversion of diffraction data. By posing the material-structure determination as a multiobjective optimization program, it has been shown that the problem can be solved accurately using a few structural constraints, but no total-energy functionals/forces, which describe the local chemistry of amorphous networks. The approach yields highly realistic models of amorphous silicon, with no or only a few coordination defects (≤1%), a narrow bond-angle distribution of width 9–11.5°, and an electronic gap of 0.8–1.4 eV. These data-driven information-based models have been found to produce electronic and vibrational properties of a -Si that match accurately with experimental data and rival that of the Wooten-Winer-Weaire models. The study confirms the effectiveness of a multiobjective optimization approach to the structural determination of complex materials, and resolves a long-standing dispute concerning the uniqueness of a model of tetrahedral amorphous semiconductors obtained via inversion of diffraction data.  more » « less
Award ID(s):
1833035
PAR ID:
10328599
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Scientific Reports
Volume:
10
Issue:
1
ISSN:
2045-2322
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We investigate 29Si nuclear magnetic resonance (NMR) chemical shifts, δiso, of silicon nitride. Our goal is to relate the local structure to the NMR signal and, thus, provide the means to extract more information from the experimental 29Si NMR spectra in this family of compounds. We apply structural modeling and the gauge-included projector augmented wave (GIPAW) method within density functional theory (DFT) calculations. Our models comprise known and hypothetical crystalline Si3N4, as well as amorphous Si3N4 structures. We find good agreement with available experimental 29Si NMR data for tetrahedral Si[4] and octahedral Si[6] in crystalline Si3N4, predict the chemical shift of a trigonal-bipyramidal Si[5] to be about −120 ppm, and quantify the impact of Si-N bond lengths on 29Si δiso. We show through computations that experimental 29Si NMR data indicates that silicon dicarbodiimide, Si(NCN)2 exhibits bent Si-N-C units with angles of about 143° in its structure. A detailed investigation of amorphous silicon nitride shows that an observed peak asymmetry relates to the proximity of a fifth N neighbor in non-bonding distance between 2.5 and 2.8 Å to Si. We reveal the impact of both Si-N(H)-Si bond angle and Si-N bond length on 29Si δiso in hydrogenated silicon nitride structure, silicon diimide Si(NH)2. 
    more » « less
  2. null (Ed.)
    Inverse multiobjective optimization provides a general framework for the unsupervised learning task of inferring parameters of a multiobjective decision making problem (DMP), based on a set of observed decisions from the human expert. However, the performance of this framework relies critically on the availability of an accurate DMP, sufficient decisions of high quality, and a parameter space that contains enough information about the DMP. To hedge against the uncertainties in the hypothetical DMP, the data, and the parameter space, we investigate in this paper the distributionally robust approach for inverse multiobjective optimization. Specifically, we leverage the Wasserstein metric to construct a ball centered at the empirical distribution of these decisions. We then formulate a Wasserstein distributionally robust inverse multiobjective optimization problem (WRO-IMOP) that minimizes a worst-case expected loss function, where the worst case is taken over all distributions in the Wasserstein ball. We show that the excess risk of the WRO-IMOP estimator has a sub-linear convergence rate. Furthermore, we propose the semi-infinite reformulations of the WRO-IMOP and develop a cutting-plane algorithm that converges to an approximate solution in finite iterations. Finally, we demonstrate the effectiveness of our method on both a synthetic multiobjective quadratic program and a real world portfolio optimization problem. 
    more » « less
  3. null (Ed.)
    Inverse multiobjective optimization provides a general framework for the unsupervised learning task of inferring parameters of a multiobjective decision making problem (DMP), based on a set of observed decisions from the human expert. However, the performance of this framework relies critically on the availability of an accurate DMP, sufficient decisions of high quality, and a parameter space that contains enough information about the DMP. To hedge against the uncertainties in the hypothetical DMP, the data, and the parameter space, we investigate in this paper the distributionally robust approach for inverse multiobjective optimization. Specifically, we leverage the Wasserstein metric to construct a ball centered at the empirical distribution of these decisions. We then formulate a Wasserstein distributionally robust inverse multiobjective optimization problem (WRO-IMOP) that minimizes a worst-case expected loss function, where the worst case is taken over all distributions in the Wasserstein ball. We show that the excess risk of the WRO-IMOP estimator has a sub-linear convergence rate. Furthermore, we propose the semi-infinite reformulations of the WRO-IMOP and develop a cutting-plane algorithm that converges to an approximate solution in finite iterations. Finally, we demonstrate the effectiveness of our method on both a synthetic multiobjective quadratic program and a real world portfolio optimization problem. 
    more » « less
  4. In this study, a machine learning based computational approach has been developed to investigate the cation distribution in spinel crystals. The computational approach integrates the construction of datasets consisting of the energies calculated from density functional theory, the training of machine learning models to derive the relationship between system energy and structural features, and atomistic Monte Carlo simulations to sample the thermodynamic equilibrium structures of spinel crystals. It is found that the support vector machine model yields excellent performance in energy predictions based on spinel crystal structures. Furthermore, the developed computational approach has been applied to predict the cation distribution in single spinel MgAl2O4 and MgFe2O4 and double spinel MgAl2-aFeaO4. Agreeing with the available experimental data, the computational approach correctly predicts that the equilibrium degree of inversion of MgAl2O4 increases with temperature, whereas the degree of inversion of MgFe2O4 decreases with temperature. Additionally, it is predicted that the equilibrium occupancy of Mg cations at the tetrahedral and octahedral sites in MgAl2-aFeaO4 could be tuned as a function of chemical composition. Therefore, this study presents a reliable computational approach that can be extended to study the variation of cation distribution with processing temperature and chemical composition in a wide range of complex multi-cation spinel oxides with numerous applications. 
    more » « less
  5. Abstract Amorphous calcium carbonate is an important precursor for biomineralization in marine organisms. Key outstanding problems include understanding the structure of amorphous calcium carbonate and rationalizing its metastability as an amorphous phase. Here we report high-quality atomistic models of amorphous calcium carbonate generated using state-of-the-art interatomic potentials to help guide fits to X-ray total scattering data. Exploiting a recently developed inversion approach, we extract from these models the effective Ca⋯Ca interaction potential governing the structure. This potential contains minima at two competing distances, corresponding to the two different ways that carbonate ions bridge Ca2+-ion pairs. We reveal an unexpected mapping to the Lennard-Jones–Gauss model normally studied in the context of computational soft matter. The empirical model parameters for amorphous calcium carbonate take values known to promote structural complexity. We thus show that both the complex structure and its resilience to crystallization are actually encoded in the geometrically frustrated effective interactions between Ca2+ions. 
    more » « less