The equation and solution of 4-point correlation function of galaxies in Gaussian approximation and its parity-odd part

Starting with the density field equation of a self-gravity fluid in a static Universe, using the Schwinger functional differentiation technique, we derive the field equation of the 4-point correlation function (4PCF) of galaxies in the Gaussian approximation, which contains hierarchically 2PCF and 3PCF. By use of the known solutions of 2PCF and 3PCF, the equation of 4PCF becomes an inhomogeneous, Helmholtz equation, and contains only two physical parameters: the mass m of galaxy and the Jeans wavenumber kJ, like the equations of the 2PCF and 3PCF. We obtain the analytical solution of 4PCF that consists of four portions, η= η⁰odd + η⁰ₑᵥₑₙ +ηFP +ηI, and has a very rich structure. η⁰odd and η⁰ₑᵥₑₙ form the homogeneous solution and depend on boundary conditions. The parity-odd η⁰odd is more interesting and qualitatively explains the observed parity-odd data of BOSS CMASS, the parity-even η⁰ₑᵥₑₙ contains the disconnected 4PCF ηᵈⁱˢᶜ (arising from a Gaussian random process), and both η⁰odd and η⁰ₑᵥₑₙ are prominent at large scales r 10Mpc, and exhibit radial oscillations determined by the Jeans wavenumber. ηFP and ηI are parity-even, and form the inhomogeneous solution. ηFP is the same as the Fry-Peebles ansatz for 4PCF, and dominates at small scales r 10Mpc. ηI is an integration of the inhomogeneous term, subdominant. We also compare the parity-even 4PCF with the observation data.