Towards a transferable nonelectrostatic model for continuum solvation: The electrostatic and nonelectrostatic energy correction model

In this work, we introduce an electrostatic and non-electrostatic (ENE) correction to the solvation energy based on the Solvent-Accessible Surface Area (SASA) of the solute and the solvent static dielectric constant. The proposed correction was developed for neutral solutes in non-aqueous solvents, considering three different implicit solvation models based on a Self-Consistent Reaction Field treatment of solute-solvent mutual polarization using an Apparent Surface Charge formalism, namely the Integral Equation Formalism of the Polarizable Continuum Model using a continuous surface charge scheme (PCM), the Solvation Model based on solute electron density (SMD), and the generalized Finite-Difference Poisson-Boltzmann (FDPB) model. The proposed correction was parametrized on a diverse training set of 4980 solvation data from the Solv@tum database of experimental solvation energies, and validated on the non-aqueous subset of the MNSOL database comprising 2140 solvation energies. The performances of the proposed ENE models with minimal and extended parameters formulations have been analyzed and the latter variant has been further compared to the widely used Cavity, Dispersion, and Solvent structural effects (CDS) non-electrostatic model originally developed for the SMx family of implicit solvation models. Overall, a very good agreement between the computed solvation energies with the ENE correction and the reference experimental data has been found on both the training and test sets for all continuum solvation models considered. Furthermore, results for the ENE correction are on par with the reference CDS non-electrostatic model for both SMD and FDPB electrostatics, but with the advantage of using a lower number of parameters and thus an improved transferability between different electrostatics treatments.

Self-trapped excitons in diamond: a Δ-SCF approach

This paper reports the first variationally-based predictions of the lowest excited state in diamond (Γ25՛ → Γ15) in the optical and thermal configurations, from direct Δ-SCF calculations based on B3LYP, PBE0, HSE06 and GGA functionals. For the B3LYP functional the energy of the optical state, 7.27 eV, is within the observed range of (7.2-7.4) eV and is predicted to be insulating, with indirect band gaps of (5.6-5.8) eV. Mulliken analyses of the excited state wavefunction indicate extensive redistributions of charge and spin resulting in a strongly excitonic state with a central charge of -0.8 ǀeǀ surrounded by charges of +0.12 ǀeǀ at the four nearest neighbour sites. Calculations of the ground and excited state relaxations lead to a Stokes' shift of 0.47 eV and predicted Γ-point luminescence energy of 6.89 eV. Assuming a similar shift at the band edge (X1), an estimate of 5.29 eV is predicted for the luminescence energy, which compares with the observed value of 5.27 eV. Excited state vibrational spectra show marked differences from the ground state, with the introduction of an infra-red peak at 1150 cm-1 and a modest shift of 2 cm -1 in the TO(X) Raman mode at 1340 cm-1. Similar calculations of the lowest energy bi- and tri-excitons predict these to be bound states in both optical and thermal configurations and plausible precursors to exciton condensation. Estimates of bi- and triexciton luminescence energies predict red shifts with respect to the single exciton line which are compared to the recently reported values.

Modeling the spectral properties of poly(x-phenylenediamine) conducting polymers using a combined TD-DFT and electrostatic embedding approach

The absorption spectra of polymers derived from ortho, meta and para phenylenediamines (o-PDA, m-PDA and p-PDA) have been simulated combining periodic density functional theory (DFT) calculations with time-dependent DFT simulations. These latter have been carried out on finite clusters embedded in a set of point charges devised to exactly reproduce the electrostatic potential of the periodic chains. The results are compared with those obtained for solvated o-PDA, m-PDA and p-PDA oligomers of increasing sizes extracted from the periodic structures. The electronic transitions involved have been investigated by a qualitative analysis based on isodensity maps completed by a quantitative analysis based on the density-based index (D_CT ). For poly-(o)- and poly-(p)- phenylenediamines the agreement with the experimental data is achieved already by modeling solvated dimers whereas the inclusion of long-range electrostatic effects is mandatory for poly-(m)-phenylenediamine highlighting the importance of an accurate treatment of the electrostatic environment when a finite cluster approach is considered.

Shry: Application of Canonical Augmentation to the Atomic Substitution Problem

A common approach for studying a solid solution or disordered system within a periodic ab initio framework is to create a supercell in which certain amounts of target elements are substituted with other elements. The key to generating supercells is determining how to eliminate symmetry-equivalent structures from many substitution patterns. Although the total number of substitutions is on the order of trillions, only symmetry-inequivalent atomic substitution patterns need to be identified, and their number is far smaller than the total. Our developed Python software package, which is called Shry (Suite for High-throughput generation of models with atomic substitutions implemented by Python), allows the selection of only symmetry-inequivalent structures from the vast number of candidates based on the canonical augmentation algorithm. Shry is implemented in Python 3 and uses the CIF format as the standard for both reading and writing the reference and generated sets of substituted structures. Shry can be integrated into another Python program as a module or can be used as a stand-alone program. The implementation was verified through a comparison with other codes with the same functionality, based on the total numbers of symmetry-inequivalent structures, and also on the equivalencies of the output structures themselves. The provided crystal structure data used for the verification are expected to be useful for benchmarking other codes and also developing new algorithms in the future.

Fast Quantum Approach for Evaluating the Energy of Non-Covalent Interactions in Molecular Crystals: The Case Study of Intermolecular H-Bonds in Crystalline Peroxosolvates

Energy/enthalpy of intermolecular hydrogen bonds (H-bonds) in crystals have been calculated in many papers. Most of the theoretical works used non-periodic models. Their applicability for describing intermolecular H-bonds in solids is not obvious since the crystal environment can strongly change H-bond geometry and energy in comparison with non-periodic models. Periodic DFT computations provide a reasonable description of a number of relevant properties of molecular crystals. However, these methods are quite cumbersome and time-consuming compared to non-periodic calculations. Here, we present a fast quantum approach for estimating the energy/enthalpy of intermolecular H-bonds in crystals. It has been tested on a family of crystalline peroxosolvates in which the H∙∙∙O bond set fills evenly (i.e., without significant gaps) the range of H∙∙∙O distances from ~1.5 to ~2.1 Å typical for strong, moderate, and weak H-bonds. Four of these two-component crystals (peroxosolvates of macrocyclic ethers and creatine) were obtained and structurally characterized for the first time. A critical comparison of the approaches for estimating the energy of intermolecular H-bonds in organic crystals is carried out, and various sources of errors are clarified.

Anharmonic Coupling of Stretching Vibrations in Ice: A Periodic VSCF and VCI Description

The anharmonicity of O-H stretching vibrations of water ice is characterized by use of a periodic implementation of the vibrational self-consistent field (VSCF) and vibrational configuration interaction (VCI) methods, which take phonon-phonon couplings explicitly into account through numerical evaluation of high-order terms of the nuclear potential. The low-temperature, proton-ordered phase of water ice (namely, ice XI) is investigated. The net effect of a coupled anharmonic treatment of stretching modes is not just a rigid blue-shift of the respective harmonic spectral frequencies but rather a complex change of their relative spectral positions, which cannot be captured by simple scaling strategies based on harmonic calculations. The adopted techniques allow for a hierarchical treatment of anharmonic terms of the nuclear potential, which is key to an effective identification of leading factors. We show that the anharmonic independent-mode approximation─only describing the "intrinsic anharmonicity" of the O-H stretches─is unable to capture the correct physics, and that couplings among O-H stretches must be described. Inspection of harmonic normal coordinates allows identification of specific features of the O-H stretching motions which most likely enable strong mode-mode couplings. Finally, by coupling O-H stretches to all other possible modes of ice XI (THz collective vibrations, molecular librations, bendings), we identify specific types of motion which significantly affect O-H stretching states: in particular, molecular librations are found to affect the stretching states more than molecular bendings.

Water Splitting on Multifaceted SrTiO3 Nanocrystals: Calculations of Raman Vibrational Spectrum

Various photocatalysts are being currently studied with the aim of increasing the photocatalytic efficiency of water splitting for production of hydrogen as a fuel and oxygen as a medical gas. A noticeable increase of hydrogen production was found recently experimentally on the anisotropic faces (facets) of strontium titanate (SrTiO₃, STO) nanoparticles. In order to identify optimal sites for water splitting, the first principles calculations of the Raman vibrational spectrum of the bulk and stepped (facet) surface of a thin STO film with adsorbed water derivatives were performed. According to our calculations, the Raman spectrum of a stepped STO surface differs from the bulk spectrum, which agrees with the experimental data. The characteristic vibrational frequencies for the chemisorption of water derivatives on the surface were identified. Moreover, it is also possible to distinguish between differently adsorbed hydrogen atoms of a split water molecule. Our approach helps to select the most efficient (size and shape) perovskite nanoparticles for efficient hydrogen/oxygen photocatalytic production.

Spin splitting in monoperiodic systems described by magnetic line groups

In this paper we report the classification of all the 81 magnetic line group families into seven spin splitting prototypes, in analogy to the similar classification previously reported for the 1651 magnetic space groups, 528 magnetic layer groups, and 394 magnetic rod groups. According to this classification, electrically induced (Pekar-Rashba) spin splitting is possible in the antiferromagnetic structures described by magnetic line groups of type I (no anti-unitary operations) and III, both in the presence and in the absence of the space inversion operation. As a specific example, a group theoretical analysis of spin splitting in CoO (8, 8) nanotube is carried out and its predictions are confirmed byab initiodensity functional theory calculations.

The structural origin of the efficient photochromism in natural minerals

SignificanceNatural photochromic minerals have been reported by geologists for decades. However, the understanding of the photochromism mechanism has a key question still unanswered: What in their structure gives rise to the photochromism's reversibility? By combining experimental and computational methods specifically developed to investigate this photochromism, this work provides the answer to this fundamental question. The specific crystal structure of these minerals allows an unusual motion of the sodium atoms stabilizing the electronic states associated to the colored forms. With a complete understanding of the photochromism mechanism in hand, it is now possible to design new families of stable and tunable photochromic inorganic materials-based devices.

Effect of Doping on Rutile TiO2 Surface Stability and Crystal Shapes

Transition-metal-(TM-)doped TiO₂ has been considered as promising electrode material for the oxygen evolution reaction (OER). OER activity is expected to depend on the coordination of the surface atoms. In this study, we theoretically investigate the stability of low-index surfaces of TM-doped rutile, (110), (100), (101) and (001), with 50 % of the Ti atoms substituted by Sc, Y, V, Nb or Ta. For Sc and Y, we also consider models with O vacancies providing the most stable oxidation state of Sc and Y. Surface energies are calculated with DFT(+U). Based on the Gibbs-Wulff theorem, the shape of the single crystals is predicted. It is observed that p-doping leads to spontaneous oxygen loss and O vacancies cause surface reconstruction. The Wulff shapes of n-doped TiO₂ have smaller contributions of the (110) facet and, for Nb and Ta, larger contributions of other facets. Given the higher coordinative unsaturation of the TM atoms in the latter, a higher catalytic activity is expected.

Disordered Rock-Salt Type Li2TiS3 as Novel Cathode for LIBs: A Computational Point of View

The development of high-energy cathode materials for lithium-ion batteries with low content of critical raw materials, such as cobalt and nickel, plays a key role in the progress of lithium-ion batteries technology. In recent works, a novel and promising family of lithium-rich sulfides has received attention. Among the possible structures and arrangement, cubic disordered Li2TiS3 has shown interesting properties, also for the formulation of new cell for all-solid-state batteries. In this work, a computational approach based on DFT hybrid Hamiltonian, localized basis functions and the use of the periodic CRYSTAL code, has been set up. The main goal of the present study is to determine accurate structural, electronic, and spectroscopic properties for this class of materials. Li2TiS3 precursors as Li2S, TiS2, and TiS3 alongside other formulations and structures such as LiTiS2 and monoclinic Li2TiS3 have been selected as benchmark systems and used to build up a consistent and robust predictive scheme. Raman spectra, XRD patterns, electronic band structures, and density of states have been simulated and compared to available literature data. Disordered rock-salt type Li2TiS3 structures have been derived via a solid solution method as implemented into the CRYSTAL code. Representative structures were extensively characterized through the calculations of their electronic and vibrational properties. Furthermore, the correlation between structure and Raman fingerprint was established.

Real Temperature Model of Dynamic Disorder in Molecular Crystals

Charge carrier mobilities in ordered organic semiconductors are limited by inherent vibrational phonons that scatter carriers. To improve a material's intrinsic mobility, restricting particularly detrimental modes with molecular substitutions may be a viable strategy. Here, we develop a probabilistic temperature-dependent displacement model that we couple with the density functional dimer projection protocol to predict effective electronic coupling fluctuations. The phonon-induced deviations from the equilibrium electronic couplings are used to infer the detriment of low-frequency phonons on charge carrier mobilities in a set of organic single crystals. We show that asymmetric sliding motions in pentacene and 2,6-diphenylanthracene induce large electronic coupling fluctuations, whereas seesawlike motions cause large fluctuations in rubrene, 9,10-diphenylanthracene, and, 2,6-diphenylanthracene. Vibrational analyses revealed that the asymmetric sliding phonon in rubrene persists only in the low-mobility direction of the crystal. Therefore, rubrene's intrinsic high mobility is likely due to the absence of this source of disorder in its high-mobility conduction channels. This model can be used to identify particularly harmful or helpful phonons in crystalline materials and may provide design rules for developing materials with intrinsically low disorder.

Unraveling the MnMoO4 polymorphism: a comprehensive DFT investigation of α, β, and ω phases

The MnMoO₄ is an environmentally friendly semiconductor material widely employed in technological devices. This material can be obtained on three different polymorphs, and although such phases were reported decades ago, some obscurity over their structure and properties is still perceived. Thus, this work provides a comprehensive DFT investigation of the α, β, and ω phases of MnMoO_4, analyzing their crystalline structure, stability, and electronic and magnetic properties. The results show that all phases of MnMoO₄ are stable at room conditions connected by pressure application or long-time high-temperature treatment. The MnMoO₄ phases are G-type antiferromagnetic with semiconductor bandgap and have enormous potential to develop magnetic, optical, and electronic devices and photocatalytic-based processes. The results also evidence potential antiviral and antibacterial activities of the three MnMoO₄ polymorphs.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10853-022-07277-7.

Study on Hydration and Dehydration of Ezetimibe by Terahertz Spectroscopy with Humidity-Controlled Measurements and Theoretical Analysis

Understanding the solid-state transitions of active pharmaceutical ingredients (APIs) is essential for quality control since differences in their forms affect the bioavailability of APIs. Terahertz (THz) frequency-domain spectroscopy is suitable for such an application since it can sensitively probe the lattice phonon modes originating in the crystal structures. THz absorption spectra were obtained for ezetimibe (EZT) and ezetimibe monohydrate (EZT-MH), which have similar crystalline structures and belong to the same space group. The observed absorption spectrum of EZT matched well with the solid-state density functional theory (ss-DFT)-simulated spectrum for the structures at 0 K and room temperature (modeled using constrained unit cell volumes). However for EZT-MH, the ss-DFT spectrum of the room-temperature structure showed better correlation with the experimental THz spectrum than that of the simulated spectrum of the 0 K structures, suggesting that the EZT-MH crystal has greater anharmonic character. Gibbs free-energy curves were calculated, and EZT-MH was found to be more stable than pure EZT and water in a broad temperature range. The hydrate stability may be influenced by the existence of more hydrogen bonds in EZT-MH. The hydration and dehydration of EZT in a pure API tablet and formulation tablets were monitored using a THz spectrometer with a humidity-controlled sample chamber. The effect of the excipient in the formulation tablet on hydration was successfully confirmed by showing that the solid-state transition of the API with excipients is significantly slower than that without it. Under a relative humidity of 60%, hydration of EZT in a pure EZT tablet occurred in 200 min, while the hydration of EZT in a formulation tablet was 50 times slower.

Computer Generated Realistic Interstellar Icy Grain Models: Physicochemical Properties and Interaction with NH3

Interstellar grains are composed by a rocky core (usually amorphous silicates) covered by an icy mantle, the most abundant molecule being H₂O followed by CO, CO₂, NH₃, and also radicals in minor quantities. In dense molecular clouds, gas-phase chemical species freeze onto the grain surface, making it an important reservoir of molecular diversity/complexity whose evolution leads to interstellar complex organic molecules (iCOMs). Many different models of water clusters have appeared in the literature, but without a systematic study on the properties of the grain (such as the H-bonds features, the oxygen radial distribution function, the dangling species present on the mantle surface, the surface electrostatic potential, etc.). In this work, we present a computer procedure (ACO-FROST) grounded on the newly developed semiempirical GFN2 tight-binding quantum mechanical method and the GFN-FF force field method to build-up structures of amorphous ice of large size. These methods show a very favorable accuracy/cost ratio as they are ideally designed to take noncovalent interactions into account. ACO-FROST program can be tuned to build grains of different composition mimicking dirty icy grains. These icy grain models allow studying the adsorption features (structure, binding energy, vibrational frequencies, etc.) of relevant species on a large variety of adsorption sites so to obtain a statistically meaningful distribution of the physicochemical properties of interest to be transferred in numerical models. As a test case, we computed the binding energy of ammonia adsorbed at the different sites of the icy grain surface, showing a broad distribution not easily accounted for by other more size limited icy grain models. Our method is also the base for further refinements, adopting the present grain in a more rigorous QM:MM treatment, capable of giving binding energies within the chemical accuracy.

Regression Machine Learning Models Used to Predict DFT-Computed NMR Parameters of Zeolites

Machine learning approaches can drastically decrease the computational time for the predictions of spectroscopic properties in materials, while preserving the quality of the computational approaches. We studied the performance of kernel-ridge regression (KRR) and gradient boosting regressor (GBR) models trained on the isotropic shielding values, computed with density-functional theory (DFT), in a series of different known zeolites containing out-of-frame metal cations or fluorine anion and organic structure-directing cations. The smooth overlap of atomic position descriptors were computed from the DFT-optimised Cartesian coordinates of each atoms in the zeolite crystal cells. The use of these descriptors as inputs in both machine learning regression methods led to the prediction of the DFT isotropic shielding values with mean errors within 0.6 ppm. The results showed that the GBR model scales better than the KRR model.

Quantum mechanical simulation of various phases of KVF3perovskite

The relative stability ΔEof the cubicPm3¯m(C), of the two tetragonalP4mbm(T1) andI4mcm(T2), and of the orthorhombicPbnm(O) phases of KVF₃has been computed both for the ferromagnetic (FM) and antiferromagnetic (AFM) solutions, by using the B3LYP full range hybrid functional and the Hartree-Fock (HF) Hamiltonian, an all-electron Gaussian type basis set and the CRYSTAL code. The stabilization of the T2 phase with respect to the C one (152μHa for B3LYP, 180μHa for HF, per 2 formula units) is due to the rotation of the VF₆octahedra with respect to thecaxis, by 4.1-4.6 degrees. The O phase is slightly less stable than the T2 phase (by 6 and 20μHa for B3LYP and HF); it is, however, a stable structure as the dynamical analysis confirms. The mechanism of the stabilization of the AFM solution with respect to the FM one is discussed through the spin density maps, and is related to the key role of the exact exchange term (20% in B3LYP, 100% in HF). The G-AFM phase (the first six neighbors of the reference V ion with spin reversed) is more stable than the FM one by about 500 (HF) and 1800 (B3LYP)μHa per two formula units. A volume reduction is observed in the C to T passage, and in the FM to AFM one, both being of the order of 0.3-0.5A˚3at the B3LYP level. Atomic charges, magnetic moments and bond populations, evaluated according to a Mulliken partition of the charge a spin density functions, complete the analysis. The IR and Raman spectra of the FM and AFM C, T2 and O cells are discussed; the only noticeable difference between the various space groups appears in the modes with wavenumbers lower than 100 cm^-1.

Lattice dynamics of BaFe2Se3

This paper presents a study of the lattice dynamics in BaFe₂Se₃. We combined first-principle calculations, infrared measurements and a thorough symmetry analysis. Our study confirms thatPnmacannot be the space group of BaFe₂Se₃, even at room temperature. The phonons assignment requiresPmto be the BaFe₂Se₃space group, not only in the magnetic phase, but also in the paramagnetic phase at room temperature. This is due to a strong coupling between a short-range spin-order along the ladders, and the lattice degrees of freedom associated with the Fe-Fe bond length. This coupling induces a change in the bond-length pattern from an alternated trapezoidal one (as inPnma) to an alternated small/large rectangular one. Out of the two patterns, only the latter is fully compatible with the observed block-type magnetic structure. Finally, we propose a complete symmetry analysis of the BaFe₂Se₃phase diagram in the 0-600 K range.

Chiral Lanthanum Metal-Organic Framework with Gated CO2 Sorption and Concerted Framework Flexibility

A metal–organic framework (MOF) CTH-17 based on lanthanum(III) and the conformationally chiral linker 1,2,3,4,5,6-hexakis(4-carboxyphenyl)benzene, cpb^6–: [La₂(cpb)]·1.5dmf was prepared by the solvothermal method in dimethylformamide (dmf) and characterized by variable-temperature X-ray powder diffraction (VTPXRD), variable-temperature X-ray single-crystal diffraction (SCXRD), and thermogravimetric analysis (TGA). CTH-17 is a rod-MOF with new topology och. It has high-temperature stability with Sohncke space groups P6₁22/P6₅22 at 90 K and P622 at 300 and 500 K, all phases characterized with SCXRD and at 293 K also with three-dimensional (3D) electron diffraction. VTPXRD indicates a third phase appearing after 620 K and stable up to 770 K. Gas sorption isotherms with N₂ indicate a modest surface area of 231 m² g^–1 for CTH-17, roughly in agreement with the crystal structure. Carbon dioxide sorption reveals a gate-opening effect of CTH-17 where the structure opens up when the loading of CO₂ reaches approximately ∼0.45 mmol g^–1 or 1 molecule per unit cell. Based on the SCXRD data, this is interpreted as flexibility based on the concerted movements of the propeller-like hexatopic cpb linkers, the movement intramolecularly transmitted by the π–π stacking of the cpb linkers and helped by the fluidity of the LaO₆ coordination sphere. This was corroborated by density functional theory (DFT) calculations yielding the chiral phase (P622) as the energy minimum and a completely racemic phase (P6/mmm), with symmetric cpb linkers representing a saddle point in a racemization process.

Partial substitution of the Mn atoms in CaMnO3 by first row transition metal atoms: effect on oxygen vacancy formation

Sustainable hydrogen production is one of the most challenging topics in modern energy economics. Electrochemical and thermal splitting of water are promising techniques, but are highly energy demanding. Efficient hydrogen evolution reaction catalysts can play a key role to lower the electrolysis potential and to make water splitting more feasible. Among other perovskites, CaMnO3 has been identified as efficient electrode material due to its relatively high oxygen vacancy concentrations at elevated temperatures. But this compound needs to be further improved for technical use. In this study, the effect of Mn substitution in CaMnO3 by 3d metals M = Sc, Ti, V, Cr, Fe, Co, Ni, Cu and Zn on the oxygen vacancy formation energy is investigated theoretically at DFT level. Vacancy formation energies, enthalpies and free enthalpies are calculated with a combination of hybrid and GGA density functionals. Configuration entropy is taken into account by calculating all possible configurations of M and oxygen vacancy site in supercell models. The calculated oxygen vacancy formation energies are strongly affected by Mn/M substitution, the most promising heteroelement being Cu. CaMn0.875Cu0.125O2.875 and CaMn0.875Cu0.125O2.75 are in equilibrium at 537 K, compared to 1231 K for CaMnO2.875/CaMnO2.75 and 1350 K for CaMnO3/CaMnO2.875.

Photocatalytic Degradation of Rhodamine B Dye and Hydrogen Evolution by Hydrothermally Synthesized NaBH4—Spiked ZnS Nanostructures

Metal sulphides, including zinc sulphide (ZnS), are semiconductor photocatalysts that have been investigated for the photocatalytic degradation of organic pollutants as well as their activity during the hydrogen evolution reaction and water splitting. However, devising ZnS photocatalysts with a high overall quantum efficiency has been a challenge due to the rapid recombination rates of charge carriers. Various strategies, including the control of size and morphology of ZnS nanoparticles, have been proposed to overcome these drawbacks. In this work, ZnS samples with different morphologies were prepared from zinc and sulphur powders via a facile hydrothermal method by varying the amount of sodium borohydride used as a reducing agent. The structural properties of the ZnS nanoparticles were analysed by X-ray diffraction (XRD), scanning electron microscopy (SEM), and X-ray photoelectron spectroscopy (XPS) techniques. All-electron hybrid density functional theory calculations were employed to elucidate the effect of sulphur and zinc vacancies occurring in the bulk as well as (220) surface on the overall electronic properties and absorption of ZnS. Considerable differences in the defect level positions were observed between the bulk and surface of ZnS while the adsorption of NaBH₄ was found to be highly favourable but without any significant effect on the band gap of ZnS. The photocatalytic activity of ZnS was evaluated for the degradation of rhodamine B dye under UV irradiation and hydrogen generation from water. The ZnS nanoparticles photo-catalytically degraded Rhodamine B dye effectively, with the sample containing 0.01 mol NaBH₄ being the most efficient. The samples also showed activity for hydrogen evolution, but with less H₂ produced compared to when untreated samples of ZnS were used. These findings suggest that ZnS nanoparticles are effective photocatalysts for the degradation of rhodamine B dyes as well as the hydrogen evolution, but rapid recombination of charge carriers remains a factor that needs future optimization.

The binding of atomic hydrogen on graphene from density functional theory and diffusion Monte Carlo calculations

In this work, density functional theory (DFT) and diffusion Monte Carlo (DMC) methods are used to calculate the binding energy of a H atom chemisorbed on the graphene surface. The DMC value of the binding energy is about 16% smaller in magnitude than the Perdew–Burke–Ernzerhof (PBE) result. The inclusion of exact exchange through the use of the Heyd–Scuseria–Ernzerhof functional brings the DFT value of the binding energy closer in line with the DMC result. It is also found that there are significant differences in the charge distributions determined using PBE and DMC approaches.

Rotamers in Crystal Structures of Xylitol, D-Arabitol and L-Arabitol

Rotamers are stereoisomers produced by rotation (twisting) about σ bonds and are often rapidly interconverting at room temperature. Xylitol-massively produced sweetener-(2R,3r,4S)-pentane-1,2,3,4,5-pentol) forms rotamers from the linear conformer by rotation of a xylitol fragment around the C2-C3 bond (rotamer 1) or the C3-C4 bond (rotamer 2). The rotamers form two distinguishable structures. Small differences in geometry of rotamers of the main carbon chain were confirmed by theoretical calculations; however, they were beyond the capabilities of the X-ray powder diffraction technique due to the almost identical unit cell parameters. In the case of rotamers of similar compounds, the rotations occurred mostly within hydroxyl groups likewise rotations in L-arabitol and D-arabitol, which are discussed in this work. Our results, supported by theoretical calculations, showed that energetic differences are slightly higher for rotamers with rotations within hydroxyl groups instead of a carbon chain.

Effect of cation configuration and solvation on the band positions of zinc ferrite (100)

Zinc ferrite ZnFe[Formula: see text]O[Formula: see text] belongs to the spinel-type ferrites that have been proposed as photocatalysts for water splitting. The electronic band gap and the band edge positions are of utmost importance for the efficiency of the photocatalytic processes. We, therefore, calculated the absolute band energies of the most stable surface of ZnFe[Formula: see text]O[Formula: see text], the Zn-terminated (100) surface at self-consistent hybrid density functional theory level. The effect of Fe- and Zn-rich environments, cation exchange as antisite defects and implicit solvation on the band positions is investigated. Calculated flat band potentials of the pristine surface model ranges from [Formula: see text] to [Formula: see text] V against SHE in vacuum. For Zn-rich (Fe-rich) models this changes 0.3-0.9 (0.0-0.7) V against SHE. Fe-rich models are closest to the experimental range of reported flat band potentials. Solvent effects lower the calculated flat band potentials by up to 1.8 eV. The calculated band gaps range from 1.5 to 2.9 eV in agreement with previous theoretical work and experiment. Overall, our calculations confirm the experimentally observed low activity of ZnFe[Formula: see text]O[Formula: see text] and its dependence on preparation conditions.

Full Configuration Interaction Quantum Monte Carlo treatment of fragments embedded in a periodic mean field

We present an embedded fragment approach for high-level quantum chemical calculations on local features in periodic systems. The fragment is defined as a set of localized orbitals (occupied and virtual) corresponding to a converged periodic Hartree-Fock solution. These orbitals serve as the basis for the in-fragment post-Hartree Fock treatment. The embedding field for the fragment, consisting of the Coulomb and exchange potential from the rest of the crystal, is included in the fragment's one-electron Hamiltonian. As an application of the embedded fragment approach we investigate the performanceof full configuration interaction quantum Monte Carlo (FCIQMC) with the adaptive shift. As the orbital choice we use the natural orbitals from the distinguishable cluster method with singles and doubles. FCIQMC is a stochastic approximation to the full CI method and can be routinely applied to much larger active spaces than the latter. This makes this method especially attractive in the context of open shell defects in crystals, where fragments of adequate size can be ratherlarge. As a test case we consider dissociation of a fluorine atom from a fluorographane surface. This process poses a challenge for high-level electronic structure models as both the static and dynamic correlations are essential here. Furthermore the active space for an adequate fragment (32 electrons in 173 orbitals) is already quite large even for FCIQMC. Despite this, FCIQMC delivers accurate dissociation and total energies.

Detailed Structural Features of the Perovskite-Related Halide RbPbI3 for Solar Cell Applications

All-inorganic lead halide perovskites like CsPbBr₃, CsPbI₃, or RbPbI₃ are good replacements for the classical hybrid organic-inorganic perovskites like CH₃NH₃PbI₃, susceptible to fast degradation in the presence of humid air. They also exhibit outstanding light absorption properties suitable for solar energy applications. Here, we describe the synthesis of RbPbI₃ by mechanochemical procedures with green credentials, avoiding toxic or expensive organic solvents; this specimen exhibits excellent crystallinity. We report neutron powder diffraction data, essential to revisit some subtle structural features around room temperature (200-400 K). In all these regimes, the orthorhombic Pnma crystal structure is characterized by the presence along the b direction of the crystal of double rows of edge-sharing PbI₆ octahedra. The lone electron pairs of Pb²⁺ ions have a strong stereochemical effect on the PbI₆ octahedral distortion. The relative covalency of Rb-I versus Pb-I bonds shows that the Pb-I-related motions are more rigid than Rb-I-related vibrations, as seen in the Debye temperatures from the evolution of the anisotropic displacements. The optical gap, measured by diffuse reflectance UV-vis spectroscopy, is ∼2.51 eV and agrees well with ab initio calculations. The thermoelectric Seebeck coefficient is 3 orders of magnitude larger than that of other halide perovskites, with a value of ∼117,000 μV·K^-1 at 460 K.

A density-functional benchmark of vibrational free-energy corrections for molecular crystal polymorphism

Many crystal structure prediction protocols only concern themselves with the electronic energy of molecular crystals. However, vibrational contributions to the free energy ( Fvib) can be significant in determining accurate stability rankings for crystal candidates. While force-field studies have been conducted to gauge the magnitude of these free-energy corrections, highly accurate results from quantum mechanical methods, such as density-functional theory (DFT), are desirable. Here, we introduce the PV17 set of 17 polymorphic pairs of organic molecular crystals, for which plane wave DFT is used to calculate the vibrational free energies and free-energy differences (Δ Fvib) between each pair. Our DFT results confirm that the vibrational free-energy corrections are small, having a mean value of 1.0 kJ/mol and a maximum value of 2.3 kJ/mol for the PV17 set. Furthermore, we assess the accuracy of a series of lower-cost DFT, semi-empirical, and force-field models for computing Δ Fvib that have been proposed in the literature. It is found that calculating Fvib using the Γ-point frequencies does not provide Δ Fvib values of sufficiently high quality. In addition, Δ Fvib values calculated using various approximate methods have mean absolute errors relative to our converged DFT results of equivalent or larger magnitude than the vibrational free-energy corrections themselves. Thus, we conclude that, in a crystal structure prediction protocol, it is preferable to forego the inclusion of vibrational free-energy corrections than to estimate them with any of the approximate methods considered here.

BSSE-corrected consistent Gaussian basis sets of triple-zeta valence with polarization quality of the fifth period for solid-state calculations

Revised versions of our published pob-TZVP basis sets [Laun, J.; Vilela Oliveira, D. and Bredow, T., J. Comput. Chem., 2018, 39 (19), 1285-1290] have been derived for periodic quantum-chemical solid-state calculations. They complete our pob-TZVP-rev2 series [Vilela Oliveira, D.; Laun, J.; Peintinger, M. F. and Bredow, T., J. Comput. Chem., 2019, 40 (27), 2364-2376 and Laun, J. and Bredow, J. Comput. Chem. 2021; 42 (15), 1064-1072] for the elements of the fifth period and are based on the fully relativistic effective core potentials (ECPs) of the Stuttgart/Cologne group and the def2-TZVP valence basis of the Ahlrichs group. The pob-TZVP-rev2 basis sets are developed to minimize the basis set superposition error (BSSE) in crystalline systems. For the applied PW1PW hybrid functional, the overall performance, transferability, and SCF stability of the resulting pob-TZVP-rev2 basis sets are significantly improved compared to the original pob-TZVP basis sets. After augmentation with single diffuse s- and p-functions, reference plane-wave band structures of metals can be accurately reproduced.

Density Functional Theory Investigation of the Biocatalytic Mechanisms of pH-Driven Biomimetic Behavior in CeO2

There is considerable interest in the pH-dependent, switchable, biocatalytic properties of cerium oxide (CeO₂) nanoparticles in biomedicine, where these materials exhibit beneficial antioxidant activity against reactive oxygen species (ROS) at a basic physiological pH but cytotoxic prooxidant activity in an acidic cancer cell pH microenvironment. While the general characteristics of the role of oxygen vacancies are known, the mechanism of their action at the atomic scale under different pH conditions has yet to be elucidated. The present work applies density functional theory (DFT) calculations to interpret, at the atomic scale, the pH-induced behavior of the stable {111} surface of CeO₂ containing oxygen vacancies. Analysis of the surface-adsorbed media species reveals the critical role of pH on the interaction between ROS (^•O₂^- and H₂O₂) and the defective CeO₂ {111} surface. Under basic conditions, the superoxide dismutase (SOD) and catalase (CAT) biomimetic reactions can be performed cyclically, scavenging and decomposing ROS to harmless products, making CeO₂ an excellent antioxidant. However, under acidic conditions, the CAT biomimetic reaction is hindered owing to the limited reversibility of Ce³⁺ ↔ Ce⁴⁺ and formation ↔ annihilation of oxygen vacancies. A Fenton biomimetic reaction (H₂O₂ + Ce³⁺ → Ce⁴⁺ + OH^- + ^•OH) is predicted to occur simultaneously with the SOD and CAT biomimetic reactions, resulting in the formation of hydroxyl radicals, making CeO₂ a cytotoxic prooxidant.

Two-Dimensional Anisotropic Flexibility of Mechanically Responsive Crystalline Cadmium(II) Coordination Polymers

Crystals of a family of six one-dimensional (1D) coordination polymers of cadmium(II) with cyanopyridines [[CdX₂L₂] _n , where X = Cl, Br, or I and L = 3-cyanopyridine (3-CNpy) or 4-cyanopyridine (4-CNpy)] presented a variety of morphologies and mechanical responses with dominant two-dimensional (2D) anisotropic flexibility, which has not been previously reported. All mechanically adaptable crystals were 2D flexible and displayed a variety of direction-dependent responses; in addition to 2D isotropic flexibility observed for solely elastic materials, 2D anisotropic flexibility was noticed for both elastic and elastic → plastic crystals. The consequences of fine and controlled structural variations on mechanical behavior were additionally explored via microfocus single-crystal X-ray diffraction and complementary theoretical studies, revealing that the relative strength and direction of the hydrogen bonding interactions were the key parameters in delivering a specific mechanical response.

High-pressure transformation of dithiazolylidene-dithiazolium polyiodide with N-H…N hydrogen bond: A Raman Spectroscopy study

The insight into the behavior of polyiodides under non-ambient conditions can enrich the practical applications due to obtaining materials with adjustable and tunable conducting properties. In this work Raman spectroscopy study in the range 0 - 6.5 GPa has been performed for dithiazolylidene-dithiazolium zigzag polyiodide with N-H…N hydrogen bond. Variations of band positions in the low-wavenumber region of Raman spectra have been attributed to the changes in the anionic part of the unit cell. The association of different interacting polyiodide subunits at pressure above 1.5 GPa leads to the emergence of the bands at 98 and 115 cm^-1 due to consistent vibrations in the I₈^2- zigzag. The emergence of the band at 159 cm^-1 above 0.28 GPa can be assigned to N-H…N hydrogen bond vibration. The obtained results depict undergoing structure transformations: the stronger are cation…cation and anion…anion interactions the more likely is the observation of conducting properties due to the formation of polyiodide zigzag and strongly-bound hydrogen-bonded fragment.

Characteristic Spectral Features of Terra Preta (TP) in the 5-15 Terahertz Range

Terra preta is a fertile anthropogenic soil found in the Amazon basin. One of the most significant differences between the terra preta and surrounding soils is that terra preta is rich in aromatic carbons. Previous infrared investigations of terra preta were reported at energies above 1000 cm^-1 where many other forms of carbon also have absorption lines. No measurements have been reported below 800 cm^-1, where many absorptions associated with aromatic carbons occur in the absence of aliphatic carbon lines. We employ Fourier transform infrared spectroscopy between 150 cm^-1 and 500 cm^-1. A comparison was made between the spectra of terra preta, several pure aromatic compounds, organic fertilizers developed to replicate terra preta and several Australian soils, some of which containing char from bushfires. The spectra in the 150-500 cm^-1 range were very similar between terra preta and the organic fertilizers, while they were very different for the natural soils. These findings indicate that the content of aromatic carbons in terra preta and organic fertilizers is different than in natural soils containing the bushfire chars, but also soils produced entirely by bacterial and fungal activities. This point to the importance of the preparation conditions of the biochars, which are essential ingredients of terra preta and organic fertilizers used in this study.

H2O2adsorption and dissociation on various CeO2(111) surface models: a first-principles study

Periodic density functional theory (DFT) calculations using the hybrid PBE0 functional and atom-centered Gaussian functions as basis sets were carried out to investigate the absorption and the first steps involved in the decomposition of hydrogen peroxide (H₂O₂) on three different models of the ceria (111) surface. One of the models is a clean surface, and the others are defective and partially hydroxylated ceria surfaces. On the clean surface, we found that the minimum energy path of hydrogen peroxide decomposition involves a three-step process, i.e., adsorption, deprotonation, and formation of the peroxide anion, stabilized through its interaction with the surface at a Ce (IV) site, with activation barriers of less than about 0.5 eV. The subsequent formation of superoxide anions and molecular oxygen species is attributed to electron transfer from the reactants to the Ce (IV) ions underneath. On the defective surface, H₂O₂dissociation is an energetically downhill reaction thermodynamically driven by the healing of the O vacancies, after the reduction and decomposition of H₂O₂into oxygen and water. On the hydroxylated surface, H₂O₂is first adsorbed by forming a favorable H-bond and then undergoes heterolytic dissociation, forming two hydroxyl groups at two vicinal Ce sites.

Two New Organic Co-Crystals Based on Acetamidophenol Molecules

Herein we present two new organic co-crystals obtained through a simple solution growth process based on an acetamidophenol molecule, either paracetamol or metacetamol, and on 7,7,8,8-tetracyanoquinodimethane (TCNQ). These co-crystals are part of a family of potential organic charge transfer complexes, where the acetamidophenol molecule behaves as an electron donor and TCNQ behaves as an electron acceptor. Due to the sub-micron size of the crystalline domains, 3D electron diffraction was employed for the structure characterization of both systems. Paracetamol-TCNQ structure was solved by standard direct methods, while the analysis of metacetamol-TCNQ was complicated by the low resolution of the available diffraction data and by the low symmetry of the system. The structure determination of metacetamol-TCNQ was eventually achieved after merging two data sets and combining direct methods with simulated annealing. Our study reveals that both paracetamol-TCNQ and metacetamol-TCNQ systems crystallize in a 1:1 stoichiometry, assembling in a mixed-stack configuration and adopting a non-centrosymmetric P1 symmetry. It appears that paracetamol and metacetamol do not form a strong structural scaffold based on hydrogen bonding, as previously observed for orthocetamol-TCNQ and orthocetamol-TCNB (1,2,4,5-tetracyanobenzene) co-crystals.

New functionalisation reactions of graphitic carbon nitrides: Computational and experimental studies

The functionalisation of two-dimensional materials is key to modify their properties and facilitate assembly into functional devices. Here, new reactions have been proposed to modify crystalline two-dimensional carbon nitrides of polytriazine imide structure. Both amine alkylation and aryl-nitrene-based reactions have been explored computationally and with exploratory synthetic trials. The approach illustrates that alkylation is unfavourable, particularly at basal-plane sites. In contrast, while initial trial reactions were inconclusive, the radical-addition of nitrenes is shown to be energetically favourable, with a preference for functionalising sheet edges to minimise steric effects.

Strength and Nature of Host-Guest Interactions in Metal-Organic Frameworks from a Quantum-Chemical Perspective

Metal‐organic frameworks (MOFs) offer a convenient means for capturing, transporting, and releasing small molecules. Their rational design requires an in‐depth understanding of the underlying non‐covalent host‐guest interactions, and the ability to easily and rapidly pre‐screen candidate architectures in silico. In this work, we devised a recipe for computing the strength and analysing the nature of the host‐guest interactions in MOFs. By assessing a range of density functional theory methods across periodic and finite supramolecular cluster scale we find that appropriately constructed clusters readily reproduce the key interactions occurring in periodic models at a fraction of the computational cost. Host‐guest interaction energies can be reliably computed with dispersion‐corrected density functional theory methods; however, decoding their precise nature demands insights from energy decomposition schemes and quantum‐chemical tools for bonding analysis such as the quantum theory of atoms in molecules, the non‐covalent interactions index or the density overlap regions indicator.

Can Mesoporous Silica Speed Up Degradation of Benzodiazepines? Hints from Quantum Mechanical Investigations

This work reports for the first time a quantum mechanical study of the interactions of a model benzodiazepine drug, i.e., nitrazepam, with various models of amorphous silica surfaces, differing in structural and interface properties. The interest in these systems is related to the use of mesoporous silica as carrier in drug delivery. The adopted computational procedure has been chosen to investigate whether silica–drug interactions favor the drug degradation mechanism or not, hindering the beneficial pharmaceutical effect. Computed structural, energetics, and vibrational properties represent a relevant comparison for future experiments. Our simulations demonstrate that adsorption of nitrazepam on amorphous silica is a strongly exothermic process in which a partial proton transfer from the surface to the drug is observed, highlighting a possible catalytic role of silica in the degradation reaction of benzodiazepines.

QUANTAS: a Python software for the analysis of thermodynamics and elastic behavior of solids from ab initio quantum mechanical simulations and experimental data

This paper presents QUANTAS, an open-source Python-based software aimed at providing a fast, flexible and easy-to-use framework to calculate the thermodynamics and elastic properties of crystalline solids. QUANTAS could be of use for researchers involved in various fields of solid-state chemistry, physics and mineralogy.

Mineralogy, petrology and materials science are fundamental disciplines not only for the basic knowledge and classification of solid phases but also for their technological applications, which are becoming increasingly demanding and challenging. Characterization and design of materials are of utmost importance and usually need knowledge of the thermodynamics and mechanical stability of solids. Alongside well known experimental approaches, in recent years the advances in both quantum mechanical methods and computational power have placed theoretical investigations as a complementary useful and powerful tool in this kind of study. In order to aid both theoreticians and experimentalists, an open-source Python-based software, QUANTAS, has been developed. QUANTAS provides a fast, flexible, easy-to-use and extensible platform for calculating the thermodynamics and elastic behavior of crystalline solid phases, starting from both experimental and ab initio data.