Computational Physicist

Extracting structure from high-dimensional, noise-dominated complex systems

I build theoretical frameworks and computational strategies for problems where signals hide in noise and structure must be inferred from incomplete data. Analytical solutions when they exist, large-scale simulations when they don't — and increasingly, machine learning that respects the underlying physics. Neural networks gain power when you encode the right constraints. The most interesting problems are the ones where brute force fails and elegance is required.

10⁹ Atoms simulated
130K+ Processor cores
$1M+ Competitive grants
20 Publications

3D Chromosome Reconstruction

Structure recovered from 2D contact frequency data
01

Expertise

Inverse Problems

Reconstructing 3D configurations from sparse 2D measurements using physics-informed constraints.

Molecular Dynamics

Large-scale simulations — billions of particles on massively parallel architectures.

Analytical Methods

Closed-form solutions for stochastic systems, polymer physics, and statistical mechanics.

Phase Transitions

Dynamical arrest, glass transitions, and the signatures of regime change.

Polymer Physics

Topology, entanglement, and collective dynamics in confined and knotted systems.

High-Performance Computing

130,000+ cores, distributed algorithms, petabyte-scale data pipelines.

02

Selected Work

I built polymer models and ran molecular dynamics simulations to reconstruct the 3D chromosomal environment where Xist RNA spreads. The modeling showed how chromatin architecture constrains the spreading — basically, the physical landscape determines the biology.

I developed a method to recover full spatiotemporal dynamics — 3D structure plus how it evolves — from Hi-C contact frequency matrices. These are sparse, noisy measurements of which parts of the genome touch each other. Polymer physics plus constrained optimization lets you pull continuous structure out of discrete contacts.

J. Comp. Chem.

Scaling MD beyond 100,000 cores

2019

The first billion-atom biomolecular simulation. One billion particles, Newtonian dynamics, 130,000+ cores on Trinity. We had to figure out how to distribute force calculations, manage state across thousands of nodes, and make sense of petabytes of trajectory data.

I derived the analytical conditions for when a polymer system arrests — transitions from flowing to stuck. Turns out connectivity alone determines the critical point. It happens exactly where constraint counting says rigidity should emerge.

Scientific Reports

Ratcheted diffusion in nanochannels

2013

How do you get directional motion from random thermal noise? Break the symmetry. I worked out the analytics and simulations for Brownian ratchets in confined geometries — the same mechanism molecular motors use.

03

Background

My path hasn't been linear. I trained as a physicist at Cambridge, spent three years at Los Alamos in the Center for Nonlinear Studies, and am now a Principal Investigator at Harvard and Mass General. Along the way: Brownian ratchets, viral proteins, chromosome dynamics.

What connects it all is the problem of finding structure when noise dominates. I go for analytical solutions when they exist — there's something satisfying about an exact result. When that's not possible, I build simulations big enough that the structure can emerge on its own.

The problems I like best sit right at the edge of tractable. The ones where the right representation, the right constraints, the right level of description turns something that looks impossible into something you can actually solve.

Tools

Python, C++, CUDA, MPI. Molecular dynamics, statistical mechanics, stochastic processes, polymer physics, HPC, signal processing.

2019 – Present

Principal Investigator

Harvard University & Mass General Hospital

2016 – 2019

CNLS Fellow

Los Alamos National Laboratory

2011 – 2015

PhD Physics

University of Cambridge

2007 – 2011

BSc Natural Sciences

University of Edinburgh

Get in touch

Always interested in hard problems.