Plane-Wave DFT Methods for Chemistry

Metal cations recognition by bowl-shaped N-pyrrolic polycyclic aromatic hydrocarbons

Bowl-shaped, nitrogen-doped polycyclic aromatic hydrocarbons were examined for the first time as molecular receptors for the recognition of metal cations. Potentiometric and spectrofluorimetric assays, supported with density functional theory computations, revealed that the title compounds recognise metal cations with a special focus on caesium (Cs+) cations.

Experimental and Computational Study of Pyrogenic Carbonaceous Matter Facilitated Hydrolysis of 2,4-Dinitroanisole (DNAN)

This study investigated the reaction pathway of 2,4-dinitroanisole (DNAN) on the pyrogenic carbonaceous matter (PCM) to assess the scope and mechanism of PCM-facilitated surface hydrolysis. DNAN degradation was observed at pH 11.5 and 25 °C with a model PCM, graphite, whereas no significant decay occurred without graphite. Experiments were performed at pH 11.5 due to the lack of DNAN decay at pH below 11.0, which was consistent with previous studies. Graphite exhibited a 1.78-fold enhancement toward DNAN decay at 65 °C and pH 11.5 relative to homogeneous solution by lowering the activation energy for DNAN hydrolysis by 54.3 ± 3.9%. This is supported by our results from the computational modeling using Car-Parrinello simulations by ab initio molecular dynamics/molecular mechanics (AIMD/MM) and DFT free energy simulations, which suggest that PCM effectively lowered the reaction barriers by approximately 8 kcal mol-1 compared to a homogeneous solution. Quaternary ammonium (QA)-modified activated carbon performed the best among several PCMs by reducing DNAN half-life from 185 to 2.5 days at pH 11.5 and 25 °C while maintaining its reactivity over 10 consecutive additions of DNAN. We propose that PCM can affect the thermodynamics and kinetics of hydrolysis reactions by confining the reaction species near PCM surfaces, thus making them less accessible to solvent molecules and creating an environment with a weaker dielectric constant that favors nucleophilic substitution reactions. Nitrite formation during DNAN decay confirmed a denitration pathway, whereas demethylation, the preferred pathway in homogeneous solution, produces 2,4-dinitrophenol (DNP). Denitration catalyzed by PCM is advantageous to demethylation because nitrite is less toxic than DNAN and DNP. These findings provide critical insights for reactive adsorbent design that has broad implications for catalyst design and pollutant abatement.

NWChem: Recent and Ongoing Developments

This paper summarizes developments in the NWChem computational chemistry suite since the last major release (NWChem 7.0.0). Specifically, we focus on functionality, along with input blocks, that is accessible in the current stable release (NWChem 7.2.0) and in the "master" development branch, interfaces to quantum computing simulators, interfaces to external libraries, the NWChem github repository, and containerization of NWChem executable images. Some ongoing developments that will be available in the near future are also discussed.

Finite Difference Interpolation for Reduction of Grid-Related Errors in Real-Space Pseudopotential Density Functional Theory

The real-space pseudopotential approach is a well-known method for large-scale density functional theory (DFT) calculations. One of its main limitations, however, is the introduction of errors associated with the positioning of the underlying real-space grid, a phenomenon usually known as the "egg-box" effect. The effect can be controlled by using a finer grid, but this raises the cost of the calculations or even undermines their feasibility altogether. Therefore, there is ongoing interest in the reduction of the effect per a given real-space grid. Here, we present a finite difference interpolation of electron orbitals as a means of exploiting the high resolution of the pseudopotential to reduce egg-box effects systematically. We implement the method in PARSEC, a finite difference real-space pseudopotential DFT code, and demonstrate error mitigation and improved convergence at a low additional computational cost.

Computational Predictions of the Hydrolysis of 2,4,6-Trinitrotoluene (TNT) and 2,4-Dinitroanisole (DNAN)

Hydrolysis is a common transformation reaction that can affect the environmental fate of many organic compounds. In this study, three proposed mechanisms of alkaline hydrolysis of 2,4,6-trinitrotoluene (TNT) and 2,4-dinitroaniline (DNAN) were investigated with plane-wave density functional theory (DFT) combined with ab initio and classical molecular dynamics (AIMD/MM) free energy simulations, Gaussian basis set DFT calculations, and correlated molecular orbital theory calculations. Most of the computations in this study were carried out using the Arrows web-based tools. For each mechanism, Meisenheimer complex formation, nucleophilic aromatic substitution, and proton abstraction reaction energies and activation barriers were calculated for the reaction at each relevant site. For TNT, it was found that the most kinetically favorable first hydrolysis steps involve Meisenheimer complex formation by attachment of OH^- at the C1 and C3 arene carbons and proton abstraction from the methyl group. The nucleophilic aromatic substitution reactions at the C2 and C4 arene carbons were found to be thermodynamically favorable. However, the calculated activation barriers were slightly lower than in previous studies, but still found to be ΔG^‡ ≈ 18 kcal/mol using PBE0 AIMD/MM free energy simulations, suggesting that the reactions are not kinetically significant. For DNAN, the barriers of nucleophilic aromatic substitution were even greater (ΔG^‡ > 29 kcal/mol PBE0 AIMD/MM). The most favorable hydrolysis reaction for DNAN was found to be a two-step process in which the hydroxyl first attacks the C1 carbon to form a Meisenheimer complex at the C1 arene carbon C1-(OCH₃)OH^-, and subsequently, the methoxy anion (-OCH₃) at the C1 arene carbon dissociates and the proton shuttles from the C1-OH to the dissociated methoxy group, resulting in methanol and an aryloxy anion.

From NWChem to NWChemEx: Evolving with the Computational Chemistry Landscape

Since the advent of the first computers, chemists have been at the forefront of using computers to understand and solve complex chemical problems. As the hardware and software have evolved, so have the theoretical and computational chemistry methods and algorithms. Parallel computers clearly changed the common computing paradigm in the late 1970s and 80s, and the field has again seen a paradigm shift with the advent of graphical processing units. This review explores the challenges and some of the solutions in transforming software from the terascale to the petascale and now to the upcoming exascale computers. While discussing the field in general, NWChem and its redesign, NWChemEx, will be highlighted as one of the early codesign projects to take advantage of massively parallel computers and emerging software standards to enable large scientific challenges to be tackled.

NWChem: Past, present, and future

Specialized computational chemistry packages have permanently reshaped the landscape of chemical and materials science by providing tools to support and guide experimental efforts and for the prediction of atomistic and electronic properties. In this regard, electronic structure packages have played a special role by using first-principle-driven methodologies to model complex chemical and materials processes. Over the past few decades, the rapid development of computing technologies and the tremendous increase in computational power have offered a unique chance to study complex transformations using sophisticated and predictive many-body techniques that describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, novel parallel algorithms have been able to take advantage of computational resources to address the polynomial scaling of electronic structure methods. In this paper, we briefly review the NWChem computational chemistry suite, including its history, design principles, parallel tools, current capabilities, outreach, and outlook.

Association of Defects and Zinc in Hematite

Zn is an essential micronutrient that is often limited in tropical, lateritic soils in part because it is sequestered in nominally refractory iron oxide phases. Stable phases such as goethite and hematite, however, can undergo reductive recrystallization without a phase change under circumneutral pH conditions and release metal impurities such as Zn into aqueous solutions. Further, the process appears to be driven by Fe vacancies. In this contribution, we used ab initio molecular dynamics informed extended X-ray absorption fine structure spectra to show that Zn incorporated in the structure of hematite is associated with coupled O-Fe and protonated Fe vacancies, providing a potential link between crystal chemistry and the bioavailability of Zn.

Corresponding Orbitals Derived from Periodic Bloch States for Electron Transfer Calculations of Transition Metal Oxides

An approach for modeling electron transfer in solids and at surfaces of iron-(oxyhydr)oxides and other redox active solids has been developed for electronic structure methods (i.e., plane-wave density functional theory) capable of performing calculations with periodic cells and large system sizes efficiently while at the same time being accurate enough to be used in the estimation of the electron-transfer coupling matrix element, V _AB, and the electron transfer transmission factor, κ_el. This method is an extension of the valence bond theory electron transfer method for molecules and clusters implemented by Dupuis and others and used extensively by Rosso and co-workers in which scaled corresponding orbitals derived from the Bloch states are used to calculate the off-diagonal matrix elements H _AB and S _AB. A key development of the present work is the formulation of algorithms to improve the accuracy of the integration of the exact exchange integral in periodic boundary conditions. This method is demonstrated on model systems for electron small polaron transfer in iron-(oxyhydr)oxides, including bare Fe²⁺-Fe³⁺ ions, and in [Fe³⁺(OH₂)₂ (OH^-)₂)] _nⁿ⁺ chains representing the common edge-sharing Fe octahedral motif in these materials.