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!

Where are the stars in chemical space?

The variability of chemical bonding in open-shell transition-metal complexes not only motivates their study as functional materials and catalysts but also challenges conventional computational modeling tools. Here, tailoring ligand chemistry can alter preferred spin or oxidation states as well as electronic structure properties and reactivity, creating vast regions of chemical space to explore when designing new materials atom by atom.

Why do conventional catalyst design rules fail?

The design of selective and active C–H activation catalysts for direct methane-to-methanol conversion is challenging. Bioinspired complexes that form high-valent metal–oxo intermediates capable of hydrogen abstraction and rebound hydroxylation are promising candidates. This promise has made them a target for computational high-throughput screening, typically simplified through the use of linear free energy relationships (LFERs). However, their mid-row transition-metal centers have numerous accessible spin and oxidation states that increase the combinatorial scale of design efforts.

When are two hydrogen bonds better than one?

Hydrogen bonds (HBs) play an essential role in the structure and catalytic action of enzymes, but a complete understanding of HBs in proteins challenges the resolution of modern structural (i.e., X-ray diffraction) techniques and mandates computationally demanding electronic structure methods from correlated wavefunction theory for predictive accuracy. Numerous amino acid sidechains contain functional groups (e.g., hydroxyls in Ser/Thr or Tyr and amides in Asn/Gln) that can act as either HB acceptors or donors (HBA/HBD) and even form simultaneous, ambifunctional HB interactions.

Detecting strong correlation with ML

Despite its widespread use in chemical discovery, approximate density functional theory (DFT) is poorly suited to many targets, such as those containing open-shell, 3d transition metals that can be expected to have strong multi-reference (MR) character. For discovery workflows to be predictive, we need automated, low-cost methods that can distinguish the regions of chemical space where DFT should be applied from those where it should not.

Moving on up: 3d vs 4d TMCs in DFT

Density functional theory (DFT) is widely used in transition-metal chemistry, yet essential properties such as spin-state energetics in transition-metal complexes (TMCs) are well known to be sensitive to the choice of the exchange–correlation functional. Increasing the amount of exchange in a functional typically shifts the preferred ground state in first-row TMCs from low-spin to high-spin by penalizing delocalization error, but the effect on properties of second-row complexes is less well known.

RACs shed light on metal-organic frameworks

Millions of distinct metal-organic frameworks (MOFs) can be made by combining metal nodes and organic linkers. At present, over 90,000 MOFs have been synthesized and over 500,000 predicted. This raises the question whether a new experimental or predicted structure adds new information. For MOF chemists, the chemical design space is a combination of pore geometry, metal nodes, organic linkers, and functional groups, but at present we do not have a formalism to quantify optimal coverage of chemical design space.

Semi-supervised learning for MR detection

Multireference (MR) diagnostics are common tools for identifying strongly correlated electronic structure that makes single reference (SR) methods (e.g., density functional theory or DFT) insufficient for accurate property prediction. However, MR diagnostics typically require computationally demanding correlated wavefunction theory (WFT) calculations, and diagnostics often disagree or fail to predict MR effects on properties.

New insights into how Nafion breaks down

Polymer electrolyte membrane fuel cells (PEMFCs) represent promising energy storage solutions, but challenges remain to maximize their utility. Nafion is frequently employed as the PEMFC membrane material, but degradation of Nafion can limit the life of PEMFCs. Using hybrid density functional theory (DFT), we carry out reaction pathway analysis on a range of candidate degradation pathways on both pristine and defect-containing models of Nafion. Degradation of pristine Nafion initiated by hydrogen radicals involves moderate (ca.

Diagnosing strong correlation with help from ML

High-throughput computational screening typically employs methods (i.e., density functional theory or DFT) that can fail to describe challenging molecules, such as those with strongly correlated electronic structure. In such cases, multireference (MR) correlated wavefunction theory (WFT) would be the appropriate choice but remains more challenging to carry out and automate than single-reference (SR) WFT or DFT.

Sampling and QM region are equally important

Quantum-mechanical/molecular-mechanical (QM/MM) methods are essential to the study of metalloproteins, but the relative importance of sampling and degree of QM treatment in achieving quantitative predictions is poorly understood. We study the relative magnitude of configurational and QM-region sensitivity of energetic and electronic properties in a representative Zn2+ metal binding site of a DNA methyltransferase. To quantify property variations, we analyze snapshots extracted from 250 ns of molecular dynamics simulation.


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: