Open Access Publications
Permanent URI for this collection
Open access publications by faculty, staff, postdocs, and graduate students from the Delaware Energy Institute.
Browse
Browsing Open Access Publications by Author "Medasani, Bharat"
Now showing 1 - 2 of 2
Results Per Page
Sort Options
Item OpenMKM: An Open-Source C++ Multiscale Modeling Simulator for Homogeneous and Heterogeneous Catalytic Reactions(Journal of Chemical Information and Modeling, 2023-06-12) Medasani, Bharat; Kasiraju, Sashank; Vlachos, Dionisios G.Microkinetic modeling is invaluable for coupling “microscale” atomistic data with “macroscale” reactor observables. We introduce an Open-source Microkinetic Modeling (OpenMKM) multiscale mean-field microkinetics modeling toolkit targeting mainly heterogeneous catalytic reactions but applies equally to homogeneous reactions. OpenMKM is a modular, object-oriented, C++ software, built on top of the robust open-source Cantera built mainly for homogeneous reactions. Reaction mechanisms can be input from human-readable files or automatic reaction generators, avoiding tedious work and errors. The governing equations are also built automatically, unlike Matlab and Python manual implementations, providing speed and error-free models. OpenMKM has built-in interfaces with numerical software, SUNDIALS, for solving ordinary differential equations and differential-algebraic equations. Users can choose various ideal reactors and energy balance options, such as isothermal, adiabatic, temperature ramp, and an experimentally measured temperature profile. OpenMKM is tightly integrated with pMuTT for thermochemistry input file generation from density functional theory (DFT), streamlining the workflow from DFT to MKM and eliminating tedious work and human errors. It is also seamlessly integrated with the RenView software for visualizing the reaction pathways and performing the reaction path or flux analysis (RPA). OpenMKM includes local sensitivity analysis (LSA) by solving the augmented system of equations or using the one-at-a-time finite difference (first or second order) method. LSA can identify not only kinetically influential reactions but also species. The software provides two techniques for large reaction mechanisms for which LSA is too expensive to run. One is the Fischer Information Matrix, which is approximate but comes at nearly zero cost. The other is a new method that we term RPA-guided LSA, which is a finite difference-based method but uses RPA to select kinetically relevant reactions instead of exploiting the entire reaction network. Users can quickly set up and conduct microkinetic simulations without writing code. The user inputs are conveniently divided into reactor setup files and thermodynamic and kinetic definition files to set up different reactors. The source code and documentation are openly available at https://github.com/VlachosGroup/openmkm.Item Predicting band gaps and band-edge positions of oxide perovskites using density functional theory and machine learning(Physical Review B, 2022-10-28) Li, Wei; Wang, Zigeng; Xiao, Xia; Zhang, Zhiqiang; Janotti, Anderson; Rajasekaran, Sanguthevar; Medasani, BharatDensity functional theory (DFT) within the local or semilocal density approximations, i.e., the local density approximation (LDA) or generalized gradient approximation (GGA), has become a workhorse in the electronic structure theory of solids, being extremely fast and reliable for energetics and structural properties, yet remaining highly inaccurate for predicting band gaps of semiconductors and insulators. The accurate prediction of band gaps using first-principles methods is time consuming, requiring hybrid functionals, quasiparticle GW, or quantum Monte Carlo methods. Efficiently correcting DFT-LDA/GGA band gaps and unveiling the main chemical and structural factors involved in this correction is desirable for discovering novel materials in high-throughput calculations. In this direction, we use DFT and machine learning techniques to correct band gaps and band-edge positions of a representative subset of ABO3 perovskite oxides. Relying on the results of HSE06 hybrid functional calculations as target values of band gaps, we find a systematic band-gap correction of ∼1.5 eV for this class of materials, where ∼1eV comes from downward shifting the valence band and ∼0.5eV from uplifting the conduction band. The main chemical and structural factors determining the band-gap correction are determined through a feature selection procedure.