A specific MNDO parameterization for water.

A modified neglect of differential overlap has been parameterized specifically for water and its oligomers with the addition of polarization functions on both hydrogen and oxygen, Feynman dispersion, and a slight modification of the treatment of the hydrogen nucleus. The results show that it is possible to easily obtain good geometries and energies for hydrogen-bonded water aggregates. Data from the Benchmark Energy and Geometry Database water-cluster database were used to parameterize the new Hamiltonian for water clusters from the dimer to the decamer using MP2/aug-cc-pVDZ optimized geometries and CCSD(T)/CBS oligomerization energies. Seventy five oligomerization and rearrangement energies derived from the parameterization data are reproduced with a root mean-square error (RMSE) of 0.79 kcal mol-1 and the geometries of 38 oligomers with an RMSE of 0.17 Å. Interestingly, the Feynman dispersion term adopts a role different from that intended and tunes the atomic polarizability. The implications of these results in terms of future dedicated neglect of diatomic differential overlap Hamiltonians and those that use force-field-like atom types are discussed.

An Electrically Conducting Three-Dimensional Iron-Catecholate Porous Framework.

We report the synthesis of a unique cubic metal-organic framework (MOF), Fe-HHTP-MOF, comprising hexahydroxytriphenylene (HHTP) supertetrahedral units and FeIII ions, arranged in a diamond topology. The MOF is synthesized under solvothermal conditions, yielding a highly crystalline, deep black powder, with crystallites of 300-500 nm size and tetrahedral morphology. Nitrogen sorption analysis indicates a highly porous material with a surface area exceeding 1400 m2 g-1 . Furthermore, Fe-HHTP-MOF shows broadband absorption from 475 up to 1900 nm with excellent absorption capability of 98.5 % of the incoming light over the visible spectral region. Electrical conductivity measurements of pressed pellets reveal a high intrinsic electrical conductivity of up to 10-3  S cm-1 . Quantum mechanical calculations predict Fe-HHTP-MOF to be an efficient electron conductor, exhibiting continuous charge-carrier pathways throughout the structure.

Energy Efficient Ultrahigh Flux Separation of Oily Pollutants from Water with Superhydrophilic Nanoscale Metal-Organic Framework Architectures.

The rising demand for clean water for a growing and increasingly urban global population is one of the most urgent issues of our time. Here, we introduce the synthesis of a unique nanoscale architecture of pillar-like Co-CAT-1 metal-organic framework (MOF) crystallites on gold-coated woven stainless steel meshes with large, 50 µm apertures. These nanostructured mesh surfaces feature superhydrophilic and underwater superoleophobic wetting properties, allowing for gravity-driven, highly efficient oil-water separation featuring water fluxes of up to nearly one million L m-2 h-1 . Water physisorption experiments reveal the hydrophilic nature of Co-CAT-1 with a total water vapor uptake at room temperature of 470 cm3 g-1 . Semiempirical molecular orbital calculations shed light on water affinity of the inner and outer pore surfaces. The MOF-based membranes enable high separation efficiencies for a number of liquids tested, including the notorious water pollutant, crude oil, affording chemical oxygen demand (COD) concentrations below 25 mg L-1 of the effluent. Our results demonstrate the great impact of suitable nanoscale surface architectures as a means of encoding on-surface extreme wetting properties, yielding energy-efficient water-selective large-aperture membranes.

EMPIRE: a highly parallel semiempirical molecular orbital program: 3: Born-Oppenheimer molecular dynamics.

Direct NDDO-based Born-Oppenheimer molecular dynamics (MD) have been implemented in the semiempirical molecular orbital program EMPIRE. Fully quantum mechanical MD simulations on unprecedented time and length scales are possible, since the calculation of self-consistent wavefunctions and gradients is performed in a massively parallel manner. MD simulations can be performed in the NVE and NVT ensembles, using either deterministic (Berendsen) or stochastic (Langevin) thermostats. Furthermore, dynamics for condensed-phase systems can be performed under periodic boundary conditions. We show three exemplary applications: the dynamics of molecular reorganization upon ionization, long timescale dynamics of an endohedral fullerene, and calculation of the vibrational spectrum of a nanoparticle consisting of more than eight hundred atoms. Graphical AbstractA snapshot from an MNDO-H simulation of NH4+@C60 at 4000 K shortly before a proton crosses the fullerene wall to give NH3@C60H+.

Propagation of Holes and Electrons in Metal-Organic Frameworks.

Charge transport in two zinc metal-organic frameworks (MOFs) has been investigated using periodic semiempirical molecular orbital calculations with the AM1* Hamiltonian. Restricted Hartree-Fock calculations underestimate the band gap using Koopmans theorem (ca. 2 eV compared to the experimental value of 2.8 eV). However, it almost doubles when the constraint on the wave function to remain spin-restricted is removed and the energies of the UHF Natural Orbitals are used. Charge-transport simulations using propagation of the electron- or hole-density in imaginary time allow charge-transport paths and mechanisms to be determined. The calculated relative mobilities in the directions of the three crystal axes agree with experimental expectations, but the absolute values are not reliable using the current technique. Hole-mobility along the crystal c-axis (along the metal stacks) is found to be 13 times higher in the zinc MOF with anthracene linker (Zn-ANMOF-74) than in the other directions, whereas the factor is far smaller (1.7) for electron mobility. Directional preferences are far less distinct in the equivalent structure with phenyl linkers (Zn-MOF-74). The imaginary-time simulation technique does not give quantitative mobilities. The simulations reveal a change in mechanism between the different directions: Coherent polaron migration is observed along the stacks but tunneling hops between them.

Size-Dependent Local Ordering in Melanin Aggregates and Its Implication on Optical Properties.

We present atomic scale models of differently sized eumelanin nanoaggregates from molecular dynamics simulations combined with a simulated annealing procedure. The analysis reveals the formation of secondary structures due to π-stacking on one hand, but on the other hand a broad distribution of stack geometries in terms of stack size, horizontal displacement angles, and relative torsion angles. The displacement angle distribution, which is a measure of the occurrence of zigzag and linear stacking motives, respectively, strongly depends on the aggregate size-and is hence controlled by the interplay of surface and bulk energy terms. Semiempirical spectra calculations of small stacks (up to five protomolecules) reveal a strong dependence on the precise stack structure and allow for a direct structure-property correlation. The observed spectral shifts result in an overall spectral broadening and, hence, further support the geometric disorder model, which complements the chemical disorder model in the interpretation of eumelanin's monotonically increasing broad-band absorption.

Correction to: The Feynman dispersion correction for MNDO extended to F, Cl, Br and I.

A small coding error in the development version of EMPIRE led to some inconsistencies in the above article. They are corrected in this erratum.

The Feynman dispersion correction for MNDO extended to F, Cl, Br and I.

The recently introduced "Feynman" dispersion correction for MNDO (MNDO-F) has been extended to include the elements fluorine, chlorine, bromine and iodine and the original parameterization for hydrogen, carbon, nitrogen and oxygen improved by allowing individual damping radii for the elements. MNDO-F gives a root-mean-square deviation to reference interaction energies of 0.35 kcal mol-1 for the complete parameterization dataset of H, C, N, O, F, Cl, Br and I containing compounds. Graphical Abstract The electrostatic potential at the 0.001 a.u. isodensity surface of the π-complex between benzene and 1,3,5-triodobenzene calculated at the MNDO-F optimized geometry.

Solvatochromic covalent organic frameworks.

Covalent organic frameworks (COFs) are an emerging class of highly tuneable crystalline, porous materials. Here we report the first COFs that change their electronic structure reversibly depending on the surrounding atmosphere. These COFs can act as solid-state supramolecular solvatochromic sensors that show a strong colour change when exposed to humidity or solvent vapours, dependent on vapour concentration and solvent polarity. The excellent accessibility of the pores in vertically oriented films results in ultrafast response times below 200 ms, outperforming commercially available humidity sensors by more than an order of magnitude. Employing a solvatochromic COF film as a vapour-sensitive light filter, we demonstrate a fast humidity sensor with full reversibility and stability over at least 4000 cycles. Considering their immense chemical diversity and modular design, COFs with fine-tuned solvatochromic properties could broaden the range of possible applications for these materials in sensing and optoelectronics.

The hpCADD NDDO Hamiltonian: Parametrization.

A neglect of diatomic differential overlap (NDDO) Hamiltonian has been parametrized as an electronic component of a polarizable force field. Coulomb and exchange potentials derived directly from the NDDO Hamiltonian in principle can be used with classical potentials, thus forming the basis for a new generation of efficiently applicable multipolar polarizable force fields. The new hpCADD Hamiltonian uses force-field-like atom types and reproduces the electrostatic properties (dipole moment, molecular electrostatic potential) and Koopmans' theorem ionization potentials closely, as demonstrated for a large training set and an independent test set of small molecules. The Hamiltonian is not intended to reproduce geometries or total energies well, as these will be controlled by the classical force-field potentials. In order to establish the hpCADD Hamiltonian as an electronic component in force-field-based calculations, we tested its performance in combination with the 3D reference interaction site model (3D RISM) for aqueous solutions. Comparison of the resulting solvation free energies for the training and test sets to atomic charges derived from standard procedures, exact solute-solvent electrostatics based on high-level quantum-chemical reference data, and established semiempirical Hamiltonians demonstrates the advantages of the hpCADD parametrization.

Self-consistent field convergence for proteins: a comparison of full and localized-molecular-orbital schemes.

Proteins in the gas phase present an extreme (and unrealistic) challenge for self-consistent-field iteration schemes because their ionized groups are very strong electron donors or acceptors, depending on their formal charge. This means that gas-phase proteins have a very small band gap but that their frontier orbitals are localized compared to "normal" conjugated semiconductors. The frontier orbitals are thus likely to be separated in space so that they are close to, but not quite, orthogonal during the SCF iterations. We report full SCF calculations using the massively parallel EMPIRE code and linear scaling localized-molecular-orbital (LMO) calculations using Mopac2009. The LMO procedure can lead to artificially over-polarized wavefunctions in gas-phase proteins. The full SCF iteration procedure can be very slow to converge because many cycles are needed to overcome the over-polarization by inductive charge shifts. Example molecules have been constructed to demonstrate this behavior. The two approaches give identical results if solvent effects are included.

Quantum mechanics-based properties for 3D-QSAR.

We have used a set of four local properties based on semiempirical molecular orbital calculations (electron density (ρ), hydrogen bond donor field (HDF), hydrogen bond acceptor field (HAF), and molecular lipophilicity potential (MLP)) for 3D-QSAR studies to overcome the limitations of the current force field-based molecular interaction fields (MIFs). These properties can be calculated rapidly and are thus amenable to high-throughput industrial applications. Their statistical performance was compared with that of conventional 3D-QSAR approaches using nine data sets (angiotensin converting enzyme inhibitors (ACE), acetylcholinesterase inhibitors (AchE), benzodiazepine receptor ligands (BZR), cyclooxygenase-2 inhibitors (COX2), dihydrofolate reductase inhibitors (DHFR), glycogen phosphorylase b inhibitors (GPB), thermolysin inhibitors (THER), thrombin inhibitors (THR), and serine protease factor Xa inhibitors (fXa)). The 3D-QSAR models generated were tested thoroughly for robustness and predictive ability. The average performance of the quantum mechanical molecular interaction field (QM-MIF) models for the nine data sets is better than that of the conventional force field-based MIFs. In the individual data sets, the QM-MIF models always perform better than, or as well as, the conventional approaches. It is particularly encouraging that the relative performance of the QM-MIF models improves in the external validation. In addition, the models generated showed statistical stability with respect to model building procedure variations such as grid spacing size and grid orientation. QM-MIF contour maps reproduce the features important for ligand binding for the example data set (factor Xa inhibitors), demonstrating the intuitive chemical interpretability of QM-MIFs.

Improving the charge transport in self-assembled monolayer field-effect transistors: from theory to devices.

A three-pronged approach has been used to design rational improvements in self-assembled monolayer field-effect transistors: classical molecular dynamics (MD) simulations to investigate atomistic structure, large-scale quantum mechanical (QM) calculations for electronic properties, and device fabrication and characterization as the ultimate goal. The MD simulations reveal the effect of using two-component monolayers to achieve intact dielectric insulating layers and a well-defined semiconductor channel. The QM calculations identify improved conduction paths in the monolayers that consist of an optimum mixing ratio of the components. These results have been used both to confirm the predictions of the calculations and to optimize real devices. Monolayers were characterized with X-ray reflectivity measurements and by electronic characterization of complete devices.

Predicting the sites and energies of noncovalent intermolecular interactions using local properties.

Feed-forward artificial neural nets have been used to recognize H-bond donor and acceptor sites on drug-like molecules based on local properties (electron density, molecular electrostatic potential and local ionization energy, electron affinity, and polarizability) calculated at grid points around the molecule. Interaction energies for training were obtained from B97-D and ωB97X-D/aug-cc-pVDZ density-functional theory calculations on a series of model central molecules and H-bond acceptor and donor probes constrained to the grid points used for training. The resulting models provide maps of both classical and unusual H- and halogen-bonding sites. Note that these reactions result even though only classical H-bond donors and acceptors were used as probes around the central molecules. Some examples demonstrate the ability of the models to take the electronics of the central molecule into consideration and to provide semiquantitative estimates of interaction energies at low computational cost.

Polarization-induced σ-holes and hydrogen bonding.

The strong collinear polarizability of the A-H bond in A-H···B hydrogen bonds is shown to lead to an enhanced σ-hole on the donor hydrogen atom and hence to stronger hydrogen bonding. This effect helps to explain the directionality of hydrogen bonds, the well known cooperative effect in hydrogen bonding, and the occurrence of blue-shifting. The latter results when significant additional electron density is shifted into the A-H bonding region by the polarization effect. The shift in the A-H stretching frequency is shown to depend essentially linearly on the calculated atomic charge on the donor hydrogen for all donors in which A belongs to the same row of the periodic table. A further result of the polarization effect, which is also expected for other σ-hole bonds, is that the strength of the non-covalent interaction depends strongly on external electric fields.

CypScore: Quantitative prediction of reactivity toward cytochromes P450 based on semiempirical molecular orbital theory.

CypScore is an in silico approach for predicting the likely sites of cytochrome P450-mediated metabolism of druglike organic molecules. It consists of multiple models for the most important P450 oxidation reactions such as aliphatic hydroxylation, N-dealkylation, O-dealkylation, aromatic hydroxylation, double-bond oxidation, N-oxidation, and S-oxidation. Each of these models is based on atomic reactivity descriptors derived from surface-based properties calculated with ParaSurf and based on AM1 semiempirical molecular orbital theory. The models were trained with data derived from Bayer Schering Pharma's in-house MajorMetabolite Database with more than 2300 transformations and more than 800 molecules collected from the primary literature. The models have been balanced to allow the treatment of relative intramolecular, intra-chemotype, and inter-chemotype reactivities of the labile sites toward oxidation. The models were evaluated with promising hit rates on three public datasets of varying quality in the annotation of the experimental positions. For 39 well-characterized compounds from 14 in-house lead optimization programs, we could detect at least one major metabolite for the three highest-ranked positions in 87 % of the compounds and overall more than 62 % of all major metabolites, with promising true- to false-positive ratios of 0.9.

The effect of a complexed lithium cation on a norcarane-based radical clock.

Density-functional theory (DFT) and ab initio calculations have been used to investigate the effect of a complexed lithium cation on the radical-clock rearrangement of the 2-norcaranyl radical to the 3-cyclohexenylmethyl radical. As found earlier for ring-closing radical clocks, complexation with a metal ion leads to a significant lowering of the barrier to rearrangement. DFT calculations on a model for the norcaranyl clock in cytochrome P450 confirm the two-state reactivity proposal of Shaik et al. and indicate that the porphyrin exerts little or no electrostatic effect on the rearrangement barrier.