Tuning the quasi-harmonic treatment of crystalline ionic liquids within the density functional theory

First-Principles Models of Polymorphism of Pharmaceuticals: Maximizing the Accuracy-to-Cost Ratio

Accuracy and sophistication of in silico models of structure, internal dynamics, and cohesion of molecular materials at finite temperatures increase over time. Applicability limits of ab initio polymorph ranking that would be feasible at reasonable costs currently represent crystals of moderately sized molecules (less than 20 nonhydrogen atoms) and simple unit cells (containing rather only one symmetry-irreducible molecule). Extending the applicability range of the underlying first-principles methods to larger systems with a real-life significance, and enabling to perform such computations in a high-throughput regime represent additional challenges to be tackled in computational chemistry. This work presents a novel composite method that combines the computational efficiency of density-functional tight-binding (DFTB) methods with the accuracy of density-functional theory (DFT). Being rooted in the quasi-harmonic approximation, it uses a cheap method to perform all of the costly scans of how static and dynamic characteristics of the crystal vary with respect to its volume. Such data are subsequently corrected to agree with a higher-level model, which must be evaluated only at a single volume of the crystal. It thus enables predictions of structural, cohesive, and thermodynamic properties of complex molecular materials, such as pharmaceuticals or organic semiconductors, at a fraction of the original computational cost. As the composite model retains the solid physical background, it suffers from a minimum accuracy deterioration compared to the full treatment with the costly approach. The novel methodology is demonstrated to provide consistent results for the structural and thermodynamic properties of real-life molecular crystals and their polymorph ranking.

Anisotropy, segmental dynamics and polymorphism of crystalline biogenic carboxylic acids

Carboxylic acids of the Krebs cycle possess invaluable biochemical significance. Still, there are severe gaps in the availability of thermodynamic and crystallographic data, as well as ambiguities prevailing in the literature on the thermodynamic characterization and polymorph ranking. Providing an unambiguous description of the structure, thermodynamics and polymorphism of their neat crystalline phases requires a complex multidisciplinary approach. This work presents results of an extensive investigation of the structural anisotropy of the thermal expansion and local dynamics within these crystals, obtained from a beneficial cooperation of NMR crystallography and ab initio calculations of non-covalent interactions. The observed structural anisotropy and spin-lattice relaxation times are traced to large spatial variations in the strength of molecular interactions in the crystal lattice, especially in the orientation of the hydrogen bonds. A completely resolved crystal structure for oxaloacetic acid is reported for the first time. Thanks to multi-instrumental calorimetric effort, this work clarifies phase behavior, determines third-law entropies of the crystals, and states definitive polymorph ranking for succinic and fumaric acids. These thermodynamic observations are then interpreted in terms of first-principles quasi-harmonic calculations of cohesive properties. A sophisticated model capturing electronic, thermal, and configurational-entropic effects on the crystal structure approaches captures the subtle Gibbs energy differences governing polymorph ranking for succinic and fumaric acids, representing another success story of computational chemistry.

Regulating the Electronic Configuration of Spinel Zinc Manganate Derived from Metal-Organic Frameworks: Controlled Synthesis and Application in Anode Materials for Lithium-Ion Batteries

In recent years, transition metal oxides have been considered as the most promising anode materials due to their high theoretical capacity, low price, and abundant natural reserves. Among them, zinc manganate is used as an electrode material for anodes, whose application is mostly hindered due to its poor ionic/electronic conductivity. In this work, a series of ZnMn₂O₄ (ZMO) are synthesized by a hydrothermal technique coordinated with a metal-organic framework-based high-temperature calcination process for their application as an anode in lithium-ion batteries (LIBs). Meanwhile, this study systematically explores the influence of carbon doping and the types of organic ligands and oxygen vacancies on the electrochemical properties of the synthesized ZMO. Density functional theory (DFT) calculations and experimental investigations reveal that the introduction of carbon and oxygen vacancies can enhance electronic conductivity, more active sites and faster Li⁺ adsorption, resulting in better electrochemical performances. As expected, all ZMOs with carbon doping (PMA-ZMO, MI-ZMO, and BDC-ZMO) derived from 1,2,4,5-benzenetetracarboxylic acid, 2-methylimidazole, and 1,4-dicarboxybenzene achieve outstanding electrochemical performance. Meanwhile, the introduction of oxygen vacancies can enhance the electronic conductivity and can significantly reduce the activation energy of Li⁺ transport, thereby accelerating the Li⁺ diffusion kinetics in the lithiation/delithiation process. Furthermore, an optimal ZMO anode material synthesized by 2-methylimidazole delivers a high reversible capacity of 1174.7 mA h g^-1 after 300 cycles at 0.1 A g^-1 and 600 mA h g^-1 at 0.5 A g^-1 after 300 cycles. After high-rate charge and discharge cycles, the specific capacity rapidly recovers to a value greater than the initial value, which proves the unusual activation and thereby an excellent rate property of the electrode. Hence, we conclude that ZMO provides potential application prospects as an anode electrode material for LIBs.

Molecular Dynamics Calculations of Partial Molar Volumes of Amino Acids in Aqueous Solutions

Partial molar volumes of amino acids in their zwitterionic and molecular forms have been calculated using molecular dynamics simulations of these systems in aqueous solutions. Calculations performed with the TIP4P, SPC (rigid and flexible), SPC/E, and polarizable water models show that the choice of water model can have a significant impact on the calculated volumes. The effect of treatment of long-range electrostatic interactions on the calculated results was also investigated. Volumes obtained in simulations with a properly chosen water model fit well the experimental data for both molecular and zwitterionic forms of amino acids.