Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userExploration of potential and free energy surfaces of the neutral Be$_4$B$_{8}$ chiral clusters and their stabilities at finite temperatures
137
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The lowest-energy structure, distribution of isomers, and their molecular properties depend significantly on the geometry and temperature. The total energy computations under DFT methodology are typically carried out at zero temperature; thereby, entropic contributions to total energy are neglected, even though functional materials work at finite temperature. In the present study, the probability of occurrence of one particular Be$_4$B$_8$ isomer at temperature T is estimated within the framework of quantum statistical mechanics and nanothermodynamics. To locate a list of all possible low-energy chiral and achiral structures, an exhaustive and efficient exploration of the potential/free energy surface is done by employing a multilevel multistep global genetic algorithm search coupled to DFT. Moreover, 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 computed VCD/IR spectrum of each isomer is multiplied by their corresponding Boltzmann weight at temperature T; then, they are summed together to produce a final Boltzmann weighted 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-enantiomer 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.
rate research
Read More
Interatomic potential models based on machine learning (ML) are rapidly developing as tools for materials simulations. However, because of their flexibility, they require large fitting databases that are normally created with substantial manual selection and tuning of reference configurations. Here, we show that ML potentials can be built in a largely automated fashion, exploring and fitting potential-energy surfaces from the beginning (de novo) within one and the same protocol. The key enabling step is the use of a configuration-averaged kernel metric that allows one to select the few most relevant structures at each step. The resulting potentials are accurate and robust for the wide range of configurations that occur during structure searching, despite only requiring a relatively small number of single-point DFT calculations on small unit cells. We apply the method to materials with diverse chemical nature and coordination environments, marking a milestone toward the more routine application of ML potentials in physics, chemistry, and materials science.
The starting point to understanding cluster properties is the putative global minimum and all the nearby local energy minima; however, locating them is computationally expensive and challenging due to a combinatorial explosion problem. The relative populations and spectroscopic properties of a molecule that are a function of temperature can be approximately computed by employing statistical thermodynamics. Here, we investigate temperature-entropy-driven isomers distribution on Be$_6$B$_{11}^{-}$ fluxional cluster and the effect of temperature on their infrared spectroscopy and relative populations. We identify the vibration modes possessed by the cluster that significantly contribute to the zero-point energy. A couple of steps are considered for computing the temperature-dependent relative population: First, using a genetic algorithm coupled to density functional theory, we performed an extensive and systematic exploration of the potential/free energy surface of Be$_6$B$_{11}^{-}$ cluster to locate the putative global minimum and elucidate the low-energy structures. Second, the relative populations temperature effects are determined by considering the thermodynamic properties and Boltzmann factors. The temperature-dependent relative populations show that the entropies and temperature are essential for determining the global minimum. We compute the temperature-dependent total infrared spectra employing the Boltzmann factor weighted sums of each isomers infrared spectrum and find that at finite temperature, the total infrared spectrum is composed of an admixture of infrared spectra that corresponds to the spectrum of the lowest energy structure and its isomers located at higher energies. The methodology and results describe the thermal effects in the relative population and the infrared spectra.
This study reports the lowest energy structure of bare Cu$_{13}$ nanoclusters as a pair of enantiomers for temperatures ranging from 20 to 1200 K. Moreover, we compute the enantiomerization energy for the interconversion from $mathcal{M}$ to $mathcal{P}$ structures in the chiral putative global minimum for temperatures ranging from 20 to 1300 K. Additionally, employing statistical thermodynamics and nanothermodynamics, we compute the Boltzmann Probability for each particular isomer as a function of temperature. To achieve that, we explore the free energy surface of the Cu$_{13}$ cluster, employing a genetic algorithm coupled with density functional theory and statistical thermodynamics. Moreover, we discuss the energetic ordering of structures computed at the DFT level and compared to high level DLPNO-CCSD(T1) reference energies, and we present chemical bonding analysis using the AdNDP method in the chiral putative global minimum. Based on the computed relative abundances, our results show that the chiral putative global minimum strongly dominates for temperatures ranging from 20 to 1100 K.
We show that the transition origins of electronic excitations identified by quantified natural transition orbital (QNTO) analysis can be employed to connect potential energy surfaces (PESs) according to their character across a widerange of molecular geometries. This is achieved by locating the switching of transition origins of adiabatic potential surfaces as the geometry changes. The transition vectors for analysing transition origins are provided by linear response time-dependent density functional theory (TDDFT) calculations under the Tamm-Dancoff approximation. We study the photochemical CO ring opening of oxirane as an example and show that the results corroborate the traditional Gomer-Noyes mechanism derived experimentally. The knowledge of specific states for the reaction also agrees well with that given by previous theoretical work using TDDFT surface-hopping dynamics that was validated by high-quality quantum Monte Carlo calculations. We also show that QNTO can be useful for considerably larger and more complex systems: by projecting the excitations to those of a reference oxirane molecule, the approach is able to identify and analyse specific excitations of a trans-2,3-diphenyloxirane molecule.
The adsorption and dissociation of O$_{2}$ molecules at the Be(0001) surface is studied by using density-functional theory within the generalized gradient approximation and a supercell approach. The physi- and chemisorbed molecular precursor states are identified to be along the parallel and vertical channels, respectively. It is shown that the HH-Z (see the text for definition) channel is the most stable channel for the molecular chemisorption of O$_{2}$. The electronic and magnetic properties of this most stable chemisorbed molecular state are studied, which shows that the electrons transfer forth and back between the spin-resolved antibonding $pi^{ast}$ molecular orbitals and the surface Be $sp$ states. A distinct covalent weight in the molecule-metal bond is also shown. The dissociation of O$_{2}$ is determined by calculating the adiabatic potential energy surfaces, wherein the T-Y channel is found to be the most stable and favorable for the dissociative adsorption of O$_{2}$. Remarkably, we predict that unlike the other simple $sp$ metal surfaces such as Al(111) and Mg(0001), the textit{adiabatic} dissociation process of O$_{2}$ at Be(0001) is an activated type with a sizeable energy barrier.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
