Quantum Mechanical Electronic Structure Research Articles

MolTwister – a molecular systems construction, manipulation and statistical mechanical calculation tool

To perform molecular dynamics (MD) simulations, Monte Carlo (MC) simulations, quantum mechanical (QM) electronic structure calculations, or similar atomistic calculations, it is necessary to first construct and define the molecular system of interest. This involves creating an initial configuration of atoms, where MD and MC simulations require forcefield assignments, for example in the form of non-bonded, bonded, angular and dihedral potentials. Once simulations or calculations have been performed, large sets of data are available (often several gigabytes). These contain atomic trajectories and other relevant static or dynamic information, from which static and dynamic properties (such as density profiles, vibrational density of states and velocity autocorrelation functions) can be obtained through statistical mechanical calculations. MolTwister is an open source software platform that addresses the construction of molecular systems, basic 3D visualization of these, the generation of input files for selected MD packages, as well as calculation of properties from atomistic simulation data. It also contains a GPU accelerated MD simulator suited for smaller tasks such as molecular thermalization. The software package is written in C++14 and can be used as a basis for further development, where efforts have been made to make access to underlying functionality easy. Moreover, it supports Python, where scripts have access to the majority of program functionality. Program summaryProgram Title: MolTwisterCPC Library link to program files:https://doi.org/10.17632/ph6pwst9mn.1Developer's repository link:http://moltwister.comLicensing provisions: GPLv3 or laterProgramming language: C++14Nature of problem: Lack of SW tools for easy construction of complex molecular systems, as well as tools enabling analysis of long molecular trajectories.Solution method: Developed command line tool with associated 3D view, where molecular construction commands can be executed in complex manners through a Python interface. The program also has commands for analysis of molecular data that process one molecular time-frame at a time, thus enabling analysis of long molecular trajectories. Moreover, the program structure is designed for easy implementation of new commands by other researchers.Additional comments including restrictions and unusual features: Supported file formats: xyz, pdb, mtt (MolTwister trajectory), dcd, xtc, qepos, python, script (MolTwister scripts)Supported MD engines: HOOMD-blue, LAMMPS, MolTwister internal (moldyn)

Mixed Precision Fermi-Operator Expansion on Tensor Cores from a Machine Learning Perspective.

We present a second-order recursive Fermi-operator expansion scheme using mixed precision floating point operations to perform electronic structure calculations using tensor core units. A performance of over 100 teraFLOPs is achieved for half-precision floating point operations on Nvidia's A100 tensor core units. The second-order recursive Fermi-operator scheme is formulated in terms of a generalized, differentiable deep neural network structure, which solves the quantum mechanical electronic structure problem. We demonstrate how this network can be accelerated by optimizing the weight and bias values to substantially reduce the number of layers required for convergence. We also show how this machine learning approach can be used to optimize the coefficients of the recursive Fermi-operator expansion to accurately represent the fractional occupation numbers of the electronic states at finite temperatures.

Correcting Rate Constants from Anharmonic Molecular Dynamics for Quantum Effects.

Anharmonicity can greatly affect rate constants. One or even several orders of magnitude of deviation are found for obtaining rate constants using the standard rigid-rotor harmonic-oscillator model. In turn, reactive molecular dynamics (MD) simulations are a powerful way to explore chemical reaction networks and calculate rate constants from the fully anharmonic potential energy surface. However, the classical nature of the dynamics and the required numerical efficiency of the force field limit the accuracy of the resulting kinetics. We combine the best of both worlds by presenting an approximation that pairs anharmonic information intrinsic to classical MD with high-accuracy energies and frequencies from quantum-mechanical electronic structure calculations. The proposed scheme is applied to hydrogen abstractions in the methane system, which allows for the benchmarking of rate constants corrected by our approach against experimental rate constants. This comparison reveals a standard deviation of factor 2.6. Two archetypes of possible failure are identified in the course of a detailed investigation of the CH3• + H• → CH22• + H2 reaction. From this follows the application range of the method, within which the method shows a standard deviation of factor 2.1. The computational efficiency and beneficial scaling of the method allow for application to larger systems, as shown for hydrogen abstraction from 2-butanone by HO2•.

SAMPL6 logP challenge: machine learning and quantum mechanical approaches.

Two different types of approaches: (a) approaches that combine quantitative structure activity relationships, quantum mechanical electronic structure methods, and machine-learning and, (b) electronic structure vertical solvation approaches, were used to predict the logP coefficients of 11 molecules as part of the SAMPL6 logP blind prediction challenge. Using electronic structures optimized with density functional theory (DFT), several molecular descriptors were calculated for each molecule, including van der Waals areas and volumes, HOMO/LUMO energies, dipole moments, polarizabilities, and electrophilic and nucleophilic superdelocalizabilities. A multilinear regression model and a partial least squares model were used to train a set of 97 molecules. As well, descriptors were generated using the molecular operating environment and used to create additional machine learning models. Electronic structure vertical solvation approaches considered include DFT and the domain-based local pair natural orbital methods combined with the solvated variant of the correlation consistent composite approach.

Proton-transfer in 1,1,3,3 tetramethyl guanidine by means of ultrasonic relaxation and Raman spectroscopies and molecular orbital calculations.

In this work, we report on the structure and dynamics of the 1,1,3,3 tetramethyl guanidine (TMG) aqueous solutions in a wide concentration and temperature range by combining vibrational and ultrasonic spectroscopies. The experimental Raman spectra have been compared with the corresponding spectra obtained by ab initio quantum mechanical and density functional theory electronic structure calculations. This comparison indicated that only a single mechanism occurs when dissolving TMG in water and this is the proton transfer reaction, while the formation of byproducts during hydrolysis of TMG is dubious. This observation is further supported by the concentration dependence of the Raman spectra. The analysis of the ultrasonic relaxation data also revealed that the system exhibits a single relaxation process associated with this proton transfer reaction. It has been also observed that both relaxation amplitude and frequency exhibit a clear monotonous increase with increasing amine concentration in the solutions supporting the concept of the proton transfer reaction. The corresponding activation enthalpy was estimated directly from the temperature dependence of the acoustic data and found equal to ΔH*=5.56±0.34kcal/mol, which seems to be reasonable for hydrogen-bond formation. Furthermore, the concentration dependence of the acoustic parameters and kinematic viscosity data has been used as a probe for the molecular association in these solutions. The results have been discussed in relation to the ability or inability of water molecules to form stable clathrates after the addition of amine molecules in the solutions.

Linear programming applied to polarized Raman spectroscopy for elucidating molecular structure at surfaces

We present a framework for using linear programming to solve a challenging problem in surface science, the elucidation of the structure and composition of adsorbed molecules from a mixture, using simulated data from polarized Raman experiments. In the past, methods applied in order to interpret such spectroscopic information were combinatorial approaches that are limited in scalability or accuracy. Quantum mechanical electronic structure calculations yield the optical response of a single molecule, from which spectra of a mixture can be determined by appropriate weighting. Furthermore, spectral obtained in different beam polarizations provide projections of the signal in the laboratory frame. We demonstrate that linear programming is an ideal tool for utilizing all of this information in order to provide the sought structural picture.

Revealing the Sequence-Structure-Electronic Property Relation of Self-Assembling π-Conjugated Oligopeptides by Molecular and Quantum Mechanical Modeling.

Self-assembled nanoaggregates of π-conjugated synthetic peptides present a biocompatible and highly tunable alternative to silicon-based optical and electronic materials. Understanding the relationship between structural morphology and electronic properties of these assemblies is critical for understanding and controlling their mechanical, optical, and electronic responses. In this work, we combine all-atom classical molecular simulations with quantum mechanical electronic structure calculations to ascertain the sequence-structure-electronic property relationship within a family of Asp-X-X-quaterthiophene-X-X-Asp (DXX-OT4-XXD) oligopeptides in which X is one of the five amino acids {Ala, Phe, Gly, Ile, Val} ({A, F, G, I, V}). Molecular dynamics simulations reveal that smaller amino acid substituents (A, G) favor linear stacking within a peptide dimer, whereas larger groups (F, I, V) induce larger twist angles between the peptides. Density functional theory calculations on the dimer show the absorption spectrum to be dominated by transitions between carbon and sulfur p orbitals. Although the absorption spectrum is largely insensitive to the relative twist angle, the highest occupied molecular orbital strongly localizes onto one molecule within the dimer at large twist angles, impeding the efficiency of transport between molecules. Our results provide a fundamental understanding of the relation between peptide orientation and electronic structure and offer design precepts for rational engineering of these systems.

2-body:Many-body QM:QM study of structures, energetics, and vibrational frequencies for microhydrated halide ions

ABSTRACTThe optimised structures, electronic dissociation energies and harmonic vibrational frequencies of microhydrated halide ions have been probed with the MP2 and CCSD(T) quantum mechanical (QM) electronic structure techniques as well as with the 2-body:Many-body QM:QM method (2b:Mb CCSD(T):MP2) in which the 1- and 2-body interactions are computed at the CCSD(T) level of theory while all higher-order contributions are captured at the MP2 level of theory. In general, the accuracy and efficiency of the 2b:Mb approach previously demonstrated for neutral, homogeneous HF and HO clusters extends to these negatively charged heterogeneous systems. Results indicate the MP2 method struggles to provide a reliable description of the vibrational frequencies for hydrated fluoride clusters even though it performs well for the other hydrated halides (chloride, bromide and iodide). Consequently, the accuracy of the 2b:Mb CCSD(T):MP2 approach for fluoride clusters is poorer than the other halides, and MP2 is not a recommended low-level method for high-accuracy frequency analyses of hydrated fluoride complexes with the 2b:Mb method or related QM:QM approaches.

Development of reactive force fields using ab initio molecular dynamics simulation minimally biased to experimental data.

Incorporation of quantum mechanical electronic structure data is necessary to properly capture the physics of many chemical processes. Proton hopping in water, which involves rearrangement of chemical and hydrogen bonds, is one such example of an inherently quantum mechanical process. Standard ab initio molecular dynamics (AIMD) methods, however, do not yet accurately predict the structure of water and are therefore less than optimal for developing force fields. We have instead utilized a recently developed method which minimally biases AIMD simulations to match limited experimental data to develop novel multiscale reactive molecular dynamics (MS-RMD) force fields by using relative entropy minimization. In this paper, we present two new MS-RMD models using such a parameterization: one which employs water with harmonic internal vibrations and another which uses anharmonic water. We show that the newly developed MS-RMD models very closely reproduce the solvation structure of the hydrated excess proton in the target AIMD data. We also find that the use of anharmonic water increases proton hopping, thereby increasing the proton diffusion constant.

Quantum mechanical electronic structure calculation reveals orientation dependence of hydrogen bond energy in proteins.

Hydrogen bond plays a unique role in governing macromolecular interactions with exquisite specificity. These interactions govern the fundamental biological processes like protein folding, enzymatic catalysis, molecular recognition. Despite extensive research work, till date there is no proper report available about the hydrogen bond's energy surface with respect to its geometric parameters, directly derived from proteins. Herein, we have deciphered the potential energy landscape of hydrogen bond directly from the macromolecular coordinates obtained from Protein Data Bank using quantum mechanical electronic structure calculations. The findings unravel the hydrogen bonding energies of proteins in parametric space. These data can be used to understand the energies of such directional interactions involved in biological molecules. Quantitative characterization has also been performed using Shannon entropic calculations for atoms participating in hydrogen bond. Collectively, our results constitute an improved way of understanding hydrogen bond energies in case of proteins and complement the knowledge-based potential. Proteins 2017; 85:1046-1055. © 2017 Wiley Periodicals, Inc.

Validation of Density Functionals for Adsorption Energies on Transition Metal Surfaces.

The quantitative prediction of adsorption energies of radicals and molecules on surfaces is essential for the design and understanding of heterogeneous catalytic processes. A recent paper by Wellendorff et al. collected an experimental database of 39 reaction energies involving adsorption energies on transition metal surfaces that can be used as benchmarks for testing quantum mechanical electronic structure methods, and we compared the experimental data to Kohn-Sham density functional calculations with six exchange-correlation functionals. In this paper, we rearranged the data into two categories: open-shell radical adsorption reactions and closed-shell molecular adsorption reactions. We recalculated the adsorption energies with PBE, and we also calculated them with three functionals, M06-L, GAM, and MN15-L, that were not studied in the Wellendorff et al. paper; then we compared our results to the benchmark data. Of the nine functionals that have been compared to the databases, we find that BEEF-vdW, GAM, and RPBE perform best for the open-shell radical adsorption reactions, and MN15-L performs best for the closed-shell molecular adsorption, followed by BEEF-vdW and M06-L.

Adsorbate-Induced Changes in Magnetic Interactions in Fe2(dobdc) with Adsorbed Hydrocarbon Molecules

Metal–organic frameworks, which are a special case of coordination polymers, form a class of materials with numerous applications due to their high porosities and large internal surface areas and in some cases also due to paramagnetic metal ions. A family of such materials of general formula M2(dobdc) (where M is a divalent metal ion, and dobdc4– is 2,5-dioxido-1,4-benzenedicarboxylate) has attracted considerable attention for gas separation, catalysis, and magnetism. In this work, we explore the magnetic properties of a member of this family, Fe2(dobdc), both in the activated form (bare MOF after solvent removal) and when hydrocarbons are bound to the open coordination sites of the metal. We report quantum mechanical electronic structure calculations using both cluster models and periodic models, and we compare our results to previously reported theoretical studies. We find that hydrocarbon adsorption only mildly affects the isotropic couplings but that the isotropic magnetic couplings obtained with hybrid exchange–correlation functionals for the cluster models were found to be in good agreement with the experimental values, and the ones obtained using local exchange–correlation functionals with empirical Coulomb and exchange integrals for the periodic models were also found to agree well with experiments. Furthermore, local density functionals with empirical Coulomb and exchange integrals for periodic models are found to give good agreement with the experimental result that the adsorption of ethylene changes the magnetic ordering. We also used second-order n-electron valence state perturbation theory and contracted spin–orbit configuration interaction to study the role of adsorption on the crystal-field splitting of the quintet manifold and on the single-ion anisotropy of the iron center.

Theoretical prediction of the 5 V rechargeable Li ion battery using Li2CoMn3O8 as a cathode

Due to the rapid development of computers and cutting edge research methods, it is now possible to use a quantum mechanical electronic structure theory of solids to obtain, completely from first principles, the average voltage of a battery based on intercalation reaction energetics. The results of ab initio calculations by means of the full potential linearized augmented plane wave (FP-LAPW) method, using the computer code WIEN2k, for spin-polarized mixtures of (x = 0, 1, 2, 3 and 4), treating exchange and correlation effects within the generalized gradient approximation (GGA), are presented. The calculated average battery voltage for Li2CoMn3O8 cathode material is around 5 volts. The ab initio calculation result for a 5 V average battery voltage perfectly describes recently experimentally-synthesized LiCo0.5Mn1.5O4 battery cathode material, which showed a discharge plateau starting at around 5 V. The calculated average battery voltages for other x values in the mixture Li2CoxMn4−xO8 (x = 0, 2, 3 and 4) (3.95 V; 4.47 V; 4.19 V and 3.99 V) are considerably below 5 V.

Atomic Oxygen Recombination at Surface Defects on Reconstructed (0001) α-Quartz Exposed to Atomic and Molecular Oxygen

The surface chemistry of silica is strongly affected by the nature of chemically active sites (or defects) occurring on the surface. Here, we employ quantum mechanical electronic structure calculations to study an uncoordinated silicon defect, a non-bridging oxygen defect, and a peroxyl defect on the reconstructed (0001) surface of α-quartz. We characterized the spin states and energies of the defects, and calculated the reaction profiles for atomic oxygen recombination at the defects. We elucidated the diradical character by analyzing the low-lying excited states using multireference wave function methods. We show that the diradical defects consist of weakly coupled doublet radicals, and the atomic oxygen recombination can take place through a barrierless process at defects. We have delineated the recombination mechanism and computed the formation energy of the peroxyl and non-bridging oxygen defects. We found that key recombination reaction paths are barrierless. In addition, we characterize the electro...

Dynamical, spectroscopic and computational imaging of bond breaking in photodissociation: roaming and role of conical intersections.

Recent experimental and theoretical advances in the study of the dissociation of excited molecules are revealing unexpected mechanisms, when their outcomes are tackled by combining (i) space-time ion imaging of translational features, with (ii) spectroscopic probing of rotational and vibrational distributions; crucial is the assistance of (iii) the quantum chemistry of structural investigations of rearrangements of chemical bonds, and of (iv) the simulations of molecular dynamics to follow the evolution of selective bond stretching and breaking. Here we present results of such an integrated approach to methyl formate, HCOOCH3, the simplest of esters; the main focus is on the rotovibrationally excited CO (v=1) product and in general on the energy distribution in the fragments. Previous laser studies of dissociation into CO and CH3OH at a sequence of various wavelengths discovered signatures of a roaming mechanism by the late arrival of CO (v=0) products in time-of-flight ion imaging. Subsequent detailed investigations as a function of excitation energy provided the assessment of the threshold, which opens for triple breakdown into CO and further fragments H and CH3O, as spectroscopically characterized by ion imaging and FTIR respectively. Accompanying quantum mechanical electronic structure calculations and classical molecular dynamics simulations clarify the origin of these fragments through "roaming" pathways involving incipient radical intermediates at energies below the triple fragmentation threshold: a specific role is played by nonadiabatic transitions at a conical intersection between ground and excited states; alternative pathways focalize our attention to regions of the potential energy surfaces other than those in the neighbourhoods of saddle points along minimum energy paths: eventually this leads us to look for avenues in reaction kinetics beyond those of venerable transition state theories.

Peptide bond formation via glycine condensation in the gas phase.

Four unique gas phase mechanisms for peptide bond formation between two glycine molecules have been mapped out with quantum mechanical electronic structure methods. Both concerted and stepwise mechanisms, each leading to a cis and trans glycylglycine product (four mechanisms total), were examined with the B3LYP and MP2 methods and Gaussian atomic orbital basis sets as large as aug-cc-pVTZ. Electronic energies of the stationary points along the reaction pathways were also computed with explicitly correlated MP2-F12 and CCSD(T)-F12 methods. The CCSD(T)-F12 computations indicate that the electronic barriers to peptide bond formation are similar for all four mechanisms (ca. 32-39 kcal mol(-1) relative to two isolated glycine fragments). The smallest barrier (32 kcal mol(-1)) is associated with the lone transition state for the concerted mechanism leading to the formation of a trans peptide bond, whereas the largest barrier (39 kcal mol(-1)) was encountered along the concerted pathway leading to the cis configuration of the glycylglycine dipeptide. Two significant barriers are encountered for the stepwise mechanisms. For both the cis and trans pathways, the early electronic barrier is 36 kcal mol(-1) and the subsequent barrier is approximately 1 kcal mol(-1) lower. A host of intermediates and transition states lie between these two barriers, but they all have very small relative electronic energies (ca. ± 4 kcal mol(-1)). The isolated cis products (glycylglycine + H2O) are virtually isoenergetic with the isolated reactants (within -1 kcal mol(-1)), whereas the trans products are about 5 kcal mol(-1) lower in energy. In both products, however, the water can hydrogen bond to the dipeptide and lower the energy by roughly 5-9 kcal mol(-1). This study indicates that the concerted process leading to a trans configuration about the peptide bond is marginally favored both thermodynamically (exothermic by ca. 5 kcal mol(-1)) and kinetically (barrier height ≈ 32 kcal mol(-1)) according to the CCSD(T)-F12/haTZ electronic energies. The other pathways have slightly larger barrier heights (by 4-8 kcal mol(-1)).

Photoelectron spectroscopic and computational study of hydrated pyrimidine anions.

The stabilization of the pyrimidine anion by the addition of water molecules is studied experimentally using photoelectron spectroscopy of mass-selected hydrated pyrimidine clusters and computationally using quantum-mechanical electronic structure theory. Although the pyrimidine molecular anion is not observed experimentally, the addition of a single water molecule is sufficient to impart a positive electron affinity. The sequential hydration data have been used to extrapolate to -0.22 eV for the electron affinity of neutral pyrimidine, which agrees very well with previous observations. These results for pyrimidine are consistent with previous studies of the hydrated cluster anions of uridine, cytidine, thymine, adenine, uracil, and naphthalene. This commonality suggests a universal effect of sequential hydration on the electron affinity of similar molecules.

Dehydroergosterol as an analogue for cholesterol: why it mimics cholesterol so well-or does it?

Although dehydroergosterol (DHE) is one of the most commonly used cholesterol (CHOL) reporters, it has remained unclear why it performs well compared with most other CHOL analogues and what its possible limitations are. We present a comprehensive study of the properties of DHE using a combination of time-resolved fluorescence spectroscopy, quantum-mechanical electronic structure computations, and classical atomistic molecular dynamics simulations. We first establish that DHE mimics CHOL behavior, as previous studies have suggested, and then move on to elucidate and discuss the particular properties that render DHE so superior. We found that the main reason why DHE mimics CHOL so well is due to its ability to stand upright in a membrane in a manner that is almost identical to that of CHOL. The minor difference in how DHE and CHOL tilt with respect to membrane normal has only faint effects on structural membrane properties, and even the lateral pressure profiles of model membranes with CHOL or DHE are almost identical. These results suggest that the mechanical/elastic effects of DHE on the function of mechanically sensitive membrane proteins are not substantially different from those of CHOL. Our study highlights similar dynamical behavior of CHOL and DHE, which implies that DHE can mimic CHOL in processes with free energies close to the thermal energy.

Quantum mechanical fragment methods based on partitioning atoms or partitioning coordinates.

Conspectus The development of more efficient and more accurate ways to represent reactive potential energy surfaces is a requirement for extending the simulation of large systems to more complex systems, longer-time dynamical processes, and more complete statistical mechanical sampling. One way to treat large systems is by direct dynamics fragment methods. Another way is by fitting system-specific analytic potential energy functions with methods adapted to large systems. Here we consider both approaches. First we consider three fragment methods that allow a given monomer to appear in more than one fragment. The first two approaches are the electrostatically embedded many-body (EE-MB) expansion and the electrostatically embedded many-body expansion of the correlation energy (EE-MB-CE), which we have shown to yield quite accurate results even when one restricts the calculations to include only electrostatically embedded dimers. The third fragment method is the electrostatically embedded molecular tailoring approach (EE-MTA), which is more flexible than EE-MB and EE-MB-CE. We show that electrostatic embedding greatly improves the accuracy of these approaches compared with the original unembedded approaches. Quantum mechanical fragment methods share with combined quantum mechanical/molecular mechanical (QM/MM) methods the need to treat a quantum mechanical fragment in the presence of the rest of the system, which is especially challenging for those parts of the rest of the system that are close to the boundary of the quantum mechanical fragment. This is a delicate matter even for fragments that are not covalently bonded to the rest of the system, but it becomes even more difficult when the boundary of the quantum mechanical fragment cuts a bond. We have developed a suite of methods for more realistically treating interactions across such boundaries. These methods include redistributing and balancing the external partial atomic charges and the use of tuned fluorine atoms for capping dangling bonds, and we have shown that they can greatly improve the accuracy. Finally we present a new approach that goes beyond QM/MM by combining the convenience of molecular mechanics with the accuracy of fitting a potential function to electronic structure calculations on a specific system. To make the latter practical for systems with a large number of degrees of freedom, we developed a method to interpolate between local internal-coordinate fits to the potential energy. A key issue for the application to large systems is that rather than assigning the atoms or monomers to fragments, we assign the internal coordinates to reaction, secondary, and tertiary sets. Thus, we make a partition in coordinate space rather than atom space. Fits to the local dependence of the potential energy on tertiary coordinates are arrayed along a preselected reaction coordinate at a sequence of geometries called anchor points; the potential energy function is called an anchor points reactive potential. Electrostatically embedded fragment methods and the anchor points reactive potential, because they are based on treating an entire system by quantum mechanical electronic structure methods but are affordable for large and complex systems, have the potential to open new areas for accurate simulations where combined QM/MM methods are inadequate.

Quest for a universal density functional: the accuracy of density functionals across a broad spectrum of databases in chemistry and physics

Kohn-Sham density functional theory is in principle an exact formulation of quantum mechanical electronic structure theory, but in practice we have to rely on approximate exchange-correlation (xc) functionals. The objective of our work has been to design an xc functional with broad accuracy across as wide an expanse of chemistry and physics as possible, leading--as a long-range goal--to a functional with good accuracy for all problems, i.e. a universal functional. To guide our path towards that goal and to measure our progress, we have developed-building on earlier work of our group-a set of databases of reference data for a variety of energetic and structural properties in chemistry and physics. These databases include energies of molecular processes, such as atomization, complexation, proton addition and ionization; they also include molecular geometries and solid-state lattice constants, chemical reaction barrier heights, and cohesive energies and band gaps of solids. For this paper, we gather many of these databases into four comprehensive databases, two with 384 energetic data for chemistry and solid-state physics and another two with 68 structural data for chemistry and solid-state physics, and we test two wave function methods and 77 density functionals (12 Minnesota meta functionals and 65 others) in a consistent way across this same broad set of data. We especially highlight the Minnesota density functionals, but the results have broader implications in that one may see the successes and failures of many kinds of density functionals when they are all applied to the same data. Therefore, the results provide a status report on the quest for a universal functional.