We have performed cascade genetic algorithm and ab initio atomistic thermodynamics under the framework of first-principles density functional theory to study the (meta-)stability of a wide range of LixNy clusters. We found that hybrid xc-functional i
s essential to address this problem as a local/semi-local functional simply fails even to predict a qualitative prediction. Most importantly, we find that though in bulk Lithium Nitride, Li rich phase, i.e. Li3N, is the stable stoichiometry, in small LixNy clusters N-rich phases are more stable at thermodynamic equilibrium. We further show a that these N-rich clusters are promising hydrogen storage material because of their easy adsorption and desorption ability at respectively low (< 300K) and moderately high temperature (> 600K).
In order to design clusters with desired properties, we have implemented a suite of genetic algorithms tailored to optimize for low total energy, high vertical electron affinity (VEA), and low vertical ionization potential (VIP). Applied to (TiO$_2$)
$_n$ clusters, the property-based optimization reveals the underlying structure-property relations and the structural features that may serve as active sites for catalysis. High VEA and low VIP are correlated with the presence of several dangling-O atoms and their proximity, respectively. We show that the electronic properties of (TiO$_2$)$_n$ up to n=20 correlate more strongly with the presence of these structural features than with size.
A first-principles based methodology for efficiently and accurately finding thermodynamically stable and metastable atomic structures is introduced and benchmarked. The approach is demonstrated for gas-phase metal-oxide clusters in thermodynamic equi
librium with a reactive (oxygen) atmosphere at finite pressure and temperature. It consists of two steps. At first, the potential-energy surface is scanned by means of a global-optimization technique, i.e., a massive-parallel first-principles cascade genetic algorithm for which the choice of all parameters is validated against higher-level methods. In particular, we validate a) the criteria for selection and combination of structures used for the assemblage of new candidate structures, and b) the choice of the exchange-correlation functional. The selection criteria are validated against a fully unbiased method: replica-exchange molecular dynamics. Our choice of the exchange-correlation functional, the van-der-Waals-corrected PBE0 hybrid functional, is justified by comparisons up to highest level currently achievable within density-functional theory, i.e., the renormalized second-order perturbation theory, rPT2. In the second step, the low-energy structures are analyzed by means of ab initio atomistic thermodynamics in order to determine compositions and structures that minimize the Gibbs free energy at given temperature and pressure of the reactive atmosphere.
We have carried out first-principles based DFT calculation on electronic properties of silicene monolayer on various (111) semi-conducting surfaces. We find that the relative stability and other properties of the silicene overlayer depends sensitivel
y on whether the interacting top layer of the substrate is metal or non-metal terminated. The nature of silicene-monolayer on the metal termi- nated surface can be metallic or even magnetic, depending upon the choice of the substrate. The silicene overlayer undergoes n-type doping on metal terminated surface while it undergoes p-type doping on non metal terminated surfaces of the semiconductor substrates.
By applying a genetic algorithm in a cascade approach of increasing accuracy, we calculate the composition and structure of MgMOx clusters at realistic temperatures and oxygen pressures. The stable and metastable systems are identified by ab initio a
tomistic thermodynamics. We find that small clusters (M <= 5) are in thermodynamic equilibrium when x > M. The non-stoichiometric clusters exhibit peculiar magnetic behavior, suggesting the possibility of tuning magnetic properties by changing environmental pressure and temperature conditions. Furthermore, we show that density-functional theory (DFT) with a hybrid exchange-correlation (xc) functional is needed for predicting accurate phase diagrams of metal-oxide clusters. Neither a (sophisticated) force field nor DFT with (semi)local xc functionals are sufficient for even a qualitative prediction.
