Radical perfluoroalkylation – easy access to 2-perfluoroalkylindol-3-imines via electron catalysis

Arylisonitriles react with alkyl and perfluoroalkyl radicals to form 2-alkylated indole-3-iminesviatwo sequential additions to the isonitrile moiety followed by homolytic aromatic substitution.

Photochemically induced dynamic nuclear polarization NMR on photosystem II: donor cofactor observed in entire plant

The solid-state photo-CIDNP (photochemically induced dynamic nuclear polarization) effect allows for increase of signal and sensitivity in magic-angle spinning (MAS) NMR experiments. The effect occurs in photosynthetic reaction centers (RC) proteins upon illumination and induction of cyclic electron transfer. Here we show that the strength of the effect allows for observation of the cofactors forming the spin-correlated radical pair (SCRP) in isolated proteins, in natural photosynthetic membranes as well as in entire plants. To this end, we measured entire selectively 13C isotope enriched duckweed plants (Spirodela oligorrhiza) directly in the MAS rotor. Comparison of 13C photo-CIDNP MAS NMR spectra of photosystem II (PS2) obtained from different levels of RC isolation, from entire plant to isolated RC complex, demonstrates the intactness of the photochemical machinery upon isolation. The SCRP in PS2 is structurally and functionally very similar in duckweed and spinach (Spinacia oleracea). The analysis of the photo-CIDNP MAS NMR spectra reveals a monomeric Chl a donor. There is an experimental evidence for matrix involvement, most likely due to the axial donor histidine, in the formation of the SCRP. Data do not suggest a chemical modification of C-131 carbonyl position of the donor cofactor.

Ab Initio Characterization of the Electrostatic Complexes Formed by H2 Molecule and Cr(+), Mn(+), Cu(+), and Zn(+) Cations

Equilibrium structures, dissociation energies, and rovibrational energy levels of the electrostatic complexes formed by molecular hydrogen and first-row S-state transition metal cations Cr(+), Mn(+), Cu(+), and Zn(+) are investigated ab initio. Extensive testing of the CCSD(T)-based approaches for equilibrium structures provides an optimal scheme for the potential energy surface calculations. These surfaces are calculated in two dimensions by keeping the H-H internuclear distance fixed at its equilibrium value in the complex. Subsequent variational calculations of the rovibrational energy levels permits direct comparison with data obtained from equilibrium thermochemical and spectroscopic measurements. Overall accuracy within 2-3% is achieved. Theoretical results are used to examine trends in hydrogen activation, vibrational anharmonicity, and rotational structure along the sequence of four electrostatic complexes covering the range from a relatively floppy van der Waals system (Mn(+)···H2) to an almost a rigid molecular ion (Cu(+)···H2).

Theoretical Assessment of Hinge-Type Models for Electron Donors in Reaction Centers of Photosystems I and II as well as of Purple Bacteria

Hinge-type molecular models for electron donors in reaction centers of Photosystems I and II and purple bacteria were investigated using a two-state computational approach based on frozen-density embedding (FDE). This methodology, dubbed FDE-diab, is known to avoid consequences of the self-interaction error as far as intermolecular phenomena are concerned, which allows a prediction of qualitatively correct spin densities for large biomolecular systems. The calculated spin density distributions are in a good agreement with available experimental results and demonstrated a very high sensitivity to changes in the relative orientation of cofactors and amino acid protonation states. This allows a validation of the previously proposed hinge-type models providing hints on possible protonation states of axial histidine molecules.

Adaptive density-guided approach to double incremental potential energy surface construction

We present a combination of the recently developed double incremental expansion of potential energy surfaces with the well-established adaptive density-guided approach to grid construction. This unique methodology is based on the use of an incremental expansion for potential energy surfaces, known as n-mode expansion; an incremental many-body representation of the electronic energy; and an efficient vibrational density-guided approach to automated determination of grid dimensions and granularity. The reliability of the method is validated calculating potential energy surfaces and obtaining fundamental excitation energies for three moderate-size chain-like molecular systems. The use of our methodology leads to considerable computational savings for potential energy surface construction compared to standard approaches while maintaining a high level of accuracy in the resulting potential energy surfaces. Additional investigations indicate that our method can be applied to covalently bound and strongly interacting molecular systems, even though these cases are known to be very unfavorable for fragmentation schemes. We therefore conclude that the presented methodology is a robust and flexible approach to potential energy surface construction, which introduces considerable computational savings without compromising the accuracy of vibrational spectra calculations.

Multi-state formulation of the frozen-density embedding quasi-diabatization approach

We present a multi-state implementation of the recently developed frozen-density embedding diabatization (FDE-diab) methodology [D. G. Artiukhin and J. Neugebauer, J. Chem. Phys. 148, 214104 (2018)] in the Serenity program. The new framework extends the original approach such that any number of charge-localized quasi-diabatic states can be coupled, giving an access to calculations of ground and excited state spin-density distributions as well as to excitation energies. We show that it is possible to obtain results similar to those from correlated wave function approaches such as the complete active space self-consistent field method at much lower computational effort. Additionally, we present a series of approximate computational schemes, which further decrease the overall computational cost and systematically converge to the full FDE-diab solution. The proposed methodology enables computational studies on spin-density distributions and related properties for large molecular systems of biochemical interest.

Developing orbital-dependent corrections for the non-additive kinetic energy in subsystem density functional theory

We present a novel route to constructing cost-efficient semi-empirical approximations for the non-additive kinetic energy in subsystem density functional theory. The developed methodology is based on the use of Slater determinants composed of non-orthogonal Kohn-Sham-like orbitals for the evaluation of kinetic energy expectation values and the expansion of the inverse molecular-orbital overlap matrix into a Neumann series. By applying these techniques, we derived and implemented a series of orbital-dependent approximations for the non-additive kinetic energy, which are employed self-consistently. Our proof-of-principle computations demonstrated quantitatively correct results for potential energy curves and electron densities and hinted on the applicability of the introduced empirical parameters to different types of molecular systems and intermolecular interactions. Therefore, we conclude that the presented study is an important step toward constructing accurate and efficient orbital-dependent approximations for the non-additive kinetic energy applicable to large molecular systems.

Gaussian process regression adaptive density-guided approach: Toward calculations of potential energy surfaces for larger molecules

We present a new program implementation of the Gaussian process regression adaptive density-guided approach [Schmitz et al., J. Chem. Phys. 153, 064105 (2020)] for automatic and cost-efficient potential energy surface construction in the MidasCpp program. A number of technical and methodological improvements made allowed us to extend this approach toward calculations of larger molecular systems than those previously accessible and maintain the very high accuracy of constructed potential energy surfaces. On the methodological side, improvements were made by using a Δ-learning approach, predicting the difference against a fully harmonic potential, and employing a computationally more efficient hyperparameter optimization procedure. We demonstrate the performance of this method on a test set of molecules of growing size and show that up to 80% of single point calculations could be avoided, introducing a root mean square deviation in fundamental excitations of about 3 cm-1. A much higher accuracy with errors below 1 cm-1 could be achieved with tighter convergence thresholds still reducing the number of single point computations by up to 68%. We further support our findings with a detailed analysis of wall times measured while employing different electronic structure methods. Our results demonstrate that GPR-ADGA is an effective tool, which could be applied for cost-efficient calculations of potential energy surfaces suitable for highly accurate vibrational spectra simulations.