Perspective: Fifty years of density-functional theory in chemical physics

Orthogonal effect of exchange and correlation parameters in density functional theory to compute geometric and spectroscopic quantities

The influence of the composition of the functional used for density functional theory computations on one structural parameter (a dihedral angle) and a spectroscopic parameter (absorption wavelength) is assessed in this study on the basis of two molecules (flavonols). In this kind of molecules, these two parameters should be correlated according to the nature of the electronic transition involved. However, it is shown herein that by varying the proportion of true exchange and correlation while building a functional, it is possible to obtain independently a large range of values for these parameters without any relation with the underlying real values. Therefore, it is concluded that the choice of a functional after a benchmark, especially using user-defined functionals, should be carried out with great care to avoid such effects.

Machine learning for accuracy in density functional approximations

Machine learning techniques have found their way into computational chemistry as indispensable tools to accelerate atomistic simulations and materials design. In addition, machine learning approaches hold the potential to boost the predictive power of computationally efficient electronic structure methods, such as density functional theory, to chemical accuracy and to correct for fundamental errors in density functional approaches. Here, recent progress in applying machine learning to improve the accuracy of density functional and related approximations is reviewed. Promises and challenges in devising machine learning models transferable between different chemistries and materials classes are discussed with the help of examples applying promising models to systems far outside their training sets.

Accuracy limit of non-polarizable four-point water models: TIP4P/2005 vs OPC. Should water models reproduce the experimental dielectric constant?

The last generation of four center non-polarizable models of water can be divided into two groups: those reproducing the dielectric constant of water, as OPC, and those significantly underestimating its value, as TIP4P/2005. To evaluate the global performance of OPC and TIP4P/2005, we shall follow the test proposed by Vega and Abascal in 2011 evaluating about 40 properties to fairly address this comparison. The liquid-vapor and liquid-solid equilibria are computed, as well as the heat capacities, isothermal compressibilities, surface tensions, densities of different ice polymorphs, the density maximum, equations of state at high pressures, and transport properties. General aspects of the phase diagram are considered by comparing the ratios of different temperatures (namely, the temperature of maximum density, the melting temperature of hexagonal ice, and the critical temperature). The final scores are 7.2 for TIP4P/2005 and 6.3 for OPC. The results of this work strongly suggest that we have reached the limit of what can be achieved with non-polarizable models of water and that the attempt to reproduce the experimental dielectric constant deteriorates the global performance of the water force field. The reason is that the dielectric constant depends on two surfaces (potential energy and dipole moment surfaces), whereas in the absence of an electric field, all properties can be determined simply from just one surface (the potential energy surface). The consequences of the choice of the water model in the modeling of electrolytes in water are also discussed.

Advancements and Challenges in the Integration of Indium Arsenide and Van der Waals Heterostructures

The strategic integration of low-dimensional InAs-based materials and emerging van der Waals systems is advancing in various scientific fields, including electronics, optics, and magnetics. With their unique properties, these InAs-based van der Waals materials and devices promise further miniaturization of semiconductor devices in line with Moore's Law. However, progress in this area lags behind other 2D materials like graphene and boron nitride. Challenges include synthesizing pure crystalline phase InAs nanostructures and single-atomic-layer 2D InAs films, both vital for advanced van der Waals heterostructures. Also, diverse surface state effects on InAs-based van der Waals devices complicate their performance evaluation. This review discusses the experimental advances in the van der Waals epitaxy of InAs-based materials and the working principles of InAs-based van der Waals devices. Theoretical achievements in understanding and guiding the design of InAs-based van der Waals systems are highlighted. Focusing on advancing novel selective area growth and remote epitaxy, exploring multi-functional applications, and incorporating deep learning into first-principles calculations are proposed. These initiatives aim to overcome existing bottlenecks and accelerate transformative advancements in integrating InAs and van der Waals heterostructures.

Theoretical study of the structural and thermodynamic properties of U-He compounds under high pressure

Uranium is considered as a very important nuclear energy material because of the huge amount of energy it releases. As the main product of the spontaneous decay of uranium, it is difficult for helium to react with uranium because of its chemical inertness. Therefore, bubbles will be formed inside uranium, which could greatly reduce the performance of uranium or cause safety problems. Additionally, nuclear materials are usually operated in an environment of high-temperature and high-pressure, so it is necessary to figure out the exact state of helium inside uranium under extreme conditions. Here, we explored the structural stability of the U-He system under high pressure and high temperature by using density functional theory calculations. Two metastable phases are found between 50 and 400 GPa: U4He with space group Fmmm and U6He with space group P1̄. Both are metallic and adopt layered structures. Electron localization function calculation combined with charge density difference analysis indicates that there are covalent bonds between U and U atoms in both Fmmm-U4He and P1̄-U6He. Regarding the elastic modulus of α-U, the addition of helium has certain influence on the mechanical properties of uranium. Besides, first-principles molecular dynamics simulations were carried out to study the dynamical behavior of Fmmm-U4He and P1̄-U6He at high-temperature. It was found that Fmmm-U4He and P1̄-U6He undergo one-dimensional superionic phase transitions at 150 GPa. Our study revealed the exotic structure of U-He compounds beyond the formation of bubbles under high-pressure and high-temperature, which might be relevant to the performance and safety issues of nuclear materials under extreme conditions.

Efficient exact exchange using Wannier functions and other related developments in planewave-pseudopotential implementation of RT-TDDFT

The plane-wave pseudopotential (PW-PP) formalism is widely used for the first-principles electronic structure calculation of extended periodic systems. The PW-PP approach has also been adapted for real-time time-dependent density functional theory (RT-TDDFT) to investigate time-dependent electronic dynamical phenomena. In this work, we detail recent advances in the PW-PP formalism for RT-TDDFT, particularly how maximally localized Wannier functions (MLWFs) are used to accelerate simulations using the exact exchange. We also discuss several related developments, including an anti-Hermitian correction for the time-dependent MLWFs (TD-MLWFs) when a time-dependent electric field is applied, the refinement procedure for TD-MLWFs, comparison of the velocity and length gauge approaches for applying an electric field, and elimination of long-range electrostatic interaction, as well as usage of a complex absorbing potential for modeling isolated systems when using the PW-PP formalism.

Insights into Antisite Defect Complex Induced High Ferro-Piezoelectric Properties in KNbO3 Perovskite: First-Principles Study

Improving ferro-piezoelectric properties of niobate-based perovskites is highly desirable for developing eco-friendly high-performance sensors and actuators. Although electro-strain coupling is usually obtained by constructing multiphase boundaries via complex chemical compositions, defect engineering can also create opportunities for novel property and functionality advancements. In this work, a representative tetragonal niobate-based perovskite, i.e., KNbO3, is studied by using first-principles calculations. Two intrinsic types of Nb antisite defect complexes are selected to mimic alkali-deficiency induced excess Nb antisites in experiments. The formation energy, electronic profiles, polarization, and piezoelectric constants are systematically analyzed. It is shown that the structural distortion and chemical heterogeneity around the energetically favorable antisite pair defects, i.e., (NbK4·+KNb4'), lower the crystal symmetry of KNbO3 from tetragonal to triclinic phase, and facilitate polarization emergence and reorientation to substantially enhance intrinsic ferro-piezoelectricity (i.e., spontaneous polarization Ps of 68.2 μC/cm2 and piezoelectric strain constant d33 of 228.3 pC/N) without complicated doping and alloying.

Condensed Matter Systems Exposed to Radiation: Multiscale Theory, Simulations, and Experiment

This roadmap reviews the new, highly interdisciplinary research field studying the behavior of condensed matter systems exposed to radiation. The Review highlights several recent advances in the field and provides a roadmap for the development of the field over the next decade. Condensed matter systems exposed to radiation can be inorganic, organic, or biological, finite or infinite, composed of different molecular species or materials, exist in different phases, and operate under different thermodynamic conditions. Many of the key phenomena related to the behavior of irradiated systems are very similar and can be understood based on the same fundamental theoretical principles and computational approaches. The multiscale nature of such phenomena requires the quantitative description of the radiation-induced effects occurring at different spatial and temporal scales, ranging from the atomic to the macroscopic, and the interlinks between such descriptions. The multiscale nature of the effects and the similarity of their manifestation in systems of different origins necessarily bring together different disciplines, such as physics, chemistry, biology, materials science, nanoscience, and biomedical research, demonstrating the numerous interlinks and commonalities between them. This research field is highly relevant to many novel and emerging technologies and medical applications.

Exploring crystal structure-physical property relationships with pressure

From its conception, X-ray crystallography has provided a unique understanding of the structure, bonding and electronic state of materials, which, in turn, unlocks a means of examining the properties and function of crystalline systems. Using state-of-the-art single-crystal X-ray diffraction, along with UV-Vis spectroscopy and DFT calculations, Zwolenik et al. [(2024). IUCrJ, 11, 519-527] have provided a comprehensive study of the structure-optical property relationship of 1,3-diacetylpyrene with methodologies that are increasingly accessible to non-specialist laboratories.

Propagated (fragment) Pipek-Mezey Wannier functions in real-time time-dependent density functional theory

Localization procedures are an important tool for analysis of complex systems in quantum chemistry, since canonical molecular orbitals are delocalized and can, therefore, be difficult to align with chemical intuition and obscure information at the local level of the system. This especially applies to calculations obeying periodic boundary conditions. The most commonly used approach to localization is Foster-Boys Wannier functions, which use a unitary transformation to jointly minimize the second moment of the orbitals. This procedure has proven to be robust and fast but has a side effect of often mixing σ- and π-type orbitals. σ/π-separation is achieved by the Pipek-Mezey Wannier function (PMWF) approach [Lehtola and Jónsson, J. Chem. Theory Comput. 10, 642 (2014) and Jónsson et al., J. Chem. Theory Comput. 13, 460 (2017)], which defines the spread functional in terms of partial charges instead. We have implemented a PMWF algorithm in the CP2K software package using the Cardoso-Souloumiac algorithm to enable their application to real-time time-dependent density functional theory. The method is demonstrated on stacked CO2 molecules, linear acetylenic carbon, boron and nitrogen co-doped graphene, and nitrogen-vacancy doped diamond. Finally, we discuss its computational scaling and recent efforts to improve it with fragment approaches.

Triple Synergism Effect of Ammonium Nitrilotriacetate on the Chemical Mechanical Polishing Performance of Ruthenium Barrier Layers

As the feature size of integrated circuits continues to decrease, ruthenium (Ru) has been suggested as the successor to traditional Ta/TaN bilayers for barrier layer materials due to its unique properties. This research delves into the effects of ammonium nitrilotriacetate (NTA(NH4)3) on the chemical mechanical polishing (CMP) performance of Ru in H2O2-based slurry. The removal rate (RR) of Ru surged from 47 to 890 Å min-1, marking an increase of about 17 times. The essence of this mechanism lies in the triple synergistic effects of NTA(NH4)3 in promoting ruthenium (Ru) removal: 1) The interaction betweenNH 4 + ${\mathrm{NH}}_{\mathrm{4}}^{\mathrm{+}}$ from NTA(NH4)3 and SiO2 abrasives; 2) The chelating action of [(NH4)N(CH2COO)3]2- from NTA(NH4)3 on Ru and its oxides; 3) The ammoniation and chelation of Ru and its oxides byNH 4 + ${\mathrm{NH}}_{\mathrm{4}}^{\mathrm{+}}$ from NTA(NH4)3, which enhance the dissolution and corrosion of oxidized Ru, making its removal during the barrier layer CMP process more efficient through mechanical means. This research introduces a synergistic approach for the effective removal of Ru, shedding light on potential applications of CMP in the field of the integrated circuits.

Tucker Tensor Approach for Accelerating Fock Exchange Computations in a Real-Space Finite-Element Discretization of Generalized Kohn-Sham Density Functional Theory

The evaluation of Fock exchange is often the computationally most expensive part of hybrid functional density functional theory calculations in a systematically improvable, complete basis. In this work, we employ a Tucker tensor based approach that substantially accelerates the evaluation of the action of Fock exchange by transforming three-dimensional convolutional integrals into a tensor product of one-dimensional convolution integrals. Our numerical implementation uses a parallelization strategy that balances the memory and communication bottlenecks, alongside overlapping compute and communication operations to enhance computational efficiency and parallel scalability. The accuracy and computational efficiency are demonstrated on various systems, including Pt clusters of various sizes and a TiO2 cluster with 3684 electrons.

X-rays, electrons, and neutrons as probes of atomic matter

X-rays, electrons, and neutrons probe different properties of matter. X-rays feel electron density (ED). Electrons sense the electrostatic potential (ESP) of electrons and nuclei. Neutrons are sensitive to nuclear coherent scattering length (NCSL). While NCSL maps are widely understood to be different, ED and ESP maps are tacitly assumed to be similar. Here, I show that the belief in ED and ESP map equivalence is mistaken, but contains a grain of truth. Using density functional theory (DFT), the Bethe-Mott (BM) relation, and the Thomas-Fermi (TF) and Cromer-Mann (CM) atomic models, I show that ED and ESP maps are indeed more similar to each other than to NCSL maps. Nevertheless, peak and integrated map values depend differently on the atomic order number and on the contributions from electrons in the inner and outer CM shells. ED and ESP maps also differ in the sign and relative magnitude of excess charge effects.

Kohn-Sham Density in a Slater Orbital Basis Set

Finite, atom-centered Slater basis sets are used to determine approximate Kohn-Sham molecular orbitals. This is achieved by minimizing the kinetic energy plus the sum-squared difference between the Kohn-Sham density and the full configuration interaction density. As a result of the finite basis, a weight factor is introduced to balance the two minimization components. Results herein show that this can be done systematically, without sensitive dependence on the choice of scaling factor. In addition, the algorithm is applied to the LiH diatomic for fractional electron counts, where stretching the bond introduces significant reorganization of the electron density. The analysis will show the correct KS orbital structure and reveal the effects of correlation and electron locality on the KS solutions.

Spurious Oscillations Caused by Density Functional Approximations: Who is to Blame? Exchange or Correlation?

We analyze the varying susceptibilities of different density functional approximations (DFAs) to present spurious oscillations on the profiles of several vibrational properties. Among other problems, these spurious oscillations cause significant errors in harmonic and anharmonic IR and Raman frequencies and intensities. This work hinges on a judicious strategy to dissect the exchange and correlation components of DFAs and pinpoint the origins of these oscillations. We identify spurious oscillations in derivatives of all energy components with respect to nuclear displacements, including those energy terms that do not involve numerical integrations. These indirect spurious oscillations are attributed to suboptimal electron densities resulting from a self-consistent field procedure using a DFA that exhibits direct spurious oscillations. Direct oscillations stem from inaccurate numerical integration of the exchange and correlation energy density functionals. A thorough analysis of direct spurious oscillations reveals that only a handful of exchange and correlation components are insensitive to spurious oscillations, giving rise to three families of functionals, BH&H, LSDA, and BLYP. Among the functionals in these families, we encounter four widespread DFAs: BLYP, B3LYP, LC-BLYP, and CAM-B3LYP. Certain DFAs like PBE appear less sensitive to spurious oscillations due to compensatory cancellations between their energy components. Additionally, we found non-negligible but small oscillations in PBE and TPSS, which could be safely employed provided a sufficiently large integration grid is used in the calculations. These findings hint at the key components of current approximations to be improved and emphasize the necessity to develop accurate DFAs suitable for studying molecular spectroscopies.

Autobench V1.0: Benchmarking Automation for Electronic Structure Calculations

This work reports on new software for automatic conformer energy benchmarking calculations for flexible molecules. The software workflow consists of four parts: conformational search, preoptimization, optimization, and frequency calculations at a higher level and last calculations using several theoretical levels. The software was written to be user-friendly and versatile to be used by nonexperts in computational chemistry. Any theoretical levels available in either Gaussian 16 or ORCA 5 may be applied in the benchmarking study. The workflow will automatically run conformational search calculations and deal with conformers that converge to the same minimum and those that show a negative frequency. At the end of the workflow, the user will have the mean absolute deviations and the most accurate method/DFT functional and basis set in comparison to the benchmark to be applied for the molecular system of interest. Case examples are given at the end of the paper that may help users to get insight into the software's main features.

Ferroelectric Polarization and Single-Atom Catalyst Synergistically Promoting CO2 Photoreduction of CuBiP2Se6

Further improving the activity and selectivity of photocatalytic CO2 reduction remains a challenge. Herein, we propose a new strategy for synergistically promoting photocatalytic CO2 reduction by combining two-dimensional (2D) ferroelectric polarization and single-atom catalysis. Our calculations showed that the ferroelectric polarization of CuBiP2Se6 provides the internal driving force for the separation and migration of photogenerated carriers, which provides a prerequisite for enhancing the photocatalytic efficiency. In addition, the introduction of single Ag atoms can act as an electron reservoir to significantly modify the bonding configurations on the surface through proper static electron transfer, thus effectively promoting the adsorption and activation of CO2 molecules. More importantly, we found that switching the ferroelectric polarization can synergistically optimize the limiting potential as well as control the final products. This study provides a new approach for enhancing the catalytic activity and selectivity of photocatalytic CO2 reduction.

Comparative first principles investigation on the structural, optoelectronic and vibrational properties of strain-engineered graphene-like AlC3, BC3and C3N monolayers

Three cardinal two-dimensional semiconductorsviz., AlC3, BC3and C3N, closely resembling the graphene structure, are intriguing contenders for emerging optoelectronic and thermomechanical applications. Starting from a critical stability analysis, this density functional theory study delves into a quantitative assessment of structural, mechanical, electronic, optical, vibrational and thermodynamical properties of these monolayers as a function of biaxial strain(ε)in a sublinear regime(-2%⩽ε⩽4%)of elastic deformation. The structures with cohesive energies slightly smaller than graphene, manifest exceptional mechanical stiffness, flexibility and breaking stress. The mechanical parameters have been deployed to further cultivate acoustic attributes and thermal conductivity. The hexagonal structures with mixed ionic-covalent molecular bonds have indirect electronic band-gap and work-function acutely sensitive toε. Dispersions of optical dielectric function, energy loss, refractive index, extinction coefficient, reflectivity, absorption coefficient and conductivity are deciphered in the UV-Vis-NIR regime against strain, where particular frequency bands featuring high polarization, dissipation, absorbance or reflectance are identified. Phonon band-structure and density of states testify dynamic stability in the ground state for all systems except the compressed ones. A comprehensive group theoretical analysis is performed to cultivate rotational; infrared and Raman-active modes, and the nature of molecular vibrations is delineated. The red-shifting of phonon bands andE2g/A1gRaman peaks with increasingε, associates estimation of Grüneisen parameter. Finally, strain-induced alterations of thermodynamic quantities such as entropy, enthalpy, free energy, heat capacity and Debye temperature are studied, followed by a molecular dynamics-based stability assessment under canonical ensemble.

Charge density waves in two-dimensional transition metal dichalcogenides

Charge density wave (CDW is one of the most ubiquitous electronic orders in quantum materials. While the essential ingredients of CDW order have been extensively studied, a comprehensive microscopic understanding is yet to be reached. Recent research efforts on the CDW phenomena in two-dimensional (2D) materials provide a new pathway toward a deeper understanding of its complexity. This review provides an overview of the CDW orders in 2D with atomically thin transition metal dichalcogenides (TMDCs) as the materials platform. We mainly focus on the electronic structure investigations on the epitaxially grown TMDC samples with angle-resolved photoemission spectroscopy and scanning tunneling microscopy/spectroscopy as complementary experimental tools. We discuss the possible origins of the 2D CDW, novel quantum states coexisting with them, and exotic types of charge orders that can only be realized in the 2D limit.

Magnetic Ionic Liquids: Current Achievements and Future Perspectives with a Focus on Computational Approaches

Magnetic ionic liquids (MILs) stand out as a remarkable subclass of ionic liquids (ILs), combining the desirable features of traditional ILs with the unique ability to respond to external magnetic fields. The incorporation of paramagnetic species into their structures endows them with additional attractive features, including thermochromic behavior and luminescence. These exceptional properties position MILs as highly promising materials for diverse applications, such as gas capture, DNA extractions, and sensing technologies. The present Review synthesizes key experimental findings, offering insights into the structural, thermal, magnetic, and optical properties across various MIL families. Special emphasis is placed on unraveling the influence of different paramagnetic species on MILs' behavior and functionality. Additionally, the Review highlights recent advancements in computational approaches applied to MIL research. By leveraging molecular dynamics (MD) simulations and density functional theory (DFT) calculations, these computational techniques have provided invaluable insights into the underlying mechanisms governing MILs' behavior, facilitating accurate property predictions. In conclusion, this Review provides a comprehensive overview of the current state of research on MILs, showcasing their special properties and potential applications while highlighting the indispensable role of computational methods in unraveling the complexities of these intriguing materials. The Review concludes with a forward-looking perspective on the future directions of research in the field of magnetic ionic liquids.

Computational chemistry for water-splitting electrocatalysis

Electrocatalytic water splitting driven by renewable electricity has attracted great interest in recent years for producing hydrogen with high-purity. However, the practical applications of this technology are limited by the development of electrocatalysts with high activity, low cost, and long durability. In the search for new electrocatalysts, computational chemistry has made outstanding contributions by providing fundamental laws that govern the electron behavior and enabling predictions of electrocatalyst performance. This review delves into theoretical studies on electrochemical water-splitting processes. Firstly, we introduce the fundamentals of electrochemical water electrolysis and subsequently discuss the current advancements in computational methods and models for electrocatalytic water splitting. Additionally, a comprehensive overview of benchmark descriptors is provided to aid in understanding intrinsic catalytic performance for water-splitting electrocatalysts. Finally, we critically evaluate the remaining challenges within this field.

Decoding energy decomposition analysis: Machine-learned Insights on the impact of the density functional on the bonding analysis

The concept of chemical bonding is a crucial aspect of chemistry that aids in understanding the complexity and reactivity of molecules and materials. However, the interpretation of chemical bonds can be hindered by the choice of the theoretical approach and the specific method utilized. This study aims to investigate the effect of choosing different density functionals on the interpretation of bonding achieved through energy decomposition analysis (EDA). To achieve this goal, a data set was created, representing four bonding groups and various combinations of functionals and dispersion correction schemes. The calculations showed significant variation among the different functionals for the EDA terms, with the dispersion correction terms exhibiting the highest variability. More information was extracted by using machine learning in combination with dimensionality reduction on the data set. Results indicate that, despite the differences in the EDA terms obtained from different functionals, the functional has the least significant impact, suggesting minimal influence on the bonding interpretation.

Improved Carrier Separation and Recombination by Ferroelectric Polarization in the CuBiP2Se6/C2N Heterostructure: A Nonadiabatic Molecular Dynamics Study

The rapid recombination of photogenerated carriers heavily restricts the photocatalytic efficiency. Here, we propose a new strategy to improve catalytic efficiency based on the ferroelectric van der Waals heterostructure (CuBiP2Se6/C2N). Combining density functional theory and the nonadiabatic molecular dynamics (NAMD) method, we have systematically analyzed the ground-state properties and carrier dynamics images in the CuBiP2Se6/C2N heterostructure. Our calculations showed that the ferroelectric polarization of CuBiP2Se6 provides the internal driving force for the photogenerated carriers separation. NAMD results demonstrate that the excited-state carrier transfer and recombination processes in the CuBiP2Se6/C2N are consistent with a type II mechanism. Meanwhile, constructing the ferroelectric heterostructure can effectively prolong the carrier lifetime, from ∼65.98 to ∼124.54 ps. Moreover, the high quantum efficiency and tunable band edge positions mean that the CuBiP2Se6/C2N heterostructure is an excellent potential candidate material for photocatalytic water splitting.

The energy landscape of the ribosome

The ribosome is a prototypical assembly that can be used to establish general principles and techniques for the study of biological molecular machines. Motivated by the fact that the dynamics of every biomolecule is governed by an underlying energy landscape, there has been great interest to understand and quantify ribosome energetics. In the present review, we will focus on theoretical and computational strategies for probing the interactions that shape the energy landscape of the ribosome, with an emphasis on more recent studies of the elongation cycle. These efforts include the application of quantum mechanical methods for describing chemical kinetics, as well as classical descriptions to characterize slower (microsecond to millisecond) large-scale (10-100 Å) rearrangements, where motion is described in terms of diffusion across an energy landscape. Together, these studies provide broad insights into the factors that control a diverse range of dynamical processes in this assembly.

Real-Time Time-Dependent Density Functional Theory for Simulating Nonequilibrium Electron Dynamics

The explicit real-time propagation approach for time-dependent density functional theory (RT-TDDFT) has increasingly become a popular first-principles computational method for modeling various time-dependent electronic properties of complex chemical systems. In this Perspective, we provide a nontechnical discussion of how this first-principles simulation approach has been used to gain novel physical insights into nonequilibrium electron dynamics phenomena in recent years. Following a concise overview of the RT-TDDFT methodology from a practical standpoint, we discuss our recent studies on the electronic stopping of DNA in water and the Floquet topological phase as examples. Our discussion focuses on how RT-TDDFT simulations played a unique role in deriving new scientific understandings. We then discuss existing challenges and some new advances at the frontier of RT-TDDFT method development for studying increasingly complex dynamic phenomena and systems.

Correlation consistent basis sets designed for density functional theory: Third-row atoms (Ga-Br)

The correlation consistent basis sets (cc-pVnZ with n = D, T, Q, 5) for the Ga-Br elements have been redesigned, tuning the sets for use for density functional approximations. Steps to redesign these basis sets for an improved correlation energy recovery and efficiency include truncation of higher angular momentum functions, recontraction of basis set coefficients, and reoptimization of basis set exponents. These redesigned basis sets are compared with conventional cc-pVnZ basis sets and other basis sets, which are, in principle, designed to achieve systematic improvement with respect to increasing basis set size. The convergence of atomic energies, bond lengths, bond dissociation energies, and enthalpies of formation to the Kohn-Sham limit is improved relative to other basis sets where convergence to the Kohn-Sham limit is typically not observed.

Facet-dependent dispersion and aggregation of aqueous hematite nanoparticles

Nanoparticle aggregates in solution controls surface reactivity and function. Complete dispersion often requires additive sorbents to impart a net repulsive interaction between particles. Facet engineering of nanocrystals offers an alternative approach to produce monodisperse suspensions simply based on facet-specific interaction with solvent molecules. Here, we measure the dispersion/aggregation of three morphologies of hematite (α-Fe2O3) nanoparticles in varied aqueous solutions using ex situ electron microscopy and in situ small-angle x-ray scattering. We demonstrate a unique tendency of (104) hematite nanoparticles to maintain a monodisperse state across a wide range of solution conditions not observed with (001)- and (116)-dominated particles. Density functional theory calculations reveal an inert, densely hydrogen-bonded first water layer on the (104) facet that favors interparticle dispersion. Results validate the notion that nanoparticle dispersions can be controlled through morphology for specific solvents, which may help in the development of various nanoparticle applications that rely on their interfacial area to be highly accessible in stable suspensions.

From Molecules to Devices: Insights into Electronic and Optical Properties of Pyridine-Derived Compounds Using Density Functional Theory Calculations

In this study, we delve into the electronic structure, spectroscopic, and optical properties of five benzo derivatives of pyridine, namely, 5-(4-chlorophenyl)-2-fluoropyridine (1), 2-fluoro-5-(4-fluorophenyl)pyridine (2), 4-(2-fluoropyridin-5-yl)phenol (3), 5-(2,3-dichlorophenyl)-2-fluoropyridine (4), and 5-(5-bromo-2-methoxyphenyl)-2-fluoropyridine (5). Utilizing quantum chemical density functional theory calculations at the B3LYP and Perdew-Burke-Ernzerhof levels of theory combined with the 6-311G(d,p) and 6-311++G(d,p) basis sets, we investigated the electronic and optical characteristics of these compounds. Band structure calculations were conducted for their crystalline structures, revealing a direct band gap varying from 3.018 to 3.558 eV, with the valence band maximum and conduction band minimum located at the G point in the Brillouin zone. The optical properties were analyzed, including the dielectric functions, reflectivity, and refractive index. Notably, reflectivity was found to be minimal in the photon energy range of 0.0-3.0 eV, and the static refractive index, n(0), ranged from 1.55 to 1.70. The research also involved assessing the reactivity of the compounds through calculation of the frontier orbital energy gaps (ΔE), indicating a significant charge transfer and high reactivity. Additionally, we performed frequency analysis to unveil the Fourier-transform infrared spectra of compounds 1-5 at room temperature. Molecular electrostatic potential surfaces of the optimized structures were employed to map the electrophilic and nucleophilic regions of the compounds. This investigation provides a comprehensive understanding of the electronic and optical properties of these pyridine derivatives, shedding light on their potential applications in optoelectronics.

Photolysis of Dissolved Organic Matter over Hematite Nanoplatelets

Solar photoexcitation of chromophoric groups in dissolved organic matter (DOM), when coupled to photoreduction of ubiquitous Fe(III)-oxide nanoparticles, can significantly accelerate DOM degradation in near-surface terrestrial systems, but the mechanisms of these reactions remain elusive. We examined the photolysis of chromophoric soil DOM coated onto hematite nanoplatelets featuring (001) exposed facets using a combination of molecular spectroscopies and density functional theory (DFT) computations. Reactive oxygen species (ROS) probed by electron paramagnetic resonance (EPR) spectroscopy revealed that both singlet oxygen and superoxide are the predominant ROS responsible for DOM degradation. DFT calculations confirmed that Fe(II) on the hematite (001) surface, created by interfacial electron transfer from photoexcited chromophores in DOM, can reduce dioxygen molecules to superoxide radicals (•O2-) through a one-electron transfer process. 1H nuclear magnetic resonance (NMR) and electrospray ionization Fourier-transform ion cyclotron resonance mass spectrometry (ESI-FTICR-MS) spectroscopies show that the association of DOM with hematite enhances the cleavage of aromatic groups during photodegradation. The findings point to a pivotal role for organic matter at the interface that guides specific ROS generation and the subsequent photodegradation process, as well as the prospect of using ROS signatures as a forensic tool to help interpret more complicated field-relevant systems.

Stretching Bonds without Breaking Symmetries in Density Functional Theory

Kohn-Sham density functional theory (KS-DFT) stands out among electronic structure methods due to its balance of accuracy and computational efficiency. However, to achieve chemically accurate energies, standard density functional approximations in KS-DFT often need to break underlying symmetries, a long-standing "symmetry dilemma". By employing fragment spin densities as the main variables in calculations (rather than total molecular densities, as in KS-DFT), we present an embedding framework in which this symmetry dilemma is understood and partially resolved. The spatial overlap between fragment densities is used as the main ingredient to construct a simple, physically motivated approximation to a universal functional of the fragment densities. This "overlap approximation" is shown to significantly improve semilocal KS-DFT binding energies of molecules without artificially breaking either charge or spin symmetries. The approach is shown to be applicable to covalently bonded molecules and to systems of the "strongly correlated" type.

Solvent-Mediated Modulation of the Au-S Bond in Dithiol Molecular Junctions

Gold-dithiol molecular junctions have been studied both experimentally and theoretically. However, the nature of the gold-thiolate bond as it relates to the solvent has seldom been investigated. It is known that solvents can impact the electronic structure of single-molecule junctions, but the correlation between the solvent and dithiol-linked single-molecule junction conductance is not well understood. We study molecular junctions formed with thiol-terminated phenylenes from both 1-chloronaphthalene and 1-bromonaphthalene solutions. We find that the most probable conductance and the distribution of conductances are both affected by the solvent. First-principles calculations show that junction conductance depends on the binding configurations (adatom, atop, and bridge) of the thiolate on the Au surface, as has been shown previously. More importantly, we find that brominated solvents can restrict the binding of thiols to specific Au sites. This mechanism offers new insight into the effects of the solvent environment on covalent bonding in molecular junctions.

Exploring HCl-HCl interactions: QZVPP calculations, improved Lennard-Jones potential, and second virial coefficient analysis for thermodynamics and industrial applications

In this paper, we present a comprehensive analysis of HCl-HCl interactions, including QZVPP calculations, energy fitting, conformation validation, and the determination of the second virial coefficient B using improved Lennard-Jones (ILJ) potential parameters. To acquire accurate interaction energies, initial QZVPP calculations are performed on approximately 1851 randomly generated HCl-HCl conformations. Then, these energies are used to fit an improved Lennard-Jones potential energy surface, allowing for a robust description of HCl-HCl interactions. The ILJ potential parameters are then used to validate particular HCl dimer conformations, ensuring their stability and consistency with experimental observations. The correlation between calculated and experimental conformations strengthens the validity of the ILJ potential parameters. In addition, the second viral coefficient B is calculated at various temperatures using the ILJ potential. The obtained B values are compared to experimental data, demonstrating close agreement, and validating the ILJ potential's ability to accurately capture the intermolecular interactions and gas-phase behavior of the HCl-HCl system. The results of this study demonstrate the effective implementation of QZVPP calculations, energy fitting, and ILJ potential parameters in validating HCl-HCl conformations and accurately determining the second virial coefficient B. The high degree of concordance between calculated B values and experimental data demonstrates the validity of the ILJ potential and its suitability for modeling HCl-HCl interactions. This research contributes to a greater comprehension of HCl-HCl interactions and their implications for numerous chemical and atmospheric processes. The validated conformations, energy fitting method, and calculated second virial coefficients provide valuable instruments for future research and pave the way for more accurate modeling and simulations of HCl-HCl systems.

Doubling down on density-functional theory

In a recent paper, Becke et al. [J. Chem. Phys. 158, 151103 (2023)] presented a novel double hybrid density functional, "DH23," whose terms are based on good physics. Its 12 coefficients were trained on the GMTKN55 (general main-group thermochemistry, kinetics, and noncovalent interactions) chemical database of Goerigk et al. [Phys. Chem. Chem. Phys. 19, 32184 (2017)]. The lowest GMTKN55 "WTMAD2" error to date for any hybrid or double hybrid density functional was obtained (1.76 kcal/mol). Here, we make some revisions to DH23 and test its efficacy on reference data beyond GMTKN55, namely, organometallic reaction energies and barrier heights. The results confirm that DH23 is robust outside its training set. In the process, a slightly smaller GMTKN55 WTMAD2 of 1.73 kcal/mol is achieved.

New Quinoid Bio-Inspired Materials Using Para-Azaquinodimethane Moiety

Quinoid single molecules are regarded as promising materials for electronic applications due to their tunable chemical structure-driven properties. A series of three single bio-inspired quinoid materials containing para-azaquinodimethane (p-AQM) moiety were designed, synthesized and characterized. AQM1, AQM2 and AQM3, prepared using aldehydes derived from almonds, corncobs and cinnamon, respectively, were studied as promising quinoid materials for optoelectronic applications. The significance of facile synthetic procedures is highlighted through a straightforward two-step synthesis, using Knoevenagel condensation. The synthesized molecules showed molar extinction coefficients of 22,000, 32,000 and 61,000 L mol-1 cm-1, respectively, for AQM1, AQM2 and AQM3. The HOMO-LUMO energy gaps were calculated experimentally, theoretically showing the same trends: AQM3 < AQM2 < AQM1. The role of the aryl substituent was studied and showed an impact on the electronic properties. DFT calculations show planar structures with quinoidal bond length alternation, in agreement with the experimental results. Finally, these bio-based materials showed high thermal stabilities between 290 °C and 340 °C and a glassy behavior after the first heating-cooling scan. These results highlight these bio-based single molecules as potential candidates for electronic or biomedical applications.

First-Principles Calculations of Magnetite (Fe3O4) above the Verwey Temperature by Using Self-Consistent DFT + U + V

In this report, we have used the DFT + U + V approach, an extension of the DFT + U approach that takes into account both on-site and intersite interactions, to simulate structural, magnetic, and electronic properties together with the Fe and O K-edge XAS spectra of Fe3O4 above the Verwey temperature (Tv). Moreover, we compared the simulated XAS spectra with experimental XAS data. We examined both orthogonalized and nonorthogonalized atomic orbital projectors and compared DFT + U + V to DFT, DFT + U, and HSE as a hybrid functional. It is noteworthy that, despite the widespread use of the same Hubbard U value for Feoct and Fetet at the DFT + U level in the literature, the HP code identified two distinct values for them using the Hubbard approaches (DFT + U and DFT + U + V). The resulting Hubbard U and V parameters are strongly dependent on the chosen orbital projectors. This study demonstrates how DFT + U + V can improve the structural, magnetic, and electronic properties of Fe3O4 compared to approximate DFT and DFT + U. In this context, DFT + U + V supports the half-metallic character of the bulk crystal Fe3O4 above Tv, since the Fermi level is found in the t2g band with a Feoct down-spin. Thus, the observations in the current study emphasize the significance of intersite interactions in the theoretical analysis of Fe3O4 above the Tv.

Long-range Through-space Charge Transfer in a pH-responsive Mixed Cyclophane of Pyridine and Teropyrene

A large, strained (SE=44.2 kcal/mol) and conformationally flexible mixed cyclophane of pyridine and teropyrene was synthesized using two intramolecular Wurtz coupling reactions and an unprecedented Scholl reaction between the unreactive 2 positions of the pyrene systems in a triply bridged pyrenophane. Protonation of the pyridine unit results in a greatly enhanced preference for nesting in the cavity of the highly bent teropyrene system (θcalc =162.6°) and emergence of a charge transfer absorption band (λmax =592 nm) due to a long range (5.0-5.5 Å), through-space intramolecular transition between the teropyrene and pyridinium units, which does not exist in the neutral cyclophane.

Do Optimally Tuned Range-Separated Hybrid Functionals Require a Reparametrization of the Dispersion Correction? It Depends

For ground- and excited-state studies of large molecules, it is the state of the art to combine (time-dependent) DFT with dispersion-corrected range-separated hybrid functionals (RSHs), which ensures an asymptotically correct description of exchange effects and London dispersion. Specifically for studying excited states, it is common practice to tune the range-separation parameter ω (optimal tuning), which can further improve the accuracy. However, since optimal tuning essentially changes the functional, it is unclear if and how much the parameters used for the dispersion correction depend on the chosen ω value. To answer this question, we explore this interdependency by refitting the DFT-D4 dispersion model for six established RSHs over a wide range of ω values (0.05-0.45 a0-1) using a set of noncovalently bound molecular complexes. The results reveal some surprising differences among the investigated functionals: While PBE-based RSHs and ωB97M-D4 generally exhibit a weak interdependency and robust performance over a wide range of ω values, B88-based RSHs, specifically LC-BLYP, are strongly affected. For these, even a minor reduction of ω from the default value manifests in strong systematic overbinding and poor performance in the typical range of optimally tuned ω values. Finally, we discuss strategies to mitigate these issues and reflect the results in the context of the employed D4 parameter optimization algorithm and fit set, outlining strategies for future improvements.

Basis Set Limit of CCSD(T) Energies: Explicit Correlation Versus Density-Based Basis-Set Correction

A thorough comparison is carried out for explicitly correlated and density-based basis-set correction approaches, which were primarily developed to mitigate the basis-set incompleteness error of wave function methods. An efficient implementation of the density-based scheme is also presented, utilizing the density-fitting approximation. The performance of these approaches is comprehensively tested for the second-order Møller-Plesset (MP2), coupled-cluster singles and doubles (CCSD), and CCSD with perturbative triples [CCSD(T)] methods with respect to the corresponding complete basis set references. It is demonstrated that the density-based correction together with complementary auxiliary basis set (CABS)-corrected Hartree-Fock energies is highly robust and effectively reduces the error of the standard approaches; however, it does not outperform the corresponding explicitly correlated methods. Nevertheless, what still makes the density-corrected CCSD and CCSD(T) methods competitive is that their computational costs are roughly half of those of the corresponding explicitly correlated variants. Additionally, an incremental approach for standard CCSD and CCSD(T) is introduced. In this simple scheme, the total energies are corrected with the CABS correction and explicitly correlated MP2 contributions. As demonstrated, the resulting methods yield surprisingly good results, below 1 kcal/mol for thermochemical properties even with a double-ζ basis, while their computational expenses are practically identical to those of the density-based basis-set correction approaches.

Range-Separated Local Hybrid Functionals with Small Fractional-Charge and Fractional-Spin Errors: Escaping the Zero-Sum Game of DFT Functionals

Extending recent developments on strong-correlation (sc) corrections to local hybrid functionals to the recent accurate ωLH22t range-separated local hybrid, a series of highly flexible strong-correlation-corrected range-separated local hybrids (scRSLHs) has been constructed and evaluated. This has required the position-dependent reduction of both short- and long-range exact-exchange admixtures in regions of space characterized by strong static correlations. Using damping procedures provides scRSLHs that retain largely the excellent performance of ωLH22t for weakly correlated situations and, in particular, for accurate quasiparticle energies of a wide variety of systems while reducing dramatically static-correlation errors, e.g., in stretched-bond situations. An additional correction to the local mixing function to reduce delocalization errors in abnormal open-shell situations provides further improvements in thermochemical and kinetic parameters, making scRSLH functionals such as ωLH23tdE or ωLH23tdP promising tools for complex molecular or condensed-phase systems, where low fractional-charge and fractional-spin errors are simultaneously important. The proposed rung 4 functionals thereby largely escape the usual zero-sum game between these two quantities and are expected to open new areas of accurate computations by Kohn-Sham DFT. At the same time, they require essentially no extra computational effort over the underlying ωLH22t functional, which means that their use is only moderately more demanding than that of global, local, or range-separated hybrid functionals.

Electron-rich Au nanocrystals/Co3O4 interface for enhanced electrochemical nitrate reduction into ammonia

Solar-driven electrochemical NO3- reduction reaction (NO3-RR) is a clean and sustainable strategy that can convert pollutant NO3- in wastewater to value-added NH3. In recent years, cobalt oxides-based catalysts have shown their intrinsic catalytic properties toward NO3-RR but still have room for improvement through catalyst design. Coupling metal oxides with noble metal has been demonstrated to improve electrochemical catalytic efficiency. Here, we use Au species to tune the surface structure of Co3O4 and improve the efficiency of NO3-RR to NH3. The obtained Au nanocrystals-Co3O4 catalyst exhibited an onset potential of 0.54 V vs RHE, NH3 yield rate of 27.86 µg/h·cm2, and Faradaic efficiency (FE) of 83.1% at 0.437 V vs RHE in an H-cell, which is much higher than Au small species (Au clusters or single atoms)-Co3O4 (15.12 µg/h·cm2) and pure Co3O4 (11.38 µg/h·cm2), respectively. Combined experiments with theory calculations, we attributed the enhanced performance of Au nanocrystals-Co3O4 to the reduced energy barrier of *NO hydrogenation to the *NHO and suppression of HER, which originated from the charge transfer from Au to Co3O4. Using an amorphous silicon triple-junction (a-Si TJ) as the solar cell and an anion exchange membrane electrolyzer (AME), an unassisted solar-driven NO3-RR to NH3 prototype was realized with a yield rate of 4.65 mg/h and FE of 92.1%.

A Resolution of Identity Technique to Speed up TDDFT with Hybrid Functionals: Implementation and Application to the Magic Cluster Series Au8n+4(SC6H5)4n+8 (n = 3-6)

The Resolution of Identity (RI) technique has been employed to speed up the use of hybrid exchange-correlation (xc) functionals at the TDDFT level using the Hybrid Diagonal Approximation. The RI has been implemented within the polTDDFT algorithm (a complex damped polarization method) in the AMS/ADF suite of programs. A speedup factor of 30 has been obtained with respect to a previous numerical implementation, albeit with the same level of accuracy. Comparison of TDDFT simulations with the experimental photoabsorption spectra of the cluster series Au8n+4(SR)4n+8(n = 3-6; R = C6H5) showed the excellent accuracy and efficiency of the method. Results were compared with those obtained via the more simplified and computationally cheaper TDDFT+TB and sTDDFT methods. The present method represents an accurate as well as computationally affordable approach to predict photoabsorption spectra of complex species, realizing an optimal compromise between accuracy and computational efficiency, and is suitable for applications to large metal clusters with sizes up to several hundreds of atoms.

Short-range screened density matrix functional for proper descriptions of thermochemistry, thermochemical kinetics, nonbonded interactions, and singlet diradicals

Common one-electron reduced density matrix (1-RDM) functionals that depend on Coulomb and exchange-only integrals tend to underestimate dynamic correlation, preventing reduced density matrix functional theory (RDMFT) from achieving comparable accuracy to density functional theory in main-group thermochemistry and thermochemical kinetics. The recently developed ωP22 functional introduces a semi-local density functional to screen the erroneous short-range portion of 1-RDM functionals without double-counting correlation, potentially providing a better treatment of dynamic correlation around equilibrium geometries. Herein, we systematically evaluate the performance of this functional model, which consists of two parameters, on main-group thermochemistry, thermochemical kinetics, nonbonded interactions, and more. Tests on atomization energies, vibrational frequencies, and reaction barriers reveal that the ωP22 functional model can reliably predict properties at equilibrium and slightly away from equilibrium geometries. In particular, it outperforms commonly used density functionals in the prediction of reaction barriers, nonbonded interactions, and singlet diradicals, thus enhancing the predictive power of RDMFT for routine calculations of thermochemistry and thermochemical kinetics around equilibrium geometries. Further development is needed in the future to refine short- and long-range approximations in the functional model in order to achieve an excellent description of properties both near and far from equilibrium geometries.

Recent Advancements in Mechanistic Studies of Palladium- and Nickel-Catalyzed Ethylene Copolymerization with Polar Monomers

The introduction of polar functional groups into polyolefin chain structures creates opportunities to enhance specific properties, such as adhesion, dyeability, printability, compatibility, thermal stability, and electrical conductivity, which widen the range of potential applications for these modified materials. Transition metal catalysts, especially late transition metals, have proven to be highly effective in copolymerization processes due to their reduced Lewis acidity and electrophilicity. However, when compared to the significant progress and summary of synthetic methods, there is a distinct lack of a comprehensive summary of mechanistic studies pertaining to the catalytic systems involved in ethylene copolymerization catalyzed by palladium and nickel catalysts. In this review, we have provided a comprehensive summary of the latest developments in mechanistic studies of ethylene copolymerization with polar monomers catalyzed by late-transition-metal complexes. Experimental and computational methods were employed to conduct a detailed investigation of these organic and organometallic systems. It is mainly focused on ligand substitution, changes in binding modes, ethylene/polar monomer insertion, chelate opening, and β-H elimination. Factors that control the catalytic activity, molecular weight, comonomer incorporation ratios, and branch content are analyzed, these include steric repulsions between ligands and monomers, electronic effects arising from both ligands and monomers, and so on.

Interactions of arylhydroxylamines and alkylaldoximes with a rhodium porphyrin

Heme enzymes are involved in the binding and metabolism of hydroxylamine (RNHOH) and aldoxime (RCH=NOH) compounds (R = H, alkyl, aryl). We report the synthesis and X-ray crystal structure of a metalloporphyrin in complex with an arylhydroxylamine, namely that of (TPP)Rh(PhNHOH)(C6H4Cl) (TPP = tetraphenylpophryinato dianion). The crystal structure reveals, in addition to N-binding of PhNHOH to Rh, the presence of an intramolecular H-bond between the hydroxylamine -OH proton and a porphyrin N-atom. Results from density functional theory (DFT) calculations support the presence of this intramolecular H-bond in this global minimum structure, and a natural bond order (NBO) analysis reveals that this H-bond comprises a donor π N=C (porphyrin) to acceptor σ* O-H (hydroxylamine) interaction of 2.32 kcal/mol. While DFT calculations predict the presence of similar intramolecular H-bond interactions in the related aldoxime complexes (TPP)Rh(RCH=NOH)(C6H4Cl) in their global minima structures, the X-ray crystal structure obtained for the (TPP)Rh(CH3(CH2)2CH=NOH)(C6H4Cl) complex is consistent with the local (non-global) minima conformation that does not have this intramolecular H-bond interaction.

Jahn-Teller distortion in Sr2FeO4: group-theoretical analysis and hybrid DFT calculations

We present theoretical justification for distorted Ruddlesden-Popper (RP) phases of the first-order by using hybrid density functional theory (DFT) calculations and group-theoretical analysis. We, thus, demonstrate the existence of the Jahn-Teller effect around an Fe[Formula: see text] ion in Sr[Formula: see text]FeO[Formula: see text]. On the calculation side, we have established a combination of Wu-Cohen (WC) exchange and Perdew-Wang (PW) correlation in a three-parameter functional WC3PW, giving the most accurate description of Sr[Formula: see text]FeO[Formula: see text] from the comparison of three hybrid DFT functionals. Self-consistently obtained Hartree-Fock exact exchange of 0.16 demonstrates consistent results with the experimental literature data. Importantly, we explain conditions for co-existing proper and pseudo-Jahn-Teller effects from the crystalline orbitals, symmetry-mode analysis and irreps products. Moreover, phonon frequency calculations support and confirm the results of symmetry-mode analysis. In particular, the symmetry-mode analysis identifies a dominating irreducible representation of the Jahn-Teller mode (X2+) and corresponding space group (SG) of ground state structure (SG Cmce model). Therefore, the usually suggested high-symmetry tetragonal crystal structure (SG I4/mmm model) is higher in energy by 121 meV/f.u. (equivalent to the Jahn-Teller stabilization energy) compared with the distorted low-symmetry structure (SG Cmce model). We also present diffraction patterns for the two crystal symmetries to discuss the differences. Therefore, our results shed light on the existence of low-symmetry RP phases and make possible direct comparisons with future experiments.

Accurate Potential Energy Surfaces Using Atom-Centered Potentials and Minimal High-Level Data

We demonstrate that a Δ-density functional theory (Δ-DFT) approach based on atom-centered potentials (ACPs) represents a computationally inexpensive and accurate method for representing potential energy surfaces (PESs) for the HONO and HFCO molecules and vibrational frequencies derived therefrom. Using as few as 100 CCSD(T)-F12a reference energies, ACPs developed for use with B3LYP/def2-TZVPP are shown to produce PESs for HONO and HFCO with mean absolute errors of 27.7 and 5.8 cm-1, respectively. Application of the multiconfigurational time-dependent Hartree (MCTDH) method with ACP-corrected B3LYP/def2-TZVPP PESs produces vibrational frequencies for cis- and trans-HONO with mean absolute percent errors (MAPEs) of 0.8 and 1.1, compared to 0.8 obtained for the two isomers with CCSD(T)-F12a/cc-pVTZ-F12/MCTDH. For HFCO, the vibrational frequencies obtained using the present (Δ-DFT)/MCTDH approach give a MAPE of 0.1, which is the error obtained with CCSD(T)-F12a/cc-pVTZ-F12/MCTDH. The ACP approach is therefore successful in representing a PES calculated at a high level of theory (CCSD(T)-F12a) and a promising method for the development of a general protocol for the representation of accurate molecular PESs and the calculation of molecular properties from them.

PubChemQC B3LYP/6-31G*//PM6 Data Set: The Electronic Structures of 86 Million Molecules Using B3LYP/6-31G* Calculations

The presented "PubChemQC B3LYP/6-31G*//PM6" data set is composed of the electronic properties of 85,938,443 molecules, encompassing a broad spectrum of molecules from essential compounds to biomolecules with a molecular weight up to 1000. These molecules account for 94.0% of the original PubChem Compound catalog as of August 29, 2016. The electronic properties, including orbitals, orbital energies, total energies, dipole moments, and other pertinent properties, were computed by using the B3LYP/6-31G* and PM6 methods. The data set, available in three formats, namely, GAMESS quantum chemistry program files, selected JSON output files, and a PostgreSQL database, provides researchers with the ability to query molecular properties. It is further subdivided into five subdata sets for more specific data. The first two subsets encompass molecules with carbon, hydrogen, oxygen, and nitrogen with molecular weights under 300 and 500, respectively. The third and fourth subsets incorporate molecules with carbon, hydrogen, nitrogen, oxygen, phosphorus, sulfur, fluorine, and chlorine, with molecular weights under 300 and 500, respectively. The fifth subset comprises molecules with carbon, hydrogen, nitrogen, oxygen, phosphorus, sulfur, fluorine, chlorine, sodium, potassium, magnesium, and calcium, with a molecular weight of under 500. The coefficients of determination for the highest occupied molecular orbital-lowest unoccupied molecular orbital energy gap range from 0.892 (for CHON500) to 0.803 (for the whole data set). These comprehensive results pave the way for applications in drug discovery and materials science, among others. The data sets can be accessed under the Creative Commons Attribution 4.0 International license at the following web address: https://nakatamaho.riken.jp/pubchemqc.riken.jp/b3lyp_pm6_datasets.html.

Reproducibility of density functional approximations: How new functionals should be reported

Density functional theory is the workhorse of chemistry and materials science, and novel density functional approximations are published every year. To become available in program packages, the novel density functional approximations (DFAs) need to be (re)implemented. However, according to our experience as developers of Libxc [Lehtola et al., SoftwareX 7, 1 (2018)], a constant problem in this task is verification due to the lack of reliable reference data. As we discuss in this work, this lack has led to several non-equivalent implementations of functionals such as Becke-Perdew 1986, Perdew-Wang 1991, Perdew-Burke-Ernzerhof, and Becke's three-parameter hybrid functional with Lee-Yang-Parr correlation across various program packages, yielding different total energies. Through careful verification, we have also found many issues with incorrect functional forms in recent DFAs. The goal of this work is to ensure the reproducibility of DFAs. DFAs must be verifiable in order to prevent the reappearance of the above-mentioned errors and incompatibilities. A common framework for verification and testing is, therefore, needed. We suggest several ways in which reference energies can be produced with free and open source software, either with non-self-consistent calculations with tabulated atomic densities or via self-consistent calculations with various program packages. The employed numerical parameters-especially the quadrature grid-need to be converged to guarantee a ≲0.1 μEh precision in the total energy, which is nowadays routinely achievable in fully numerical calculations. Moreover, as such sub-μEh level agreement can only be achieved when fully equivalent implementations of the DFA are used, the source code of the reference implementation should also be made available in any publication describing a new DFA.

Efficient immobilization and detoxification of gaseous elemental mercury by nanoflower/rod WSe2/halloysite composite: Performance and mechanisms

Gaseous mercury pollution control technologies with low stability and high releasing risks always face with great challenges. Herein, we developed one halloysite nanotubes (HNTs)-supported tungsten diselenide (WSe2) composite (WSe2/HNTs) by one-pot solvothermal approach, curing Hg0 from complicated flue gas (CFG) and reducing second environment risks. WSe2 as a monolayer with nano-flower structure and HNTs with rod shapes in the as-prepared sorbent exhibited outstanding synergy efficiency, resulting in exceptional performance for Hg0 removal with high capture capacity of 30.6 mg·g-1 and rate of 9.09 μg·g-1·min-1, which benefited from the high affinity of selenium and mercury (1 ×1045) and the adequate exposure of Se-terminated. The adsorbent showed beneficial tolerance to high amount of NOx and SOx. An online lab-built thermal decomposition system (TPD-AFS) was employed to explore Hg species on the used-sorbent, finding that the adsorbed-mercury species were principally mercury selenide (HgSe). Density functional theory calculations indicated that the hollow-sites were the major adsorption sites and exhibited excellent selectivity for Hg0, as well as HgSe generation needed to overcome the 0.32 eV energy barrier. The adsorbed mercury displayed high environmental stability after the leaching toxicity test, which significantly decreased its secondary environmental risks. With these advantages, WSe2/HNTs possess enormous potential to achieve the effective and permanent immobilization of gaseous mercury from CFG in the future.

Better performance of Hartree-Fock over DFT: a quantum mechanical investigation on pyridinium benzimidazolate types of zwitterions in the light of localization/delocalization issues

CONTEXT

With the advent of fast computing facilities, combined with rapid emerges of many new and intricate quantum mechanical functionals, computations with pure Hartree-Fock (HF) theory are now-a-days regarded as trivial or obsolete, or even considered as not reliable by many researchers. Consequently, current trends in computational chemistry show extensive use of post-HF theories for smaller molecular systems and various DFT methods for organic and inorganic chemistry related problems (larger molecules/systems). In this contribution, I have tried to show that sometimes, HF might be more suitable over DFT methodologies in addressing structure-property correlations. Molecules studied here were previously synthesized by Boyd in 1966 and important experimental data were produced by Alcalde and co-workers in 1987. Comparison of computed and experimental results clearly shows that HF method was more effective in reproducing the experimental data compared to especially the DFT methodologies. Reliability of HF method was further assured from the very similar results shown by the CCSD, CASSCF, CISD and QCISD methods. Current study also indicates that the localization issue associated with HF proved to be advantageous over delocalization issue of DFT based methodologies, in correctly describing the structure-property correlation for zwitterion systems.

METHODS

All computations were performed with Gaussian 09. A wide-range of quantum mechanical methodologies, HF, B3LYP, CAM-B3LYP, BMK, B3PW91, TPSSh, LC-ωPBE, M06-2X, M06-HF, ωB97xD, MP2, CASSCF, CCSD, QCISD, CISD and semi-empirical methods like, Huckel, CNDO, AM1, PM3MM and PM6, were used for investigations.