Most applications of Bayesian Inference for parameter estimation and model selection in astrophysics involve the use of Markov Chain Monte Carlo (MCMC) techniques. In this work, we introduce Variational Inference as an alternative to solve these problems, and compare how the results hold up to MCMC methods. Variational Inference converts the inference problem into an optimization problem by approximating the posterior from a known family of distributions and using Kullback-Leibler divergence to measure closeness. Variational Inference takes advantage of fast optimization techniques, which make it ideal to deal with large datasets and also makes it trivial to parallelize. As a proof of principle, we apply Variational Inference to four different problems in astrophysics, where MCMC techniques were previously used. These include measuring exoplanet orbital parameters from radial velocity data, tests of periodicities in measurements of Newton's constant $G$, assessing the significance of a turnover in the spectral lag data of GRB 160625B , and estimating the mass of a galaxy cluster using weak gravitational lensing. We find that Variational Inference is much faster than MCMC for these problems.

Gamma-Ray Bursts (GRBs) have been conventionally bifurcated into two distinct categories: ``short'' and ``long'' with durations less than and greater than two seconds respectively. However, there is a lot of literature (although with conflicting results) regarding the existence of a third intermediate class. To investigate this issue, we extend a recent study (arXiv:1612.08235) on classification of GRBs to two dimensions by incorporating the GRB hardness in addition to the observed durations. We carry out this unified analysis on GRB datasets from four detectors, viz. BATSE, RHESSI, Swift (observed and intrinsic frame), and Fermi-GBM. We consider the duration and hardness features in log-scale for each of these datasets and determine the best-fit parameters using Gaussian Mixture Model. This is followed by information theoretic criterion (AIC and BIC) to determine if a three-component fit is favored compared to a two-component one or vice-versa. For BATSE, we find that both AIC and BIC show preference for three components with decisive significance. For Fermi and RHESSI, both AIC and BIC show preference for two components, although the significance is marginal from AIC, but decisive using BIC. For Swift dataset in both the observed and rest frame, we find that three components are favored according to AIC with decisive significance, and two are preferred with BIC with marginal to strong significance.