Energy-based model explained

An energy-based model (EBM) (also called a Canonical Ensemble Learning(CEL) or Learning via Canonical Ensemble (LCE)) is an application of canonical ensemble formulation of statistical physics for learning from data problems. The approach prominently appears in generative models (GMs).

EBMs provide a unified framework for many probabilistic and non-probabilistic approaches to such learning, particularly for training graphical and other structured models.

An EBM learns the characteristics of a target dataset and generates a similar but larger dataset. EBMs detect the latent variables of a dataset and generate new datasets with a similar distribution.

Energy-based generative neural networks ^[1] ^[2] is a class of generative models, which aim to learn explicit probability distributions of data in the form of energy-based models whose energy functions are parameterized by modern deep neural networks.

Boltzmann machines are a special form of energy-based models with a specific parametrization of the energy.^[3]

Description

For a given input

, the model describes an energy

E_\theta(x)

such that the Boltzmann distribution

P_{\theta(x)=\exp(-\beta}E_{\theta(x))/Z(\theta)}

is a probability (density) and typically

\beta=1

Since the normalization constant

Z(\theta):=\int_xdx\exp(-\betaE_\theta(x))

, also known as partition function, depends on all the Boltzmann factors of all possible inputs

it cannot be easily computed or reliably estimated during training simply using standard maximum likelihood estimation.

However for maximizing the likelihood during training, the gradient of the log likelihood of a single training example

is given by using the chain rule

\partial_\thetalog\left(P_{\theta(x)\right)=E}


	x'\simP_\theta

[\partial_\thetaE_{\theta(x')]-\partial}_\thetaE_\theta(x)(*)

The expectation in the above formula for the gradient can be approximately estimated by drawing samples

from the distribution

P_\theta

using Markov chain Monte Carlo (MCMC)^[4]

Early energy-based models like the 2003 Boltzmann machine by Hinton estimated this expectation using block Gibbs sampler. Newer approaches make use of more efficient Stochastic Gradient Langevin Dynamics (LD) drawing samples using:^[5]

x_0'\simP_0,x_i+1'=x_i'-

	\alpha
	2

	\partialE_\theta(x_i')
	\partialx_i'

+\epsilon,

and

\epsilon\siml{N}(0,\alpha)

. A replay buffer of past values

x_i'

is used with LD to initialize the optimization module.

The parameters

\theta

of the neural network are, therefore, trained in a generative manner by MCMC-based maximum likelihood estimation:^[6] The learning process follows an "analysis by synthesis" scheme, where within each learning iteration, the algorithm samples the synthesized examples from the current model by a gradient-based MCMC method, e.g., Langevin dynamics or Hybrid Monte Carlo, and then updates the model parameters

\theta

based on the difference between the training examples and the synthesized ones, see equation

(*)

. This process can be interpreted as an alternating mode seeking and mode shifting process, and also has an adversarial interpretation.^[7]

In the end, the model learns a function

E_\theta

that associates low energies to correct values, and higher energies to incorrect values.

After training, given a converged energy model

E_\theta

, the Metropolis–Hastings algorithm can be used to draw new samples.The acceptance probability is given by:

P_acc(x_i\tox^{*)=min\left(1,}

	*)
P
	\theta(x

P_\theta(x_i)

\right).

History

The term "energy-based models" was first coined in a 2003 JMLR paper^[8] where the authors defined a generalisation of independent components analysis to the overcomplete setting using EBMs.Other early work on EBMs proposed models that represented energy as a composition of latent and observable variables.

Characteristics

EBMs demonstrate useful properties:

Simplicity and stability–The EBM is the only object that needs to be designed and trained. Separate networks need not be trained to ensure balance.
Adaptive computation time–An EBM can generate sharp, diverse samples or (more quickly) coarse, less diverse samples. Given infinite time, this procedure produces true samples.
Flexibility–In Variational Autoencoders (VAE) and flow-based models, the generator learns a map from a continuous space to a (possibly) discontinuous space containing different data modes. EBMs can learn to assign low energies to disjoint regions (multiple modes).
Adaptive generation–EBM generators are implicitly defined by the probability distribution, and automatically adapt as the distribution changes (without training), allowing EBMs to address domains where generator training is impractical, as well as minimizing mode collapse and avoiding spurious modes from out-of-distribution samples.
Compositionality–Individual models are unnormalized probability distributions, allowing models to be combined through product of experts or other hierarchical techniques.

Experimental results

On image datasets such as CIFAR-10 and ImageNet 32x32, an EBM model generated high-quality images relatively quickly. It supported combining features learned from one type of image for generating other types of images. It was able to generalize using out-of-distribution datasets, outperforming flow-based and autoregressive models. EBM was relatively resistant to adversarial perturbations, behaving better than models explicitly trained against them with training for classification.

Applications

Target applications include natural language processing, robotics and computer vision.

The first energy-based generative neural network is the generative ConvNet proposed in 2016 for image patterns, where the neural network is a convolutional neural network.^[9] ^[10] The model has been generalized to various domains to learn distributions of videos,^[11] and 3D voxels.^[12] They are made more effective in their variants.^[13] ^[14] ^[15] ^[16] ^[17] ^[18] They have proven useful for data generation (e.g., image synthesis, video synthesis, 3D shape synthesis, etc.), data recovery (e.g., recovering videos with missing pixels or image frames, 3D super-resolution, etc), data reconstruction (e.g., image reconstruction and linear interpolation).

Alternatives

EBMs compete with techniques such as variational autoencoders (VAEs), generative adversarial networks (GANs) or normalizing flows.

Extensions

Joint energy-based models

Joint energy-based models (JEM), proposed in 2020 by Grathwohl et al., allow any classifier with softmax output to be interpreted as energy-based model. The key observation is that such a classifier is trained to predict the conditional probability

p_\theta(y|x)=

	\vec{f
e
	\theta(x)[y]

} \ \ \text y = 1, \dotsc, K \text \vec_\theta = (f_1, \dotsc, f_K) \in \R^K,where

\vec{f}_\theta(x)[y]

is the y-th index of the logits

\vec{f}

corresponding to class y.Without any change to the logits it was proposed to reinterpret the logits to describe a joint probability density:

\theta(y,x)=

	\vec{f
e
	\theta(x)[y]

},with unknown partition function

Z(\theta)

and energy

E_\theta(x,y)=-f_\theta(x)[y]

.By marginalization, we obtain the unnormalized density

p_{\theta(x)=\sum}_yp_\theta(y,x)=\sum_y

	\vec{f
e
	\theta(x)[y]

}=:\exp(-E_\theta(x)),therefore,

E_{\theta(x)=-log\left(\sum}_y

	\vec{f
e
	\theta(x)[y]

}\right),so that any classifier can be used to define an energy function