How Water Interacts with the NOH Group: The Rotational Spectrum of the 1:1 N,N-diethylhydroxylamine·Water Complex

The rotational spectrum of the 1:1 N,N-diethylhydroxylamine-water complex has been investigated using pulsed jet Fourier transform microwave spectroscopy in the 6.5–18.5 GHz frequency region. The most stable conformer has been detected as well as the 13C monosubstituted isotopologues in natural abundance and the 18O enriched water species, allowing to determine the nitrogen nuclear quadrupole coupling constants and the molecular structure in the vibrational ground state. The molecule has a Cs symmetry and the water lies in the bc symmetry plane forming two hydrogen bonds with the NOH frame with length: dHOH·NOH = 1.974 Å and dH2O·HON = 2.096 Å. From symmetry-adapted perturbation theory calculations coupled to atoms in molecule approach, the corresponding interaction energy values are estimated to be 24 and 13 kJ·mol−1, respectively. The great strength of the intermolecular interaction involving the nitrogen atom is in agreement with the high reactivity of hydroxylamine compounds at the nitrogen site.     

Data-Driven and Multiscale Modeling of DNA-Templated Dye Aggregates

Dye aggregates are of interest for excitonic applications, including biomedical imaging, organic photovoltaics, and quantum information systems. Dyes with large transition dipole moments (μ) are necessary to optimize coupling within dye aggregates. Extinction coefficients (ε) can be used to determine the μ of dyes, and so dyes with a large ε (>150,000 M⁻¹cm⁻¹) should be engineered or identified. However, dye properties leading to a large ε are not fully understood, and low-throughput methods of dye screening, such as experimental measurements or density functional theory (DFT) calculations, can be time-consuming. In order to screen large datasets of molecules for desirable properties (i.e., large ε and μ), a computational workflow was established using machine learning (ML), DFT, time-dependent (TD-) DFT, and molecular dynamics (MD). ML models were developed through training and validation on a dataset of 8802 dyes using structural features. A Classifier was developed with an accuracy of 97% and a Regressor was constructed with an R2 of above 0.9, comparing between experiment and ML prediction. Using the Regressor, the ε values of over 18,000 dyes were predicted. The top 100 dyes were further screened using DFT and TD-DFT to identify 15 dyes with a μ relative to a reference dye, pentamethine indocyanine dye Cy5. Two benchmark MD simulations were performed on Cy5 and Cy5.5 dimers, and it was found that MD could accurately capture experimental results. The results of this study exhibit that our computational workflow for identifying dyes with a large μ for excitonic applications is effective and can be used as a tool to develop new dyes for excitonic applications.     

Thermodynamic Properties of Crystalline Cellulose Allomorphs Studied with Dispersion-Corrected Density Functional Methods

The phonon properties and thermodynamics of four crystalline cellulose allomorphs, Iα, Iβ, II, and III1, have been investigated using dispersion-corrected density functional theory (DFT). In line with experimental findings, the free energy differences between the studied cellulose allomorphs are small, less than 1 kJ/mol per atom. The calculated specific heat at constant volume (Cv) has been compared with the available experimental data in the temperature range 10–300 K. Quasiharmonic approximation has been employed to study thermodynamics and specific heat at constant pressure (Cp). For the studied temperature range of 10–400 K, the specific heat of all cellulose allomorphs shows very similar behavior. The calculated and experimental specific heat agree well at low temperatures below 100 K, but the deviation between theory and experiment increases with temperature. This may be due to increasing phonon anharmonicity as the temperature increases.     

Influence of Liquid Crystallinity and Mechanical Deformation on the Molecular Relaxations of an Auxetic Liquid Crystal Elastomer

Liquid Crystal Elastomers (LCEs) combine the anisotropic ordering of liquid crystals with the elastic properties of elastomers, providing unique physical properties, such as stimuli responsiveness and a recently discovered molecular auxetic response. Here, we determine how the molecular relaxation dynamics in an acrylate LCE are affected by its phase using broadband dielectric relaxation spectroscopy, calorimetry and rheology. Our LCE is an excellent model system since it exhibits a molecular auxetic response in its nematic state, and chemically identical nematic or isotropic samples can be prepared by cross-linking. We find that the glass transition temperatures (Tg) and dynamic fragilities are similar in both phases, and the T-dependence of the α relaxation shows a crossover at the same T* for both phases. However, for T>T*, the behavior becomes Arrhenius for the nematic LCE, but only more Arrhenius-like for the isotropic sample. We provide evidence that the latter behavior is related to the existence of pre-transitional nematic fluctuations in the isotropic LCE, which are locked in by polymerization. The role of applied strain on the relaxation dynamics and mechanical response of the LCE is investigated; this is particularly important since the molecular auxetic response is linked to a mechanical Fréedericksz transition that is not fully understood. We demonstrate that the complex Young’s modulus and the α relaxation time remain relatively unchanged for small deformations, whereas for strains for which the auxetic response is achieved, significant increases are observed. We suggest that the observed molecular auxetic response is coupled to the strain-induced out-of-plane rotation of the mesogen units, in turn driven by the increasing constraints on polymer configurations, as reflected in increasing elastic moduli and α relaxation times; this is consistent with our recent results showing that the auxetic response coincides with the emergence of biaxial order.     

Non-Peptide Opioids Differ in Effects on Mu-Opioid (MOP) and Serotonin 1A (5-HT1A) Receptors Heterodimerization and Cellular Effectors (Ca2+, ERK1/2 and p38) Activation

The importance of the dynamic interplay between the opioid and the serotonin neuromodulatory systems in chronic pain is well recognized. In this study, we investigated whether these two signalling pathways can be integrated at the single-cell level via direct interactions between the mu-opioid (MOP) and the serotonin 1A (5-HT_1A) receptors. Using fluorescence cross-correlation spectroscopy (FCCS), a quantitative method with single-molecule sensitivity, we characterized in live cells MOP and 5-HT_1A interactions and the effects of prolonged (18 h) exposure to selected non-peptide opioids: morphine, codeine, oxycodone and fentanyl, on the extent of these interactions. The results indicate that in the plasma membrane, MOP and 5-HT_1A receptors form heterodimers that are characterized with an apparent dissociation constant Kdapp = (440 ± 70) nM). Prolonged exposure to all non-peptide opioids tested facilitated MOP and 5-HT_1A heterodimerization and stabilized the heterodimer complexes, albeit to a different extent: Kd, Fentanylapp = (80 ± 70) nM), Kd,Morphineapp = (200 ± 70) nM, Kd, Codeineapp = (100 ± 70) nM and Kd, Oxycodoneapp = (200 ± 70) nM. The non-peptide opioids differed also in the extent to which they affected the mitogen-activated protein kinases (MAPKs) p38 and the extracellular signal-regulated kinase (Erk1/2), with morphine, codeine and fentanyl activating both pathways, whereas oxycodone activated p38 but not ERK1/2. Acute stimulation with different non-peptide opioids differently affected the intracellular Ca²⁺ levels and signalling dynamics. Hypothetically, targeting MOP–5-HT_1A heterodimer formation could become a new strategy to counteract opioid induced hyperalgesia and help to preserve the analgesic effects of opioids in chronic pain.     

The Structural Rule Distinguishing a Superfold: A Case Study of Ferredoxin Fold and the Reverse Ferredoxin Fold

Superfolds are folds commonly observed among evolutionarily unrelated multiple superfamilies of proteins. Since discovering superfolds almost two decades ago, structural rules distinguishing superfolds from the other ordinary folds have been explored but remained elusive. Here, we analyzed a typical superfold, the ferredoxin fold, and the fold which reverses the N to C terminus direction from the ferredoxin fold as a case study to find the rule to distinguish superfolds from the other folds. Though all the known structural characteristics for superfolds apply to both the ferredoxin fold and the reverse ferredoxin fold, the reverse fold has been found only in a single superfamily. The database analyses in the present study revealed the structural preferences of αβ- and βα-units; the preferences separate two α-helices in the ferredoxin fold, preventing their collision and stabilizing the fold. In contrast, in the reverse ferredoxin fold, the preferences bring two helices near each other, inducing structural conflict. The Rosetta folding simulations suggested that the ferredoxin fold is physically much more realizable than the reverse ferredoxin fold. Therefore, we propose that minimal structural conflict or minimal frustration among secondary structures is the rule to distinguish a superfold from ordinary folds. Intriguingly, the database analyses revealed that a most stringent structural rule in proteins, the right-handedness of the βαβ-unit, is broken in a set of structures to prevent the frustration, suggesting the proposed rule of minimum frustration among secondary structural units is comparably strong as the right-handedness rule of the βαβ-unit.     

Anharmonic Vibrational Frequencies of Water Borane and Associated Molecules

Water borane (BH3OH2) and borinic acid (BH2OH) have been proposed as intermediates along the pathway of hydrogen generation from simple reactants: water and borane. However, the vibrational spectra for neither water borane nor borinic acid has been investigaged experimentally due to the difficulty of isolating them in the gas phase, making accurate quantum chemical predictions for such properties the most viable means of their determination. This work presents theoretical predictions of the full rotational and fundamental vibrational spectra of these two potentially application-rich molecules using quartic force fields at the CCSD(T)-F12b/cc-pCVTZ-F12 level with additional corrections included for the effects of scalar relativity. This computational scheme is further benchmarked against the available gas-phase experimental data for the related borane and HBO molecules. The differences are found to be within 3 cm−1 for the fundamental vibrational frequencies and as close as 15 MHz in the B0 and C0 principal rotational constants. Both BH2OH and BH3OH2 have multiple vibrational modes with intensities greater than 100 km mol−1, namely ν2 and ν4 in BH2OH, and ν1, ν3, ν4, ν9, and ν13 in BH3OH2. Finally, BH3OH2 has a large dipole moment of 4.24 D, which should enable it to be observable by rotational spectroscopy, as well.     

Metastable Kitaev Magnets

Nearly two decades ago, Alexei Kitaev proposed a model for spin-1/2 particles with bond-directional interactions on a two-dimensional honeycomb lattice which had the potential to host a quantum spin-liquid ground state. This work initiated numerous investigations to design and synthesize materials that would physically realize the Kitaev Hamiltonian. The first generation of such materials, such as Na2IrO3, α-Li2IrO3, and α-RuCl3, revealed the presence of non-Kitaev interactions such as the Heisenberg and off-diagonal exchange. Both physical pressure and chemical doping were used to tune the relative strength of the Kitaev and competing interactions; however, little progress was made towards achieving a purely Kitaev system. Here, we review the recent breakthrough in modifying Kitaev magnets via topochemical methods that has led to the second generation of Kitaev materials. We show how structural modifications due to the topotactic exchange reactions can alter the magnetic interactions in favor of a quantum spin-liquid phase.     

Removal of Carbamazepine in Aqueous Solution by CoS2/Fe2+/PMS Process

Carbamazepine (CBZ), as a typical pharmaceutical and personal care product (PPCP), cannot be efficiently removed by the conventional drinking water and wastewater treatment process. In this work, the CoS₂/Fe²⁺/PMS process was applied for efficient elimination of CBZ. The CBZ removal efficiency of CoS₂/Fe²⁺/PMS was 2.5 times and 23 times higher than that of CoS₂/PMS and Fe²⁺/PMS, respectively. The intensity of DMPO-HO• and DMPO-SO4•− followed the order of Fe²⁺/PMS < CoS₂/PMS < CoS₂/Fe²⁺/PMS, also suggesting the CoS₂/Fe²⁺/PMS process has the highest oxidation activity. The effects of reaction conditions (e.g., CoS₂ dosage, Fe²⁺ concentration, PMS concentration, initial CBZ concentration, pH, temperature) and water quality parameters (e.g., SO42−, NO3−, H2PO4−, Cl−, NH4+, humic acid) on the degradation of CBZ were also studied. Response surface methodology analysis was carried out to obtain the best conditions for the removal of CBZ, which are: Fe²⁺ = 70 µmol/L, PMS = 240 µmol/L, CoS₂ = 0.59 g/L. The sustainability test demonstrated that the repeated use of CoS₂ for 8 successive cycles resulted in little function decrease (<10%). These findings suggest that CoS₂/Fe²⁺/PMS may be a promising method for advanced treatment of tailwater from sewage treatment plant.     

Galvano-Fenton Engineering Solution with Spontaneous Catalyst’s Generation from Waste: Experimental Efficiency,Parametric Analysis and Modeling Interpretation Applied to a Clean Technology for Dyes Degradation in Water

In this paper, the degradation of the diazo dye naphthol blue black (NBB) using the Galvano-Fenton process is studied experimentally and numerically. The simulations are carried out based on the anodic, cathodic, and 34 elementary reactions evolving in the electrolyte, in addition to the oxidative attack of NBB by HO• at a constant rate of 3.35×107 mol−1·m3·s−1 during the initiation stage of the chain reactions. The selection of the operating conditions including the pH of the electrolyte, the stirring speed, and the electrodes disposition is performed by assessing the kinetics of NBB degradation; these parameters are set to 3, 350 rpm and a parallel disposition with a 3 cm inter-electrode distance, respectively. The kinetics of Fe(III) in the electrolyte were monitored using the principles of Fricke dosimetry and simulated numerically. The model showed more than a 96% correlation with the experimental results in both the blank test and the presence of the dye. The effects of H2O2 and NBB concentrations on the degradation of the dye were examined jointly with the evolution of the simulated H2O2, Fe2+, and HO• concentrations in the electrolyte. The model demonstrated a good correlation with the experimental results in terms of the initial degradation rates, with correlation coefficients exceeding 98%.     

On the Coexistence of the Carbene⋯H-D Hydrogen Bond and Other Accompanying Interactions in Forty Dimers of N-Heterocyclic-Carbenes (I,IMe2, IiPr2, ItBu2, IMes2, IDipp2, IAd2; I = imidazol-2-ylidene) and Some Fundamental Proton Donors (HF,HCN, H2O,MeOH, NH3)

The subject of research is forty dimers formed by imidazol-2-ylidene (I) or its derivative (IR2) obtained by replacing the hydrogen atoms in both N-H bonds with larger important and popular substituents of increasing complexity (methyl = Me, iso-propyl = iPr, tert-butyl = tBu, phenyl = Ph, mesityl = Mes, 2,6-diisopropylphenyl = Dipp, 1-adamantyl = Ad) and fundamental proton donor (HD) molecules (HF, HCN, H2O, MeOH, NH3). While the main goal is to characterize the generally dominant C⋯H-D hydrogen bond engaging a carbene carbon atom, an equally important issue is the often omitted analysis of the role of accompanying secondary interactions. Despite the often completely different binding possibilities of the considered carbenes, and especially HD molecules, several general trends are found. Namely, for a given carbene, the dissociation energy values of the IR2⋯HD dimers increase in the following order: NH3< H2O < HCN ≤ MeOH ≪ HF. Importantly, it is found that, for a given HD molecule, IDipp2 forms the strongest dimers. This is attributed to the multiplicity of various interactions accompanying the dominant C⋯H-D hydrogen bond. It is shown that substitution of hydrogen atoms in both N-H bonds of the imidazol-2-ylidene molecule by the investigated groups leads to stronger dimers with HF, HCN, H2O or MeOH. The presented results should contribute to increasing the knowledge about the carbene chemistry and the role of intermolecular interactions, including secondary ones.     

Spectral Signatures of Hydrogen Thioperoxide (HOSH) and Hydrogen Persulfide (HSSH): Possible Molecular Sulfur Sinks in the Dense ISM

For decades, sulfur has remained underdetected in molecular form within the dense interstellar medium (ISM), and somewhere a molecular sulfur sink exists where it may be hiding. With the discovery of hydrogen peroxide (HOOH) in the ISM in 2011, a natural starting point may be found in sulfur-bearing analogs that are chemically similar to HOOH: hydrogen thioperoxide (HOSH) and hydrogen persulfide (HSSH). The present theoretical study couples the accuracy in the anharmonic fundamental vibrational frequencies from the explicitly correlated coupled cluster theory with the accurate rotational constants provided by canonical high-level coupled cluster theory to produce rovibrational spectra for use in the potential observation of HOSH and HSSH. The ν6 mode for HSSH at 886.1 cm−1 is within 0.2 cm−1 of the gas-phase experiment, and the B0 rotational constant for HSSH of 6979.5 MHz is within 9.0 MHz of the experimental benchmarks, implying that the unknown spectral features (such as the first overtones and combination bands) provided herein are similarly accurate. Notably, a previous experimentally-attributed 2ν1 mode, at 7041.8 cm−1, has been reassigned to the ν1+ν5 combination band based on the present work’s ν1+ν5 value at 7034.3 cm−1. The most intense vibrational transitions for each molecule are the torsions, with HOSH having a more intense transition of 72 km/mol compared to HSSH’s intensity of 14 km/mol. Furthermore, HOSH has a larger net dipole moment of 1.60 D compared to HSSH’s 1.15 D. While HOSH may be the more likely candidate of the two for possible astronomical observation via vibrational spectroscopy due to the notable difference in their intensities, both HSSH and HOSH have large enough net dipole moments to be detectable by rotational spectroscopy to discover the role these molecules may have as possible molecular sulfur sinks in the dense ISM.     

Columnar Aggregates of Azobenzene Stars: Exploring Intermolecular Interactions,Structure, and Stability in Atomistic Simulations

We present a simulation study of supramolecular aggregates formed by three-arm azobenzene (Azo) stars with a benzene-1,3,5-tricarboxamide (BTA) core in water. Previous experimental works by other research groups demonstrate that such Azo stars assemble into needle-like structures with light-responsive properties. Disregarding the response to light, we intend to characterize the equilibrium state of this system on the molecular scale. In particular, we aim to develop a thorough understanding of the binding mechanism between the molecules and analyze the structural properties of columnar stacks of Azo stars. Our study employs fully atomistic molecular dynamics (MD) simulations to model pre-assembled aggregates with various sizes and arrangements in water. In our detailed approach, we decompose the binding energies of the aggregates into the contributions due to the different types of non-covalent interactions and the contributions of the functional groups in the Azo stars. Initially, we investigate the origin and strength of the non-covalent interactions within a stacked dimer. Based on these findings, three arrangements of longer columnar stacks are prepared and equilibrated. We confirm that the binding energies of the stacks are mainly composed of π–π interactions between the conjugated parts of the molecules and hydrogen bonds formed between the stacked BTA cores. Our study quantifies the strength of these interactions and shows that the π–π interactions, especially between the Azo moieties, dominate the binding energies. We clarify that hydrogen bonds, which are predominant in BTA stacks, have only secondary energetic contributions in stacks of Azo stars but remain necessary stabilizers. Both types of interactions, π–π stacking and H-bonds, are required to maintain the columnar arrangement of the aggregates.     

Studies of the Formation of Inclusion Complexes Derivatives of Cinnamon Acid with α-Cyclodextrin in a Wide Range of Temperatures Using Conductometric Methods

The electrical conductivities of aqueous solutions of sodium salts of trans-4-hydroxycinnamic acid (trans-p-coumaric acid), trans-3,4-dihydroxycinnamic acid (trans-caffeic acid), trans-4-hydroxy-3-methoxycinnamic acid, (trans-ferulic acid) and trans-3-phenylacrylic acid (trans-cinnamic acid) with α-cyclodextrin were measured in the temperature range of 288.15 K–318.15 K. For the first time in the literature, using the limiting molar conductivity (Λmo) obtained from conductivity measurements, the values of the complexation constants (K_f) of the salts of phenolic acid derivatives with α-cyclodextrin were determined using a modified low concentration chemical model (IcCM). An attempt was also made to analyze the individual thermodynamic functions ΔGo, ΔHo and ΔSo describing the complexation process as a function of temperature changes. The obtained results show that the process of formation of inclusion complexes is exothermic and is spontaneous.     

The Relationship between the IC50 Values and the Apparent Inhibition Constant in the Study of Inhibitors of Tyrosinase Diphenolase Activity Helps Confirm the Mechanism of Inhibition

Tyrosinase is the enzyme involved in melanization and is also responsible for the browning of fruits and vegetables. Control of its activity can be carried out using inhibitors, which is interesting in terms of quantitatively understanding the action of these regulators. In the study of the inhibition of the diphenolase activity of tyrosinase, it is intriguing to know the strength and type of inhibition. The strength is indicated by the value of the inhibition constant(s), and the type can be, in a first approximation: competitive, non-competitive, uncompetitive and mixed. In this work, it is proposed to calculate the degree of inhibition (iD), varying the concentration of inhibitor to a fixed concentration of substrate, L-dopa (D). The non-linear regression adjustment of iD with respect to the initial inhibitor concentration [I]0 allows for the calculation of the inhibitor concentration necessary to inhibit the activity by 50%, at a given substrate concentration (IC₅₀), thus avoiding making interpolations between different values of iD. The analytical expression of the IC₅₀, for the different types of inhibition, are related to the apparent inhibition constant (KIapp). Therefore, this parameter can be used: (a) To classify a series of inhibitors of an enzyme by their power. Determining these values at a fixed substrate concentration, the lower IC₅₀, the more potent the inhibitor. (b) Checking an inhibitor for which the type and the inhibition constant have been determined (using the usual methods), must confirm the IC₅₀ value according to the corresponding analytical expression. (c) The type and strength of an inhibitor can be analysed from the study of the variation in iD and IC₅₀ with substrate concentration. The dependence of IC₅₀ on the substrate concentration allows us to distinguish between non-competitive inhibition (iD does not depend on [D]0) and the rest. In the case of competitive inhibition, this dependence of iD on [D]0 leads to an ambiguity between competitive inhibition and type 1 mixed inhibition. This is solved by adjusting the data to the possible equations; in the case of a competitive inhibitor, the calculation of KI1app is carried out from the IC₅₀ expression. The same occurs with uncompetitive inhibition and type 2 mixed inhibition. The representation of iD vs. n, with n=[D]0/KmD, allows us to distinguish between them. A hyperbolic iD vs. n representation that passes through the origin of coordinates is a characteristic of uncompetitive inhibition; the calculation of KI2app is immediate from the IC₅₀ value. In the case of mixed inhibitors, the values of the apparent inhibition constant of meta-tyrosinase (Em) and oxy-tyrosinase (Eox), KI1app and the apparent inhibition constant of metatyrosinase/Dopa complexes (EmD) and oxytyrosinase/Dopa (EoxD), KI2app are obtained from the dependence of iD vs. n, and the results obtained must comply with the IC₅₀ value.     

Does Hamiltonian Replica Exchange via Lambda-Hopping Enhance the Sampling in Alchemical Free Energy Calculations?

In the context of computational drug design, we examine the effectiveness of the enhanced sampling techniques in state-of-the-art free energy calculations based on alchemical molecular dynamics simulations. In a paradigmatic molecule with competition between conformationally restrained E and Z isomers whose probability ratio is strongly affected by the coupling with the environment, we compare the so-called λ-hopping technique to the Hamiltonian replica exchange methods assessing their convergence behavior as a function of the enhanced sampling protocols (number of replicas, scaling factors, simulation times). We found that the pure λ-hopping, commonly used in solvation and binding free energy calculations via alchemical free energy perturbation techniques, is ineffective in enhancing the sampling of the isomeric states, exhibiting a pathological dependence on the initial conditions. Correct sampling can be restored in λ-hopping simulation by the addition of a “hot-zone” scaling factor to the λ-stratification (FEP⁺ approach), provided that the additive hot-zone scaling factors are tuned and optimized using preliminary ordinary replica-exchange simulation of the end-states.     

Give or Take: Effects of Electron-Accepting/-Withdrawing Groups in Red-Fluorescent BODIPY Molecular Rotors

Mapping microviscosity, temperature, and polarity in biosystems is an important capability that can aid in disease detection. This can be achieved using fluorescent sensors based on a green-emitting BODIPY group. However, red fluorescent sensors are desired for convenient imaging of biological samples. It is known that phenyl substituents in the β position of the BODIPY core can shift the fluorescence spectra to longer wavelengths. In this research, we report how electron-withdrawing (EWG) and -donating (EDG) groups can change the spectral and sensory properties of β-phenyl-substituted BODIPYs. We present a trifluoromethyl-substituted (EWG) conjugate with moderate temperature sensing properties and a methoxy-substituted (EDG) molecule that could be used as a lifetime-based polarity probe. In this study, we utilise experimental results of steady-state and time-resolved fluorescence, as well as quantum chemical calculations using density functional theory (DFT). We also explain how the energy barrier height (Ea) for non-radiative relaxation affects the probe’s sensitivity to temperature and viscosity and provide appropriate Ea ranges for the best possible sensitivity to viscosity and temperature.     

Naphthazarin Derivatives in the Light of Intra- and Intermolecular Forces

Our long-term investigations have been devoted the characterization of intramolecular hydrogen bonds in cyclic compounds. Our previous work covers naphthazarin, the parent compound of two systems discussed in the current work: 2,3-dimethylnaphthazarin (1) and 2,3-dimethoxy-6-methylnaphthazarin (2). Intramolecular hydrogen bonds and substituent effects in these compounds were analyzed on the basis of Density Functional Theory (DFT), Møller–Plesset second-order perturbation theory (MP2), Coupled Clusters with Singles and Doubles (CCSD) and Car-Parrinello Molecular Dynamics (CPMD). The simulations were carried out in the gas and crystalline phases. The nuclear quantum effects were incorporated a posteriori using the snapshots taken from ab initio trajectories. Further, they were used to solve a vibrational Schrödinger equation. The proton reaction path was studied using B3LYP, ωB97XD and PBE functionals with a 6-311++G(2d,2p) basis set. Two energy minima (deep and shallow) were found, indicating that the proton transfer phenomena could occur in the electronic ground state. Next, the electronic structure and topology were examined in the molecular and proton transferred (PT) forms. The Atoms In Molecules (AIM) theory was employed for this purpose. It was found that the hydrogen bond is stronger in the proton transferred (PT) forms. In order to estimate the dimers’ stabilization and forces responsible for it, the Symmetry-Adapted Perturbation Theory (SAPT) was applied. The energy decomposition revealed that dispersion is the primary factor stabilizing the dimeric forms and crystal structure of both compounds. The CPMD results showed that the proton transfer phenomena occurred in both studied compounds, as well as in both phases. In the case of compound 2, the proton transfer events are more frequent in the solid state, indicating an influence of the environmental effects on the bridged proton dynamics. Finally, the vibrational signatures were computed for both compounds using the CPMD trajectories. The Fourier transformation of the autocorrelation function of atomic velocity was applied to obtain the power spectra. The IR spectra show very broad absorption regions between 700 cm−1–1700 cm−1 and 2300 cm−1–3400 cm−1 in the gas phase and 600 cm−1–1800 cm−1 and 2200 cm−1–3400 cm−1 in the solid state for compound 1. The absorption regions for compound 2 were found as follows: 700 cm−1–1700 cm−1 and 2300 cm−1–3300 cm−1 for the gas phase and one broad absorption region in the solid state between 700 cm−1 and 3100 cm−1. The obtained spectroscopic features confirmed a strong mobility of the bridged protons. The inclusion of nuclear quantum effects showed a stronger delocalization of the bridged protons.     

Fluorescent Liquid Tetrazines

Tetrazines with branched alkoxy substituents are liquids at ambient temperature that despite the high chromophore density retain the bright orange fluorescence that is characteristic of this exceptional fluorophore. Here, we study the photophysical properties of a series of alkoxy-tetrazines in solution and as neat liquids. We also correlate the size of the alkoxy substituents with the viscosity of the liquids. We show using time-resolved spectroscopy that intersystem crossing is an important decay pathway competing with fluorescence, and that its rate is higher for 3,6-dialkoxy derivatives than for 3-chloro-6-alkoxytetrazines, explaining the higher fluorescence quantum yields for the latter. Quantum chemical calculations suggest that the difference in rate is due to the activation energy required to distort the tetrazine core such that the nπ*S1 and the higher-lying ππ*T2 states cross, at which point the spin-orbit coupling exceeding 10 cm−1 allows for efficient intersystem crossing to occur. Femtosecond time-resolved anisotropy studies in solution allow us to measure a positive relationship between the alkoxy chain lengths and their rotational correlation times, and studies in the neat liquids show a fast decay of the anisotropy consistent with fast exciton migration in the neat liquid films.