Galoisian approach to complex oscillation theory of some Hill equations

We apply Kovacic's algorithm from differential Galois theory to show that all complex non-oscillatory solutions (finite exponent of convergence of zeros) of certain Hill equations considered by Bank and Laine using Nevanlinna theory must be Liouvillian solutions. These are solutions obtainable by suitable differential field extension constructions. In particular, we establish a full correspondence between solutions of non-oscillatory type and Liouvillian solutions for a particular Hill equation. Explicit closed-form solutions are obtained via both methods for this Hill equation whose potential is a combination of four exponential functions in the Bank-Laine theory. The differential equation is a periodic form of a biconfluent Heun equation. We further show that these Liouvillian solutions exhibit novel single and double orthogonality, and satisfy Fredholm integral equations over suitable integration regions in $\mathbb{C}$ that mimic single/double orthogonality for the corresponding Liouvillian solutions of the Lamé and Whittaker-Hill equations, discovered by Whittaker and Ince almost a century ago.