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: A Unified Analysis of Penalty-Based Collision Energies
We analyze a wide class of penalty energies used for contact response through the lens of a reduced frame. Applying our analysis to both spring-based and barrier-based energies, we show that we can obtain closed-form, analytic eigensystems that can be used to guarantee positive semidefiniteness in implicit solvers. Our approach is both faster than direct numerical methods, and more robust than approximate methods such as Gauss-Newton. Over the course of our analysis, we investigate physical interpretations for two separate notions of length. Finally, we showcase the stability of our analysis on challenging strand, cloth, and volume scenarios with large timesteps on the order of 1/40 s.  more » « less
Award ID(s):
2132280
PAR ID:
10431370
Author(s) / Creator(s):
;
Editor(s):
Ye, Yuting; Wang Huamin
Date Published:
Journal Name:
Proceedings of the ACM on computer graphics and interactive techniques
ISSN:
2577-6193
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Ye, Yuting; Wang Huamin (Ed.)
    Angle-based energies appear in numerous physics-based simulation models, including thin-shell bending and isotropic elastic strands. We present a generic analysis of these energies that allows us to analytically filter the negative eigenvalues of the second derivative (Hessian), which is critical for stable, implicit time integration. While these energies are usually formulated in terms of angles and positions, we propose an abstract edge stencil that succinctly parameterizes the edge deformation, and allows us to derive generic, closed-form analytical expressions for the energy eigensystems. The resultant eigenvectors have straightforward geometric interpretations. We demonstrate that our method is readily applicable to a variety of 2D and 3D angle-based elastic energies, including both cloth and strands, and is up to 7x faster than numerical eigendecomposition. 
    more » « less
  2. In this work, we study vacancy energetics in the equiatomic Nb-Mo-Ta-W alloy, especially vacancy formation and migration energies, using molecular statics calculations based on a spectral neighbor analysis potential specifically developed for Nb-Mo-Ta-W. We consider vacancy properties in bulk environments as well as near edge dislocation cores, including the effect of short-range order (SRO) by preparing supercells through Metropolis Monte-Carlo relaxations and temperature on the calculation. The nudged elastic band (NEB) method is applied to study vacancy migration energies. Our results show that both vacancy formation energies and vacancy migration energies are statistically distributed with a wide spread, on the order of 1.0 eV in some cases, and display a noticeable dependence on SRO. We find that, in some cases, vacancies can form with very low energies at edge dislocation cores, from which we hypothesize the formation of stable ‘superjogs’ on edge dislocation lines. Moreover, the large spread in vacancy formation energies results in an asymmetric thermal sampling of the formation energy distribution towards lower values. This gives rise to effective vacancy formation energies that are noticeably lower than the distribution averages. We study the effect that this phenomenon has on the vacancy diffusivity in the alloy and discuss the implications of our findings on the structural features of Nb-Mo-Ta-W. 
    more » « less
  3. This paper presents a new approach for predicting thermodynamic properties of perovskites that harnesses deep learning and crystal structure fingerprinting based on Hirshfeld surface analysis. It is demonstrated that convolutional neural network methods capture critical features embedded in two-dimensional Hirshfeld surface fingerprints that enable a quantitative assessment of the formation energy of perovskites. Building on our recent work on lattice parameter prediction from Hirshfeld surface calculations, we show how transfer learning can be used to speed up the training of the neural network, allowing multiple properties to be trained using the same feature extraction layers. We also predict formation energies for various perovskite polymorphs, and our predictions are found to give generally improved performance over a well-established graph network method, but with the methods better suited to different types of datasets. Analysis of the structure types within the dataset reveals the Hirshfeld surface-based method to excel for the less symmetric and similar structures, while the graph network performs better for very symmetric and similar structures. 
    more » « less
  4. The molecules 1,4-cyclohexadiene (unconjugated 1,4-CHD) and 1,3-cyclohexadiene (conjugated 1,3-CHD) both have two double bonds, but these bonds interact in different ways. These molecules have long served as examples of through-bond and through-space interactions, respectively, and their electronic structures have been studied in detail both experimentally and theoretically, with the experimental assignments being especially complete. The existence of Rydberg states interspersed with the valence states makes the quantum mechanical calculation of their spectra a challenging task. In this work, we explore the electronic excitation energies of 1,4-CHD and 1,3-CHD for both valence and Rydberg states by means of complete active space second-order perturbation theory (CASPT2), extended multi-state CASPT2 (XMS-CASPT2), and multiconfiguration pair-density functional theory (MC-PDFT); it is shown by comparison to experiment that MC-PDFT yields the most accurate results. We found that the inclusion of Rydberg orbitals in the active space not only enables the calculation of Rydberg excitation energies but also improves the accuracy of the valence ones. A special characteristic of the present analysis is the calculation of the second moments of the excited-state orbitals. Because we find that the CASPT2 densities agree well with the CASSCF ones and since the MC-PDFT methods gets accurate excitation energies based on the CASSCF densities, we believe that we can trust these moments as far as giving a more accurate picture of the diffuseness of the excited-state orbitals in these prototype molecules than has previously been available. 
    more » « less
  5. Pursuing ever-smaller feature size in laser-based lithography is a research topic of vital importance to keep this technique competitive with other micro-/nano-fabrication methods. Features smaller than the diffraction-limited spot size can be obtained by “thresholding”, which utilizes the deterministic nature of damage threshold with ultrashort laser pulses and is achieved by precisely tuning pulse energies so that only the central portion of the focal spot produces permanent modification. In this paper, we examine the formulation commonly used to describe thresholding and show that the relationship between feature size (r) and laser fluence (F) is invariant with respect to the nature of laser absorption. Verified by our experiments performed on metal, semiconductor, and dielectric samples, such invariance is used to predict the smallest feature size that can be achieved for different materials in a real-world system. 
    more » « less