Electron and proton acceleration in three-dimensional electric and magnetic fields is studied through test particle simulations. The fields are obtained by a three-dimensional magnetohydrodynamic simulation of magnetic reconnection in slab geometry. The nonlinear evolution of the system is characterized by the growth of many unstable modes and the initial current sheet is fragmented with formation of small scale structures. We inject at random points inside the evolving current sheet a Maxwellian distribution of particles. In relatively short time (less than a millisecond) the particles develop a power law tail. The acceleration is extremely efficient and the electrons absorb a large percentage of the available energy in a small fraction of the characteristic time of the MHD simulation, suggesting that resistive MHD codes, used extensively in the current literature, are unable to represent the full extent of particle acceleration in 3D reconnection.It is widely accepted that magnetic reconnection plays a significant role in converting magnetic energy to thermal energy and kinetic energy of electrons and protons in laboratory plasmas, the Earth's magnetosphere, the solar corona, and in extragalactic jets [1,2].In resistive magnetohydrodynamic models, resistivity breaks the frozen-in law in a boundary layer, allowing reconnection to occur. A current sheet can be spontaneously unstable to resistive instabilities, like the tearing modes, which lead to magnetic reconnection [3,4]. Many numerical codes have been developed to study the nonlinear evolution of tearing modes in two-dimensional approximations [5]. However, three-dimensional effects may become important in modifying the spatial structure of the current sheets and the reconnection rate [6,7]. Different numerical studies have been performed to investigate collisionless magnetic reconnection using fluid models, where magnetic reconnection is made possible by electron inertia, and kinetic simulations, in two and three dimensional configurations [8,9]. It has been shown that in the three-dimensional kinetic reconnection the characteristic time scale of the instability is much faster than that of the two-dimensional tearing mode instability.The change of the topology of the magnetic field due to magnetic reconnection allows the release of magnetic energy, which can be responsible for the acceleration of particles. In two-dimensional reconnection configurations particle acceleration has been extensively studied both analytically and numerically [10,11]. The acceleration is caused by the motion of particles along the electric field in the current sheet, but the magnetic field plays a significant role since it influences the trajectory and therefore the energy gain of the particles. Recently, it has become clear that it is essential to include in the model the longitudinal component of the magnetic field, which is parallel to the electric field in the current sheet [12]- [14]. Studies of particle acceleration with a longitudinal magnetic field component have also be...