M. Holst

C. T. Lee, J. B. Moody, R. E. Amaro et al.

We introduce CASC: a new, modern, and header-only C++ library which provides a data structure to represent arbitrary dimension abstract simplicial complexes (ASC) with user-defined classes stored directly on the simplices at each dimension. This is accomplished by using the latest C++ language features including variadic template parameters introduced in C++11 and automatic function return type deduction from C++14. Effectively CASC decouples the representation of the topology from the interactions of user data. We present the innovations and design principles of the data structure and related algorithms. This includes a metadata aware decimation algorithm which is general for collapsing simplices of any dimension. We also present an example application of this library to represent an orientable surface mesh.

F. Grogan, M. Holst, L. Lindblom et al.

Molecular dynamics (MD) simulations are used in biochemistry, physics, and other fields to study the motions, thermodynamic properties, and the interactions between molecules. Computational limitations and the complexity of these problems, however, create the need for approximations to the standard MD methods and for uncertainty quantification and reliability assessment of those approximations. In this paper, we exploit the intrinsic two-scale nature of MD to construct a class of large-scale dynamics approximations. The reliability of these methods are evaluated here by measuring the differences between full, classical MD simulations and those based on these large-scale approximations. Molecular dynamics evolutions are non-linear and chaotic, so the complete details of molecular evolutions cannot be accurately predicted even using full, classical MD simulations. This paper provides numerical results that demonstrate the existence of computationally efficient large-scale MD approximations which accurately model certain large-scale properties of the molecules: the energy, the linear and angular momenta, and other macroscopic features of molecular motions.

M. Holst, V. Kungurtsev

The Einstein constraint equations have been the subject of study for more than fifty years. The introduction of the conformal method in the 1970's as a parameterization of initial data for the Einstein equations led to increased interest in the development of a complete solution theory for the constraints, with the theory for constant mean curvature (CMC) spatial slices and closed manifolds completely developed by 1995. The first general non-CMC existence result was establish by Holst et al. in 2008, with extensions to rough data by Holst et al. in 2009, and to vacuum spacetimes by Maxwell in 2009. The non-CMC theory remains mostly open; moreover, recent work of Maxwell on specific symmetry models sheds light on fundamental non-uniqueness problems with the conformal method as a parameterization in non-CMC settings. In parallel with these mathematical developments, computational physicists have uncovered surprising behavior in numerical solutions to the extended conformal thin sandwich formulation of the Einstein constraints. In particular, numerical evidence suggests the existence of multiple solutions with a quadratic fold, and a recent analysis of a simplified model supports this conclusion. In this article, we examine this apparent bifurcation phenomena in a methodical way, using modern techniques in bifurcation theory and in numerical homotopy methods. We first review the evidence for the presence of bifurcation in the Hamiltonian constraint in the time-symmetric case. We give a brief introduction to the mathematical framework for analyzing bifurcation phenomena, and then develop the main ideas behind the construction of numerical homotopy, or path-following, methods in the analysis of bifurcation phenomena. We then apply the continuation software package AUTO to this problem, and verify the presence of the fold with homotopy-based numerical methods.

M. Holst, M. Licht

We present a new technique to apply finite element methods to partial differential equations over curved domains. A change of variables along a coordinate transformation satisfying only low regularity assumptions can translate a Poisson problem over a curved physical domain to a Poisson problem over a polyhedral parametric domain. This greatly simplifies both the geometric setting and the practical implementation, at the cost of having globally rough non-trivial coefficients and data in the parametric Poisson problem. Our main result is that a recently developed broken Bramble-Hilbert lemma is key in harnessing regularity in the physical problem to prove higher-order finite element convergence rates for the parametric problem. Numerical experiments are given which confirm the predictions of our theory.