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.     

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.     

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.     

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%.     

Efficiency of the same neat silica column in hydrophilic interaction chromatography and per aqueous liquid chromatography

The dependencies on the mobile phase flow velocity of the efficiency of a column packed with shell particles of neat porous silica (Halo) was measured under two different sets of experimental conditions. These conditions corresponded to the retention mechanisms of per   aqueous liquid chromatography (PALC) at low acetonitrile concentrations and of hydrophilic interaction chromatography (HILIC) at high acetonitrile concentrations. The results are compared. Small amounts of a diluted solution of caffeine were injected in order to record the chromatograms under strictly linear conditions. These efficiencies were measured in both water-rich (PALC retention mechanism) and acetonitrile-rich (HILIC mechanism) mobile phases for the same retention factors, between 0.25 and 2.5. The mobile phases were mixtures of acetonitrile and water containing neither supporting salt nor buffer component. At low retention factors, the efficiency of caffeine is better in the PALC than in the HILIC mode. For k^′=0.5$k^{\prime }=0.5$, the minimum reduced height equivalent to a theoretical plate (HETP) is close to 2.5 in PALC while it exceeds 5 in HILIC. The converse is true for high retention factors. For k^′>2.5$k^{\prime }>2.5$, the HETP is lower in HILIC than in PALC, because the major contribution to band broadening and peak tailing in this latter mode originates from the heterogeneous thermodynamics of retention and eventually restricts column performance in PALC. Most interestingly, the reduced HETP measured in HILIC for caffeine never falls below 4. This suggests that the mass transfer of caffeine between the multilayer adsorbed phase (due to the interactions of the strong solvent and the silanol groups) and the acetonitrile-rich bulk eluent is slow.     

Solubility of carbon dioxide in the solvent system (2-amino-2-methyl-1-propanol + sulfolane + water)

New sets of data for the solubility of CO₂ in the amine solvent system of 2-amino-2-methyl-1-propanol (1) + sulfolane (2) + water (3) were presented in this work. The measurements were done at temperatures of 313.2, 333.2, 353.2, and 373.2 K and CO₂ partial pressures up to 193 kPa. The investigated compositions were as follows: (i) w₁=16.5%$w_1=16.5\%$, w₂=32.2%$w_2=32.2\%$; (ii) w₁=8.2%$w_1=8.2\%$, w₂=41.2%$w_2=41.2\%$; (iii) w₁=22.3%$w_1=22.3\%$, w₂=27.7%$w_2=27.7\%$; and (iv) w₁=30.6%$w_1=30.6\%$, w₂=19.4%$w_2=19.4\%$, where w$w$ is the mass percent of the component. The present solubility data was correlated by a modified Kent–Eisenberg model. The model reasonably represents the present solubility data, not only over the considered conditions, but also for a wider range of temperatures, partial pressures, and compositions.     

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.     

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.     

Ab Initio Study of Fine and Hyperfine Interactions in Triplet POH

Phosphorous-containing molecules have a great relevance in prebiotic chemistry in view of the fact that phosphorous is a fundamental constituent of biomolecules, such as RNA, DNA, and ATP. Its biogenic importance has led astrochemists to investigate the possibility that P-bearing species could have formed in the interstellar medium (ISM) and subsequently been delivered to early Earth by rocky bodies. However, only two P-bearing molecules have been detected so far in the ISM, with the chemistry of interstellar phosphorous remaining poorly understood. Here, in order to shed further light on P-carriers in space, we report a theoretical spectroscopic characterisation of the rotational spectrum of POH in its 3A″ ground electronic state. State-of-the-art coupled-cluster schemes have been employed to derive rotational constants, centrifugal distortion terms, and most of the fine and hyperfine interaction parameters, while the electron spin–spin dipolar coupling has been investigated using the multi-configuration self-consistent-field method. The computed spectroscopic parameters have been used to simulate the appearance of triplet POH rotational and ro-vibrational spectra in different conditions, from cold to warm environments, either in gas-phase experiments or in molecular clouds. Finally, we point out that the predicted hyperfine structures represent a key pattern for the recognition of POH in laboratory and interstellar spectra.     

Interaction Energy Analysis of Monovalent Inorganic Anions in Bulk Water Versus Air/Water Interface

Soft anions exhibit surface activity at the air/water interface that can be probed using surface-sensitive vibrational spectroscopy, but the structural implications of this surface activity remain a matter of debate. Here, we examine the nature of anion–water interactions at the air/water interface using a combination of molecular dynamics simulations and quantum-mechanical energy decomposition analysis based on symmetry-adapted perturbation theory. Results are presented for a set of monovalent anions, including Cl⁻, Br⁻, I⁻, CN⁻, OCN⁻, SCN⁻, NO2−, NO3−, and ClOn− (n=1,2,3,4), several of which are archetypal examples of surface-active species. In all cases, we find that average anion–water interaction energies are systematically larger in bulk water although the difference (with respect to the same quantity computed in the interfacial environment) is well within the magnitude of the instantaneous fluctuations. Specifically for the surface-active species Br⁻(aq), I⁻(aq), ClO4−(aq), and SCN⁻(aq), and also for ClO⁻(aq), the charge-transfer (CT) energy is found to be larger at the interface than it is in bulk water, by an amount that is greater than the standard deviation of the fluctuations. The Cl⁻(aq) ion has a slightly larger CT energy at the interface, but NO3−(aq) does not; these two species are borderline cases where consensus is lacking regarding their surface activity. However, CT stabilization amounts to <20% of the total induction energy for each of the ions considered here, and CT-free polarization energies are systematically larger in bulk water in all cases. As such, the role of these effects in the surface activity of soft anions remains unclear. This analysis complements our recent work suggesting that the short-range solvation structure around these ions is scarcely different at the air/water interface from what it is in bulk water. Together, these observations suggest that changes in first-shell hydration structure around soft anions cannot explain observed surface activities.     

Dehydration performance of a hydrophobic DD3R zeolite membrane

The all silica DDR membrane turns out to be well suited to separate water from organic solvents under pervaporation conditions, despite its hydrophobic character. All-silica zeolites are chemically and hydrothermally more stable than aluminum containing ones and are therefore preferred for membrane applications, including for dehydration, even though these type of membranes are hydrophobic. Permeation of water, ethanol and methanol through an all-silica DDR membrane has been measured at temperatures ranging from 344 to 398 K. The hydrophobic membrane shows high water fluxes (up to 20 kg m⁻² h⁻¹). The pure water permeance is insensitive to temperature and is well described assuming weak adsorption. Excellent performance in dewatering ethanol (N=2$N=2$ kg m⁻² h⁻¹and _αw=1500${\alpha }_{\text{w}}=1500$ at 373 K and _xw=0.18$x_{\text{w}}=0.18$) is observed and the membrane is also able to selectively remove water from methanol (N=5$N=5$ kg m⁻² h⁻¹ and _αw=9${\alpha }_{\text{w}}=9$). Water could also be removed from methanol/ethanol/water (_αwater_/EtOH=1500${\alpha }_{\text{water}/\text{EtOH}}=1500$, _αMeOH_/EtOH=70${\alpha }_{\text{MeOH}/\text{EtOH}}=70$ at 373 K) mixtures, even at water feed concentrations below 1.5 mol%.     

Hard biaxial ellipsoids revisited: Numerical results

Monte Carlo simulations are undertaken for hard ellipsoids for a number of values of its semi-axes in the range c/a∈{0.1,10}$c/a\in \{0.1\text{,}10\}$. The isotropic phase results are compared to the Vega isotropic phase equation of state for hard biaxial ellipsoids [C. Vega, Virial coefficients and equation of state of hard ellipsoids, Mol. Phys. 92 (1997) 651]. The Vega equation of state is found to perform exceedingly well for all of the models studied. The biaxial phase is seen to form only after the previous formation of a discotic phase.     

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.     

Photochemistry of Thymine in Protic Polar Nanomeric Droplets Using Electrostatic Embeding TD-DFT/MM

Thymine photochemistry is important for understanding DNA photodamage. In the gas phase, thymine undergoes a fast non-radiative decay from S 2 to S 1 . In the S 1 state, it gets trapped for several picoseconds until returning to the ground-state S 0 . Here, we explore the electrostatic effects of nanomeric droplets of methanol and water on the excited states of thymine. For this purpose, we develop and implement an electrostatic embedding TD-DFT/MM method based on a QM/MM coupling defined through electrostatic potential fitting charges. We show that both in methanol and water, the mechanism is similar to the gas phase. The solvent molecules participate in defining the branching plane of S 0 /S 1 intersection and have a negligible effect on the S 1 /S 2 intersection. Despite the wrong topology of the ground/excited state intersections, electrostatic embedding TD-DFT/MM allows for a fast exploration of the potential energy surfaces and a qualitative picture of the photophysics of thymine in solvent droplets.