Quantum computing in pharma: A multilayer embedding approach for near future applications.

Quantum computers are special purpose machines that are expected to be particularly useful in simulating strongly correlated chemical systems. The quantum computer excels at treating a moderate number of orbitals within an active space in a fully quantum mechanical manner. We present a quantum phase estimation calculation on F2 in a (2,2) active space on Rigetti's Aspen-11 QPU. While this is a promising start, it also underlines the need for carefully selecting the orbital spaces treated by the quantum computer. In this work, a scheme for selecting such an active space automatically is described and simulated results obtained using both the quantum phase estimation (QPE) and variational quantum eigensolver (VQE) algorithms are presented and combined with a subtractive method to enable accurate description of the environment. The active occupied space is selected from orbitals localized on the chemically relevant fragment of the molecule, while the corresponding virtual space is chosen based on the magnitude of interactions with the occupied space calculated from perturbation theory. This protocol is then applied to two chemical systems of pharmaceutical relevance: the enzyme [Fe] hydrogenase and the photosenzitizer temoporfin. While the sizes of the active spaces currently amenable to a quantum computational treatment are not enough to demonstrate quantum advantage, the procedure outlined here is applicable to any active space size, including those that are outside the reach of classical computation.

Unraveling individual host-guest interactions in molecular recognition from first principles quantum mechanics: Insights into the nature of nicotinic acetylcholine receptor agonist binding.

Drug binding to a protein target is governed by a complex pattern of noncovalent interactions between the ligand and the residues in the protein's binding pocket. Here we introduce a generally applicable, parameter-free, computational method that allows for the identification, quantification, and analysis of the key ligand-residue interactions responsible for molecular recognition. Our strategy relies on Local Energy Decomposition analysis at the "gold-standard" coupled cluster DLPNO-CCSD(T) level. In the study case shown in this paper, nicotine and imidacloprid binding to the nicotinic acetylcholine receptor, our approach provides new insights into how individual amino acids in the active site determine sensitivity and selectivity of the ligands, extending and refining classical pharmacophore hypotheses. By inference, the method is applicable to any kind of host/guest interactions with potential applications in industrial biocatalysis and protein engineering.

NAMD goes quantum: an integrative suite for hybrid simulations.

Hybrid methods that combine quantum mechanics (QM) and molecular mechanics (MM) can be applied to studies of reaction mechanisms in locations ranging from active sites of small enzymes to multiple sites in large bioenergetic complexes. By combining the widely used molecular dynamics and visualization programs NAMD and VMD with the quantum chemistry packages ORCA and MOPAC, we created an integrated, comprehensive, customizable, and easy-to-use suite (http://www.ks.uiuc.edu/Research/qmmm). Through the QwikMD interface, setup, execution, visualization, and analysis are streamlined for all levels of expertise.

A binding mode hypothesis for prothioconazole binding to CYP51 derived from first principles quantum chemistry.

In order to assess safety and efficacy of small molecule drugs as well as agrochemicals, it is key to understanding the nature of protein-ligand interaction on an atomistic level. Prothioconazole (PTZ), although commonly considered to be an azole-like inhibitor of sterol 14-α demethylase (CYP51), differs from classical azoles with respect to how it binds its target. The available evidence is only indirect, as crystallographic elucidation of CYP51 complexed with PTZ have not yet been successful. We derive a binding mode hypothesis for PTZ binding its target, compare to DPZ, a triazole-type metabolite of PTZ, and set our findings into context of its biochemistry and spectroscopy. Quantum Theory of Atoms in Molecules (QTAIM) analysis of computed DFT electron densities is used to qualitatively understand the topology of binding, revealing significant differences of how R- and S-enantiomers are binding and, in particular, how the thiozolinthione head of PTZ binds to heme compared to DPZ's triazole head. The difference of binding enthalpy is calculated at coupled cluster (DLPNO-CCSD(T)) level of theory, and we find that DPZ binds stronger to CYP51 than PTZ by more than ΔH ~ 11 kcal/mol.

Addressing the System-Size Dependence of the Local Approximation Error in Coupled-Cluster Calculations.

Over the last two decades, the local approximation has been successfully used to extend the range of applicability of the "gold standard" singles and doubles coupled-cluster method with perturbative triples CCSD(T) to systems with hundreds of atoms. The local approximation error grows in absolute value with the increasing system size, i.e., by increasing the number of electron pairs in the system. In this study, we demonstrate that the recently introduced two-point extrapolation scheme for approaching the complete pair natural orbital (PNOs) space limit in domain-based pair natural orbital CCSD(T) calculations drastically reduces the dependence of the error on the system size, thus opening up unprecedented opportunities for the calculation of benchmark quality relative energies for large systems.

The ORCA quantum chemistry program package.

In this contribution to the special software-centered issue, the ORCA program package is described. We start with a short historical perspective of how the project began and go on to discuss its current feature set. ORCA has grown into a rather comprehensive general-purpose package for theoretical research in all areas of chemistry and many neighboring disciplines such as materials sciences and biochemistry. ORCA features density functional theory, a range of wavefunction based correlation methods, semi-empirical methods, and even force-field methods. A range of solvation and embedding models is featured as well as a complete intrinsic to ORCA quantum mechanics/molecular mechanics engine. A specialty of ORCA always has been a focus on transition metals and spectroscopy as well as a focus on applicability of the implemented methods to "real-life" chemical applications involving systems with a few hundred atoms. In addition to being efficient, user friendly, and, to the largest extent possible, platform independent, ORCA features a number of methods that are either unique to ORCA or have been first implemented in the course of the ORCA development. Next to a range of spectroscopic and magnetic properties, the linear- or low-order single- and multi-reference local correlation methods based on pair natural orbitals (domain based local pair natural orbital methods) should be mentioned here. Consequently, ORCA is a widely used program in various areas of chemistry and spectroscopy with a current user base of over 22 000 registered users in academic research and in industry.

Linear scaling perturbative triples correction approximations for open-shell domain-based local pair natural orbital coupled cluster singles and doubles theory [DLPNO-CCSD(T0/T)].

The coupled cluster method with single-, double-, and perturbative triple excitations [CCSD(T)] is considered to be one of the most reliable quantum chemistry theories. However, the steep scaling of CCSD(T) has limited its application to small or medium-sized systems for a long time. In our previous work, the linear scaling domain based local pair natural orbital CCSD variant (DLPNO-CCSD) has been developed for closed-shell and open-shell. However, it is known from extensive benchmark studies that triple-excitation contributions are important to reach chemical accuracy. In the present work, two linear scaling (T) approximations for open-shell DLPNO-CCSD are implemented and compared: (a) an algorithm based on the semicanonical approximation, in which off-diagonal Fock matrix elements in the occupied space are neglected [referred to as DLPNO-(T0)]; and (b) an improved algorithm in which the triples amplitudes are computed iteratively [referred to as DLPNO-(T)]. This work is based on the previous open-shell DLPNO-CCSD algorithm [M. Saitow et al., J. Chem. Phys. 146, 164105 (2017)] as well as the iterative (T) correction for closed-shell systems [Y. Guo et al., J. Chem. Phys. 148, 011101 (2018)]. Our results show that the new open-shell perturbative corrections, DLPNO-(T0/T), can predict accurate absolute and relative correlation energies relative to the canonical reference calculations with the same basis set. The absolute energies from DLPNO-(T) are significantly more accurate than those of DLPNO-(T0). The additional computational effort of DLPNO-(T) relative to DLPNO-(T0) is a factor of 4 on average. We report calculations on systems with more than 4000 basis functions.

A near-linear scaling equation of motion coupled cluster method for ionized states.

In this work, a domain-based local pair natural orbital (DLPNO) version of the equation of motion coupled cluster theory with single and double excitations for ionization potentials (IP-EOM-CCSD) equations has been formulated and implemented. The method uses ground state localized occupied and pair natural virtual orbitals and applies the DLPNO machinery to arrive at a linear scaling implementation of the IP-EOM-CCSD method. The accuracy of the method is controllable using ground state truncation parameters. Using default thresholds, the method predicts ionization potential (IP) values with good accuracy (mean absolute error of 0.08 eV). We demonstrate that our code can be used to compute IP values for systems with more than 1000 atoms and 10 000 basis functions.

Communication: An improved linear scaling perturbative triples correction for the domain based local pair-natural orbital based singles and doubles coupled cluster method [DLPNO-CCSD(T)].

In this communication, an improved perturbative triples correction (T) algorithm for domain based local pair-natural orbital singles and doubles coupled cluster (DLPNO-CCSD) theory is reported. In our previous implementation, the semi-canonical approximation was used and linear scaling was achieved for both the DLPNO-CCSD and (T) parts of the calculation. In this work, we refer to this previous method as DLPNO-CCSD(T0) to emphasize the semi-canonical approximation. It is well-established that the DLPNO-CCSD method can predict very accurate absolute and relative energies with respect to the parent canonical CCSD method. However, the (T0) approximation may introduce significant errors in absolute energies as the triples correction grows up in magnitude. In the majority of cases, the relative energies from (T0) are as accurate as the canonical (T) results of themselves. Unfortunately, in rare cases and in particular for small gap systems, the (T0) approximation breaks down and relative energies show large deviations from the parent canonical CCSD(T) results. To address this problem, an iterative (T) algorithm based on the previous DLPNO-CCSD(T0) algorithm has been implemented [abbreviated here as DLPNO-CCSD(T)]. Using triples natural orbitals to represent the virtual spaces for triples amplitudes, storage bottlenecks are avoided. Various carefully designed approximations ease the computational burden such that overall, the increase in the DLPNO-(T) calculation time over DLPNO-(T0) only amounts to a factor of about two (depending on the basis set). Benchmark calculations for the GMTKN30 database show that compared to DLPNO-CCSD(T0), the errors in absolute energies are greatly reduced and relative energies are moderately improved. The particularly problematic case of cumulene chains of increasing lengths is also successfully addressed by DLPNO-CCSD(T).

Binding free energies in the SAMPL5 octa-acid host-guest challenge calculated with DFT-D3 and CCSD(T).

We have tried to calculate the free energy for the binding of six small ligands to two variants of the octa-acid deep cavitand host in the SAMPL5 blind challenge. We employed structures minimised with dispersion-corrected density-functional theory with small basis sets and energies were calculated using large basis sets. Solvation energies were calculated with continuum methods and thermostatistical corrections were obtained from frequencies calculated at the HF-3c level. Care was taken to minimise the effects of the flexibility of the host by keeping the complexes as symmetric and similar as possible. In some calculations, the large net charge of the host was reduced by removing the propionate and benzoate groups. In addition, the effect of a restricted molecular dynamics sampling of structures was tested. Finally, we tried to improve the energies by using the DLPNO-CCSD(T) approach. Unfortunately, results of quite poor quality were obtained, with no correlation to the experimental data, systematically too positive affinities (by ~50 kJ/mol) and a mean absolute error (after removal of the systematic error) of 11-16 kJ/mol. DLPNO-CCSD(T) did not improve the results, so the accuracy is not limited by the energy function. Instead, four likely sources of errors were identified: first, the minimised structures were often incorrect, owing to the omission of explicit solvent. They could be partly improved by performing the minimisations in a continuum solvent with four water molecules around the charged groups of the ligands. Second, some ligands could bind in several different conformations, requiring sampling of reasonable structures. Third, there is an indication the continuum-solvation model has problems to accurately describe the binding of both the negatively and positively charged guest molecules. Fourth, different methods to calculate the thermostatistical corrections gave results that differed by up to 30 kJ/mol and there is an indication that HF-3c overestimates the entropy term. In conclusion, it is a challenge to calculate binding affinities for this octa-acid system with quantum-mechanical methods.

Pair natural orbital and canonical coupled cluster reaction enthalpies involving light to heavy alkali and alkaline earth metals: the importance of sub-valence correlation.

In this work, we tested canonical and domain based pair natural orbital coupled cluster methods (CCSD(T) and DLPNO-CCSD(T), respectively) for a set of 32 ligand exchange and association/dissociation reaction enthalpies involving ionic complexes of Li, Be, Na, Mg, Ca, Sr, Ba and Pb(ii). Two strategies were investigated: in the former, only valence electrons were included in the correlation treatment, giving rise to the computationally very efficient FC (frozen core) approach; in the latter, all non-ECP electrons were included in the correlation treatment, giving rise to the AE (all electron) approach. Apart from reactions involving Li and Be, the FC approach resulted in non-homogeneous performance. The FC approach leads to very small errors (<2 kcal mol-1) for some reactions of Na, Mg, Ca, Sr, Ba and Pb, while for a few reactions of Ca and Ba deviations up to 40 kcal mol-1 have been obtained. Large errors are both due to artificial mixing of the core (sub-valence) orbitals of metals and the valence orbitals of oxygen and halogens in the molecular orbitals treated as core, and due to neglecting core-core and core-valence correlation effects. These large errors are reduced to a few kcal mol-1 if the AE approach is used or the sub-valence orbitals of metals are included in the correlation treatment. On the technical side, the CCSD(T) and DLPNO-CCSD(T) results differ by a fraction of kcal mol-1, indicating the latter method as the perfect choice when the CPU efficiency is essential. For completely black-box applications, as requested in catalysis or thermochemical calculations, we recommend the DLPNO-CCSD(T) method with all electrons that are not covered by effective core potentials included in the correlation treatment and correlation-consistent polarized core valence basis sets of cc-pwCVQZ(-PP) quality.

A new near-linear scaling, efficient and accurate, open-shell domain-based local pair natural orbital coupled cluster singles and doubles theory.

The Coupled-Cluster expansion, truncated after single and double excitations (CCSD), provides accurate and reliable molecular electronic wave functions and energies for many molecular systems around their equilibrium geometries. However, the high computational cost, which is well-known to scale as O(N6) with system size N, has limited its practical application to small systems consisting of not more than approximately 20-30 atoms. To overcome these limitations, low-order scaling approximations to CCSD have been intensively investigated over the past few years. In our previous work, we have shown that by combining the pair natural orbital (PNO) approach and the concept of orbital domains it is possible to achieve fully linear scaling CC implementations (DLPNO-CCSD and DLPNO-CCSD(T)) that recover around 99.9% of the total correlation energy [C. Riplinger et al., J. Chem. Phys. 144, 024109 (2016)]. The production level implementations of the DLPNO-CCSD and DLPNO-CCSD(T) methods were shown to be applicable to realistic systems composed of a few hundred atoms in a routine, black-box fashion on relatively modest hardware. In 2011, a reduced-scaling CCSD approach for high-spin open-shell unrestricted Hartree-Fock reference wave functions was proposed (UHF-LPNO-CCSD) [A. Hansen et al., J. Chem. Phys. 135, 214102 (2011)]. After a few years of experience with this method, a few shortcomings of UHF-LPNO-CCSD were noticed that required a redesign of the method, which is the subject of this paper. To this end, we employ the high-spin open-shell variant of the N-electron valence perturbation theory formalism to define the initial guess wave function, and consequently also the open-shell PNOs. The new PNO ansatz properly converges to the closed-shell limit since all truncations and approximations have been made in strict analogy to the closed-shell case. Furthermore, given the fact that the formalism uses a single set of orbitals, only a single PNO integral transformation is necessary, which offers large computational savings. We show that, with the default PNO truncation parameters, approximately 99.9% of the total CCSD correlation energy is recovered for open-shell species, which is comparable to the performance of the method for closed-shells. UHF-DLPNO-CCSD shows a linear scaling behavior for closed-shell systems, while linear to quadratic scaling is obtained for open-shell systems. The largest systems we have considered contain more than 500 atoms and feature more than 10 000 basis functions with a triple-ζ quality basis set.

SparseMaps-A systematic infrastructure for reduced scaling electronic structure methods. V. Linear scaling explicitly correlated coupled-cluster method with pair natural orbitals.

In this work, we present a linear scaling formulation of the coupled-cluster singles and doubles with perturbative inclusion of triples (CCSD(T)) and explicitly correlated geminals. The linear scaling implementation of all post-mean-field steps utilizes the SparseMaps formalism [P. Pinski et al., J. Chem. Phys. 143, 034108 (2015)]. Even for conservative truncation levels, the method rapidly reaches near-linear complexity in realistic basis sets, e.g., an effective scaling exponent of 1.49 was obtained for n-alkanes with up to 200 carbon atoms in a def2-TZVP basis set. The robustness of the method is benchmarked against the massively parallel implementation of the conventional explicitly correlated coupled-cluster for a 20-water cluster; the total dissociation energy of the cluster (∼186 kcal/mol) is affected by the reduced scaling approximations by only ∼0.4 kcal/mol. The reduced scaling explicitly correlated CCSD(T) method is used to examine the binding energies of several systems in the L7 benchmark data set of noncovalent interactions.

Stability of surface protons in pyridine-catalyzed CO2 reduction at p-GaP photoelectrodes.

Adsorbed protons that develop hydride character have been proposed to play a role in the mechanism of CO2 reduction catalyzed by pyridine on GaP photoelectrodes. Investigating their stability represents an important step towards vetting this mechanism. In this contribution, the relative stability of the adsorbed protons is determined using cluster models with dispersion-corrected density functional theory and continuum solvation. Proton acidity constants computed under typical experimental conditions are compared to the acidity constants of other relevant species. The adsorbed protons are predicted to be very stable, suggesting that they will be present on the surface and available to be reduced to surface hydrides that could possibly react with adsorbed pyridine to form adsorbed dihydropyridine, a previously proposed co-catalyst. However, the high stability of such protons also suggests that the surface does not represent a significant proton source; as a consequence, protons required in the proposed mechanism must be provided by a different source such as the acidified aqueous solution in contact with the electrode surface.

SparseMaps--A systematic infrastructure for reduced-scaling electronic structure methods. IV. Linear-scaling second-order explicitly correlated energy with pair natural orbitals.

We present a formulation of the explicitly correlated second-order Møller-Plesset (MP2-F12) energy in which all nontrivial post-mean-field steps are formulated with linear computational complexity in system size. The two key ideas are the use of pair-natural orbitals for compact representation of wave function amplitudes and the use of domain approximation to impose the block sparsity. This development utilizes the concepts for sparse representation of tensors described in the context of the domain based local pair-natural orbital-MP2 (DLPNO-MP2) method by us recently [Pinski et al., J. Chem. Phys. 143, 034108 (2015)]. Novel developments reported here include the use of domains not only for the projected atomic orbitals, but also for the complementary auxiliary basis set (CABS) used to approximate the three- and four-electron integrals of the F12 theory, and a simplification of the standard B intermediate of the F12 theory that avoids computation of four-index two-electron integrals that involve two CABS indices. For quasi-1-dimensional systems (n-alkanes), the ON DLPNO-MP2-F12 method becomes less expensive than the conventional ON(5) MP2-F12 for n between 10 and 15, for double- and triple-zeta basis sets; for the largest alkane, C200H402, in def2-TZVP basis, the observed computational complexity is N(∼1.6), largely due to the cubic cost of computing the mean-field operators. The method reproduces the canonical MP2-F12 energy with high precision: 99.9% of the canonical correlation energy is recovered with the default truncation parameters. Although its cost is significantly higher than that of DLPNO-MP2 method, the cost increase is compensated by the great reduction of the basis set error due to explicit correlation.

Sparse maps--A systematic infrastructure for reduced-scaling electronic structure methods. II. Linear scaling domain based pair natural orbital coupled cluster theory.

Domain based local pair natural orbital coupled cluster theory with single-, double-, and perturbative triple excitations (DLPNO-CCSD(T)) is a highly efficient local correlation method. It is known to be accurate and robust and can be used in a black box fashion in order to obtain coupled cluster quality total energies for large molecules with several hundred atoms. While previous implementations showed near linear scaling up to a few hundred atoms, several nonlinear scaling steps limited the applicability of the method for very large systems. In this work, these limitations are overcome and a linear scaling DLPNO-CCSD(T) method for closed shell systems is reported. The new implementation is based on the concept of sparse maps that was introduced in Part I of this series [P. Pinski, C. Riplinger, E. F. Valeev, and F. Neese, J. Chem. Phys. 143, 034108 (2015)]. Using the sparse map infrastructure, all essential computational steps (integral transformation and storage, initial guess, pair natural orbital construction, amplitude iterations, triples correction) are achieved in a linear scaling fashion. In addition, a number of additional algorithmic improvements are reported that lead to significant speedups of the method. The new, linear-scaling DLPNO-CCSD(T) implementation typically is 7 times faster than the previous implementation and consumes 4 times less disk space for large three-dimensional systems. For linear systems, the performance gains and memory savings are substantially larger. Calculations with more than 20 000 basis functions and 1000 atoms are reported in this work. In all cases, the time required for the coupled cluster step is comparable to or lower than for the preceding Hartree-Fock calculation, even if this is carried out with the efficient resolution-of-the-identity and chain-of-spheres approximations. The new implementation even reduces the error in absolute correlation energies by about a factor of two, compared to the already accurate previous implementation.

Sparse maps­A systematic infrastructure for reduced-scaling electronic structure methods. I. An efficient and simple linear scaling local MP2 method that uses an intermediate basis of pair natural orbitals.

In this work, a systematic infrastructure is described that formalizes concepts implicit in previous work and greatly simplifies computer implementation of reduced-scaling electronic structure methods. The key concept is sparse representation of tensors using chains of sparse maps between two index sets. Sparse map representation can be viewed as a generalization of compressed sparse row, a common representation of a sparse matrix, to tensor data. By combining few elementary operations on sparse maps (inversion, chaining, intersection, etc.), complex algorithms can be developed, illustrated here by a linear-scaling transformation of three-center Coulomb integrals based on our compact code library that implements sparse maps and operations on them. The sparsity of the three-center integrals arises from spatial locality of the basis functions and domain density fitting approximation. A novel feature of our approach is the use of differential overlap integrals computed in linear-scaling fashion for screening products of basis functions. Finally, a robust linear scaling domain based local pair natural orbital second-order Möller-Plesset (DLPNO-MP2) method is described based on the sparse map infrastructure that only depends on a minimal number of cutoff parameters that can be systematically tightened to approach 100% of the canonical MP2 correlation energy. With default truncation thresholds, DLPNO-MP2 recovers more than 99.9% of the canonical resolution of the identity MP2 (RI-MP2) energy while still showing a very early crossover with respect to the computational effort. Based on extensive benchmark calculations, relative energies are reproduced with an error of typically <0.2 kcal/mol. The efficiency of the local MP2 (LMP2) method can be drastically improved by carrying out the LMP2 iterations in a basis of pair natural orbitals. While the present work focuses on local electron correlation, it is of much broader applicability to computation with sparse tensors in quantum chemistry and beyond.

Mechanistic contrasts between manganese and rhenium bipyridine electrocatalysts for the reduction of carbon dioxide.

[Re(bpy)(CO)3](-) is a well-established homogeneous electrocatalyst for the reduction of CO2 to CO. Recently, substitution of the more abundant transition metal Mn for Re yielded a similarly active electrocatalyst, [Mn(bpy)(CO)3](-). Compared to the Re catalyst, this Mn catalyst operates at a lower applied reduction potential but requires the presence of a weak acid in the solution for catalytic activity. In this study, we employ quantum chemistry combined with continuum solvation and microkinetics to examine the mechanism of CO2 reduction by each catalyst. We use cyclic voltammetry experiments to determine the turnover frequencies of the Mn catalyst with phenol as the added weak acid. The computed turnover frequencies for both catalysts agree to within one order of magnitude of the experimental ones. The different operating potentials for these catalysts indicate that different reduction pathways may be favored during catalysis. We model two different pathways for both catalysts and find that, at their respective operating potentials, the Mn catalyst indeed is predicted to take a different reaction route than the Re catalyst. The Mn catalyst can access both catalytic pathways, depending on the applied potential, while the Re catalyst does not show this flexibility. Our microkinetics analysis predicts which intermediates should be observable during catalysis. These intermediates for the two catalyzed reactions have qualitatively different electronic configurations, depending on the applied potential. The observable intermediate at higher applied potentials possesses an unpaired electron and therefore should be EPR-active; however, the observable intermediate at lower applied potentials, accessible only for the Mn catalyst, is diamagnetic and therefore should be EPR-silent. The differences between both catalysts are rationalized on the basis of their electronic structure and different ligand binding affinities.

Hydrogen-bonded cyclic water clusters nucleated on an oxide surface.

We report the observation and molecular-scale scanning probe electronic structure (dI/dV) mapping of hydrogen-bonded cyclic water clusters nucleated on an oxide surface. The measurements are made on a new type of cyclic water cluster that is characterized by simultaneous and cooperative bonding interactions among molecules as well as with both metal and oxygen sites of an oxide surface. Density functional theory + U + D calculations confirm the stability of these clusters and are used to discuss other potential water-oxide bonding scenarios. The calculations show that the spatial distributions of electronic states in the system are similar in character to those of the lowest unoccupied molecular orbitals of hydrogen-bonded water molecules. On the partially oxidized Cu(111) investigated here, experiment and theory together suggest that Cu vacancies in the growing islands of cuprous oxide inhibit water adsorption in the centers of the islands (which have reached thermodynamic equilibrium). A stoichiometric, less stable cuprous oxide likely exists at island edges (the growth front) and selectively binds these water clusters.

New insights into the nature of observable reaction intermediates in cytochrome P450 NO reductase by using a combination of spectroscopy and quantum mechanics/molecular mechanics calculations.

Cytochrome P450 NO reductase is an unusual member of the cytochrome P450 superfamily. It catalyzes the reduction of nitric oxide to nitrous oxide. The reaction intermediates were studied in detail by a combination of experimental and computational methods. They have been characterized experimentally by UV/Vis, EPR, Mössbauer, and MCD spectroscopy. In conjunction with quantum mechanics/molecular mechanics (QM/MM) calculations, we sought to characterize the resting state and the two detectable intermediates in detail and to elucidate the nature of the key intermediate I of the reaction. Six possible candidates were taken into account for the unknown key intermediate in the computational study, differing in protonation state and electronic structure. Two out of the six candidates could be identified as putative intermediates I with the help of the spectroscopic data: singlet diradicals Fe(III)-NHO(·)(-) and Fe(III)-NHOH(.). In a companion publication (C. Riplinger, F. Neese, ChemPhysChem- 2011, 12, 3192) we have used QM/MM models based on these structures and performed a kinetic simulation. The combination of these two studies shows the nature of the key intermediate to be the singlet diradical, Fe(III)-NHOH(·).