CMU Mathematics Undergraduate Research, 2026
In early 2020, as the global coronavirus pandemic brought the world to a screeching halt, many of our students found that their summer programs had been cancelled. Many of these students were very interested in using this suddenly free summer to engage with open ended-projects in mathematics. In response to this demand, the Department of Mathematical Sciences curated a brand new summer program: Summer Experiences in Mathematical Sciences (SEMS). This program is offered in parallel with the traditional SURF program, which provides stipends to students who stay at Ä¢¹½Av in the summer months to work on research projects. Outlined below are the projects undertaken during the summer of 2026 as part of the SURF and SEMS programs.
▼ SURF projects
► CAILA, the Calculus-Aligned Interactive Learning Assistant
Bjorn Beck, Elliot Li, Eduardo Chen, Lucy Wilson-Sadlowski, Daniel Sium
Advisor: Jay Howell
Abstract: The CAILA project developed and tested a prototype AI learning assistant for Calculus I intended to support students through guided, office-hours-style conversations.
Using customized Google Gemini Gems, the project incorporated course-specific materials, Socratic questioning, pedagogical guardrails, and a confidence-aware scaffolding framework that adjusted the level of support during a dialogue. Multiple versions were evaluated using structured student personas, standard calculus prompts, and a rubric focused on mathematical correctness, instructional quality, scope compliance, and adaptive behavior. The project demonstrated both the potential of customized LLMs for calculus support and the limitations of commercial custom-LLM platforms, informing ongoing work toward a standalone, API-based implementation.
► Continuous equivariant learning with linearly many frames
Nathaniel Kenschaft
Advisor: Florian Frick
Abstract: Many AI systems in science predict physical quantities, such as the forces between atoms in a molecule, and those predictions must respect symmetry: rotating the molecule should rotate the forces. A popular approach first turns every input into a standard orientation, runs an ordinary neural network, and turns the answer back. But a standard orientation cannot generally vary smoothly across all inputs, and abrupt switches can introduce artificial jumps in predicted forces. Like covering the Earth with overlapping maps, one can instead use several orientations and blend their predictions with gradually shifting weights. Each orientation, however, costs another run of the network.
Using topology, the mathematics of shape, this work shows that the number needed grows in proportion to the number of particles, and gives an explicit recipe that achieves this scaling, for example, 61 candidate orientations for 24 particles, compared with more than 1,000 in previous constructions.
► Replacement principles for vector valued random sums
Steven Ma
Advisor: Tomasz Tkocz
Abstract: A nonstandard tool for establishing precise Gaussian tail bounds for sums of independent random variables is a robust replacement principle based on convex stochastic domination. The classical theory concerns only scalar Rademacher sums. Building on recent work on unimodal sums, we extend such replacement principles to the vector-valued setting. As a corollary, we obtain sharp tail bounds for spherically symmetric sums.
► Bore wave solutions to the viscous Saint-Venant system with capillarity
Keming Meng
Advisor: Ian Tice
Abstract: Bore waves are a special type of traveling wave observed on the surfaces of rivers and canals, characterized by distinct upstream and downstream limits. In this project we studied the existence of bore waves within the context of the viscous Saint-Venant system with incline, slip, and capillarity. We proved the existence of bore waves in this model using techniques from dynamical systems.
► Cycle-backbone random graphs
Wilson Pan
Advisor: Wes Pegden
Abstract: Fix a cycle $C_k$ and, for each of its vertices, associate a label set by including label $i$ from $\{1,...,n\}$ in the label set with probability $p$. This generates a random graph with correlations among edges. How does this change basic properties of the resulting graph?
We show, for example, that one can use this model to generate random graphs that are likely to be disconnected even though they are unlikely to have any isolated vertices, in contrast to the uncorrelated Erdos-Renyi model. On the other hand, if $k$ is large relative to $n$, we show that indeed this does converge in distribution to the uncorrelated case, allowing one to generate random graphs with a single model that interpolate between low and high correlation outputs.
► Gradient-flow methods for sampling and PDE approximation
Kausthub Satluri and Jerry Yin
Advisor:
Abstract: Gradient flows provide a common framework for studying evolution partial differential equations (PDEs) and developing sampling algorithms. This project investigated how regularizing the Wasserstein geometry on the space of probability measures changes the dynamics of Kullback–Leibler (KL) gradient flows, and whether these changes can lead to effective algorithms for sampling and PDE approximation. We implemented several gradient-flow algorithms and evaluated their performance, including comparisons with the “blob” method. In the examples tested, the numerical results suggest that algorithms based on the proposed flows may offer computational advantages over the ``blob" method.

► Interpolation and fractional Sobolev embeddings
Asher Wheaton
Advisors: Giovanni Leoni and Ian Tice
Abstract: Sobolev spaces measure a function's size along with the size of its derivatives, and embedding theorems say how much smoothness can be traded for integrability. Abstract interpolation takes two such spaces and builds a continuous family of intermediate ones between them, so a result proved at the two extremes yields estimates for all the intermediate spaces.
The project had two parts. In the first, we studied a classical embedding between fractional Sobolev spaces, for which Nguyen, Díaz and Nguyen recently found a simple proof using only the intrinsic Gagliardo seminorm; their argument does not reach the critical case, and we gave a short proof of that case by combining their approach with classical interpolation. In the second, we identified the intermediate spaces between the homogeneous Sobolev spaces of equal order $\dot W^{k,1}$ and $\dot W^{k,\infty}$, and showed they are the homogeneous Sobolev spaces $\dot W^{k,p}$.
► Gaussian volume of two dimensional sections of the cube
Emiliano Zetune
Advisor: Tomasz Tkocz
Abstract: The study of extremal-volume sections of convex bodies is a classical topic in convex geometry, with roots in the geometry of numbers. These problems are well understood for certain classes of convex bodies. We investigate analogous problems in Gaussian space, that is, Euclidean space equipped with the standard Gaussian measure. As a first step, we establish new sharp results on maximal planar sections.
▼ SEMS projects
► The $k$-Color Poset Ramsey Number of the Diamond is $2k$
Will Cong, Teresa Huang, John Li, and BangTam Ngo
Advisors: Tom Bohman and Quentin Dubroff
Abstract: Ramsey theory focuses on the appearance of order in a substructure given a structure of a specific size. The $k$-color Ramsey number of a poset $P$, a partially ordered set where some pairs of elements may not be comparable, is the smallest dimension $d$ such that any $k$-coloring of the $d$-dimensional Boolean lattice contains a monochromatic induced copy of $P$, our substructure. Our structure is the collection of all subsets of a set with $d$ elements, ordered by inclusion. Furthermore, an induced copy preserves exactly the ordering relations of $P$.
For $j\geq 2$, the $j$-diamond is the poset consisting of a maximum element, a minimum element, and an antichain of size $j$, a collection of $j$ elements that are pairwise incomparable, between them.

We prove that the $k$-color Ramsey number of the $j$-diamond is at most $jk$. Colorings of the Boolean lattice that assign sets of the same size the same color provide matching lower bounds for $j\in{2,3,4}$, showing that this bound is exact in these cases.
► New Bounds in Multicolor Ramsey Numbers of Hypergraphs
Will Schremmer and Lillian Xia
Advisors: Tom Bohman and Quentin Dubroff
Abstract: Traditional Ramsey Theory is the study of finding monochromatic subgraphs in an arbitrary coloring of complete graphs - graphs that have every possible edge given a set of vertices. We in particular looked at multicolor Ramsey numbers for uniform hypergraphs. An $r$-uniform hypergraph is a generalization of a graph where each edge connects $r$ vertices at once, rather than the standard two vertices per edge. Given an $r$-uniform hypergraph $H$, we define $R(k; H)$ as the Ramsey number of coloring the complete graph with $k$ colors and looking for a monochromatic $H$. Our main goal was to investigate the behavior of $R(k; H)$ for fixed $H$ and $k \to \infty$. We used Keevash's theorem on the existence of designs to find lower bounds of $R(k; H)$ for infinitely many values of $k$ for every hypergraph $H$, and developed a weighting lemma to find upper bounds. We ended up with matching upper and lower bounds for $R(k; H)$ for a large class of hypergraphs $H$.

► Random walks with jumps
Jack Xu and Noah Kim
Advisor: Gautam Iyer
Abstract: Consider a random walk alternated with deterministic jumps. Although these jumps typically accelerate mixing, they can occasionally have the opposite effect. We investigate the exceptional cases in which deterministic jumps slow down mixing, and reduce the rate of convergence to the stationary distribution.

► On the Banach–Mazur distance to the three-dimensional cube
Leo Zhuang and Clement Zhou
Advisor: Konstantin Tikhomirov
Abstract: We use computational methods that combine results from convex geometry with optimization algorithms to improve the best previously known uniform bound on the Banach–Mazur distance to the three-dimensional cube.

Specifically, we show that for every origin-symmetric 3-dimensional convex body K, its Banach–Mazur distance to the cube is at most 2.45, improving upon the Taschuk bound.

2026 SURF team
Students: Bjorn Beck, Eduardo Chen, Nathaniel Kenschaft, Elliot Li, Steven Ma, Keming Meng, Wilson Pan, Kausthub Satluri, Daniel Sium, Asher Wheaton, Lucy Wilson-Sadlowski, Jerry Yin, Emiliano Zetune
Faculty: Florian Frick, Elias Hess-Childs, Jay Howell, Giovanni Leoni, Wes Pegden, Ian Tice, Tomasz Tkocz
2026 SEMS research team
Students: Will Cong, Teresa Huang, Noah Kim, John Li, BangTam Ngo, Will Schremmer, Lillian Xia, Jack Xu, Clement Zhou, Leo Zhuang
Faculty: Tom Bohman, Quentin Dubroff, Gautam Iyer, Konstantin Tikhomirov