Efficiency and Mechanism of Catalytic Siloxane Exchange in Vitrimer Polymers: Modeling and Density Functional Theory Investigations.

Of late, siloxane-containing vitrimers have gained significant interest due to their fast dynamic characteristics over a reasonable temperature range (180-220 °C), making them well-suited for diverse applications. The exchange reaction pathway in the siloxane vitrimers is accountable for the covalent adaptive network, with the reaction's effectiveness being regulated by either organic or organometallic catalysts. However, directly studying the exchange reaction pathway in the bulk phase using experimental approaches is challenging because of the intricate and interconnected structure of these vitrimers. Here, we perform comprehensive density functional theory (DFT) and experimental investigations to discover the detailed catalytic efficacy of siloxane exchange and provide direction for the reaction process using a 1,5,7-triazabicyclo[4.4.0]dec-5-ene (TBD) catalyst. The calculated transition barrier energy and catalytic efficiency of hexamethyldisiloxane and dihydroxy-dimethylsilane exchange derived from the nudged elastic band with transition-state calculations strongly agree with the experimental findings. In addition, Fukui indices, along with partial charges, are employed to evaluate the nucleophilic and electrophilic behaviors of silanol and siloxane molecules. Our analysis revealed that by utilizing the Fukui indices of both the acid and the base, we can make an approximate estimation of the respective kinetics of the SN2 process in the siloxane exchange reaction mechanism. These findings establish a foundation for comprehending a crucial aspect of the exchange mechanism in siloxane vitrimer systems and could aid in the development of novel catalysts.

Ab initio molecular dynamics modeling of single polyethylene chains: Scission kinetics and influence of radical under mechanical strain.

Ab initio molecular dynamics was used to estimate the response to constant imposed strain on a short polyethylene (PE) chain and a radical chain with a removed hydrogen atom. Two independent types of simulations were run. In the first case, the chains were strained by expanding a periodic cell, restraining the length but allowing the internal degrees of freedom to reach equilibrium. From these simulations, the average force on the chain was computed, and the resulting force was integrated to determine the Helmholtz free energy for chain stretching. In the second set of simulations, chains were constrained to various lengths, while a bond was restrained at various bond lengths using umbrella sampling. This provided free energy of bond scission for various chain strains. The sum of the two free energy functions results in an approximation of the free energy of chain scission under various strains and gives a realistic and new picture of the effect of chain strain on bond breaking. Unimolecular scission rates for each chain type were examined as a function of chain strain. The scission rate for the radical chain is several orders of magnitude larger than that of the pristine chain at smaller strains and at equilibrium. This highlights the importance of radical formation in PE rupture and is consistent with experiments. Constant strain results were used to derive a constant-force model for the radical chain that demonstrates a roll over in rate similar to the "catch-bond" behavior observed in protein membrane detachment experiments.

The generality of the GUGA MRCI approach in COLUMBUS for treating complex quantum chemistry.

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

Analytic non-adiabatic derivative coupling terms for spin-orbit MRCI wavefunctions. I. Formalism.

Analytic gradients of electronic eigenvalues require one calculation per nuclear geometry, compared to at least 3n + 1 calculations for finite difference methods, where n is the number of nuclei. Analytic nonadiabatic derivative coupling terms (DCTs), which are calculated in a similar fashion, are used to remove nondiagonal contributions to the kinetic energy operator, leading to more accurate nuclear dynamics calculations than those that employ the Born-Oppenheimer approximation, i.e., that assume off-diagonal contributions are zero. The current methods and underpinnings for calculating both of these quantities, gradients and DCTs, for the State-Averaged MultiReference Configuration Interaction with Singles and Doubles (MRCI-SD) wavefunctions in COLUMBUS are reviewed. Before this work, these methods were not available for wavefunctions of a relativistic MRCI-SD Hamiltonian. Calculation of these terms is critical in successfully modeling the dynamics of systems that depend on transitions between potential energy surfaces split by the spin-orbit operator, such as diode-pumped alkali lasers. A formalism for calculating the transition density matrices and analytic derivative coupling terms for such systems is presented.

Analytic non-adiabatic derivative coupling terms for spin-orbit MRCI wavefunctions. II. Derivative coupling terms and coupling angle for KHeA2Π1/2⇔KHeB2Σ1/2.

A method for calculating the analytic nonadiabatic derivative coupling terms (DCTs) for spin-orbit multi-reference configuration interaction wavefunctions is reviewed. The results of a sample calculation using a Stuttgart basis for KHe are presented. Additionally, the DCTs are compared with a simple calculation based on the Nikitin's 3 × 3 description of the coupling between the Σ and Π surfaces, as well as a method based on Werner's analysis of configuration interaction coefficients. The nonadiabatic coupling angle calculated by integrating the radial analytic DCTs using these different techniques matches extremely well. The resultant nonadiabatic energy surfaces for KHe are presented.

Bond breaking in epoxy systems: A combined QM/MM approach.

A novel method to combine quantum mechanics (QM) and molecular mechanics has been developed to accurately and efficiently account for covalent bond breaking in polymer systems under high strain without the use of predetermined break locations. Use of this method will provide a better fundamental understanding of the mechano-chemical origins of fracture in thermosets. Since classical force fields cannot accurately account for bond breaking, and QM is too demanding to simulate large systems, a hybrid approach is required. In the method presented here, strain is applied to the system using a classical force field, and all bond lengths are monitored. When a bond is stretched past a threshold value, a zone surrounding the bond is used in a QM energy minimization to determine which, if any, bonds break. The QM results are then used to reconstitute the system to continue the classical simulation at progressively larger strain until another QM calculation is triggered. In this way, a QM calculation is only computed when and where needed, allowing for efficient simulations. A robust QM method for energy minimization has been determined, as well as appropriate values for the QM zone size and the threshold bond length. Compute times do not differ dramatically from classical molecular mechanical simulations.

Computation of the binding free energy of peptides to graphene in explicit water.

The characteristic properties of graphene make it useful in an assortment of applications. One particular application--the use of graphene in biosensors--requires a thorough understanding of graphene-peptide interactions. In this study, the binding of glycine (G) capped amino acid residues (termed GXG tripeptides) to trilayer graphene surfaces in aqueous solution was examined and compared to results previously obtained for peptide binding to single-layer free-standing graphene [A. N. Camden, S. A. Barr, and R. J. Berry, J. Phys. Chem. B 117, 10691-10697 (2013)]. In order to understand the interactions between the peptides and the surface, binding enthalpy and free energy values were calculated for each GXG system, where X cycled through the typical 20 amino acids. When the GXG tripeptides were bound to the surface, distinct conformations were observed, each with a different binding enthalpy. Analysis of the binding energy showed the binding of peptides to trilayer graphene was dominated by van der Waals interactions, unlike the free-standing graphene systems, where the binding was predominantly electrostatic in nature. These results demonstrate the utility of computational materials science in the mechanistic explanation of surface-biomolecule interactions which could be applied to a wide range of systems.

Accurate prediction by AlphaFold2 for ligand binding in a reductive dehalogenase and implications for PFAS (per- and polyfluoroalkyl substance) biodegradation.

Despite the success of AlphaFold2 (AF2), it is unclear how AF2 models accommodate for ligand binding. Here, we start with a protein sequence from Acidimicrobiaceae TMED77 (T7RdhA) with potential for catalyzing the degradation of per- and polyfluoroalkyl substances (PFASs). AF2 models and experiments identified T7RdhA as a corrinoid iron-sulfur protein (CoFeSP) which uses a norpseudo-cobalamin (BVQ) cofactor and two Fe4S4 iron-sulfur clusters for catalysis. Docking and molecular dynamics simulations suggest that T7RdhA uses perfluorooctanoic acetate (PFOA) as a substrate, supporting the reported defluorination activity of its homolog, A6RdhA. We showed that AF2 provides processual (dynamic) predictions for the binding pockets of ligands (cofactors and/or substrates). Because the pLDDT scores provided by AF2 reflect the protein native states in complex with ligands as the evolutionary constraints, the Evoformer network of AF2 predicts protein structures and residue flexibility in complex with the ligands, i.e., in their native states. Therefore, an apo-protein predicted by AF2 is actually a holo-protein awaiting ligands.

M + Ng potential energy curves including spin-orbit coupling for M = K, Rb, Cs and Ng = He, Ne, Ar.

The X(2)Σ(1/2)(+), A(2)Π(1∕2), A(2)Π(3∕2), and B(2)Σ(1/2)(+) potential energy curves and associated dipole matrix elements are computed for M + Ng at the spin-orbit multi-reference configuration interaction level, where M = K, Rb, Cs and Ng = He, Ne, Ar. Dissociation energies and equilibrium positions for all minima are identified and corresponding vibrational energy levels are computed. Difference potentials are used together with the quasistatic approximation to estimate the position of satellite peaks of collisionally broadened D2 lines. The comparison of potential energy curves for different alkali atom and noble gas atom combinations is facilitated by using the same level of theory for all nine M + Ng pairs.

AlphaFold2 models indicate that protein sequence determines both structure and dynamics.

AlphaFold 2 (AF2) has placed Molecular Biology in a new era where we can visualize, analyze and interpret the structures and functions of all proteins solely from their primary sequences. We performed AF2 structure predictions for various protein systems, including globular proteins, a multi-domain protein, an intrinsically disordered protein (IDP), a randomized protein, two larger proteins (> 1000 AA), a heterodimer and a homodimer protein complex. Our results show that along with the three dimensional (3D) structures, AF2 also decodes protein sequences into residue flexibilities via both the predicted local distance difference test (pLDDT) scores of the models, and the predicted aligned error (PAE) maps. We show that PAE maps from AF2 are correlated with the distance variation (DV) matrices from molecular dynamics (MD) simulations, which reveals that the PAE maps can predict the dynamical nature of protein residues. Here, we introduce the AF2-scores, which are simply derived from pLDDT scores and are in the range of [0, 1]. We found that for most protein models, including large proteins and protein complexes, the AF2-scores are highly correlated with the root mean square fluctuations (RMSF) calculated from MD simulations. However, for an IDP and a randomized protein, the AF2-scores do not correlate with the RMSF from MD, especially for the IDP. Our results indicate that the protein structures predicted by AF2 also convey information of the residue flexibility, i.e., protein dynamics.

Transesterification in Vitrimer Polymers Using Bifunctional Catalysts: Modeled with Solution-Phase Experimental Rates and Theoretical Analysis of Efficiency and Mechanisms.

Recently, thermoset vitrimer polymers have shown significant promise for structural applications because of their ability to be reshaped and remolded due to their covalent adaptive network (CAN). In these vitrimers, the transesterification reaction is responsible for the CAN, where the efficiency of the reaction is controlled either by organic or by organometallic catalysts. Understanding the mechanism of the transesterification reaction in the bulk phase using direct experimental techniques is extremely difficult due to the highly cross-linked complex structure of thermosetting vitrimers. Therefore, we use solution-phase experiments to investigate the catalytic efficiency and to guide density functional theory (DFT) simulations of the transesterification reaction mechanism with catalysts triazabicyclodecene (TBD), zinc acetate (Zn(OAc)2), 1-methylimidazole (1-MI), and dibutyltin oxide (DBTO). The estimated catalytic efficiency from the detailed DFT reaction path calculations follows the order TBD ≳ DBTO ≳ Zn(OAc)2 > 1-MI, which agrees with the experimental results. In addition to reaction path modeling, the mechanism and the relative rates of the transesterification reaction are analyzed with the assistance of Fukui indices as a measure of electrophilicity and nucleophilicity of atomic sites and with partial charges. It was found that the sum of the nucleophilicity index of the base and the electrophilicity index of the acid of the bifunctional catalysts correlates with the SN2 transition state and tetrahedral intermediate energies, which are related to the barrier of the rate-limiting step. This correlation provides a hypothesis for computational prescreening of potentially better catalysts that have an index in a range of values. These results provide a basis for understanding an important part of the mechanism of transesterification in vitrimer systems and may assist with designing new catalysts.

Observation of two discrete conductivity states in quinone-oligo(phenylene vinylene).

The single-molecule conductivity of quinone-oligo(phenylene vinylene) (Q-OPV) attached to a gold substrate was studied using electrochemical scanning tunnelling microscopy. The results show that the molecule has two discrete conductivity states: a low-conductivity state, when it is oxidized, and a high-conductivity state, when reduced. The electron transport through the molecule in both states occurs via coherent tunnelling. The molecular conductivity in either oxidation state is independent from the electrochemical gate potential; however, the gate potential can be used to switch the oxidation state of the molecule. Numerical calculations suggest that the highest occupied molecular orbital (HOMO) of Q-OPV controls tunnelling through the molecule and that the independence of conductivity from the electrochemical gate in either oxidation state originates from strong penetration of HOMO into the substrate. In addition, the greater delocalization of HOMO in the reduced state than in the oxidized state explains the greater conductivity of Q-OPV in the former than in the latter.

Structural reorganizations control intermolecular conductance and charge trapping in paraquat-tetraphenylborate inverse photochemical cell.

Miniaturization of electronic devices to the level of single molecules requires detailed understanding of the mechanisms of their operation. One of the questions here is the identification of the role of structural alterations in charge separation and stabilization in photoactive complexes. To address this question, we calculate optimized molecular and electronic structures, and optical and vibrational spectra of l,l'-dimethyl 4,4'-bipyridinium-bis tetraphenylborate PQ(BPh(4))(2) complex ab initio using density functional theory approach and compare them with the experimentally observed UV-Vis and Raman spectra of the molecules in solid-state films. The results indicate that the association of PQ and BPh(4) leads to the formation of an internally ionized structure that is accompanied by the structural reorganization of both PQ (the twisting of pyridinium rings) and BPh(4) (phenyl rings rotation) moieties. The quanta of light do not seem to be directly involved in the formation of this ionized structure, but provide energy for fast recombination of the separated charges between BPh(4)(-) and PQ(2+). The high efficiency of the dark charge separation and the stabilization of separated charges in the complex permit the using of PQ(BPh(4))(2) in various charge-transfer devices like molecular probes, photovoltaic devices or chemical memory units.

Advances in methods and algorithms in a modern quantum chemistry program package.

Advances in theory and algorithms for electronic structure calculations must be incorporated into program packages to enable them to become routinely used by the broader chemical community. This work reviews advances made over the past five years or so that constitute the major improvements contained in a new release of the Q-Chem quantum chemistry package, together with illustrative timings and applications. Specific developments discussed include fast methods for density functional theory calculations, linear scaling evaluation of energies, NMR chemical shifts and electric properties, fast auxiliary basis function methods for correlated energies and gradients, equation-of-motion coupled cluster methods for ground and excited states, geminal wavefunctions, embedding methods and techniques for exploring potential energy surfaces.