Energy and energy gradient matrix elements with N-particle explicitly correlated complex Gaussian basis functions with L=1.