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.
From first principles calculations, we investigate the stability and physical properties of single layer h-BN sheet chemically functionalized by various groups viz. H, F, OH, CH3, CHO, CN, NH2 etc. We find that full functionalization of h-BN sheet wi
th these groups lead to decrease in its electronic band gap, albeit to different magnitudes varying from 0.3 eV to 3.1 eV, depending upon the dopant group. Functionalization by CHO group, in particular, leads to a sharp decrease in the electronic band gap of the pristine BN sheet to ~ 0.3 eV, which is congenial for its usage in transistor based devices. The phonon calculations on these sheets show that frequencies corresponding to all their vibrational modes are real (positive), thereby suggesting their inherent stability. The chemisorption energies of these groups to the B and N atoms of the sheet are found to lie in the range of 1.5 -6 eV.
Strain induced band gap deformations of hydrogenated/fluorinated graphene and hexagonal BN sheet have been investigated using first principles density functional calculations. Within harmonic approximation, the deformation is found to be higher for h
ydrogenated systems than for the fluorinated systems. Interestingly, our calculated band gap deformation for hydrogenated/fluorinated graphene and BN sheets are positive, while those for pristine graphene and BN sheet are found to be negative. This is due to the strong overlap between nearest neighbor {pi} orbitals in the pristine sheets, that is absent in the passivated systems. We also estimate the intrinsic strength of these materials under harmonic uniaxial strain, and find that the in-plane stiffness of fluorinated and hydrogenated graphene are close, but larger in magnitude as compared to those of fluorinated and hydrogenated BN sheet.
We propose, on the basis of our first principles density functional based calculations, a new isomer of graphane, in which the C-H bonds of a hexagon alternate in 3-up, 3-down fashion on either side of the sheet. This 2D puckered structure called sti
rrup has got a comparable stability with the previously discovered chair and boat conformers of graphane. The physico-chemical properties of this third conformer are found to be similar to the other two conformers of graphane with an insulating direct band gap of 3.1 eV at the {Gamma} point. Any other alternative hydrogenation of the graphene sheet disrupts its symmetric puckered geometry and turns out to be energetically less favorable.
