Stability Trends in disubstituted Cobaltocenium Based on the Analysis of the Machine Learning Models.

Cobaltocenium derivatives have shown great potential as components of anion exchange membranes in fuel cells because they exhibit excellent thermal and alkaline stability under operating conditions while allowing for high anion mobility. The properties of the cobaltocenium-anion complexes can be chemically tuned through the substituent groups on the cyclopentadienyl (Cp) rings of the cation CoCp2+. However, the synthesis and characterization of the full range of possible derivatives are very challenging and time-consuming, and while the computational tools can greatly expedite this process, full screening of the electronic structure at a high level of theory is still computationally intensive. Therefore, in this work, we consider the machine learning (ML) modeling as a tool of predicting stability of disubstituted [CoCp2]OH complexes measured by their bond-dissociation energy (BDE). The relevant process here is the dissociation of the cobaltocenium-hydroxide complex into fragments [CoCpY']OH and CpY, where Y and Y' each represent one out of 42 substituent groups of experimental interest. In agreement with the previous ML study of 120 mono- and selected disubstituted species [Wetthasinghe et al. J. Chem. Phys. A (2022) 126], our analysis of the data set expanded to all possible disubstituted cobaltoceniums, points to the highest occupied and lowest unoccupied molecular orbitals, along with the Hirshfeld charge on the singly substituted benzene, to be the key features predicting the BDE of the unseen complexes. Based on the examination of the outliers, the acidity of substituents ((CO)NH2 in our case) is found to be of special significance for the cobaltocenium stability and for the model development. Moreover, we demonstrate that upon the data set refinement, the conventional ML models are capable of predicting the BDE close to 1 kcal/mol based on the properties of just the fragments, thereby greatly reducing the total number of species and of the computational time of each calculation. Such fragment-based "combinatorial" approach to the BDE modeling is noteworthy, since the geometry optimization of complexes in solution is conceptually challenging and computationally demanding, even when leveraging high-performance computing resources.

Cooperative and Orthogonal Switching in the Solid State Enabled by Metal-Organic Framework Confinement Leading to a Thermo-Photochromic Platform.

Cooperative behavior and orthogonal responses of two classes of coordinatively integrated photochromic molecules towards distinct external stimuli were demonstrated on the first example of a photo-thermo-responsive hierarchical platform. Synergetic and orthogonal responses to temperature and excitation wavelength are achieved by confining the stimuli-responsive moieties within a metal-organic framework (MOF), leading to the preparation of a novel photo-thermo-responsive spiropyran-diarylethene based material. Synergistic behavior of two photoswitches enables the study of stimuli-responsive resonance energy transfer as well as control of the photoinduced charge transfer processes, milestones required to advance optoelectronics development. Spectroscopic studies in combination with theoretical modeling revealed a nonlinear effect on the material electronic structure arising from the coordinative integration of photoresponsive molecules with distinct photoisomerization mechanisms. Thus, the reported work covers multivariable facets of not only fundamental aspects of photoswitch cooperativity, but also provides a pathway to modulate photophysics and electronics of multidimensional functional materials exhibiting thermo-photochromism.

Confinement-Driven Photophysics in Hydrazone-Based Hierarchical Materials.

Confinement-imposed photophysics was probed for novel stimuli-responsive hydrazone-based compounds demonstrating a conceptual difference in their behavior within 2D versus 3D porous matrices for the first time. The challenges associated with photoswitch isomerization arising from host interactions with photochromic compounds in 2D scaffolds could be overcome in 3D materials. Solution-like photoisomerization rate constants were realized for sterically demanding hydrazone derivatives in the solid state through their coordinative immobilization in 3D scaffolds. According to steady-state and time-resolved photophysical measurements and theoretical modeling, this approach provides access to hydrazone-based materials with fast photoisomerization kinetics in the solid state. Fast isomerization of integrated hydrazone derivatives allows for probing and tailoring resonance energy transfer (ET) processes as a function of excitation wavelength, providing a novel pathway for ET modulation.

Structure-property investigations in urea tethered iodinated triphenylamines.

Herein, we report structural, computational, and conductivity studies on urea-directed self-assembled iodinated triphenylamine (TPA) derivatives. Despite numerous reports of conductive TPAs, the challenges of correlating their solid-state assembly with charge transport properties hinder the efficient design of new materials. In this work, we compare the assembled structures of a methylene urea bridged dimer of di-iodo TPA (1) and the corresponding methylene urea di-iodo TPA monomer (2) with a di-iodo mono aldehyde (3) control. These modifications lead to needle shaped crystals for 1 and 2 that are organized by urea hydrogen bonding, π⋯π stacking, I⋯I, and I⋯π interactions as determined by SC-XRD, Hirshfeld surface analysis, and X-ray photoelectron spectroscopy (XPS). The long needle shaped crystals were robust enough to measure the conductivity by two contact probe methods with 2 exhibiting higher conductivity values (∼6 × 10-7 S cm-1) compared to 1 (1.6 × 10-8 S cm-1). Upon UV-irradiation, 1 formed low quantities of persistent radicals with the simple methylurea 2 displaying less radical formation. The electronic properties of 1 were further investigated using valence band XPS, which revealed a significant shift in the valence band upon UV irradiation (0.5-1.9 eV), indicating the potential of these materials as dopant free p-type hole transporters. The electronic structure calculations suggest that the close packing of TPA promotes their electronic coupling and allows effective charge carrier transport. Our results show that ionic additives significantly improve the conductivity up to ∼2.0 × 10-6 S cm-1 in thin films, enabling their implementation in functional devices such as perovskite or solid-state dye sensitized solar cells.

Modeling the Ligand Effect on the Structure of CYP 450 Within the Density Functional Theory.

An improved understanding of the P450 structure is relevant to the development of biomimetic catalysts and inhibitors for controlled CH-bond activation, an outstanding challenge of synthetic chemistry. Motivated by the experimental findings of an unusually short Fe-S bond of 2.18 Šfor the wild-type (WT) OleT P450 decarboxylase relative to a cysteine pocket mutant form (A369P), a computational model that captures the effect of the thiolate axial ligand on the iron-sulfur distance is presented. With the computational efficiency and streamlined analysis in mind, this model combines a cluster representation of the enzyme─40-110 atoms, depending on the heme and ligand truncation level─with a density functional theory (DFT) description of the electronic structure (ES) and is calibrated against the experimental data. The optimized Fe-S distances show a difference of 0.25 Šbetween the low and high spin states, in agreement with the crystallographic structures of the OleT WT and mutant forms. We speculate that this difference is attributable to the packing of the ligand; the mutant is bulkier due to an alanine-to-proline replacement, meaning that it is excluded from the energetically favored low-spin minimum because of steric constraints. The presence of pure spin-state pairs and the intersection of the low/high spin states for the enzyme model is indicative of the limitations of single-reference ES methods in such systems and emphasizes the significance of using the proper state when modeling the hydrogen atom transfer (HAT) reaction catalyzed by OleT. At the same time, the correct characterization of both the short and long Fe-S bonds within a small DFT-based model of 42 atoms paves the way for quantum dynamics modeling of the HAT step, which initiates the OleT decarboxylation reaction.

Correlation between the Stability of Substituted Cobaltocenium and Molecular Descriptors.

Metallocenium cations, used as a component in an anion exchange membrane of a fuel cell, demonstrate excellent thermal and alkaline stability, which can be improved by the chemical modification of the cyclopentadienyl rings with substituent groups. In this work, the relation between the bond dissociation energy (BDE) of the cobaltocenium (CoCp2+) derivatives, used as a measure of the cation stability, and chemistry-informed descriptors obtained from the electronic structural calculations is established. The analysis of 12 molecular descriptors for 118 derivatives reveals a nonlinear dependence of the BDE on the electron donating-withdrawing character of the substituent groups coupled to the energy of the frontier molecular orbitals. A chemistry-informed feed-forward neural network trained using k-fold cross-validation over the modest data set is able to predict the BDE from the molecular descriptors with the mean absolute error of about 1 kcal/mol. The theoretical analysis suggests some promising modifications of cobaltocenium for experimental research. The results demonstrate that even for modest data sets the incorporation of the chemistry knowledge into the neural network architecture, e.g., through mindful selection and screening of the descriptors and their interactions, paves the way to gain new insight into molecular properties.

A Metal-Organic Framework (MOF)-Based Multifunctional Cargo Vehicle for Reactive-Gas Delivery and Catalysis.

The efficient delivery of reactive and toxic gaseous reagents to organic reactions was studied using metal-organic frameworks (MOFs). The simultaneous cargo vehicle and catalytic capabilities of several MOFs were probed for the first time using the examples of aromatization, aminocarbonylation, and carbonylative Suzuki-Miyaura coupling reactions. These reactions highlight that MOFs can serve a dual role as a gas cargo vehicle and a catalyst, leading to product formation with yields similar to reactions employing pure gases. Furthermore, the MOFs can be recycled without sacrificing product yield, while simultaneously maintaining crystallinity. The reported findings were supported crystallographically and spectroscopically (e.g., diffuse reflectance infrared Fourier transform spectroscopy), foreshadowing a pathway for the development of multifunctional MOF-based reagent-catalyst cargo vessels for reactive gas reagents as an attractive alternative to the use of toxic pure gases or gas generators.

Assembled triphenylamine bis-urea macrocycles: exploring photodriven electron transfer from host to guests.

Absorption of electronic acceptors in the accessible channels of an assembled triphenylamine (TPA) bis-urea macrocycle 1 enabled the study of electron transfer from the walls of the TPA framework to the encapsulated guests. The TPA host is isoskeletal in all host-guest structures analyzed with guests 2,1,3-benzothiadiazole, 2,5-dichlorobenzoquinone and I2 loading in single-crystal-to-single-crystal transformations. Analysis of the crystal structures highlights how the spatial proximity and orientation of the TPA host and the entrapped guests influence their resulting photophysical properties and allow direct comparison of the different donor-acceptor complexes. Diffuse reflectance spectroscopy shows that upon complex formation 1·2,5-dichlorobenzoquinone exhibits a charge transfer (CT) transition. Whereas, the 1·2,1,3-benzothiadiazole complex undergoes a photoinduced electron transfer (PET) upon irradiation with 365 nm LEDs. The CT absorptions were also identified with the aid of time dependent density functional theory (TD-DFT) calculations. Cyclic voltammetry experiments show that 2,1,3-benzothiadiazole undergoes reversible reduction within the host-guest complex. Moreover, the optical band gaps of the host 1·2,5-dichlorobenzoquinone (1.66 eV), and host 1·2,1,3-benzothiadiazole (2.15 eV) complexes are significantly smaller as compared to the free host 1 material (3.19 eV). Overall, understanding this supramolecular electron transfer strategy should pave the way towards designing lower band gap inclusion complexes.

Local Measure of Quantum Effects in Quantum Dynamics.

The Madelung-de Broglie-Bohm formulation of the Schrödinger equation casts the time-evolution of a wave function as dynamics of an ensemble of quantum, or Bohmian, trajectories, interacting via the nonlocal quantum potential. This trajectory perspective gives insight into the quantumness (or classicality) of a given system due to clear partitioning of the energy into classical and quantum components. Here, we propose a system-independent measure of the quantumness of dynamics, based on the energy time-change, referred to as "quantum power". This measure is local in the coordinate space. Based on applications to model chemical systems, we argue that during the transition from the quantum to classical regime, defined as compression of quantization, the quantum features in dynamics do not "disappear" but are pushed forward in time. This feature may be used to gauge the validity of the semiclassical and other approximate dynamics approaches in applications to anharmonic systems.

From Incident Light to Persistent and Regenerable Radicals of Urea-Assembled Benzophenone Frameworks: A Structural Investigation.

Herein we probe the effects of crystalline structure and geometry on benzophenone photophysics, self-quenching, and the regenerable formation of persistent triplet radical pairs at room temperature. Radical pairs are not observed in solution but appear via an emergent pathway within the solid-state assembly. Single crystal X-ray diffraction (SC-XRD) of two sets of constitutional isomers, benzophenone bis-urea macrocycles, and methylene urea-tethered dibenzophenones are compared. Upon irradiation with 365 nm light-emitting diodes (LEDs), each forms photogenerated radicals as monitored by electron paramagnetic resonance (EPR). Once generated, the radicals exhibit half-lives from 2 to 60 days before returning to starting material without degradation. Re-exposure to light regenerates the radicals with similar efficiency. Subtle differences in the structure of the crystalline frameworks modulates the maximum concentration of photogenerated radicals, phosphorescence quantum efficiency (φ), and n-type self-quenching as observed using laser flash photolysis (LFP). These studies along with the electronic structure analysis based on the time-dependent density functional theory (TD-DFT) suggest the microenvironment surrounding benzophenone largely dictates the favorability of self-quenching or radical formation and affords insights into structure/function correlations. Advances in understanding how structure determines the excited state pathway solid-state materials undertake will aid in the design of new radical initiators, components of OLEDs, and NMR polarizing agents.

Confinement-Driven Photophysics in Cages, Covalent-Organic Frameworks, Metal-Organic Frameworks, and DNA.

Photophysics tunability through alteration of framework aperture (metal-organic framework (MOF) = variable; guest = constant) was probed for the first time in comparison with previously explored concepts (MOF = constant; guest = variable). In particular, analysis of the confinement effect on a photophysical response of integrated 5-(3-chlorobenzylidene)-2,3-dimethyl-3,5-dihydro-4H-imidazol-4-one (Cl-BI) chromophore allowed us to establish a photophysics-aperture relationship. To shed light on the observed correlation, the framework confined environment was replicated using a molecular cage, Pd6(TPT)4 (TPT = 2,4,6-tri(pyridin-4-yl)-1,3,5-triazine), thus allowing for utilization of crystallography, spectroscopy, and theoretical simulations to reveal the effect a confined space has on the chromophore's molecular conformation (including disruption of strong hydrogen bonding and novel conformer formation) and any associated changes on a photophysical response. Furthermore, the chosen Cl-oHBI@Pd6(TPT)4 (Cl-oHBI = 5-(5-chloro-2-hydroxybenzylidene)-2,3-dimethyl-3,5-dihydro-4H-imidazol-4-one, chromophore) system was applied as a tool for targeted cargo delivery of a chromophore to the confined space of DNA, which resulted in promotion of chromophore-DNA interactions through a well-established intercalation mechanism. Moreover, the developed principles were applied toward utilizing a HBI-based chromophore as a fluorescent probe on the example of macrophage cells. For the first time, suppression of non-radiative decay pathways of a chromophore was tested by anchoring the chromophore to a framework metal node, portending a potential avenue to develop an alternative to natural biomarkers. Overall, these studies are among the first attempts to demonstrate the unrevealed potential of a confined scaffold environment for tailoring a material's photophysical response.

Multidimensional Tunneling Dynamics Employing Quantum-Trajectory Guided Adaptable Gaussian Bases.

An efficient basis representation of time-dependent wavefunctions is essential for theoretical studies of high-dimensional molecular systems exhibiting large-amplitude motion. For fully coupled anharmonic systems, the complexity of a general wavefunction scales exponentially with the system size; therefore, for practical reasons, it is desirable to adapt the basis to the time-dependent wavefunction at hand. Often times on this quest for a minimal basis representation, time-dependent Gaussians are employed, in part because of their localization in both configuration and momentum spaces and also because of their direct connection to classical and semiclassical dynamics, guiding the evolution of the basis function parameters. In this work, the quantum-trajectory guided adaptable Gaussian (QTAG) bases method [ J. Chem. Theory Comput. 2020, 16, 18-34] is generalized to include correlated, i.e., non-factorizable, basis functions, and the performance of the QTAG dynamics is assessed on benchmark system/bath tunneling models of up to 20 dimensions. For the popular choice of initial conditions describing tunneling between the reactant/product wells, the minimal "semiclassical" description of the bath modes using essentially a single multidimensional basis function combined with the multi-Gaussian representation of the tunneling mode is shown to capture the dominant features of dynamics in a highly efficient manner.

A Dual Threat: Redox-Activity and Electronic Structures of Well-Defined Donor-Acceptor Fulleretic Covalent-Organic Materials.

The effect of donor (D)-acceptor (A) alignment on the materials electronic structure was probed for the first time using novel purely organic porous crystalline materials with covalently bound two- and three-dimensional acceptors. The first studies towards estimation of charge transfer rates as a function of acceptor stacking are in line with the experimentally observed drastic, eight-fold conductivity enhancement. The first evaluation of redox behavior of buckyball- or tetracyanoquinodimethane-integrated crystalline was conducted. In parallel with tailoring the D-A alignment responsible for "static" changes in materials properties, an external stimulus was applied for "dynamic" control of the electronic profiles. Overall, the presented D-A strategic design, with stimuli-controlled electronic behavior, redox activity, and modularity could be used as a blueprint for the development of electroactive and conductive multidimensional and multifunctional crystalline porous materials.

Connecting Wires: Photoinduced Electronic Structure Modulation in Metal-Organic Frameworks.

Electronic structure modulation of metal-organic frameworks (MOFs) through the connection of linker "wires" as a function of an external stimulus is reported for the first time. The established correlation between MOF electronic properties and photoisomerization kinetics as well as changes in an absorption profile is unprecedented for extended well-defined structures containing coordinatively integrated photoresponsive linkers. The presented studies were carried out on both single crystal and bulk powder with preservation of framework integrity. An LED-containing electric circuit, in which the switching behavior was driven by the changes in MOF electronic profile, was built for visualization of experimental findings. The demonstrated concept could be used as a blueprint for development of stimuli-responsive materials with dynamically controlled electronic behavior.

Stack the Bowls: Tailoring the Electronic Structure of Corannulene-Integrated Crystalline Materials.

We report the first examples of purely organic donor-acceptor materials with integrated π-bowls (πBs) that combine not only crystallinity and high surface areas but also exhibit tunable electronic properties, resulting in a four-orders-of-magnitude conductivity enhancement in comparison with the parent framework. In addition to the first report of alkyne-azide cycloaddition utilized for corannulene immobilization in the solid state, we also probed the charge transfer rate within the Marcus theory as a function of mutual πB orientation for the first time, as well as shed light on the density of states near the Fermi edge. These studies could foreshadow new avenues for πB utilization for the development of optoelectronic devices or a route for highly efficient porous electrodes.

Quantum Dynamics with Gaussian Bases Defined by the Quantum Trajectories.

Development of a general approach to construction of efficient high-dimensional bases is an outstanding challenge in quantum dynamics describing large amplitude motion of molecules and fragments. A number of approaches, proposed over the years, utilize Gaussian bases whose parameters are somehow-usually by propagating classical trajectories or by solving coupled variational equations-tailored to the shape of a wave function evolving in time. In this paper we define the time-dependent Gaussian bases through an ensemble of quantum or Bohmian trajectories, known to provide a very compact representation of a wave function due to conservation of the probability density associated with each trajectory. Though the exact numerical implementation of the quantum trajectory dynamics itself is, generally, impractical, the quantum trajectories can be obtained from the wave function expanded in a basis. The resulting trajectories are used to guide compact Gaussian bases, as illustrated on several model problems.

Dynamics in the quantum/classical limit based on selective use of the quantum potential.

A classical limit of quantum dynamics can be defined by compensation of the quantum potential in the time-dependent Schrödinger equation. The quantum potential is a non-local quantity, defined in the trajectory-based form of the Schrödinger equation, due to Madelung, de Broglie, and Bohm, which formally generates the quantum-mechanical features in dynamics. Selective inclusion of the quantum potential for the degrees of freedom deemed "quantum," defines a hybrid quantum/classical dynamics, appropriate for molecular systems comprised of light and heavy nuclei. The wavefunction is associated with all of the nuclei, and the Ehrenfest, or mean-field, averaging of the force acting on the classical degrees of freedom, typical of the mixed quantum/classical methods, is avoided. The hybrid approach is used to examine evolution of light/heavy systems in the harmonic and double-well potentials, using conventional grid-based and approximate quantum-trajectory time propagation. The approximate quantum force is defined on spatial domains, which removes unphysical coupling of the wavefunction fragments corresponding to distinct classical channels or configurations. The quantum potential, associated with the quantum particle, generates forces acting on both quantum and classical particles to describe the backreaction.

The Schrödinger equation with friction from the quantum trajectory perspective.

Similarity of equations of motion for the classical and quantum trajectories is used to introduce a friction term dependent on the wavefunction phase into the time-dependent Schrödinger equation. The term describes irreversible energy loss by the quantum system. The force of friction is proportional to the velocity of a quantum trajectory. The resulting Schrödinger equation is nonlinear, conserves wavefunction normalization, and evolves an arbitrary wavefunction into the ground state of the system (of appropriate symmetry if applicable). Decrease in energy is proportional to the average kinetic energy of the quantum trajectory ensemble. Dynamics in the high friction regime is suitable for simple models of reactions proceeding with energy transfer from the system to the environment. Examples of dynamics are given for single and symmetric and asymmetric double well potentials.

Factorized Electron-Nuclear Dynamics with an Effective Complex Potential.

We present a quantum dynamics approach for molecular systems based on wave function factorization into components describing the light and heavy particles, such as electrons and nuclei. The dynamics of the nuclear subsystem can be viewed as motion of the trajectories defined in the nuclear subspace, evolving according to the average nuclear momentum of the full wave function. The probability density flow between the nuclear and electronic subsystems is facilitated by the imaginary potential, derived to ensure a physically meaningful normalization of the electronic wave function for each configuration of the nuclei, and conservation of the probability density associated with each trajectory in the Lagrangian frame of reference. The imaginary potential, defined in the nuclear subspace, depends on the momentum variance in the nuclear coordinates averaged over the electronic component of the wave function. An effective real potential, driving the dynamics of the nuclear subsystem, is defined to minimize motion of the electronic wave function in the nuclear degrees of freedom. Illustration and the analysis of the formalism are given for a two-dimensional model system of vibrationally nonadiabatic dynamics.

Theoretical Examination of the Hydroxide Transport in Cobaltocenium-Containing Polyelectrolytes.

Polymers incorporating cobaltocenium groups have received attention as promising components of anion-exchange membranes (AEMs), exhibiting a good balance of chemical stability and high ionic conductivity. In this work, we analyze the hydroxide diffusion in the presence of cobaltocenium cations in an aqueous environment based on the molecular dynamics of model systems confined in one dimension to mimic the AEM channels. In order to describe the proton hopping mechanism, the forces are obtained from the electronic structure computed at the density-functional tight-binding level. We find that the hydroxide diffusion depends on the channel size, modulation of the electrostatic interactions by the solvation shell, and its rearrangement ability. Hydroxide diffusion proceeds via both the vehicular and structural diffusion mechanisms with the latter playing a larger role at low diffusion coefficients. The highest diffusion coefficient is observed under moderate water densities (around half the density of liquid water) when there are enough water molecules to form the solvation shell, reducing the electrostatic interaction between ions, yet there is enough space for the water rearrangements during the proton hopping. The effects of cobaltocenium separation, orientation, chemical modifications, and the role of nuclear quantum effects are also discussed.