Noncovalent Interactions in SIESTA Using the vdW-DF Functional: S22 Benchmark and Macrocyclic Structures

Mechanoresponsive Flexible Crystals

Flexible crystals have gained significant attention owing to their remarkable pliability, plasticity, and adaptability, making them highly popular in various research and application fields. The main challenges in developing flexible crystals lie in the rational design, preparation, and performance optimization of such crystals. Therefore, a comprehensive understanding of the fundamental origins of crystal flexibility is crucial for establishing evaluation criteria and design principles. This Perspective offers a retrospective analysis of the development of flexible crystals over the past two decades. It summarizes the elastic standards and possible plastic bending mechanisms tailored to diverse flexible crystals and analyzes the assessment of their theoretical basis and applicability. Meanwhile, the compatibility between crystal elasticity and plasticity has been discussed, unveiling the immense prospects of elastic/plastic crystals for applications in biomedicine, flexible electronic devices, and flexible optics. Furthermore, this Perspective presents state-of-the-art experimental avenues and analysis methods for investigating molecular interactions in molecular crystals, which is vital for the future exploration of the mechanisms of crystal flexibility.

Non-covalent interactions for carbonaceous materials: impacts of doping, curving and their combination

Non-covalent interactions and their influencing factors are addressed after studying the binding of ion pairs with graphene (doped, curved and combination).

Adaptively Compressed Exchange Operator for Large-Scale Hybrid Density Functional Calculations with Applications to the Adsorption of Water on Silicene

Density functional theory (DFT) calculations using hybrid exchange-correlation functionals have been shown to provide an accurate description of the electronic structures of nanosystems. However, such calculations are often limited to small system sizes due to the high computational cost associated with the construction and application of the Hartree-Fock (HF) exchange operator. In this paper, we demonstrate that the recently developed adaptively compressed exchange (ACE) operator formulation [J. Chem. Theory Comput. 2016, 12, 2242-2249] can enable hybrid functional DFT calculations for nanosystems with thousands of atoms. The cost of constructing the ACE operator is the same as that of applying the exchange operator to the occupied orbitals once, while the cost of applying the Hamiltonian operator with a hybrid functional (after construction of the ACE operator) is only marginally higher than that associated with applying a Hamiltonian constructed from local and semilocal exchange-correlation functionals. Therefore, this new development significantly lowers the computational barrier for using hybrid functionals in large-scale DFT calculations. We demonstrate that a parallel planewave implementation of this method can be used to compute the ground-state electronic structure of a 1000-atom bulk silicon system in less than 30 wall clock minutes and that this method scales beyond 8000 computational cores for a bulk silicon system containing about 4000 atoms. The efficiency of the present methodology in treating large systems enables us to investigate adsorption properties of water molecules on Ag-supported two-dimensional silicene. Our computational results show that water monomer, dimer, and trimer configurations exhibit distinct adsorption behaviors on silicene. In particular, the presence of additional water molecules in the dimer and trimer configurations induces a transition from physisorption to chemisorption, followed by dissociation on Ag-supported silicene. This is caused by the enhanced effect of hydrogen bonds on charge transfer and proton transfer processes. Such a hydrogen bond autocatalytic effect is expected to have broad applications for silicene as an efficient surface catalyst for oxygen reduction reactions and water dissociation.

Molecular modeling and infrared and Raman spectroscopy of the crystal structure of the chiral antiparasitic drug Praziquantel

Anthelmintic praziquantel (PZQ) and its molecular and crystal lattice structures were studied by means of atomistic calculations based on empirical interatomic potentials and quantum mechanical methods (DFT). This chiral drug presents several crystal polymorphs due to the enantiomers (either R or S) or the racemic crystal, and different molecular conformations. The relative configurations of the carbonyl groups in PZQ define these conformations that produce different polymorphs. The polarity of these conformers is quite different and their relative population can vary in media with different polarity. Crystal structures of PZQ were studied by infrared and Raman spectroscopy and their spectroscopical properties were calculated by quantum mechanical methods, assigning many of their bands, and finding a good agreement with experimental data.

Structural Collapse of the Hydroquinone-Formic Acid Clathrate: A Pressure-Medium-Dependent Phase Transition

The energy landscape governing a new pressure-induced phase transition in the hydroquinone-formic acid clathrate is reported in which the host structure collapses, opening up the cavity channels within which the guest molecules migrate and order. The reversible isosymmetric phase transition causes significant changes in the morphology and the birefringence of the crystal. The subtle intermolecular interaction energies in the clathrate are quantified at varying pressures using novel model energies and energy frameworks. These calculations show that the high-pressure phase forms a more stable host network at the expense of less-stable host-guest interactions. The phase transition can be kinetically hindered using a nonhydrostatic pressure-transmitting medium, enabling the comparison of intermolecular energies in two polymorphic structures in the same pressure range. Overall this study illustrates a need for accurate intermolecular energies when analyzing self-assembly structures and supramolecular aggregates.

Silver(I), gold(I) and palladium(II) complexes of a NHC-pincer ligand with an aminotriazine core: a comparison with pyridyl analogues

Dinuclear silver, di- and tetra-nuclear gold, and mononuclear palladium complexes with chelating C,N,C diethylaminotriazinyl-bridged bis(NHC) pincer ligands were prepared and characterised. The silver and gold complexes exist in a twisted, helical conformation in both the solution- and the solid state. In contrast, an analogous dinuclear gold complex with pyridyl-bridged NHCs exists in a linear conformation. Computational studies have been performed to rationalise the formation of twisted/helical vs. linear forms.

Accurate and Efficient Model Energies for Exploring Intermolecular Interactions in Molecular Crystals

The energy of interaction between molecules is commonly expressed in terms of four key components: electrostatic, polarization, dispersion, and exchange-repulsion. Using monomer wave functions to obtain accurate estimates of electrostatic, polarization, and repulsion energies along with Grimme's dispersion corrections, a series of energy models are derived by fitting to dispersion-corrected DFT energies for a large number of molecular pairs extracted from organic and inorganic molecular crystals. The best performing model reproduces B3LYP-D2/6-31G(d,p) counterpoise-corrected energies with a mean absolute deviation (MAD) of just over 1 kJ mol(-1) but in considerably less computation time. It also performs surprisingly well against benchmark CCSD(T)/CBS energies, with a MAD of 2.5 kJ mol(-1) for a combined data set including Hobza's X40, S22, A24, and S66 dimers. Two of these energy models, the most accurate and the fastest, are expected to find widespread application in investigations of molecular crystals.

A study of planar anchor groups for graphene-based single-molecule electronics

To identify families of stable planar anchor groups for use in single molecule electronics, we report detailed results for the binding energies of two families of anthracene and pyrene derivatives adsorbed onto graphene. We find that all the selected derivatives functionalized with either electron donating or electron accepting substituents bind more strongly to graphene than the parent non-functionalized anthracene or pyrene. The binding energy is sensitive to the detailed atomic alignment of substituent groups over the graphene substrate leading to larger than expected binding energies for -OH and -CN derivatives. Furthermore, the ordering of the binding energies within the anthracene and pyrene series does not simply follow the electron affinities of the substituents. Energy barriers to rotation or displacement on the graphene surface are much lower than binding energies for adsorption and therefore at room temperature, although the molecules are bound to the graphene, they are almost free to move along the graphene surface. Binding energies can be increased by incorporating electrically inert side chains and are sensitive to the conformation of such chains.