On $p$-adic multiple Barnes-Euler zeta functions and the corresponding log gamma functions

Suppose that $\omega_1,\ldots,\omega_N$ are positive real numbers and $x$ is a complex number with positive real part. The multiple Barnes-Euler zeta function $\zeta_{E,N}(s,x;\bar\omega)$ with parameter vector $\bar\omega=(\omega_1,\ldots,\omega_N)$ is defined as a deformation of the Barnes multiple zeta function as follows $$\zeta_{E,N}(s,x;\bar\omega)=\sum_{t_1=0}^\infty\cdots\sum_{t_N=0}^\infty \frac{(-1)^{t_1+\cdots+t_N}}{(x+\omega_1t_1+\cdots+\omega_Nt_N)^s}.$$ In this paper, based on the fermionic $p$-adic integral, we define the $p$-adic analogue of multiple Barnes-Euler zeta function $\zeta_{E,N}(s,x;\bar\omega)$ which we denote by $\zeta_{p,E,N}(s,x;\bar\omega).$ We prove several properties of $\zeta_{p,E,N}(s,x; \bar\omega)$, including the convergent Laurent series expansion, the distribution formula, the difference equation, the reflection functional equation and the derivative formula. By computing the values of this kind of $p$-adic zeta function at nonpositive integers, we show that it interpolates the higher order Euler polynomials $E_{N,n}(x;\bar\omega)$ $p$-adically. Furthermore, we define the corresponding multiple $p$-adic Diamond-Euler Log Gamma function. We also show that the multiple $p$-adic Diamond-Euler Log Gamma function ${\rm Log}\, \Gamma_{\! D,E,N}(x;\bar\omega)$ has an integral representation by the multiple fermionic $p$-adic integral, and it satisfies the distribution formula, the difference equation, the reflection functional equation, the derivative formula and also the Stirling's series expansions.