Providing theoretical insight into the role of symmetry in the photoisomerization mechanism of a non-symmetric dithienylethene photoswitch

Dithienylethene (DTE) molecular photoswitches have shown to be excellent candidates in the design of efficient optoelectronic devices, due to their high photoisomerization quantum yield (QY), for which symmetry is suggested to play a crucial role. Here, we present a theoretical study on the photochemistry of a non-symmetric dithienylethene photoswitch, with a special emphasis on the effect of asymmetric substitution on the photocyclization and photoreversion mechanisms. We used the Spin-Flip Time Dependent Density Functional Theory (SF-TDDFT) method to locate and characterize the main structures (conical intersections and minima) of the ground state and the first two excited states, S1 and S2, along the ring-opening/closure reaction coordinate of the photocyclization and photoreversion processes, and to identify the important coordinates governing the radiationless decay pathways. Our results suggest that while the main features that characterize the photoisomerization of symmetric DTEs are also present for the photoisomerization of the non-symmetric DTE, the lower energy barrier on S1 along the cycloreversion reaction speaks in favor of a more efficient and therefore a higher cycloreversion QY for the non-symmetric DTEs, making them a better candidate for molecular optoelectronic devices than their symmetric counterparts

Providing theoretical insight into the role of symmetry in the photoisomerization mechanism of a non-symmetric dithienylethene photoswitch

Dithienylethene (DTE) molecular photoswitches have shown to be excellent candidates in the design of efficient optoelectronic devices, due to their high photoisomerization quantum yield (QY), for which symmetry is suggested to play a crucial role. Here, we present a theoretical study on the photochemistry of a non-symmetric dithienylethene photoswitch, with a special emphasis on the effect of asymmetric substitution on the photocyclization and photoreversion mechanisms. We used the Spin-Flip Time Dependent Density Functional Theory (SF-TDDFT) method to locate and characterize the main structures (conical intersections and minima) of the ground state and the first two excited states, S1 and S2, along the ring-opening/closure reaction coordinate of the photocyclization and photoreversion processes, and to identify the important coordinates governing the radiationless decay pathways. Our results suggest that while the main features that characterize the photoisomerization of symmetric DTEs are also present for the photoisomerization of the non-symmetric DTE, the lower energy barrier on S1 along the cycloreversion reaction speaks in favor of a more efficient and therefore a higher cycloreversion QY for the non-symmetric DTEs, making them a better candidate for molecular optoelectronic devices than their symmetric counterparts.</p

Effect of stacking interactions on charge transfer states in photoswitches interacting with ion channels

The activity of ion channels can be reversibly photo-controlled via the binding of molecular photoswitches, often based on an azobenzene scaffold. Those azobenzene derivatives interact with aromatic residues of the protein via stacking interactions. In the present work, the effect of face-to-face and t-shaped stacking interactions on the excited state electronic structure of azobenzene and p-diaminoazobenzene integrated into the Na V1.4 channel is computationally investigated. The formation of a charge transfer state, caused by electron transfer from the protein to the photoswitches, is observed. This state is strongly red shifted when the interaction takes place in a face-to-face orientation and electron donating groups are present on the aromatic ring of the amino acids. The low-energy charge transfer state can interfere with the photoisomerization process after excitation to the bright state by leading to the formation of radical species. </p

The impact of anion elements on the engineering of the electronic and T optical characteristics of the two dimensional monolayer janus MoSSe for nanoelectronic device applications

Two-dimensional (2D) materials have gained prominent attention in the nano-electronics arena, owing to their tunable electronic and optical features. Here, the physical properties of a janus MoSSe monolayer are examined upon the chemically co-doping of S/Se sites by non-metallic and halogen elements (C, Si, N, P, As, and F) employing first-principles calculations. Accordingly, an alteration of both the upper valence and the lower conduction states is revealed for janus MoSSe monolayer upon the replacement of both S and Se anion host atoms by sp-elements (C, Si, N, P, As, and F). A shift in the lowest conduction band underneath the Fermi energy level (EF) occurs in janus MoSSe monolayer when both S and Se elements are replaced by (F, F) atoms. This effectively conducted to a system with an n-type character. In contrast, the highest valence bands moved upward EF owing to the co-doping effect of C, Si, N, P, and As atoms on the janus MoSSe monolayer with p-type nature. The key features of the optical spectra, such as the optical absorption, reflectivity, and electron loss functions of the co-doped janus MoSSe monolayer are inspected. Our results imply a modification in the low-energy photon regime of the co-doped janus MoSSe monolayer at S and Se host atoms by non-metallic sp-elements comparatively to the free-standing monolayer. A reduction in the optical absorption and an increase in the reflectivity at low-energy photon window are detected when the janus MoSSe monolayer is co-doped by (C, Si), (N, P), (P,As), and (F,F) elements, respectively at S and Se chalcogen atoms. The current study infers that the co-doping S and Se sites of janus MoSSe monolayers, with sp- elements, can be beneficial in the future applications of 2D materials for the field-effect transistors and nano-electronic devices

Sesquicarbene Complexes:Bonding at the Interface Between M-C Single Bonds and M=C Double Bonds

Allylic dimetalated complexes [M2C(vinyl)]+ (M = Au(IPr) and Cr(CO)5 -) incorporate a new coordination mode of carbon. Digold complexes of this type have recently been detected experimentally. The intrinsic bond orbitals, partial charges, and structural parameters of the gold complexes and of chromium analogs were studied computationally and compared to those of the respective monometalated species and hydrocarbons. This showed that such digold complexes have a carbene character at both Au-C bonds comparable to typical carbene complexes of gold. Dichromium complexes with their stronger π-backdonation compete for interaction with carbon's π-orbital; each of the chromium atoms partakes in double bonding that is significant but weaker than that in the carbene analogs. Containing two M-C bonds on the interface between single and double bonds, these bridged complexes can be conceived as "sesquicarbene complexes". The π-system acted in a very adaptive manner and employed additional stabilization of the vinyl system only where needed. Significant carbene character is found simultaneously in both M-C bonds at the same carbon center. The discovery of these complexes with relatively strong double bond character between one carbon and two metal atoms could bring unusual single-carbon-centered organometallic cascade reactions to the horizon

Influence of the crystal packing in singlet fission:One step beyond the gas phase approximation

Singlet fission (SF), a multiexciton generation process, has been proposed as an alternative to enhance the performance of solar cells. The gas phase dimer model has shown its utility to study this process, but it does not always cover all the physics and the effect of the surrounding atoms has to be included in such cases. In this contribution, we explore the influence of crystal packing on the electronic couplings, and on the so-called exciton descriptors and electron-hole correlation plots. We have studied three tetracene dimers extracted from the crystal structure, as well as several dimers and trimers of the α and β polymorphs of 1,3-diphenylisobenzofuran (DPBF). These polymorphs show different SF yields. Our results highlight that the character of the excited states of tetracene depends on both the mutual disposition of molecules and inclusion of the environment. The latter does however not change significantly the interpretation of the SF mechanism in the studied systems. For DPBF, we establish how the excited state analysis is able to pinpoint differences between the polymorphs. We observe strongly bound correlated excitons in the β polymorph which might hinder the formation of the 1TT state and, consequently, explain its low SF yield

Theoretical insights into the effect of size and substitution patterns of azobenzene derivatives on the DNA G-quadruplex

Introducing photoswitches into the DNA G-quadruplex provides excellent opportunities to control folding and unfolding of these assemblies, demonstrating their potential in the development of novel nanodevices with medical and nanotechnology applications. Using a quantum mechanics/molecular mechanics (QM/MM) scheme, we carried out a series of simulations to identify the effect of the size and substitution patterns of three azobenzene derivatives (AZ1, AZ2 and AZ3) on the excitation energies of the two lowest excited states of the smallest photoswitchable G-quadruplex reported to date. We demonstrated that the size and the substitution pattern do not affect the ultrafast cis-trans photoiomerization mechanism of the azobenzene derivatives significantly, in agreement with the experiment. However, molecular dynamics simulations revealed that while AZ2 and AZ3 G-quadruplexes are structurally stable during the simulations, the AZ1 G-quadruplex undergoes larger structural changes and shows two ground state populations that differ in the azobenzene backbone adopting two different conformations. AZ1, with para-para substitution pattern, provides more flexibility to the whole G-quadruplex structure compared to AZ2 and AZ3, and can thus facilitate the photoisomerization reaction between a nonpolymorphic, stacked, tetramolecular G-quadruplex and an unstructured state after trans-cis isomerization occurring in a longer time dynamics, in agreement with the experimental findings. The QM/MM simulations of the absorption spectra indicated that the thermal fluctuation plays a more crucial role in the main absorption band of the azobenzene derivatives than the inclusion of the G-quadruplex, implying that the influence of the G-quadruplex environment is minimal. We propose that the latter is attributed to the position of the azobenzene linkers in the G-quadruplexes, i.e. the edgewise loops containing the azobenzene moieties that are located above the G-quartets, not being fully embedded inside or involved in the stacked structure. Our theoretical findings provide support to a recent study of the photoresponsive formation of photoswitchable G-quadruplex motifs

Influence of the Environment on Shaping the Absorption of Monomeric Infrared Fluorescent Proteins

Infrared fluorescent proteins (iRFPs) are potential candidates for deep-tissue in vivo imaging. Here, we provide molecular-level insights into the role of the protein environment in the structural stability of the chromophore within the protein binding pocket through the flexible hydrogen-bonding network using molecular dynamics simulation. Furthermore, we present systematic excited-state analysis to characterize the nature of the first two excited states and the role of the environment in shaping the nature of the chromophore's excited states within the hybrid quantum mechanics/molecular mechanics framework. Our results reveal that the environment red-shifts the absorption of the chromophore by about 0.32 eV compared to the isolated counterpart, and besides the structural stability, the protein environment does not alter the nature of the excited state of the chromophore significantly. Our study contributes to the fundamental understanding of the excited-state processes of iRFPs in a complex environment and provides a design principle for developing iRFPs with desired spectral properties

PySurf:A Framework for Database Accelerated Direct Dynamics

The greatest restriction to the theoretical study of the dynamics of photoinduced processes is computationally expensive electronic structure calculations. Machine learning algorithms have the potential to reduce the number of these computations significantly. Here, PySurf is introduced as an innovative code framework, which is specifically designed for rapid prototyping and development tasks for data science applications in computational chemistry. It comes with powerful Plugin and Workflow engines, which allows intuitive customization for individual tasks. Data is automatically stored through the database framework, which enables additional interpolation of properties in previously evaluated regions of the conformational space. To illustrate the potential of the framework, a code for nonadiabatic surface hopping simulations based on the Landau-Zener algorithm is presented here. Deriving gradients from the interpolated potential energy surfaces allows for full-dimensional nonadiabatic surface hopping simulations using only adiabatic energies (energy only). Simulations of a pyrazine model and ab initio-based calculations of the SO2 molecule show that energy-only calculations with PySurf are able to correctly predict the nonadiabatic dynamics of these prototype systems. The results reveal the degree of sophistication, which can be achieved by the database accelerated energy-only surface hopping simulations being competitive to commonly used semiclassical approaches

HOAX: A Hyperparameter Optimization Algorithm Explorer for Neural Networks

Computational chemistry has become an important tool to predict and understand molecular properties and reactions. Even though recent years have seen a significant growth in new algorithms and computational methods that speed up quantum chemical calculations, the bottleneck for trajectory-based methods to study photoinduced processes is still the huge number of electronic structure calculations. In this work, we present an innovative solution, in which the amount of electronic structure calculations is drastically reduced, by employing machine learning algorithms and methods borrowed from the realm of artificial intelligence. However, applying these algorithms effectively requires finding optimal hyperparameters, which remains a challenge itself. Here we present an automated user-friendly framework, HOAX, to perform the hyperparameter optimization for neural networks, which bypasses the need for a lengthy manual process. The neural network generated potential energy surfaces (PESs) reduces the computational costs compared to the ab initio-based PESs. We perform a comparative investigation on the performance of different hyperparameter optimiziation algorithms, namely grid search, simulated annealing, genetic algorithm, and bayesian optimizer in finding the optimal hyperparameters necessary for constructing the well-performing neural network in order to fit the PESs of small organic molecules. Our results show that this automated toolkit not only facilitate a straightforward way to perform the hyperparameter optimization but also the resulting neural networks-based generated PESs are in reasonable agreement with the ab initio-based PESs.Comment: 18 page