Lowest-energy structures, the distribution of isomers, and their molecular properties depend significantly on geometry and temperature. Total energy computations using DFT methodology are typically carried out at a temperature of zero K; thereby, en-tropic contributions to the total energy are neglected, even though functional materials work at finite temperatures. In the present study, the probability of the occurrence of one particular Be4B8 isomer at temperature T is estimated by employing Gibbs free energy computed within the framework of quantum statistical mechanics and nanothermody-namics. To identify a list of all possible low-energy chiral and achiral structures, an ex-haustive and efficient exploration of the potential/free energy surfaces is carried out using a multi-level multistep global genetic algorithm search coupled with DFT. In addition, we discuss the energetic ordering of structures computed at the DFT level against single-point energy calculations at the CCSD(T) level of theory. The total VCD/IR spectra as a function of temperature are computed using each isomer’s probability of occurrence in a Boltzmann-weighted superposition of each isomer’s spectrum. Additionally, we present chemical bonding analysis using the adaptive natural density partitioning method in the chiral putative global minimum. The transition state structures and the enantiomer–en-antiomer and enantiomer–achiral activation energies as a function of temperature evidence that a change from an endergonic to an exergonic type of reaction occurs at a temperature of 739 K.
Bibliographical notePublisher Copyright:
© 2021 by the authors. Licensee MDPI, Basel, Switzerland.
- Adaptive natural density partitioning method
- Beryllium–boron cluster
- Boltzmann factors
- Density functional theory
- Genetic algorithm
- Gibbs free energy
- Global minimum
- IR spectra
- Quantum statistical mechanics
- Vibrational circular dichroism