Heather J. Kulik

Bio: Heather J. Kulik is an academic researcher from Massachusetts Institute of Technology. The author has contributed to research in topics: Density functional theory & Medicine. The author has an hindex of 37, co-authored 136 publications receiving 4404 citations. Previous affiliations of Heather J. Kulik include Lawrence Livermore National Laboratory & Stanford University.

Density functional theory in transition-metal chemistry: A self-consistent Hubbard U approach

Abstract: Transition-metal centers are the active sites for a broad variety of biological and inorganic chemical reactions. Notwithstanding this central importance, density-functional theory calculations based on generalized-gradient approximations often fail to describe energetics, multiplet structures, reaction barriers, and geometries around the active sites. We suggest here an alternative approach, derived from the Hubbard U correction to solid-state problems, that provides an excellent agreement with correlated-electron quantum chemistry calculations in test cases that range from the ground state of Fe-2 and Fe-2(-) to the addition elimination of molecular hydrogen on FeO+. The Hubbard U is determined with a novel self-consistent procedure based on a linear-response approach.

Critical Knowledge Gaps in Mass Transport through Single-Digit Nanopores: A Review and Perspective

Abstract: Not all nanopores are created equal. By definition, nanopores have characteristic diameters or conduit widths between ∼1 and 100 nm. However, the narrowest of such pores, perhaps best called Single...

Understanding the diversity of the metal-organic framework ecosystem

Abstract: 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. In this work, we develop a machine learning method to quantify similarities of MOFs to analyse their chemical diversity. This diversity analysis identifies biases in the databases, and we show that such bias can lead to incorrect conclusions. The developed formalism in this study provides a simple and practical guideline to see whether new structures will have the potential for new insights, or constitute a relatively small variation of existing structures. At present there are databases with over 500,000 predicted or synthesized MOF structures, yet a method to establish whether a new material adds new information does not exist. Here the authors propose a machine-learning based approach to quantify the structural and chemical diversity in common MOF databases.

Protection of tissue physicochemical properties using polyfunctional crosslinkers.

Abstract: Understanding complex biological systems requires the system-wide characterization of both molecular and cellular features. Existing methods for spatial mapping of biomolecules in intact tissues suffer from information loss caused by degradation and tissue damage. We report a tissue transformation strategy named stabilization under harsh conditions via intramolecular epoxide linkages to prevent degradation (SHIELD), which uses a flexible polyepoxide to form controlled intra- and intermolecular cross-link with biomolecules. SHIELD preserves protein fluorescence and antigenicity, transcripts and tissue architecture under a wide range of harsh conditions. We applied SHIELD to interrogate system-level wiring, synaptic architecture, and molecular features of virally labeled neurons and their targets in mouse at single-cell resolution. We also demonstrated rapid three-dimensional phenotyping of core needle biopsies and human brain cells. SHIELD enables rapid, multiscale, integrated molecular phenotyping of both animal and clinical tissues.

Fast parallel algorithms for short-range molecular dynamics

Abstract: Three parallel algorithms for classical molecular dynamics are presented. The first assigns each processor a fixed subset of atoms; the second assigns each a fixed subset of inter-atomic forces to compute; the third assigns each a fixed spatial region. The algorithms are suitable for molecular dynamics models which can be difficult to parallelize efficiently—those with short-range forces where the neighbors of each atom change rapidly. They can be implemented on any distributed-memory parallel machine which allows for message-passing of data between independently executing processors. The algorithms are tested on a standard Lennard-Jones benchmark problem for system sizes ranging from 500 to 100,000,000 atoms on several parallel supercomputers--the nCUBE 2, Intel iPSC/860 and Paragon, and Cray T3D. Comparing the results to the fastest reported vectorized Cray Y-MP and C90 algorithm shows that the current generation of parallel machines is competitive with conventional vector supercomputers even for small problems. For large problems, the spatial algorithm achieves parallel efficiencies of 90% and a 1840-node Intel Paragon performs up to 165 faster than a single Cray C9O processor. Trade-offs between the three algorithms and guidelines for adapting them to more complex molecular dynamics simulations are also discussed.

Maximally-localized Wannier Functions: Theory and Applications

Abstract: The electronic ground state of a periodic system is usually described in terms of extended Bloch orbitals, but an alternative representation in terms of localized "Wannier functions" was introduced by Gregory Wannier in 1937. The connection between the Bloch and Wannier representations is realized by families of transformations in a continuous space of unitary matrices, carrying a large degree of arbitrariness. Since 1997, methods have been developed that allow one to iteratively transform the extended Bloch orbitals of a first-principles calculation into a unique set of maximally localized Wannier functions, accomplishing the solid-state equivalent of constructing localized molecular orbitals, or "Boys orbitals" as previously known from the chemistry literature. These developments are reviewed here, and a survey of the applications of these methods is presented. This latter includes a description of their use in analyzing the nature of chemical bonding, or as a local probe of phenomena related to electric polarization and orbital magnetization. Wannier interpolation schemes are also reviewed, by which quantities computed on a coarse reciprocal-space mesh can be used to interpolate onto much finer meshes at low cost, and applications in which Wannier functions are used as efficient basis functions are discussed. Finally the construction and use of Wannier functions outside the context of electronic-structure theory is presented, for cases that include phonon excitations, photonic crystals, and cold-atom optical lattices.

Ab initio calculation of vibrational absorption and circular dichroism spectra using density functional force fields

Abstract: : The unpolarized absorption and circular dichroism spectra of the fundamental vibrational transitions of the chiral molecule, 4-methyl-2-oxetanone, are calculated ab initio. Harmonic force fields are obtained using Density Functional Theory (DFT), MP2, and SCF methodologies and a 5S4P2D/3S2P (TZ2P) basis set. DFT calculations use the Local Spin Density Approximation (LSDA), BLYP, and Becke3LYP (B3LYP) density functionals. Mid-IR spectra predicted using LSDA, BLYP, and B3LYP force fields are of significantly different quality, the B3LYP force field yielding spectra in clearly superior, and overall excellent, agreement with experiment. The MP2 force field yields spectra in slightly worse agreement with experiment than the B3LYP force field. The SCF force field yields spectra in poor agreement with experiment.The basis set dependence of B3LYP force fields is also explored: the 6-31G* and TZ2P basis sets give very similar results while the 3-21G basis set yields spectra in substantially worse agreements with experiment. jg

Density functional theory for transition metals and transition metal chemistry

Abstract: We introduce density functional theory and review recent progress in its application to transition metal chemistry. Topics covered include local, meta, hybrid, hybrid meta, and range-separated functionals, band theory, software, validation tests, and applications to spin states, magnetic exchange coupling, spectra, structure, reactivity, and catalysis, including molecules, clusters, nanoparticles, surfaces, and solids.