Skip to content

Projects

My research centres on developing new methods, implementing new models, and enhancing workflows to bridge the scales between theory and experiment. Research topics focus on pressing challenges in catalysis, energy, materials, and the natural world. I deepen our understanding of processes such as chemical reactions, the effects of microstructure, nucleation, and degradation.

20
Publications
344
Citations
9
h-index
Metrics from Google Scholar.
Last updated: 2026-09-24

The future of modeling experiments

◀ Past
Future ▶
Complex polycrystalline barium zirconate model
Idealized periodic barium zirconate model
‹ ›

Project Highlights

4 highlights
Simulated Experiment

Simulated Experiment

Computational material synthesis

Efficient proton transport is essential for advancing solid-state fuel cells and electrolyzers, yet the mechanisms governing proton diffusion in polycrystalline barium zirconate (BaZrO3) remain incompletely understood. Experimental and theoretical studies often report discrepancies in the activation energy of proton diffusion, primarily due to limitations in modeling complex microstructures, such as grain boundaries and defects, at relevant time and length scales. With the high dimensionality of grain boundaries, new models are required to accurately capture the microstructure of polycrystalline barium zirconate.

Simulated hot press of Barium Zirconate powder

I simulated a barium zirconate powder under high pressure and temperature to simulate pellet formation via hot press. The pseudo-pellet is then considered for its microstructure and proton transport properties. With this consideration, we can better describe the proton transport properties of polycrystalline barium zirconate.

Calculating the Activity of a Single Atom Alloy Nanoparticle

Calculating the Activity of a Single Atom Alloy Nanoparticle

Combining active site sampling and microkinetics to get nanoparticle activity

Single atom alloys are unique as the only catalytically active sites are single atoms. This highly localized activity allows "simpler" kinetics due to the reactivity of the single atom being more homogeneous-like than heterogeneous in nature. Removal of surface diffusion allows the total rate equation to be simplified to the sum of probabilities of a dopant atom being in a surface site and the reactivity of that site. With this, I sampled the active sites of the nanoparticle using a machine learning interatomic potential. I then classified the sites using machine learning. To get the reaction rates, I used density functional theory, along with microkinetic modeling within the mean field approximation. This gave a total reaction rate for the whole nanoparticle catalyst.

Site sampling for total reactivity
Natural Phenomena at Scale

Natural Phenomena at Scale

Contrails from first principles

With water being the most abundant liquid on Earth's surface, ice nucleation is a critical process that continually changes the environment that we live and experience. It dictates the outcome of life itself, from the inhibition of ice nucleation through antifreeze proteins in cryophilic life combating sub-zero conditions, to the formation of snow in a treacherous blizzard. Beyond natural scenarios for ice nucleation, there are also anthropogenic causes. An example acutely tangible to us can be observed when a jet aircraft flies overhead: contrails.

Atomistic simulation of contrail formation

With experimental EDS data as an input, I performed simulated annealing with a developed Machine Learning Interatomic Potential to generate experimentally characteristic atomistic structures of soot. I could quantify ice nucleation to classify specific surface sites that drive nucleation. With this, the nucleating ability of soot from different fuels could be considered.

Reactivity of IrO2 undergoing dissolution

Reactivity of IrO2 undergoing dissolution

Coupling surface degradation and chemical activity with developed graph kinetic Monte Carlo

Catalysts are dynamic systems; if the rate of reaction changes, the catalyst changes. Despite this, models are often built assuming that the catalyst is static. This is most notable in electrocatalysis, where degradation and dissolution of the catalyst under potential is an expected occurence. Here, I employ a graph kinetic Monte Carlo method to model the reaction of OER on IrO2 while undergoing dissolution.

Surface degradation of IrO2 under reaction conditions

Traditional kinetic Monte Carlo assumed fixed lattice points and processes. With my method, that I call graph kinetic Monte Carlo, I can model a changing lattice by doing neighborhood analysis every time a new reaction occurs. With kinetic parameters generated for each possible surface state using DFT, along with the free energies of dissolution, both the reactivity and the degradation are coupled together. My method could be applied to other processes, such as oxidation/rusting, surface restructuring, and out-of-equilibrium reactions.