Shubham Vyas

Impact of Conjugation and Hyperconjugation on the Radical Stability of Allylic and Benzylic Systems: A Theoretical Study

Resonantly stabilized radicals are some of the most investigated chemical species due to their preferential formation in a wide variety of chemical environments. Density functional theory and post-Hartree-Fock calculations were utilized to elucidate the chemical interactions that contribute to the stability of two ubiquitous, resonantly stabilized radicals, allyl and benzyl radicals. The relative stability of these radical species was quantified through bond dissociation energies and relative rotational energy barriers, with a difference of only 0.1 kcal/mol. To clarify and contextualize the energetic results, natural bond orbitals were used to evaluate the atomic spin density distribution in the given molecules. The benzyl radical was found to be ∼3 kcal/mol less stable than the allyl radical, which was attributed to the inability to efficiently delocalize the spin on a phenyl unit, starkly contrary to general chemistry knowledge. Increasing the degree of π-conjugation and hyperconjugation was shown to benefit allyl radicals to a greater degree than benzyl radicals, again due to more efficient radical delocalization in allyl radicals. This work highlights that more resonance structures do not always lead to a more stabilized radical species, and provides fundamental knowledge about how conjugation and hyperconjugation impact the stabilization of nonbonding electrons in these systems.

A charge-modified general amber force field for phospholipids: improved structural properties in the tensionless ensemble

Accurately predicting the structural properties of phospholipid with a fully atomistic molecular model is critical for the study of pure phospholipid bilayers, mixed bilayer systems and bilayers containing proteins. The general amber force field (GAFF) has traditionally required the presence of a surface tension parameter to correctly model phospholipid bilayer properties such as area per lipid and order parameters. In this work, the GAFF partial charges for 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphate (POPA), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphochiline (POPC) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol (POPG) were re-parameterised utilising high-level ab initio calculations and the restrained electrostatic potential method. Simulations of pure POPA, POPC and POPG bilayers using the charge-modified GAFF and no applied surface tension are compared with available experimental data, the original GAFF model and the recent Lipid14 variant. The results indicate a significant improvement in the accuracy of the lipid model for reproducing experimental observables without the need for a surface tension parameter. The successful application of modifying the lipid charge distributions represents an alternative to the use of a surface tension parameter within GAFF, and highlights the importance of the partial charge calculations when modelling lipid bilayers.

Quantum Mechanical Study of the Reaction Mechanism for 2π-2π Cycloaddition of Fluorinated Methylene Groups

Perfluorocyclobutyl polymers are thermally and chemically stable, may be produced without a catalyst via thermal 2π-2π cycloaddition, and can form block structures, making them suitable for commercialization of specialty polymers. Thermal 2π-2π cycloaddition is a rare reaction that begins in the singlet state and proceeds through a triplet intermediate to form an energetically stable four-membered ring in the singlet state. This reaction involves two changes in spin state and, thus, two spin-crossover transitions. Presented here are density functional theory calculations that evaluate the energetics and reaction mechanisms for the dimerizations of two different polyfluorinated precursors, 1,1,2-trifluoro-2-(trifluoromethoxy)ethane and hexafluoropropylene. The spin-crossover transition states are thoroughly investigated, revealing important kinetics steps and an activation energy for the gas-phase cycloaddition of two hexafluoropropene molecules of 36.9 kcal/mol, which is in good agreement with the experimentally determined value of 34.3 kcal/mol. It is found that the first carbon-carbon bond formation is the rate-limiting step, followed by a rotation about the newly formed bond in the triplet state that results in the formation of the second carbon-carbon bond. Targeting the rotation of the C-C bond, a set of parameters were obtained that best produce high molecular weight polymers using this chemistry.

Solvation Dynamics of HEHEHP Ligand at the Liquid-Liquid Interface

Actinide-lanthanide separation (ALSEP) has been a topic of interest in recent years as it has been shown to selectively extract problematic metals from spent nuclear fuel. However, the process suffers from slow kinetics, prohibiting it from being applied to nuclear facilities. In an effort to improve the process, many fundamental studies have been performed, but the majority have only focused on the thermodynamics of separation. Therefore, to understand the mechanism behind the ALSEP process, molecular dynamics (MD) simulations were utilized to obtain the dynamics and solvation characteristics for an organic extractant, 2-ethylhexylphosphonic acid mono-2-ethylhexyl ester (HEHEHP). Simulations were conducted with both pure and biphasic solvent systems to evaluate the complex solvent interactions within the ALSEP extraction method. The MD simulations revealed solvation and dynamical behaviors that are consistent with the experimentally observed chemical properties of HEHEHP for the pure solvent systems (e.g., hydrophobic/hydrophilic behaviors of the polar head group and alkyl chains and dimer formation between the ligands within an organic solvent). When present in a biphasic solvent system, interfacial behaviors of the ligand revealed that, at low concentrations, the alkyl side chains of HEHEHP were parallel to the interfacial plane. Upon increasing the concentration to 0.75 M, tendency for the parallel orientation decreased and a more perpendicular-like orientation was observed. Analysis of ligand solvation energies in different solvents through the thermodynamic integration method demonstrated favorability toward n-dodecane and biphasic solvents, which is in agreement with the previous experimental findings.

1,2-Fluorine Radical Rearrangements: Isomerization Events in Perfluorinated Radicals

Devising effective degradation technologies for perfluoroalkyl substances (PFASs) is an active area of research, where the molecular mechanisms involving both oxidative and reductive pathways are still elusive. One commonly neglected pathway in PFAS degradation is fluorine atom migration in perfluoroalkyl radicals, which was largely assumed to be implausible because of the high C-F bond strength. Using density functional theory calculations, it was demonstrated that 1,2-F atom migrations are thermodynamically favored when the fluorine atom migrated from a less branched carbon center to a more branched carbon center. Activation barriers for these rearrangements were within 19-29 kcal/mol, which are possible to easily overcome at elevated temperatures or in photochemically activated species in the gas or aqueous phase. It was also found that the activation barriers for the 1,2-F atom migration are lowered as much as by 10 kcal/mol when common oxidative degradation products such as HF assisted the rearrangements or if the resulting radical center was stabilized by vicinal π-bonds. Natural bond orbital analyses showed that fluorine moves as a radical in a noncharge-separated state. These findings add an important reaction to the existing knowledge of mechanisms for PFAS degradation and highlights the fact that 1,2-F atom shifts may be a small channel for isomerization of these compounds, but upon availability of mineralization products, this isomerization process could become more prominent.