No Arabic abstract
The electronic structure of the nitrogenase metal cofactors is central to nitrogen fixation. However, the P-cluster and iron molybdenum cofactor, each containing eight irons, have resisted detailed characterization of their electronic properties. Through exhaustive many-electron wavefunction simulations enabled by new theoretical methods, we report on the low-energy electronic states of the P-cluster in three oxidation states. The energy scales of orbital and spin excitations overlap, yielding a dense spectrum with features we trace to the underlying atomic states and recouplings. The clusters exist in superpositions of spin configurations with non-classical spin correlations, complicating interpretation of magnetic spectroscopies, while the charges are mostly localized from reorganization of the cluster and its surroundings. Upon oxidation, the opening of the P-cluster significantly increases the density of states, which is intriguing given its proposed role in electron transfer. These results demonstrate that many-electron simulations stand to provide new insights into the electronic structure of the nitrogenase cofactors.
We report that a recent active space model of the nitrogenase FeMo cofactor, proposed in the context of quantum simulations, is not representative of the electronic structure of the FeMo cofactor ground-state. Although quantum resource estimates, outside of the cost of adiabatic state preparation, will not be much affected, conclusions should not be drawn from the complexity of classical simulations of the electronic structure of this system in this active space. We provide a different model active space for the FeMo cofactor that contains the basic open-shell qualitative character, which may be useful as a benchmark system for making classical and quantum resource estimates.
We assess the utility of Hartree-Fock (HF) trial wavefunctions in performing phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) on the uniform electron gas (UEG) model. The combination of ph-AFQMC with spin-restricted HF (RHF+ph-AFQMC), was found to be highly accurate and efficient for systems containing up to 114 electrons in 2109 orbitals, particularly for $r_s$ $le$ 2.0. Compared to spin-restricted coupled-cluster (RCC) methods, we found that RHF+ph-AFQMC performs better than CC with singles, doubles, and triples (RCCSDT) and similarly to or slightly worse than CC with singles, doubles, triples, and quadruples (RCCSDTQ) for $r_s$ $le$ 3.0 in the 14-electron UEG model. With the 54-electron, we found RHF+ph-AFQMC to be nearly exact for $r_s$ $le$ 2.0 and pointed out potential biases in existing benchmarks. Encouraged by these, we performed RHF+ph-AFQMC on the 114-electron UEG model for $r_s$ $le$ 2.0 and provided new benchmark data for future method development. We found that the UEG models with $r_s$ = 5.0 remain to be challenging for RHF+ph-AFQMC. Employing non-orthogonal configuration expansions or unrestricted HF states as trial wavefunctions was also found to be ineffective in the case of the 14-electron UEG model with $r_s$ = 5.0. We emphasize the need for a better trial wavefunction for ph-AFQMC in simulating strongly correlated systems. With the 54-electron and 114-electron UEG models, we stress the potential utility of RHF+ph-AFQMC for simulating dense solids.
A full coupled-cluster expansion suitable for sparse algebraic operations is developed by expanding the commutators of the Baker-Campbell-Hausdorff series explicitly for cluster operators in binary representations. A full coupled-cluster reduction that is capable of providing very accurate solutions of the many-body Schrodinger equation is then initiated employing screenings to the projection manifold and commutator operations. The projection manifold is iteratively updated through the single commutators $leftlangle kappa right| [hat H,hat T]left| 0 rightrangle$ comprised of the primary clusters $hat T_{lambda}$ with substantial contribution to the connectivity. The operation of the commutators is further reduced by introducing a correction, taking into account the so-called exclusion principle violating terms, that provides fast and near-variational convergence in many cases.
We carry out a first-principles atomistic study of the electronic mechanisms of ligand binding and discrimination in the myoglobin protein. Electronic correlation effects are taken into account using one of the most advanced methods currently available, namely a linear-scaling density functional theory (DFT) approach wherein the treatment of localized iron 3d electrons is further refined using dynamical mean-field theory (DMFT). This combination of methods explicitly accounts for dynamical and multi-reference quantum physics, such as valence and spin fluctuations, of the 3d electrons, whilst treating a significant proportion of the protein (more than 1000 atoms) with density functional theory. The computed electronic structure of the myoglobin complexes and the nature of the Fe-O2 bonding are validated against experimental spectroscopic observables. We elucidate and solve a long standing problem related to the quantum-mechanical description of the respiration process, namely that DFT calculations predict a strong imbalance between O2 and CO binding, favoring the latter to an unphysically large extent. We show that the explicit inclusion of many body-effects induced by the Hunds coupling mechanism results in the correct prediction of similar binding energies for oxy- and carbonmonoxymyoglobin.
Hyperconjugation is a basic conception of chemistry. Its straightforward effect is exhibited by the spatial delocalization characteristics of the electron density distributions or wavefunctions. Such effects on the electron wavefunctions of the highest-occupied molecular orbitals (HOMO) of two ethanol conformers are demonstrated with electron momentum spectroscopy together with natural bond orbital analyses, exhibiting the distinctly different symmetries of the HOMO wavefunctions in momentum space.