We present a new software package called M-Chem that is designed from scratch in C++ and parallelized on shared-memory multi-core architectures to facilitate efficient molecular simulations. Currently, M-Chem is a fast molecular dynamics (MD) engine that supports the evaluation of energies and forces from two-body to many-body all-atom potentials, reactive force fields, coarse-grained models, combined quantum mechanics molecular mechanics (QM/MM) models, and external force drivers from machine learning, augmented by algorithms that are focused on gains in computational simulation times. M-Chem also includes a range of standard simulation capabilities including thermostats, barostats, multi-timestepping, and periodic cells, as well as newer methods such as fast extended Lagrangians and high quality electrostatic potential generation. At present M-Chem is a developer friendly environment in which we encourage new software contributors from diverse fields to build their algorithms, models, and methods in our modular framework. The long-term objective of M-Chem is to create an interdisciplinary platform for computational methods with applications ranging from biomolecular simulations, reactive chemistry, to materials research.

We report an investigation of structure and photophysics of thin layers of cibalackrot, a sturdy dye derived from indigo by double annulation at the central double bond. Evaporated layers contain up to three phases, two crystalline and one amorphous. Relative amounts of all three have been determined by a combination of X-ray diffraction and FT-IR reflectance spectroscopy. Initially, excited singlet state rapidly produces a high yield of a transient intermediate whose spectral properties are compatible with charge-transfer nature. This intermediate more slowly converts to a significant yield of triplet, which, however, does not exceed 100% and may well be produced by intersystem crossing rather than singlet fission. The yields were determined by transient absorption spectroscopy and corrected for effects of partial sample alignment by a simple generally applicable procedure. Formation of excimers was also observed. In order to obtain guidance for improving molecular packing by a minor structural modification, calculations by a simplified frontier orbital method were used to find all local maxima of singlet fission rate as a function of geometry of a molecular pair. The method was tested at 48 maxima by comparison with the ab initio Frenkel-Davydov exciton model.

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.