Fall 2026 IML Projects

*Projects may be added to this page after applications open.*

Illinois Mathematics Lab provides a framework for faculty and graduate students to engage University of Illinois undergraduates in projects related to mathematical research. Projects are one semester long, during fall and spring semesters.

Applications for Fall 2026 will open by August 1 and will be due on Friday, August 14. Applicants must be current undergraduates at the University of Illinois Urbana-Champaign. In the application, you will need to answer some questions and upload an unofficial transcript. To receive email notifications about due dates and about newly added projects, join our mailing list.


Adversarial Quantum Machine Learning

A flow chart with boxes for datasets, adversarial attacks, quantum machine learning model, evaluate robustness, and new vulnerabilities.

Faculty member: Theshani Nuradha
Level: Intermediate
Keywords: data science, linear algebra, statistics
Course prerequisites: MATH 416 or an equivalent, Course on machine learning fundamentals
Coding/software prerequisites: Python, Qiskit. Experience in machine learning projects would be useful.

The rapid development of quantum technologies has created growing interest in quantum models for data analysis. In particular, quantum machine learning combines parameterized quantum circuits with classical optimization methods. These models are being explored for potential applications in areas such as healthcare, finance, and scientific data analysis. This project investigates the robustness of quantum machine learning models against adversarial attacks. Adversarial attacks include input perturbations, interference to optimization process, among others. A central objective is to determine whether improving robustness against one type of attack can unintentionally increase vulnerability to another. We will combine theoretical analysis with numerical experiments to compare different quantum and hybrid quantum-classical models under multiple threat scenarios. More broadly, the study will contribute to a better understanding of adversarial quantum machine learning and the limitations that should be considered before deploying quantum models in real-world applications.


Decidability of Variations of the Monotiling Problem

Contiguous rectangles that appear to tile the plane, each with a colored dot in the center.

Faculty member: Nicolas Chavarria Gomez
Level: Intermediate
Course prerequisites: MATH 417 is recommended
Coding/software prerequisites: none
Keywords: Tilings, decidability, algorithms, abelian groups

Given a set of planar tiles and some rules on how to put them together, can the entirety of the plane be covered by non-overlapping copies of these tiles? This is the problem known as tiling the plane, whose geometrically intuitive nature has been used to understand certain complexity problems. An early example of this idea comes from Hao Wang, who used a specific kind of tiling to understand the complexity of ∀∃∀\forall\exists\forall formulas in first-order logic. Also known as the domino problem, the tilings used by Wang consisted of squares with colors on each side and a rule mandating that two tiles could only be put next to each other if their overlapping sides matched in color. Wang believed that, whenever there was a tiling of the plane by a tileset of this kind, there was a periodic tiling. From this he deduced that there was an algorithm for deciding, given a tileset, whether a tiling existed or not. However, his student Robert Berger proved that any Turing machine could be coded as a Wang tileset, thus disproving the decidability of the problem and implying the existence of aperiodic tilings to boot!

Similar decidability problems have interested mathematicians since Berger proved the undecidability of the domino problem. One such problem is the monotiling problem. Namely, given a finite subset TT of ℤ2\mathbb{Z} ^2 (the monotile), can all of ℤ2\mathbb{Z}^2 be covered by non-intersecting translates of TT? Siddhartha Bhattacharya proved in 2016 that the monotiling problem is decidable. In other words, there is an algorithm which, when fed a monotile, effectively computes whether a tiling by it exists or not. On the other hand, Terence Tao and Rachel Greenfeld proved in 2023 that this is no longer true in higher dimensions. That is to say, no algorithm exists for tilings of ℤd\mathbb{Z}^d where d>2.d>2. Their proof goes through tilings of groups of the form ℤ2×G,\mathbb{Z}^2\cross G,where GG is finitely generated, but says nothing of the decidability of the monotiling problem for any specific GG. In this project we will attempt to bridge this gap by focusing on the decidability of tilings of, for example, ℤ×ℤ/nℤ\mathbb{Z} \cross \mathbb{Z}/n\mathbb{Z}. We will also try to understand the decidability of other variations of this problem, such as parallel tilings (i.e. tiling the plane with two monotiles using the same tiling pattern).


Dynamic Lie Algebras in Quantum Information Theory

Faculty member: Marius Junge
Level: Intermediate
Keywords: quantum information theory, Lie algebras
Course prerequisites: Linear algebra such as Math 416 (but you have to like it beyond matrix manipulations). Quantum mechanics for physics student is a welcome addition.

Quantum Information Theory is very good in borrowing mathematical tools from many areas. In this project we will review basic concepts of quantum computation and learn what magic gates are.  Then we will focus on the special role of Pauli matrices and Ising Hamiltonians (search them up).  We will end up using tools from Lie algebra theory, a good opportunity to understand some of the basics. The goal is to find an algorithmic approach for the approximation of quantum gates-which constitute quantum computers. 


Exploring I-Functions of Calabi-Yau Orbifolds/Manifolds

Faculty member: Deniz Genlik
Level: Advanced
Keywords: Algebra, Differential Equations, Mathematical Physics
Course prerequisites: An ODE course, an analysis course, an abstract algebra course.
Coding/software prerequisites: SageMath

Calabi-Yau manifolds are fascinating geometric structures at the heart of mirror symmetry, a duality predicted by string theory that connects the algebraic and analytic worlds. In genus zero, the analytic side of mirror symmetry is captured by I-functions, special functions that encode deep geometric information. This project is a continuation of an IML project that has started in Spring 2025, with new participants for Fall 2026. Students will go first go over what the last group members accomplished and build on that to get new results for the I-functions of certain Calabi-Yau orbifolds/manifolds.


Fitting a Potential Model of Dark Energy

A diagram of big bang expansion of the universe over 13.7 billion years.

Faculty member: Kay L. Kirkpatrick
Level: Intermediate
Keywords: Mathematical physics (Dark Matter and Dark Energy), PDEs and probability, scientific computation
Course prerequisites: Exposure to probability, PDEs, dark matter, and quantum mechanics
Coding/software prerequisites: The project will likely use Python, maybe Matlab


We will use an idea about a mechanism that converts dark matter to dark energy, a nonlinear quantum equation with a fractional Laplacian, and do computations and estimates for fitting this mechanism with the Hubble tension.


The Future of Math Research

Faculty member: Jacob Beckey
Level: Introductory
Keywords: Lean, Formalization, Education
Course prerequisites: There are not formal prerequisites beyond Calc III for this project. However, an ideal applicant would have a strong interest in examining the interplay between research mathematics, LLMs, interactive theorem provers, and math education.
Coding/software prerequisites: Some coding experience would be preferred, as we will be learning and then using the Lean programming language for formalizing proofs.

Large language models (LLMs) have matured to the point of being able to resolve conjectures that have been open for decades in fields ranging from combinatorics to quantum information theory, and many in between. As this incredible technology continues to evolve, it is crucial that mathematical research evolves with it. In this project, we will explore what the future of math research may look like: proof generation, formal verification using interactive theorem provers, and then broad communication/exposition facilitated by generative AI. The goal will be to provide an example of what high-quality mathematical research projects should look like in the age of AI.


Hybrid Multinomials and Their Combinatorics

An integral equation related to hybrid monomials.

Faculty member: Deniz Genlik

**This is a continuation of the summer 2026 Hoover project and will not be accepting new members.

Recent work of Cavalieri, Gillespie, and Monin introduced a family of “asymmetric multinomials,” which arise in connection with the geometry of moduli spaces of stable rational curves with n-marked points. In this project, we will define and study a new family of multinomial-type structures that interpolates between the classical symmetric multinomials and these asymmetric versions. Possible directions include computing examples, finding recursions or generating functions, studying coefficient patterns, and seeking combinatorial interpretations. If time permits, we will also explore possible connections with parking functions and related combinatorial structures.


In-game Win Probabilities

A two dimensional graph with time in the game on the horizontal axis and win probability on the vertical axis.

Faculty member: AJ Hildebrand
Level: Intermediate( theory) / Advanced (coding)
Keywords: Sports Analytics, Probabilistic Forecasting, Regression Models
Course prerequisites: None beyond Calculus III
Coding prerequisites: Due to the scale and complexity of the data involved, a high level of proficiency in Python is essential. Prior experience with large scale coding projects in Python and/or with webscraping libraries such as BeautifulSoup and Selenium would be ideal. If you have a Github site, please mention it in your application.

In recent years, sports websites such as ESPN have been posting charts of real-time in-game win probabilities such as the one shown in the accompanying image for Game 3 of the 2026 NBA Finals between the San Antonio Spurs and the New York Knicks. The chart shows that the Spurs had a 23.4% win probability at the half-time mark, but ended up winning the game. The chart for a single basketball game can contain hundreds of such individual data points, with each data point representing a win probability at a particular time point in the game. Taken together, the available archived win probability charts for games in major sports leagues such as the NBA represent an exceptionally rich trove of data containing millions of individual probabilistic predictions.

The goal of this project is to collect this data and analyze the accuracy of the predictions implied by the win probability charts from a variety of angles. In particular, we seek to answer questions such as the following: How reliable are the in-game win probabilities in terms of predicting the ultimate outcome of a game? To what extent does the reliability change/improve over the course of a game? Are there significant differences between sports in terms of the reliability of their in-game win probabilities? How do ESPN’s in-game win probabilities compare to probabilistic forecasts based on other metrics such as margin of lead in terms of their predictive power on the outcome of the game?


Optimization and Bayesian Inference in Practice: Applications to the Environmental Sciences

Two maps of the western United States with colors indicating fire locations (1984-2020) and prior distribution.

Faculty member: Gabriel Provencher Langlois
Level: Intermediate
Keywords: optimization, Bayesian inference, data science, machine learning, environmental sciences
Course prerequisites: MATH 415/416 (Linear Algebra), MATH 285/441 (Ordinary differential equations), MATH 461/STAT 400 (Probability and Statistics I).
Coding/software prerequisites: Python, MATLAB, or any other scientific computing programming language.

In this project, the students will apply optimization and Bayesian inference methods to study, analyze, and understand real environmental and weather data. Possible examples of concrete problems include, but are not limited to: ensemble forecasting of wildfire activity in the Western United States; prediction of wildfire spread using a blend of machine learning and PDE modelling; statistical inference of temperature and wind effects on electrical grids such as ERCOT. Students will implement state-of-the-art optimization methods (e.g., accelerated first-order methods, ADAM, Adagrad…) and Bayesian inference methods (e.g., Hamiltonian and Least-Squares Monte Carlo methods, Gaussian Processes, Bootstrap…) in Python or MATLAB and test them on public datasets.


Predictability in Major Sports Leagues

A chart of various statistics for teams in the Premier League.

Faculty member: AJ Hildebrand
Level: Intermediate
Keywords: Sports Analytics, Graph Theory, Combinatorics

Course prerequisites: Aside from completion of the calculus sequence, there are no hard course prerequisites. The background in graph theory necessary for the theoretical component of the project can be acquired during the course of the project, but if you have taken courses in graph theory or combinatorics such as Math 412 or Math 413, please mention this in your application.

Coding/software prerequisites: We will use Python as primary coding language, possibly supplemented by Mathematica or a similar tool for visualizations. Participants should be comfortable with at least one of these languages. If you have a Github site, please mention it in your application.

If we are given the end-of-season ranking of all teams in a sports league, how well does this ranking “predict” the results of games played within the season? In an ideal scenario in which games are completely predictable a team should win all of its games against lower-ranked opponents and lose all of its games against higher-ranked opponents. In the real world, this situation rarely happens. For example, in the 2025-2026 English Premier League season, third-ranked Manchester United lost five games against lower-ranked opponents, but also won a game each against the first and second place teams, Arsenal and Manchester City. Thus, seven of the 38 games played by Manchester United that season resulted in upsets. The frequency of such upsets within a given league is a natural measure for the randomness and (un)predictability inherent in the league. In this project we will investigate the predictability of major professional sports leagues such as the Premier League and Major League Baseball from this angle. We will also construct appropriate graph-theoretic models of those leagues, study the occurrence of upsets within these models, and compare the theoretical results from such “model leagues” to those observed in the actual sports leagues.


Real Conic Bundles: an Experimental Approach

An image of the Cartesian plan with curves dividing it into several regions, and four straight lines on the plane.

Faculty member: Lena Ji
Level: Intermediate
Keywords: Algebra, algebraic geometry, polynomials
Course prerequisites: Abstract algebra (ring theory), Linear algebra, Basic point set topology (e.g., Math 432 or 447)
Coding/software prerequisites: Basic familiarity with coding (e.g., Python) is strongly preferred. I expect that this project will use Magma or SageMath.

Algebraic varieties are geometric objects defined by polynomial equations. Classically, algebraic varieties have been studied the most over the complex numbers, where the Fundamental Theorem of Algebra guarantees the existence of such solutions. If the polynomials defining an algebraic variety have coefficients in a subfield of ℂ\mathbb C, it is interesting to study the properties of the variety over that field. In this project, we will study a certain class of varieties (conic bundles) over the real numbers, with a focus on computing explicit examples.


Sparse Small-World Networks: Geometry and First-Passage Percolation

A poster titled "Weighted Small-World Cycles with Long-Range Shortcuts: Interplay of Geometry and Randomness."

Faculty members: Partha Dey and Neeladri Maitra
Level: Advanced
Keywords: Random networks, graph distance, branching processes
Course prerequisites: Probability Theory, Graph Theory, Linear Algebra, Algorithms, Data Structures
Coding/software prerequisites: Python required, C+/Julia if needed

We consider a one-dimensional discrete cycle of length nn, in which each vertex is connected to its two nearest neighbors. In addition, we introduce random shortcut edges between pairs of vertices whose graph distance along the cycle lies between 2 and l l. The parameter, l=lnl=l_n grows with n n but edges are rare with probability p=a/lp=a/l so the mean random degree is 2a and the graph remains sparse with typical degree 2+Poisson(2a). Such small-world graphs interpolate between a simple cycle and a random graph. In the classic Newman–Watts model where l=n/2l=n/2, it is known that the graph diameter and typical distances grow logarithmically in nn. When random exponential edge weights are added, the typical weighted distance between random vertices is asymptotically (log⁡n)/λ(\log n)/\lambda for some λ>0\lambda>0. These results illustrate that sparse cycles with random long edges exhibit small-world scaling (distance ~ log n) rather than the linear distances of a bare cycle. We assign independent random weights ωnn\omega^{nn} to nearest neighbor edges and independent random weights ωsc\omega^{sc} to the random shortcut edges. There are two competing growth mechanisms: 1) Geometric propagation along the underlying cycle, and 2) Branching propagation through random shortcut edges. This project will explore the interplay between the underlying geometry of the cycle and the added randomness. In particular, we will study how the typical graph distance and weighted first-passage distance from the origin to a uniformly chosen vertex depend on the shortcut range l, the shortcut density gamma, and the distribution of the nearest-neighbor random weights ωnn\omega^{nn} and short-cut random weights ωsc\omega^{sc}. We will focus on the small-world regime l=c∗nl=c*n and the crossover regime l=c∗n/log⁡nl=c*n/\log n for some c>0c>0, and investigate the following goals: a) literature survey, b) local limit analysis, c) computational experiments via simulation, d) regime classification, e) graph distance and weighted distance analysis. This project offers a balanced combination of rigorous mathematics, probabilistic modeling, algorithm development, and computational experimentation. It is suitable for a team of strong undergraduate students interested in probability, networks, and data-driven mathematical research.


Symmetric Fixed Points of the Ricci Flow

Mathematical images of Ricci flow evolution of a sphere, a torus, and a genus two surface.

Faculty member: Eric Chen
Level: Intermediate/Advanced
Keywords: differential geometry, ODE/differential equations, asymptotic analysis
Course prerequisites:
required: MATH 441, MATH 444
highly desirable but not strictly required: MATH 423
not required but could allow for more directions: MATH 442, MATH 450
Coding prerequisites: No prerequisites but knowledge/experience with Mathematica/MATLAB or a related program would be helpful.

The sphere, the flat plane, and the hyperbolic plane are examples of 2D surfaces. The Ricci flow takes a curved 2D surface and evolves it (changes its shape), and in many cases this evolution will converge to a limiting 2D surface. These fixed points are often highly symmetric, like the examples mentioned at the beginning. The Ricci flow can also evolve higher-dimensional surfaces and was an important tool used to prove the Poincaré conjecture in 2003. We will study properties of these fixed points (solitons) in these higher dimensions assuming a high degree of symmetry. This reduces the problem to understanding the behavior of solutions to systems of ordinary differential equations. Our goal is to understand properties of these solutions, including asymptotics and possibly existence of such solutions.


The Wythoff Array: continued fractions, congruences, fractional parts, botany(phyllotaxis), and more

An array of numbers with first column 1, 4, 6, , 12, 14, ...., with second column 2, 7, 10, 15, 20, 23, ..., and with each subsequent column equal to the sum of the two columns immediate to its left.

Faculty member: Kenneth Stolarsky
Level: Intermediate/advanced
Course prerequisites: Math 314 or 347 required. Number theory and abstract algebra are desirable. Other 400 level math courses in pure math are a plus.
Coding prerequisites: proficient in either Python or Mathematica

This continues previous IML groups that studied the Wythoff Array ( WA ). Newcomers are welcome! The WA is a particular infinite matrix of positive integers such that every row satisfies the Fibonacci recurrence, and every positive integer occurs in it exactly once. The first row is the Fibonacci sequence itself, and here ratios of consecutive entries have continued fractions in which every partial quotient is 1. Can one describe in detail the form for the other rows?

There are many questions about congruences involving the WA. For example, given the first five entries of any row, must two of them be congruent modulo 5 ?

A problem of general interest in Diophantine Approximation is to study for (1) any increasing sequence of positive integers {s(n)}\{s(n)\}, and (2) any a>0a>0 ,  the sequence of fractional parts of the sequence {s(n)a}\{s(n)a\} , denoted by {frac[{s(n)a}]}\{ frac [ \{ s ( n ) a \} ] \} . For example, when a=2a = \sqrt{2} and s(n)=ns ( n ) = n , it is well known that the sequence is uniformly distributed in the unit interval. Our interest is for {s(n)}\{s(n)\} a row of WA and a=ga = g , the golden ratio. Much is already known here, but the sequence {frac[s(n)s(n+1)g]}\{ frac [ s ( n ) s ( n + 1 ) g ] \} seems unexamined. For example, in the sixth row the 7th and 8th entries are 254 and 411 respectively, with frac[254411g]=.0402215….frac [ 254 411 g ] = .0402215….  . For the 17th and 18th entries the number is .040325217….   . The latter agrees with a root of p(x)=1−25x+5x2p ( x ) = 1 – 25 x + 5 x ^ 2  to about 9 decimal places. Also, the number for the 8th and 9th entries is .959635…. . Some conjectures suggest themselves!

The first to write about the geometry of leaf arrangements (phyllotaxis) seems to be J. Kepler (1611) . He connected these biological forms with Fibonacci numbers! In recent years, the late John Conway ( famous for, inter alia, devising the ” Game of Life”, a (very simplified) mathematical model of biological behavior) left behind some very rough notes that try to explain the above Fibonacci connection. It led him to the WA, and to some correspondence with Clark Kimberling , perhaps the leading expert on the WA. These items can be found on the internet, and are worthy of some attention.(An excellent examination of current mathematical models of phyllotaxis ( which, however, does not mention the WA) has recently appeared. The reference is Emergence of the Golden Angle in Flower Heads, Przemyslaw Prusinkiewicz, Amer Math Monthly, 2026 (133) , 303 – 314 .)

Participants would also be very welcome to ask their own questions about the WA, and to pursue them.


Illinois Mathematics Lab
Altgeld Hall
1409 W. Green Street
Urbana, IL 61801
Email: math-iml@illinois.edu
Log In