Jithin John Varghese

Two-dimensional metal organic frameworks (2D MOFs) are a class of next generation materials for adsorbents,membranes and sensors; however, many 2D MOF synthesis methods lack the scalability and precision required for translation to industrial scales. Furthermore, the engineering of 2D MOF structures is challenged by complex environment-surface interactions that arise from their high anisotropy, thinness and functionally diverse surfaces. In this work we developed new understandings and methods of engineering such structures by using accelerated, high shear synthesis, and solvent exchange. With the recently developed annular flow microreactor we synthesized 2DMOFs more efficiently than conventional batch methods, by up to 5 orders of magnitude in terms of reactor space-time-yield. To accurately characterize particle size and dynamics in various organic solvents, we used liquid cell transmission electron microscopy. This technique not only visualized the oriented attachment of nanosheets, but also showed that the rate and direction of attachment is significantly influenced by solvent-surface interactions. These techniques and understandings provide rational bases for 2D MOF engineering and process design.

Using periodic Density Functional Theory calculations, propane oxidative dehydrogenation (ODH) and overoxidation over bare and hydroxylated Ni-doped CeO2 nanorods with predominantly exposed (1 1 0) facets were studied. Ab-initio thermodynamics-based surface phase analysis and computational Raman spectroscopic analysis predicted an 8.3 % surface oxygen vacancy concentration (Ce0.83Ni0.17O1.83) at typical ODH conditions. Only one-third of the surface oxygens, adjacent to the dopant, were selective for propene formation. Moreover, activated oxygen (O22–*) favored the formation of overoxidation products over propene. A monolayer hydroxyl coverage from water dissociation was stable at typical ODH conditions. This reduced the activation energy barrier for propene formation by 0.38 eV, increased the barrier for undesired acetone formation by 0.54 eV, and increased the barrier for propene activation by 0.6 eV. These promotional effects were due to the destabilization and induced hyperconjugation effects in the C3 adsorbates due to surface hydroxylation. Hence, surface hydroxylation (Lewis Base addition) is a potential strategy to improve propene selectivity.

The liquid-phase direct partial oxidation of methane to methanol is highly desirable for methane utilization. Using a combination of density functional theory and molecular dynamics simulations, this investigation reports the impact of sulfolane solvent on the reaction mechanism of methane oxidation with H2O2 in Fe-ZSM-5 catalyst and the potential product distribution. Sulfolane favored the formation of [Fe-OH] and [Fe =O] as active oxygen species from H2O2 dissociation, unlike Fe-OOH in water. On [Fe-OH], methanol formation was via a methoxy intermediate (activation barrier, Ea = 1.79 eV), while on [Fe O], it was by the radical rebound mechanism (Ea = 1.41 eV). The formation of formaldehyde was kinetically and thermodynamically less favorable than methanol and the formation of formic acid had a higher barrier. The lower diffusivity of H2O2 than CH3OH and the lower diffusivity of these in sulfolane than in water, together with the trends in reaction barriers, indicate higher methanol selectivity in sulfolane than in water.

DFT simulations predicted the C-O bond cleavage in glycerol on Brønsted acidic Pt-WOx catalyst under a hydrogen atmosphere to be via a “protonation dehydration” mechanism. Interfacial Pt-WOx sites facilitate Brønsted acid (BA) site formation by hydrogen spillover and synergistic activity of Pt and WOx. At strong BA sites, selective formation of 1,3-PDO is likely due to the large difference of over 40 kJ/mol in primary and secondary C-O bond cleavage barriers. The secondary C-O bond cleavage has a linear relation with NH3 binding energy while there is a non-monotonic trend for the primary C-O cleavage. Hence, Brønsted acid strength is a potential descriptor for the catalyst activity and 1,3-PDO selectivity. Synthetic strategies enabling fine dispersion of WOx to trimeric units on the Pt nanoparticles, while preventing W doping and agglomeration of WOx to 3D WO3 on small Pt nanoparticles are desirable for high BA strength and thereby 1,3-PDO yield.

Due to their high anisotropy and tunable chemical composition, two-dimensional metal organic frameworks (2D MOFs) have great potential as building blocks for next-generation materials in a diverse range of applications - from electrochemical catalysis to membrane separation. However, the controllable synthesis is complicated by the environment-surface interactions that arise from the high anisotropy, thinness, and functionally diverse surfaces of 2D MOFs. Liquid cell transmission electron microscopy (LCTEM) offers a unique opportunity to study these interactions in situ. In this work, we analyzed the effects of different solvent environments on the structure and aggregation dynamics of copper benzene dicarboxylic acid (CuBDC) nanosheets, which were synthesized using a high shear annular microreactor. LCTEM revealed that 2D MOF nanosheets undergo oriented attachment and that the rate and direction of oriented attachment is controlled by solvent-surface interactions. We investigated the nature of these solvent interactions using density functional theory calculations, which suggest that the binding energy of solvents to different MOF surfaces is likely responsible for this behavior. The CuBDC nanosheets were then applied as adsorbents in organic solvents, in which we showed how solvent-mediated oriented attachment could significantly affect adsorption properties.

The aim of this work is to develop insights into adsorption of hydrogen and carbon dioxide in a zeolitic imidazolate framework, ZIF-8, using high-pressure adsorption studies, adsorption isotherm model fitting, and DFT investigation of preferential adsorption sites and binding energies. The robustness of ZIF series metal–organic frameworks has drawn interest towards its utility in large scale applications in gas storage and separation. We use room temperature synthesis of ZIF-8 using DMF as a solvent, and benchmarked it against typical solvothermal synthesis. The resulting material is characterized using XRD, SEM, TG–DSC and N2 adsorption isotherm. High-pressure volumetric adsorption of the activated materials is conducted to analyze the hydrogen and carbon dioxide storage capacities up to 50 and 40 bar, respectively. ZIF-8 shows maximum H2 storage capacity of 3.13 wt% at 50 bar and 77 K, and CO2 storage capacity of 46 wt% at 40 bar and 300 K. The parameters of Unilan adsorption isotherm are estimated from the equilibrium adsorption data and isosteric heats of adsorption for H2 and CO2 on ZIF-8 are computed. DFT calculations are used to obtain preferential adsorption sites of H2 and CO2. Adsorption enthalpy values were calculated from DFT as − 7.08 and − 25.98 kJ/mol, respectively for H2 and CO2 at the most preferred sites. We found a close agreement between isosteric heat of adsorption of hydrogen (− 4.68 kJ/mol) and the enthalpy of hydrogen adsorption from DFT (− 6.04 kJ/mol) at 77 K.

Modern heterogeneous catalysis has benefitted immensely from computational predictions of catalyst structure and its evolution under reaction conditions, first-principles mechanistic investigations, and detailed kinetic modeling, which are rungs on a multiscale workflow. Establishing connections across these rungs and integration with experiments have been challenging. Here, operando catalyst structure prediction techniques using density functional theory simulations and ab initio thermodynamics calculations, molecular dynamics, and machine learning techniques are presented. Surface structure characterization by computational spectroscopic and machine learning techniques is then discussed. Hierarchical approaches in kinetic parameter estimation involving semi-empirical, data-driven, and first-principles calculations and detailed kinetic modeling via mean-field microkinetic modeling and kinetic Monte Carlo simulations are discussed along with methods and the need for uncertainty quantification. With these as the background, this article proposes a bottom-up hierarchical and closed loop modeling framework incorporating consistency checks and iterative refinements at each level and across levels.

Cu/ZnO/ZrO2/Al2O3 catalysts are widely explored for CO2 conversion to methanol due to their higher activity and stability. However, mechanistic understanding of the performance of such catalysts is lacking due to ambiguity on the actual active sites. This study focuses on unraveling the nature of different interfaces on Cu/ZnO/ZrO2/Al2O3 catalyst by coupling experiments, Density Functional Theory (DFT) simulations and a DFT-based reactor scale multi-site microkinetic model. Although DFT calculations suggested the ZrO2/Cu interface to be the CO2 adsorption site, the validated microkinetic model predicted the ZnO/Cu interface to be the crucial reaction center. Reaction pathway analysis showed that methanol is produced through the formate pathway near the reactor entrance, whereas, the carboxyl pathway dominates in the latter zones, emphasizing the occurrence of both CO2 and CO hydrogenation. This deeper understanding of the reaction behavior of such multicomponent catalysts will aid in designing better catalysts and optimizing reaction conditions and systems.

Single crystal preparation of two 1,4-dihydropyridine (1,4-DHP) derivative molecules, namely, (i) 2,6-dimethyl-4-phenyl-1, 4-dihydro-pyridine-3, 5-dicarboxylic acid diethyl ester and (ii) 4-(4‑methoxy-phenyl)-2,6-dimethyl-1, 4-dihydro-pyridine-3, 5-dicarboxylic acid diethyl ester was done. Crystal structure solution and refinement of final R values converged to 0.0478 and 0.0466. Hirshfeld surface analysis mapped dnorm surfaces, shape index, curvedness and hydrogen bonding present in the molecules. These calculations established the relative contributions of the various intermolecular contacts to the total Hirshfeld surface area. Both the molecular structures were optimized and confirmed by Density Functional Theory (DFT) simulations in an implicit aqueous solvation environment. The molecular electrostatic potential (MEP) analysis elucidated the electrophilic and nucleophilic regions in the molecules to be the hydrogen attached to the nitrogen and the carboxyl oxygen atoms. Frontier molecular orbital analysis indicated the stability of these molecules while they appeared to be electrophilic in nature from the reactivity descriptors analysis.

Propane oxidative dehydrogenation (ODH) was studied over VOx/CeO2 nanoparticle catalysts for various vanadia nuclearities using Density Functional Theory simulations. On monomers and dimers, propene formation involved V=O, bridged (V–O–Ce, V–OV) and ceria surface oxygens but with low selectivity, while on trimers and possibly higher oligomers it was kinetically favourable over ceria surface oxygens, with higher selectivity than on monomers and dimers. On trimer, overoxidation via acetone formation was inhibited, but at the expense of overall catalyst reactivity. Doping of ceria support with Ni induced strong Ni–VOx interactions, stabilizing the vanadia monomer. The V–O–Ni bridged oxygen was less accessible for overoxidation, resulting in an increase in activation energy barrier for acetone formation by 0.69 eV on VO2/Ce0.89Ni0.11O2(1 1 1), while there was a minimal impact on V3O6/Ce0.89Ni0.11O2 (1 1 1). Hence, transition metal doping of ceria support is a potential strategy to improve propene selectivity over VOx/CeO2 catalysts for a wide range of V loadings.