In this paper, we develop a Galerkin-type approximation, with quantitative error estimates, for weak solutions to the Cauchy problem for kinetic Fokker-Planck equations in the domain (0,T)×D×Rd, where D is either Td or Rd. Our approach is based on a Hermite expansion in the velocity variable only, with a hyperbolic system that appears as the truncation of the Brinkman hierarchy, as well as ideas from $\href{arXiv:1902.04037v2}{Alb+21}$ and additional energy-type estimates that we have developed. We also establish the regularity of the solution based on the regularity of the initial data and the source term.

In this paper, we investigate weak solutions and Perron-Wiener-Brelot solutions to the linear stationary Kramers-Fokker-Planck equation in bounded domains. We establish the existence of weak solutions in product domains by applying the Lions-Lax-Milgram theorem and the vanishing viscosity method. Furthermore, we show that these solutions coincide in well-behaved domains. Building on the existence of weak solutions in product domains, we develop the foundational theory of Perron-Wiener-Brelot solutions in arbitrary bounded domains. Our results rely on recent advancements in the theory of kinetic Fokker-Planck equations with rough coefficients.

We establish the boundedness of weak subsolutions for a class of degenerate Kolmogorov equations of the hypoelliptic type, compatible with a homogeneous Lie group structure, within bounded product domains using the De Giorgi iteration. We employ the renormalization formula to handle boundary values and provide energy estimates. An L¹ - L^p type embedding estimate derived from the fundamental solution is utilized to incorporate lower-order divergence terms. This work naturally extends the boundedness theory for uniformly parabolic equations, with matching exponents for the coefficients.

We introduce a variational formulation for a general class of possibly degenerate, non-self-adjoint Fokker-Planck operators in divergence form, motivated by the work of Albritton et al. (2024), and prove that it is suitable for defining the variational capacity.  Using this framework, we establish rough estimates for the equilibrium potential in the elliptic case, providing a novel approach compared to previous methods.  Finally, we derive the Eyring-Kramers formula for non-self-adjoint elliptic Fokker-Planck operators in divergence form, extending the results of Landim et al. (2019) and Lee & Seo (2022).