TiCl4 Dissolved in Ionic Liquid Mixtures from Аb Initio Molecular Dynamics Simulations

To gain a deeper understanding of the TiCl4 solvation effects in multi-component ionic liquids, we performed ab initio molecular dynamics simulations of 1-butyl-3-methylimidazolium [C4C1Im]+, tetrafluoroborate [BF4]−, chloride [Cl]− both with and without water and titanium tetrachloride TiCl4. Complex interactions between cations and anions are observed in all investigated systems. By further addition of water and TiCl4 this complex interaction network is extended. Observations of the radial distribution functions and number integrals show that water and TiCl4 not only compete with each other to interact mainly with [Cl]−, which strongly influences the cation-[BF4]− interaction, but also interact with each other, which leads to the fact that in certain systems the cation-anion interaction is enhanced. Further investigations of the Voronoi polyhedra analysis have demonstrated that water has a greater impact on the nanosegregated system than TiCl4 which is also due to the fact of the shear amount of water relative to all other components and its higher mobility compared to TiCl4. Overall, the polar network of the IL mixture collapses by including water and TiCl4. In the case of [Cl]− chloride enters the water continuum, while [BF4]− remains largely unaffected, which deeply affects the interaction of the ionic liquid (IL) network.     

Magnetic-Field-Orientation Dependent Thermal Entanglement of a Spin-1 Heisenberg Dimer: The Case Study of Dinuclear Nickel Complex with an Uniaxial Single-Ion Anisotropy

The bipartite entanglement in pure and mixed states of a quantum spin-1 Heisenberg dimer with exchange and uniaxial single-ion anisotropies is quantified through the negativity in a presence of the external magnetic field. At zero temperature the negativity shows a marked stepwise dependence on a magnetic field with two abrupt jumps and plateaus, which can be attributed to the quantum antiferromagnetic and quantum ferrimagnetic ground states. The magnetic-field-driven phase transition between the quantum antiferromagnetic and quantum ferrimagnetic ground states manifests itself at nonzero temperatures by a local minimum of the negativity, which is followed by a peculiar field-induced rise of the negativity observable in a range of moderately strong magnetic fields. The rising temperature generally smears out abrupt jumps and plateaus of the negativity, which cannot be distinguished in the relevant dependencies above a certain temperature. It is shown that the thermal entanglement is most persistent against rising temperature at the magnetic field, for which an energy gap between a ground state and a first excited state is highest. Besides, temperature variations of the negativity of the spin-1 Heisenberg dimer with an easy-axis single-ion anisotropy may exhibit a singular point-kink, at which the negativity has discontinuity in its first derivative. The homodinuclear nickel complex [Ni2(Medpt)2(μ-ox)(H2O)2](ClO4)2·2H2O provides a suitable experimental platform of the antiferromagnetic spin-1 Heisenberg dimer, which allowed us to estimate a strength of the bipartite entanglement between two exchange-coupled Ni2+ magnetic ions on the grounds of the interaction constants reported previously from the fitting procedure of the magnetization data. It is verified that the negativity of this dinuclear compound is highly magnetic-field-orientation dependent due to presence of a relatively strong uniaxial single-ion anisotropy.     

Ca+ Ions Solvated in Helium Clusters

We present a combined experimental and theoretical investigation on Ca+ ions in helium droplets, HeNCa+. The clusters have been formed in the laboratory by means of electron-impact ionization of Ca-doped helium nanodroplets. Energies and structures of such complexes have been computed using various approaches such as path integral Monte Carlo, diffusion Monte Carlo and basin-hopping methods. The potential energy functions employed in these calculations consist of analytical expressions following an improved Lennard-Jones formula whose parameters are fine-tuned by exploiting ab initio estimations. Ion yields of HeNCa+ -obtained via high-resolution mass spectrometry- generally decrease with N with a more pronounced drop between N=17 and N=25, the computed quantum HeNCa+ evaporation energies resembling this behavior. The analysis of the energies and structures reveals that covering Ca+ with 17 He atoms leads to a cluster with one of the smallest energies per atom. As new atoms are added, they continue to fill the first shell at the expense of reducing its stability, until N=25, which corresponds to the maximum number of atoms in that shell. Behavior of the evaporation energies and radial densities suggests liquid-like cluster structures.     

Kinetic Study on Alpha-Form Crystallization of Mixed-Acid Triacylglycerols POP,PPO, and Their Mixture

The crystallization behavior of the metastable α form of triacylglycerols (TAGs) plays a critical role as a precursor for the crystallization of more stable β′ and β forms for various applications in food and pharmaceutical products. However, precise analysis of the crystallization kinetics of α has not been performed, likely due to its rapid and complex behavior. This paper presents the observation results of the initial stages of the isothermal crystallization kinetics of α forms of 1,3-dipalmitoyl-2-oleoyl-glycerol (POP), 1,2-dipalmitoyl-3-oleoyl-rac-glycerol (rac-PPO), and molecular compound (MC) crystals of a POP/rac-PPO (1/1) mixture (MCPOP/PPO) using synchrotron radiation time-resolved X-ray diffraction and polarized optical microscopy. In all the TAGs, α crystals with a worm-like morphology started to grow rapidly in the first stage. Then, the α crystals slowly transformed into more stable forms in different manners for different TAG samples. In POP, the conversion was simple, as the α-2 form transformed into γ-3, whereas in rac-PPO, the lamellar distance values of the α-2 form continuously decreased with time and changed into the α-3 form. In the MCPOP/PPO crystals, in contrast, separate crystallization of α-2 of a rac-PPO fraction initially occurred, followed by the crystallization of α-2 of POP, and the two α forms merged into α-2 of MCPOP/PPO. This separate crystallization was caused by large differences in the crystallization kinetics of the α forms of POP and rac-PPO.     

Free Energies of Hydrated Halide Anions: High Through-Put Computations on Clusters to Treat Rough Energy-Landscapes

With a longer-term goal of addressing the comparative behavior of the aqueous halides F−, Cl−, Br−, and I− on the basis of quasi-chemical theory (QCT), here we study structures and free energies of hydration clusters for those anions. We confirm that energetically optimal (H2O)nX clusters, with X = Cl−, Br−, and I−, exhibit surface hydration structures. Computed free energies, based on optimized surface hydration structures utilizing a harmonic approximation, typically (but not always) disagree with experimental free energies. To remedy the harmonic approximation, we utilize single-point electronic structure calculations on cluster geometries sampled from an AIMD (ab initio molecular dynamics) simulation stream. This rough-landscape procedure is broadly satisfactory and suggests unfavorable ligand crowding as the physical effect addressed. Nevertheless, this procedure can break down when n≳4, with the characteristic discrepancy resulting from a relaxed definition of clustering in the identification of (H2O)nX clusters, including ramified structures natural in physical cluster theories. With ramified structures, the central equation for the present rough-landscape approach can acquire some inconsistency. Extension of these physical cluster theories in the direction of QCT should remedy that issue, and should be the next step in this research direction.     

A Computational Study of the S2 State in the Oxygen-Evolving Complex of Photosystem II by Electron Paramagnetic Resonance Spectroscopy

The S₂ state produces two basic electron paramagnetic resonance signal types due to the manganese cluster in oxygen-evolving complex, which are influenced by the solvents, and cryoprotectant added to the photosystem II samples. It is presumed that a single manganese center oxidation occurs on S₁ → S₂ state transition. The S₂ state has readily visible multiline and g4.1 electron paramagnetic resonance signals and hence it has been the most studied of all the Kok cycle intermediates due to the ease of experimental preparation and stability. The S₂ state was studied using electron paramagnetic resonance spectroscopy at X-band frequencies. The aim of this study was to determine the spin states of the g4.1 signal. The multiline signal was observed to arise from a ground state spin ½ centre while the g4.1 signal generated at ≈140 K NIR illumination was proposed to arise from a spin 52 center with rhombic distortion. The ‘ground’ state g4.1 signal was generated solely or by conversion from the multiline. The data analysis methods used involved numerical simulations of the experimental spectra on relevant models of the oxygen-evolving complex cluster. A strong focus in this paper was on the ‘ground’ state g4.1 signal, whether it is a rhombic 52 spin state signal or an axial 32 spin state signal. The data supported an X-band CW-EPR-generated g4.1 signal as originating from a near rhombic spin 5/2 of the S₂ state of the PSII manganese cluster.     

Carbon Atoms Speaking Out: How the Geometric Sensitivity of 13C Chemical Shifts Leads to Understanding the Colour Tuning of Phycocyanobilin in Cph1 and AnPixJ

We present a combined quantum mechanics/molecular mechanics (QM/MM) molecular dynamics–statistical approach for the interpretation of nuclear magnetic resonance (NMR) chemical shift patterns in phycocyanobilin (PCB). These were originally associated with colour tuning upon photoproduct formation in red/green-absorbing cyanobacteriochrome AnPixJg2 and red/far-red-absorbing phytochrome Cph1Δ2. We pursue an indirect approach without computation of the absorption frequencies since the molecular geometry of cofactor and protein are not accurately known. Instead, we resort to a heuristic determination of the conjugation length in PCB through the experimental NMR chemical shift patterns, supported by quantum chemical calculations. We have found a characteristic correlation pattern of 13C chemical shifts to specific bond orders within the π-conjugated system, which rests on the relative position of carbon atoms with respect to electron-withdrawing groups and the polarisation of covalent bonds. We propose the inversion of this regioselective relationship using multivariate statistics and to apply it to the known experimental NMR chemical shifts in order to predict changes in the bond alternation pattern. Therefrom the extent of electronic conjugation, and eventually the change in absorption frequency, can be derived. In the process, the consultation of explicit mesomeric formulae plays an important role to qualitatively account for possible conjugation scenarios of the chromophore. While we are able to consistently associate the NMR chemical shifts with hypsochromic and bathochromic shifts in the Pg and Pfr, our approach represents an alternative method to increase the explanatory power of NMR spectroscopic data in proteins.     

An Exchange Mechanism for the Magnetic Behavior of Er3+ Complexes

We study the magnetic properties of the erbium based compounds, Na9[Er(W5O18)2] and [(Pc)Er{Pc{N(C4H9)2}8}]·/−, in the framework of an effective spin exchange model involving delocalized electrons occupying molecular orbitals. The calculations successfully reproduce the experimental data available in the literature for the magnetic spectrum, magnetization and molar susceptibility in dc and ac fields. Owing to their similar molecular geometry, the compounds’ magnetic behaviors are interpreted in terms of the same set of active orbitals and thus the same effective spin coupling scheme. For all three complexes, the model predicts a prompt change in the ground state from a Kramer’s doublet at zero fields to a fully polarized quartet one brought about by the action of an external magnetic field without Zeeman splitting. This alteration is attributed to the enhancement of the effect of orbital interactions over the spin exchange as the magnitude of the external magnetic field increases.     

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.     

Operando Raman Shift Replaces Current in Electrochemical Analysis of Li-ion Batteries: A Comparative Study

Li-rich and catalytically active γ-LixV2O5 (x = 1.48) was investigated as a cathode for its heterogeneous charge transfer kinetics. Using a specially designed two-electrode system lithium half cell, Butler–Volmer analysis was performed, and Raman spectra were acquired in 18 mV intervals. A direct correlation was observed between the Raman shift of the active modes Ag,Bg, Au, and Bu, and the development of the Faraday current at the working electrode. The Raman intensity and the Raman shift were implemented to replace the current in a Tafel plot used for the analysis of Butler–Volmer kinetics. Striking similarities in the charge transfer proportionality constants α were found for current and Raman-based analysis. The potential of this new method of Raman-aided electrochemical detection at the diffraction limit is discussed.     

Zhang–Zhang Polynomials of Multiple Zigzag Chains Revisited: A Connection with the John–Sachs Theorem

Multiple zigzag chains Zm,n of length n and width m constitute an important class of regular graphene flakes of rectangular shape. The physical and chemical properties of these basic pericondensed benzenoids can be related to their various topological invariants, conveniently encoded as the coefficients of a combinatorial polynomial, usually referred to as the ZZ polynomial of multiple zigzag chains Zm,n. The current study reports a novel method for determination of these ZZ polynomials based on a hypothesized extension to John–Sachs theorem, used previously to enumerate Kekulé structures of various benzenoid hydrocarbons. We show that the ZZ polynomial of the Zm,n multiple zigzag chain can be conveniently expressed as a determinant of a Toeplitz (or almost Toeplitz) matrix of size m2×m2 consisting of simple hypergeometric polynomials. The presented analysis can be extended to generalized multiple zigzag chains Zkm,n, i.e., derivatives of Zm,n with a single attached polyacene chain of length k. All presented formulas are accompanied by formal proofs. The developed theoretical machinery is applied for predicting aromaticity distribution patterns in large and infinite multiple zigzag chains Zm,n and for computing the distribution of spin densities in biradical states of finite multiple zigzag chains Zm,n.     

Fractional Dynamics of HIV with Source Term for the Supply of New CD4+ T-Cells Depending on the Viral Load via Caputo–Fabrizio Derivative

Human immunodeficiency virus (HIV) is a life life-threatening and serious infection caused by a virus that attacks CD4+ T-cells, which fight against infections and make a person susceptible to other diseases. It is a global public health problem with no cure; therefore, it is highly important to study and understand the intricate phenomena of HIV. In this article, we focus on the numerical study of the path-tracking damped oscillatory behavior of a model for the HIV infection of CD4+ T-cells. We formulate fractional dynamics of HIV with a source term for the supply of new CD4⁺ T-cells depending on the viral load via the Caputo–Fabrizio derivative. In the formulation of fractional HIV dynamics, we replaced the constant source term for the supply of new CD4+ T-cells from the thymus with a variable source term depending on the concentration of the viral load, and introduced a term that describes the incidence of the HIV infection of CD4+ T-cells. We present a novel numerical scheme for fractional view analysis of the proposed model to highlight the solution pathway of HIV. We inspect the periodic and chaotic behavior of HIV for the given values of input factors using numerical simulations.     

Cross-Talk between Overlap Interactions in Biomolecules: A Case Study of the β-Turn Motif

Noncovalent interactions play a pivotal role in regulating protein conformation, stability and dynamics. Among the quantum mechanical (QM) overlap-based noncovalent interactions, n→π* is the best understood with studies ranging from small molecules to β-turns of model proteins such as GB1. However, these investigations do not explore the interplay between multiple overlap interactions in contributing to local structure and stability. In this work, we identify and characterize all noncovalent overlap interactions in the β-turn, an important secondary structural element that facilitates the folding of a polypeptide chain. Invoking a QM framework of natural bond orbitals, we demonstrate the role of several additional interactions such as n→σ* and π→π* that are energetically comparable to or larger than n→π*. We find that these interactions are sensitive to changes in the side chain of the residues in the β-turn of GB1, suggesting that the n→π* may not be the only component in dictating β-turn conformation and stability. Furthermore, a database search of n→σ* and π→π* in the PDB reveals that they are prevalent in most proteins and have significant interaction energies (∼1 kcal/mol). This indicates that all overlap interactions must be taken into account to obtain a comprehensive picture of their contributions to protein structure and energetics. Lastly, based on the extent of QM overlaps and interaction energies, we propose geometric criteria using which these additional interactions can be efficiently tracked in broad database searches.     

Significance of Shape Factor in Heat Transfer Performance of Molybdenum-Disulfide Nanofluid in Multiple Flow Situations; A Comparative Fractional Study

In this modern era, nanofluids are considered one of the advanced kinds of heat transferring fluids due to their enhanced thermal features. The present study is conducted to investigate that how the suspension of molybdenum-disulfide (MoS2) nanoparticles boosts the thermal performance of a Casson-type fluid. Sodium alginate (NaAlg) based nanofluid is contained inside a vertical channel of width d and it exhibits a flow due to the movement of the left wall. The walls are nested in a permeable medium, and a uniform magnetic field and radiation flux are also involved in determining flow patterns and thermal behavior of the nanofluid. Depending on velocity boundary conditions, the flow phenomenon is examined for three different situations. To evaluate the influence of shape factor, MoS2 nanoparticles of blade, cylinder, platelet, and brick shapes are considered. The mathematical modeling is performed in the form of non-integer order operators, and a double fractional analysis is carried out by separately solving Caputo-Fabrizio and Atangana-Baleanu operators based fractional models. The system of coupled PDEs is converted to ODEs by operating the Laplace transformation, and Zakian’s algorithm is applied to approximate the Laplace inversion numerically. The solutions of flow and energy equations are presented in terms of graphical illustrations and tables to discuss important physical aspects of the observed problem. Moreover, a detailed inspection on shear stress and Nusselt number is carried out to get a deep insight into skin friction and heat transfer mechanisms. It is analyzed that the suspension of MoS2 nanoparticles leads to ameliorating the heat transfer rate up to 9.5%. To serve the purpose of achieving maximum heat transfer rate and reduced skin friction, the Atangana-Baleanu operator based fractional model is more effective. Furthermore, it is perceived that velocity and energy functions of the nanofluid exhibit significant variations because of the different shapes of nanoparticles.     

Photophysical Deactivation Mechanisms of the Pyrimidine Analogue 1-Cyclohexyluracil

The photophysical relaxation mechanisms of 1-cyclohexyluracil, in vacuum and water, were investigated by employing the Multi-State CASPT2 (MS-CASPT2, Multi-State Complete Active-Space Second-Order Perturbation Theory) quantum chemical method and Dunning’s cc-pVDZ basis sets. In both environments, our results suggest that the primary photophysical event is the population of the S11(ππ*) bright state. Afterwards, two likely deactivation pathways can take place, which is sustained by linear interpolation in internal coordinates defined via Z-Matrix scans connecting the most important characteristic points. The first one (Route 1) is the same relaxation mechanism observed for uracil, its canonical analogue, i.e., internal conversion to the ground state through an ethylenic-like conical intersection. The other route (Route 2) is the direct population transfer from the S11(ππ*) bright state to the T23(nπ*) triplet state via an intersystem crossing process involving the (S11(ππ*)/T23(nπ*))STCP singlet-triplet crossing point. As the spin-orbit coupling is not too large in either environment, we propose that most of the electronic population initially on the S11(ππ*) state returns to the ground following the same ultrafast deactivation mechanism observed in uracil (Route 1), while a smaller percentage goes to the triplet manifold. The presence of a minimum on the S11(ππ*) potential energy hypersurface in water can help to understand why experimentally it is noticed suppression of the triplet states population in polar protic solvent.     

Hydrogen-Bonded and Halogen-Bonded: Orthogonal Interactions for the Chloride Anion of a Pyrazolium Salt

In the hydrochloride of a pyrazolyl-substituted acetylacetone, the chloride anion is hydrogen-bonded to the protonated pyrazolyl moiety. Equimolar co-crystallization with tetrafluorodiiodobenzene (TFDIB) leads to a supramolecular aggregate in which TFDIB is situated on a crystallographic center of inversion. The iodine atom in the asymmetric unit acts as halogen bond donor, and the chloride acceptor approaches the σ-hole of this TFDIB iodine subtending an almost linear halogen bond, with Cl···I = 3.1653(11) Å and Cl···I–C = 179.32(6)°. This contact is roughly orthogonal to the N–H···Cl hydrogen bond. An analysis of the electron density according to Bader’s Quantum Theory of Atoms in Molecules confirms bond critical points (bcps) for both short contacts, with ρbcp = 0.129 for the halogen and 0.321 eÅ−3 for the hydrogen bond. Our halogen-bonded adduct represents the prototype for a future class of co-crystals with tunable electron density distribution about the σ-hole contact.     

Catechin Photolysis Suppression by Aluminum Chloride under Alkaline Conditions and Assessment with Liquid Chromatography–Mass Spectrometry

Tea is rich in catechins and aluminum. In this study, the process of catechin photolysis was applied as a model for examining the effects of aluminum chloride (AlCl₃) on the structural changes of catechin and the alteration of aluminum complexes under blue light irradiation (BLI) at pH 8 using liquid chromatography and mass spectrometry techniques. Additionally, the effects of anions on catechin upon the addition of AlCl₃ and treatment with BLI were also studied. In this study, when 1 mM catechin was treated with BLI, a superoxide anion radical (O2•−) was generated in an air-saturated aqueous solution, in addition to forming a dimeric catechin (proanthocyanidin) via a photon-induced redox reaction. The relative percentage of catechin was found to be 59.0 and 95.7 for catechin treated with BLI and catechin upon the addition of 1 mM AlCl₃ treated with BLI, respectively. It suggested that catechin treated with BLI could be suppressed by AlCl₃, while AlCl₃ did not form a complex with catechin in the photolytic system. However, under the same conditions, it was also found that the addition of AlCl₃ inhibited the photolytic formation of O2•−, and reduced the generation of proanthocyanidin, suggesting that the disconnection of proanthocyanidin was achieved by AlCl₃ acting as a catalyst under treatment with BLI. The influence of 1 mM fluoride (F−) and 1 mM oxalate (C2O42−) ions on the photolysis of 1 mM catechin upon the addition of 1 mM AlCl₃ and treatment with BLI was found to be insignificant, implying that, during the photolysis of catechin, the Al species were either neutral or negatively charged and the aluminum species did not form a complex with anions in the photolytic system. Therefore, aluminum, which is an amphoteric species, has an inherent potential to stabilize the photolysis of catechin in an alkaline conditions, while suppressing the O2•− and proanthocyanidin generation via aluminum ion catalysis in the catechin/Al system under treatment with BLI.     

Understanding Topological Insulators in Real Space

A real space understanding of the Su–Schrieffer–Heeger model of polyacetylene is introduced thanks to delocalization indices defined within the quantum theory of atoms in molecules. This approach enables to go beyond the analysis of electron localization usually enabled by topological insulator indices—such as IPR—enabling to differentiate between trivial and topological insulator phases. The approach is based on analyzing the electron delocalization between second neighbors, thus highlighting the relevance of the sublattices induced by chiral symmetry. Moreover, the second neighbor delocalization index, δi,i+2, also enables to identify the presence of chirality and when it is broken by doping or by eliminating atom pairs (as in the case of odd number of atoms chains). Hints to identify bulk behavior thanks to δ1,3 are also provided. Overall, we present a very simple, orbital invariant visualization tool that should help the analysis of chirality (independently of the crystallinity of the system) as well as spreading the concepts of topological behavior thanks to its relationship with well-known chemical concepts.     

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