Research in the Kulik group

Research in the Kulik group leverages computational modeling to aid the discovery of new materials and mechanisms. Our group uses first-principles modeling to unearth fundamental aspects of structure-property relationships in catalysts and materials. By taking a computational approach, we carry out studies that allow us to make connections across a wide range of catalytic systems from biological enzymes to emerging heterogeneous single-atom catalysts. We develop computational software and machine learning models that accelerate the discovery of new materials and design rules. This approach enables the prediction of new materials properties in seconds, the exploration of million-compound design spaces, and the identification of design rules and exceptions that go beyond intuition. To ensure the predictive power of our approach, our group develops new methods to increase the accuracy of density functional theory especially for materials with challenging electronic structure such as transition metal complexes and solids.

Read more about our group’s work in the recently published papers below!

DFT for molecule-surface interactions

First-principles simulation has played an ever-increasing role in the discovery and interpretation of the chemical properties of surface–adsorbate interactions. Nevertheless, key challenges remain for the computational chemist wishing to study surface chemistry: modelling the full extent of experimental conditions, managing computational cost, minimizing human effort in simulation set-up and maximizing accuracy. Our recent work introduces new tools for streamlining surface chemistry simulation set-up and reviews some of the challenges in first-principles, density functional theory (DFT) simulation of surface phenomena. 

Density delocalization in DFT

Approximate DFT is well-known to suffer from self-interaction error, which is expected to particularly plague the localized 3d and 4f electrons of transition metal complexes. In order to diagnose SIE, energetic delocalization error, i.e. deviation from piecewise linearity, is frequently used, but errors in the density are less well-understood.

Large-scale QM/MM in enzymology: COMT

Hybrid quantum mechanical–molecular mechanical (QM/MM) simulations are widely used in studies of enzymatic catalysis. Until recently, it has been cost prohibitive to determine the asymptotic limit of key energetic and structural properties with respect to increasingly large QM regions.

Hydrogen bond design for ion separation

Selective ion separation is a major challenge with far-ranging impact from wastewater treatment to product separation in catalysis. The Hatton group here at MIT has recently pioneered the synthesis of Ferrocenium (Fc+)/ferrocene (Fc) polymeric electrode materials for catalysis and ion separation. In earlier collaborative work (Xiao Su et al Adv. Funct.

Substrate positioning in COMT reactivity

Catechol O-methyltransferase (COMT) is a SAM- and Mg2+-dependent methyltransferase that regulates neurotransmitters through methylation. Simulations and experiments have identified divergent catecholamine substrate orientations in the COMT active site: molecular dynamics simulations have favored a monodentate coordination of catecholate substrates to the active site Mg2+, and crystal structures instead preserve bidentate coordination along with short (2.65 Å) methyl donor-acceptor distances. We carry out longer dynamics (up to 350 ns) to quantify interconversion between bidentate and monodentate binding poses.

The molSimplify code automates discovery

We recently developed a code to automate the discovery and generation of inorganic complex structures. The code is available both as a user-friendly GUI and through commandline or input file interface on OS X and Linux. You can download the latest version of the code, the user guide and compiled executables by visiting

Hubbard U for correcting minimal basis sets

We recently introduced an approach to get more mileage out of minimal basis set calculations for geometry optimization on graphical processing units. In this approach, we take the Hubbard U +U correction normally applied for correcting self-interaction errrors in DFT to correct basis set incompleteness error in formally-SIE free Hartree-Fock theory.

Fullerene allotropes throughout the periodic table

Thirty years after the discovery of Buckminsterfullerene (C60), the excellent properties and potential applications of this unusual carbon allotrope continue to drive considerable scientific inquiry.

Watching QDs form with AIMD

Colloidal quantum dots (QDs), such as Indium Phosphide QDs, exhibit highly desirable size- and shape-dependent properties for applications from electronic devices to imaging. Production of InP QDs with the desired properties has lagged behind other QD materials due to a poor understanding of how to tune the growth process. Using high-temperature ab initio molecular dynamics (AIMD) simulations, we report the first direct observation of early-stage intermediates and subsequent formation of an InP cluster from indium and phosphorus precursors. In our simulations, indium agglomeration precedes formation of In–P bonds. We observe a predominantly intercomplex pathway in which In–P bonds form between one set of precursor copies, and the carboxylate ligand of a second indium precursor in the agglomerated indium abstracts a ligand from the phosphorus precursor. This process produces an indium-rich cluster with structural properties comparable to those in bulk zinc-blende InP crystals.

New paradigms for enzyme catalysis

We use a combination of long-time molecular dynamics and ab initio electronic structure theory to probe the underlying characteristic motions and interactions that facilitate catalysis of reactions in the environment of an enzyme. Recent work on a ubiquitous methyltransferase has been published in collaboration with an experimental group in the Proceedings of the National Academy of Sciences.


About Us

The Kulik group focuses on the development and application of new electronic structure methods and atomistic simulations tools in the broad area of catalysis.

Our Interests

We are interested in transition metal chemistry, with applications from biological systems (i.e. enzymes) to nonbiological applications in surface science and molecular catalysis.

Our Focus

A key focus of our group is to understand mechanistic features of complex catalysts and to facilitate and develop tools for computationally driven design.

Contact Us

Questions or comments? Let us know! Contact Dr. Kulik: