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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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

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.     

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.     

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.     

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.     

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.     

Unveiling Carbon Dioxide and Ethanol Diffusion in Carbonated Water-Ethanol Mixtures by Molecular Dynamics Simulations

The diffusion of carbon dioxide (CO2) and ethanol (EtOH) is a fundamental transport process behind the formation and growth of CO2 bubbles in sparkling beverages and the release of organoleptic compounds at the liquid free surface. In the present study, CO2 and EtOH diffusion coefficients are computed from molecular dynamics (MD) simulations and compared with experimental values derived from the Stokes-Einstein (SE) relation on the basis of viscometry experiments and hydrodynamic radii deduced from former nuclear magnetic resonance (NMR) measurements. These diffusion coefficients steadily increase with temperature and decrease as the concentration of ethanol rises. The agreement between theory and experiment is suitable for CO2. Theoretical EtOH diffusion coefficients tend to overestimate slightly experimental values, although the agreement can be improved by changing the hydrodynamic radius used to evaluate experimental diffusion coefficients. This apparent disagreement should not rely on limitations of the MD simulations nor on the approximations made to evaluate theoretical diffusion coefficients. Improvement of the molecular models, as well as additional NMR measurements on sparkling beverages at several temperatures and ethanol concentrations, would help solve this issue.     

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.