More Like this

1. DFTpy : An efficient and object‐oriented platform for orbital‐free DFT simulations

   https://doi.org/10.1002/wcms.1482  

   Shao, Xuecheng ; Jiang, Kaili ; Mi, Wenhui ; Genova, Alessandro ; Pavanello, Michele ( June 2020 , WIREs Computational Molecular Science)

   In silico materials design is hampered by the computational complexity of Kohn–Sham DFT, which scales cubically with the system size. Owing to the development of new‐generation kinetic energy density functionals (KEDFs), orbital‐free DFT (OFDFT) can now be successfully applied to a large class of semiconductors and such finite systems as quantum dots and metal clusters. In this work, we present DFTpy, an open‐source software implementing OFDFT written entirely in Python 3 and outsourcing the computationally expensive operations to third‐party modules, such as NumPy and SciPy. When fast simulations are in order, DFTpy exploits the fast Fourier transforms from PyFFTW. New‐generation, nonlocal and density‐dependent‐kernel KEDFs are made computationally efficient by employing linear splines and other methods for fast kernel builds. We showcase DFTpy by solving for the electronic structure of a million‐atom system of aluminum metal which was computed on a single CPU. The Python 3 implementation is object‐oriented, opening the door to easy implementation of new features. As an example, we present a time‐dependent OFDFT implementation (hydrodynamic DFT) which we use to compute the spectra of small metal clusters recovering qualitatively the time‐dependent Kohn–Sham DFT result. The Python codebase allows for easy implementation of application programming interfaces. We showcase the combination of DFTpy and ASE for molecular dynamics simulations of liquid metals. DFTpy is released under the MIT license.

   This article is categorized under:

   Software > Quantum Chemistry

   Electronic Structure Theory > Density Functional Theory

   Data Science > Computer Algorithms and Programming

    

   more » « less
2. The Landscape of Gold Nanocrystal Surface Chemistry

   https://doi.org/10.1021/acs.accounts.3c00109  

   Greskovich, Katherine M. ; Powderly, Kelly M. ; Kincanon, Maegen M. ; Forney, Nathan B. ; Jalomo, Catherine A. ; Wo, Anita ; Murphy, Catherine J. ( January 2023 , Accounts of Chemical Research)

   ConspectusGold nanoparticles (AuNPs) exhibit unique size- and shape-dependent properties not obtainable at the macroscale. Gold nanorods (AuNRs), with their morphology-dependent optical properties, ability to convert light to heat, and high surface-to-volume ratios, are of great interest for biosensing, medicine, and catalysis. While the gold core provides many fascinating properties, this Account focuses on AuNP soft surface coatings, which govern the interactions of nanoparticles with the local environments. Postmodification of AuNP surface chemistry can greatly alter NP colloidal stability, nano-bio interactions, and functionality. Polyelectrolyte coatings provide controllable surface-coating thickness and charge, which impact the composition of the acquired corona in biological settings. Covalent modification, in which covalently bound ligands replace the original capping layer, is often performed with thiols and disulfides due to their ability to replace native coatings. N-heterocyclic carbenes and looped peptides expand the possible functionalities of the ligand layer.The characterization of surface ligands bound to AuNPs, in terms of ligand density and dynamics, remains a challenge. Nuclear magnetic resonance (NMR) spectroscopy is a powerful tool for understanding molecular structures and dynamics. Our recent NMR work on AuNPs demonstrated that NMR data were obtainable for ligands on NPs with diameters up to 25 nm for the first time. This was facilitated by the strong proton NMR signals of the trimethylammonium headgroup, which are present in a distinct regime from other ligand protons’ signals. Ligand density analyses showed that the smallest AuNPs (below 4 nm) had the largest ligand densities, yet spin–spin T2 measurements revealed that these smallest NPs also had the most mobile ligand headgroups. Molecular dynamics simulations were able to reconcile these seemingly contradictory results.While NMR spectroscopy provides ligand information averaged over many NPs, the ligand distribution on individual particles’ surfaces must also be probed to fully understand the surface coating. Taking advantage of improvements in electron energy loss spectroscopy (EELS) detectors employed with scanning transmission electron microscopy (STEM), a single-layer graphene substrate was used to calibrate the carbon K-edge EELS signal, allowing quantitative imaging of the carbon atom densities on AuNRs with sub-nanometer spatial resolution. In collaboration with others, we revealed that the mean value for surfactant-bilayer-coated AuNRs had 10–30% reduced ligand density at the ends of the rods compared to the sides, confirming prior indirect evidence for spatially distinct ligand densities.Recent work has found that surface ligands on nanoparticles can, somewhat surprisingly, enhance the selectivity and efficiency of the electrocatalytic reduction of CO2 by controlling access to the active site, tuning its electronic and chemical environment, or denying entry to impurities that poison the nanoparticle surface to facilitate reduction. Looking to the future, while NMR and EELS are powerful and complementary techniques for investigating surface coatings on AuNPs, the frontier of this field includes the development of methods to probe the surface ligands of individual NPs in a high-throughput manner, to monitor nano-bio interactions within complex matrices, and to study structure–property relationships of AuNPs in biological systems. 

   more » « less
3. Toward a microscopic model of light absorbing dissolved organic compounds in aqueous environments: theoretical and experimental study

   https://doi.org/10.1039/D0CP06554D  

   Karimova, Natalia V. ; Alves, Michael R. ; Luo, Man ; Grassian, Vicki H. ; Gerber, R. Benny ( May 2021 , Physical Chemistry Chemical Physics)

   null (Ed.)

   Water systems often contain complex macromolecular systems that absorb light. In marine environments, these light absorbing components are often at the air–water interface and can participate in the chemistry of the atmosphere in ways that are poorly understood. Understanding the photochemistry and photophysics of these systems represents a major challenge since their composition and structures are not unique. In this study, we present a successful microscopic model of this light absorbing macromolecular species termed “marine derived chromophoric dissolved organic matter” or “m-CDOM” in water. The approach taken involves molecular dynamics simulations in the ground state using on the fly Density Functional Tight-Binding (DFTB) electronic structure theory; Time Dependent DFTB (TD-DFTB) calculations of excited states, and experimental measurements of the optical absorption spectra in aqueous solution. The theoretical hydrated model shows key features seen in the experimental data for a collected m-CDOM sample. As will be discussed, insights from the model are: (i) the low-energy A-band (at 410 nm) is due to the carbon chains combined with the diol- and the oxy-groups present in the structure; (ii) the weak B-band (at 320–360 nm) appears due to the contribution of the ionized speciated form of m-CDOM; and (iii) the higher-energy C-band (at 280 nm) is due to the two fused ring system. Thus, this is a two-speciated formed model. Although a relatively simple system, these calculations represent an important step in understanding light absorbing compounds found in nature and the search for other microscopic models of related materials remains of major interest. 

   more » « less
4. Modeling Spin-Crossover Dynamics

   https://doi.org/10.1146/annurev-physchem-101419-012625  

   Mukherjee, Saikat ; Fedorov, Dmitry A. ; Varganov, Sergey A. ( April 2021 , Annual Review of Physical Chemistry)

   In this article, we review nonadiabatic molecular dynamics (NAMD) methods for modeling spin-crossover transitions. First, we discuss different representations of electronic states employed in the grid-based and direct NAMD simulations. The nature of interstate couplings in different representations is highlighted, with the main focus on nonadiabatic and spin-orbit couplings. Second, we describe three NAMD methods that have been used to simulate spin-crossover dynamics, including trajectory surface hopping, ab initio multiple spawning, and multiconfiguration time-dependent Hartree. Some aspects of employing different electronic structure methods to obtain information about potential energy surfaces and interstate couplings for NAMD simulations are also discussed. Third, representative applications of NAMD to spin crossovers in molecular systems of different sizes and complexities are highlighted. Finally, we pose several fundamental questions related to spin-dependent processes. These questions should be possible to address with future methodological developments in NAMD. 

   more » « less
5. A chemical dynamics study on the gas-phase formation of triplet and singlet C 5 H 2 carbenes

   https://doi.org/10.1073/pnas.2019257117  

   He, Chao ; Galimova, Galiya R. ; Luo, Yuheng ; Zhao, Long ; Eckhardt, André K. ; Sun, Rui ; Mebel, Alexander M. ; Kaiser, Ralf I. ( December 2020 , Proceedings of the National Academy of Sciences)

   null (Ed.)

   Since the postulation of carbenes by Buchner (1903) and Staudinger (1912) as electron-deficient transient species carrying a divalent carbon atom, carbenes have emerged as key reactive intermediates in organic synthesis and in molecular mass growth processes leading eventually to carbonaceous nanostructures in the interstellar medium and in combustion systems. Contemplating the short lifetimes of these transient molecules and their tendency for dimerization, free carbenes represent one of the foremost obscured classes of organic reactive intermediates. Here, we afford an exceptional glance into the fundamentally unknown gas-phase chemistry of preparing two prototype carbenes with distinct multiplicities—triplet pentadiynylidene (HCCCCCH) and singlet ethynylcyclopropenylidene (c-C 5 H 2 ) carbene—via the elementary reaction of the simplest organic radical—methylidyne (CH)—with diacetylene (HCCCCH) under single-collision conditions. Our combination of crossed molecular beam data with electronic structure calculations and quasi-classical trajectory simulations reveals fundamental reaction mechanisms and facilitates an intimate understanding of bond-breaking processes and isomerization processes of highly reactive hydrocarbon intermediates. The agreement between experimental chemical dynamics studies under single-collision conditions and the outcome of trajectory simulations discloses that molecular beam studies merged with dynamics simulations have advanced to such a level that polyatomic reactions with relevance to extreme astrochemical and combustion chemistry conditions can be elucidated at the molecular level and expanded to higher-order homolog carbenes such as butadiynylcyclopropenylidene and triplet heptatriynylidene, thus offering a versatile strategy to explore the exotic chemistry of novel higher-order carbenes in the gas phase. 

   more » « less