Zuyu Liu
All work

UIUC Research · Computational Chemistry

Screening catalyst surfaces at 73,000 CPU-hour scale

Catalyst behaviour under real operating conditions can't be read off a structure — it has to be computed. I built and screened 40+ aMOC surface configurations with ASE and VASP across roughly 73,000 CPU-hours of density functional theory, assessing structural stability under reducing conditions relevant to hydrogen fuel cell applications.

Organisation
University of Illinois Urbana-Champaign
Role
Undergraduate Researcher
Location
Urbana, IL
Period
Sep 2024 – Present
Computational screening diagram. Left: a crystalline slab drawn side-on, with a distinct surface layer sitting above the bulk lattice and its bonds. Centre: a grid of 40+ candidate surface configurations, one of them marked as the configuration under screen. Right: the compute behind the screen — ~73,000 CPU-hours of queued density functional theory jobs, and a convergence trace settling onto its asymptote.SURFACE MODEL · ASECONFIGURATION SCREENDFT COMPUTE · VASPSURFACE40+AMOC SURFACE CONFIGURATIONS~73,000CPU-HOURSDFT JOBS
System diagram

Results

  • ~73,000

    CPU-hours of DFT computation

  • 40+

    aMOC surface configurations built and screened

  • ASE / VASP

    Computational stack for construction and calculation

Situation

The problem

Assessing whether a catalyst surface holds up under reducing conditions requires modelling many candidate configurations at the atomic scale. Each configuration is an expensive DFT calculation, the space of plausible surfaces is large, and stability is not obvious from structure alone — so the screening has to be systematic rather than intuition-led.

Scope

What I owned

  • Built and screened 40+ aMOC surface configurations for computational catalysis research.

  • Ran approximately 73,000 CPU-hours of density functional theory calculations using ASE and VASP.

  • Assessed structural stability of candidate surfaces under reducing conditions.

  • Connected atomic-scale modelling to potential hydrogen fuel cell applications.

  • Separately built Sankey diagrams mapping high-performance polymer supply chains (see Explorations).

Method

The approach

  1. Constructed candidate surface configurations systematically with ASE so the screen covers the space rather than sampling it arbitrarily.

  2. Ran DFT via VASP at a scale where job orchestration, convergence behaviour, and failure handling matter as much as the physics.

  3. Evaluated structural stability under reducing conditions — the environment the material would actually face.

  4. Tied results back to the application question: whether these surfaces are viable for hydrogen fuel cell contexts.

Constraints

Why it was hard

  • Compute at this scale is unforgiving: a misconfigured calculation is not a fast failure, it's thousands of wasted CPU-hours.

  • Convergence is not guaranteed. A meaningful share of the work is diagnosing why a calculation didn't settle rather than interpreting one that did.

  • 40+ configurations only produce a usable comparison if the setup is consistent across every one of them.

  • The result is a stability assessment, not a product — long-horizon work where the payoff is a credible answer, including a negative one.

Close

What this demonstrates

Comfort with long-horizon, compute-heavy technical work where the answer isn't guaranteed and the discipline is in the setup. This is the layer beneath the product — the part that has to be right before anything gets built on top of it.

Stack & methods

  • VASP
  • ASE
  • Python
  • DFT
  • HPC job orchestration

Let's build what comes after the prototype.

Open to 2027 internships and project work.