Source author record

Maxwell Hutchinson

Maxwell Hutchinson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.comp-ph cond-mat.mtrl-sci Discrete Mathematics Machine Learning math-ph math.CO math.MP physics.flu-dyn

Catalog footprint

What is connected

5works

8topics

3close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Multivariate Prediction Intervals for Random Forests

Accurate uncertainty estimates can significantly improve the performance of iterative design of experiments, as in Sequential and Reinforcement learning. For many such problems in engineering and the physical sciences, the design task depends on multiple correlated model outputs as objectives and/or constraints. To better solve these problems, we propose a recalibrated bootstrap method to generate multivariate prediction intervals for bagged models and show that it is well-calibrated. We apply the recalibrated bootstrap to a simulated sequential learning problem with multiple objectives and show that it leads to a marked decrease in the number of iterations required to find a satisfactory candidate. This indicates that the recalibrated bootstrap could be a valuable tool for practitioners using machine learning to optimize systems with multiple competing targets.

preprint2016arXiv

Direct numerical simulation of single mode three-dimensional Rayleigh-Taylor experiments

The single-mode Rayleigh-Taylor instability (smRTI) is well defined, poorly understood, and applicable to many fluid flows directly and through its relationship to multi-mode Rayleigh-Taylor models. This study reproduces three low-Atwood smRTI experimental runs (Wilkinson and Jacobs, 2007) in a specialized version of the Nek5000 spectral element code. The simulations use the initial amplitude, wavelength, acceleration, Atwood number, and viscosity from the three specific experiments and impose no-slip and no-flux boundaries on the velocity and scalar, respectively. The simulations are shown to reproduce the linear, saturation, stagnation, and re-acceleration phases of the smRTI seen in the experiments. Additionally, access to the full velocity and scalar fields demonstrates three different finite size effects: wall drag, wall lift, and a long wavelength mode along the diagonal. One of the simulations is extended by a factor of two in the vertical direction and the resulting late-time dynamics reach Froude numbers around 1.8, higher than previously reported. Finally, inspection of the span-wise flow reveals secondary flows of the first kind that transport the scalar from the bubble-spike interfaces into the bubble and spike centers. The agreement between simulations and experiments inspires confidence in the spectral element method for studying the Rayleigh-Taylor instability.

preprint2015arXiv

Enumeration of octagonal tilings

Random tilings are interesting as idealizations of atomistic models of quasicrystals and for their connection to problems in combinatorics and algorithms. Of particular interest is the tiling entropy density, which measures the relation of the number of distinct tilings to the number of constituent tiles. Tilings by squares and 45 degree rhombi receive special attention as presumably the simplest model that has not yet been solved exactly in the thermodynamic limit. However, an exact enumeration formula can be evaluated for tilings in finite regions with fixed boundaries. We implement this algorithm in an efficient manner, enabling the investigation of larger regions of parameter space than previously were possible. Our new results appear to yield monotone increasing and decreasing lower and upper bounds on the fixed boundary entropy density that converge toward S = 0.36021(3).

preprint2014arXiv

The Shirley reduced basis: a reduced order model for plane-wave DFT

The Shirley reduced basis (SRB) represents the periodic parts of Bloch functions as linear combi- nations of eigenvectors taken from a coarse sample of the Brillouin zone, orthogonalized and reduced through proper orthogonal decomposition. We describe a novel transformation of the self-consistent density functional theory eigenproblem from a plane-wave basis with ultra-soft pseudopotentials to the SRB that is independent of the k-point. In particular, the number of operations over the space of plane-waves is independent of the number of k-points. The parameter space of the transformation is explored and suitable defaults are proposed. The SRB is shown to converge to the plane-wave solution. For reduced dimensional systems, reductions in computational cost, compared to the plane-wave calculations, exceed 5x. Performance on bulk systems improves by 1.67x in molecular dynamics-like contexts. This robust technique is well-suited to efficient study of systems with strin- gent requirements on numerical accuracy related to subtle details in the electronic band structure, such as topological insulators, Dirac semi-metals, metal surfaces and nanostructures, and charge transfer at interfaces with any of these systems. The techniques used to achieve a k-independent transformation could be applied to other computationally expensive matrix elements, such as those found in density functional perturbation theory and many-body perturbation theory.

preprint2011arXiv

VASP on a GPU: application to exact-exchange calculations of the stability of elemental boron

General purpose graphical processing units (GPU's) offer high processing speeds for certain classes of highly parallelizable computations, such as matrix operations and Fourier transforms, that lie at the heart of first-principles electronic structure calculations. Inclusion of exact-exchange increases the cost of density functional theory by orders of magnitude, motivating the use of GPU's. Porting the widely used electronic density functional code VASP to run on a GPU results in a 5-20 fold performance boost of exact-exchange compared with a traditional CPU. We analyze performance bottlenecks and discuss classes of problems that will benefit from the GPU. As an illustration of the capabilities of this implementation, we calculate the lattice stability α- and β-rhombohedral boron structures utilizing exact-exchange. Our results confirm the energetic preference for symmetry-breaking partial occupation of the β-rhombohedral structure at low temperatures, but does not resolve the stability of α relative to β.