An official website of the United States government
We consider a common problem occurring after using a statistical process control (SPC) method based on three-dimensional measurements: locate where on the surface of the part that triggered an out-of-control alarm there is a significant shape difference with respect to either an in-control part or its nominal (computer-aided design (CAD)) design. In the past, only registration-based solutions existed for this problem, which first orient and locate the part and its nominal design under the same frame of reference. Recently, spectral Laplacian methods have been proposed for the SPC of discrete parts and their measured surface meshes. These techniques provide an intrinsic solution to the SPC problem: that is, a solution exclusively based on data whose coordinates lie on the surfaces without making reference to their ambient space, thus avoiding registration. Registration-free methods avoid the computationally expensive, nonconvex registration step needed to align the parts as required by previous methods, eliminating registration errors, and they are important in industry because of the increasing use of portable noncontact scanners. In this paper, we first present a new registration-free solution to the post-SPC part defect localization problem. The approach uses a spectral decomposition of the Laplace–Beltrami operator in order to construct a functional map between the CAD and measured manifolds to locate defects on the suspected part. A computational complexity analysis demonstrates the approach scales better with the mesh size and is more stable than a registration-based approach. To reduce computational expense, a new mesh partitioning algorithm is presented to find a region of interest on the surface of the part where defects are more likely to exist. The functional map method involves a large number of point-to-point comparisons based on noisy measurements, and a new statistical thresholding method used to filter the false positives in the underlying massive multiple comparisons problem is also provided. Funding: This research was partially funded by the National Science Foundation [Grant CMMI 2121625]. Data Ethics & Reproducibility Note: There are no data ethics considerations. The code capsule is available on Code Ocean at https://codeocean.com/capsule/4615101/tree/v1 and in the e-Companion to this article (available https://doi.org/10.1287/ijds.2023.0030 ).
more »
« less
Adaptive Exploration and Optimization of Materials Crystal Structures
A central problem of materials science is to determine whether a hypothetical material is stable without being synthesized, which is mathematically equivalent to a global optimization problem on a highly nonlinear and multimodal potential energy surface (PES). This optimization problem poses multiple outstanding challenges, including the exceedingly high dimensionality of the PES, and that PES must be constructed from a reliable, sophisticated, parameters-free, and thus very expensive computational method, for which density functional theory (DFT) is an example. DFT is a quantum mechanics-based method that can predict, among other things, the total potential energy of a given configuration of atoms. DFT, although accurate, is computationally expensive. In this work, we propose a novel expansion-exploration-exploitation framework to find the global minimum of the PES. Starting from a few atomic configurations, this “known” space is expanded to construct a big candidate set. The expansion begins in a nonadaptive manner, where new configurations are added without their potential energy being considered. A novel feature of this step is that it tends to generate a space-filling design without the knowledge of the boundaries of the domain space. If needed, the nonadaptive expansion of the space of configurations is followed by adaptive expansion, where “promising regions” of the domain space (those with low-energy configurations) are further expanded. Once a candidate set of configurations is obtained, it is simultaneously explored and exploited using Bayesian optimization to find the global minimum. The methodology is demonstrated using a problem of finding the most stable crystal structure of aluminum. History: Kwok Tsui served as the senior editor for this article. Funding: The authors acknowledge a U.S. National Science Foundation Grant DMREF-1921873 and XSEDE through Grant DMR170031. Data Ethics & Reproducibility Note: The code capsule is available on Code Ocean at https://codeocean.com/capsule/3366149/tree and in the e-Companion to this article (available at https://doi.org/10.1287/ijds.2023.0028 ).
more »
« less
- Award ID(s):
- 1921873
- PAR ID:
- 10443720
- Date Published:
- Journal Name:
- INFORMS Journal on Data Science
- ISSN:
- 2694-4022
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract In the 1980s, Nelson, Fraser, and Klemperer (NFK) published an experimentally derived structure of the ammonia dimer dramatically different from the structure determined computationally, which led these authors to the question “Does ammonia hydrogen bond?". This question has not yet been answered satisfactorily. To answer it, we have developed an ab initio potential energy surface (PES) for this dimer at the limits of the current computational capabilities and performed essentially exact six-dimensional calculations of the vibration-rotation-tunneling (VRT) spectra of NH 3 -NH 3 and ND 3 -ND 3 , obtaining an unprecedented agreement with experimental spectra. In agreement with other recent electronic structure calculations, the global minimum on the PES is in a substantially bent hydrogen-bonded configuration. Since the bottom of the PES is exceptionally flat, the dimer is extremely fluxional and the probability of finding it in configurations that are not hydrogen bonded is high. Nevertheless, the probability of hydrogen-bonded configurations is large enough to consider the ammonia dimer to be hydrogen bonded. We also show that NFK’s inference that the ammonia dimer is nearly rigid actually results from unusual cancellations between quantum effects that generate differences in spectra of different isotopologues.more » « less
-
We investigate the anisotropic thermal expansion behavior of a co- crystalline system composed of 4,40-azopyridine and trimesic acid (TMA-azo). Using variable-temperature single-crystal X-ray diffrac- tion (SC-XRD), low-frequency Raman spectroscopy, and terahertz time-domain spectroscopy (THz-TDS), we observe significant temperature-induced shifting and broadening of the vibrational absorption features, indicating changes in the intermolecular potential. Our findings reveal that thermal expansion is driven by anharmonic interactions and the potential energy topography, rather than increased molecular dynamics. Density functional the- ory (DFT) simulations support these results, highlighting significant softening of the potential energy surface (PES) with temperature. This comprehensive approach offers valuable insights into the relationship between structural dynamics and thermal properties, providing a robust framework for designing materials with tailored thermal expansion characteristics.more » « less
-
The quantum chemistry community has developed analytic forces for approximate electronic excited states to enable walking on excited state potential energy surfaces (PES). One can thereby computationally characterize excited state minima and saddle points. Always implicit in using this machinery is the fact that an excited state PES only exists within the realm of the Born-Oppenheimer approximation, where the nuclear and electronic degrees of freedom separate. This work demonstrates through ab initio calculations and simple nonadiabatic dynamics that some excited state minimum structures are fantastical: they appear to exist as stable configurations only as a consequence of the PES construct, rather than being physically observable. Each fantastical structure exhibits an unphysically high predicted harmonic frequency and associated force constant. This fact can serve as a valuable diagnostic of when an optimized excited state structure is non-observable. The origin of this phenomenon can be attributed to the coupling between different electronic states. As PESs approach one another, the upper surface can form a minimum that is very close to a near-touching point. The force constant, evaluated at this minimum, relates to the strength of the electronic coupling rather than to any characteristic excited state vibration. Nonadiabatic dynamics results using a Landau-Zener model illustrate that fantastical excited state structures have extremely short lifetimes on the order of a few femtoseconds. Their appearance in a calculation signals the presence of a nearby conical intersection through which the system will rapidly cross to a lower surface.more » « less
-
Metal-organic frameworks (MOFs) are nanoporous compounds composed of metal ions and organic linkers. MOFs play an important role in industrial applications such as gas separation, gas purification, and electrolytic catalysis. Important MOF properties such a potential energy are currently computed via techniques such as density functional theory (DFT). Although DFT provides accurate results, it is computationally costly. We propose a machine learning approach for estimating the potential energy of candidate MOFs, decomposing it into separate pair-wise atomic interactions using a graph neural network. Such a technique will allow high-throughput screening of candidates MOFs. We also generate a database of 50,000 spatial configurations and high quality potential energy values using DFT.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1287/ijds.2023.0028
