Catalysis of dynamical symmetry breaking by a constant magnetic field in (3 + 1) dimensions is considered. We use the three flavour NambuJona-Lasinio type model with 't Hooft and eight-quark interaction terms. It is shown that the multi-quark interactions introduce new additional features to this phenomenon: (a) the local minimum of the effective potential catalyzed by the constant magnetic field is smoothed out with increasing strength of the field at the characteristic scale H ∼ 10 19 G, (b) the multi-quark forces generate independently another local minimum associated with a larger dynamical fermion mass. This state may exist even for multi-quark interactions with a subcritical set of couplings, and is globally stable with respect to a further increase of the magnetic field. © 2007 Elsevier B.V. All rights reserved. PACS: 11.30.Rd; 11.30.Qc It has been shown in a series of papers [1][2][3] that in (2 + 1) and (3 + 1) dimensions a constant magnetic field H = 0 catalyzes the dynamical symmetry breaking leading to a fermion mass even at the weakest attractive four-fermion interaction between particles, and the symmetry is not restored at any arbitrarily large H . Soon thereafter it became also clear [4][5][6] that the zero-energy surface of the lowest Landau level (LLL) plays a crucial role in the dynamics of such fermion pairing. It has been found that the dynamics of the fermion pairing in the homogeneous magnetic field is essentially (1 + 1)-dimensional, and a deep analogy of this phenomenon with the dynamics of electron pairing in BCS [7] has been stressed. The generated fermion mass, M dyn , turned out to be much smaller than the Landau gap ∼ √ |eH |. The existence of a zero-energy surface in the spectrum of a Dirac particle is ensured for any homogeneous magnetic field with a fixed direction by a quantum mechanical supersymme-* Corresponding author.E-mail address: osipov@nu.jinr.ru (A.A. Osipov).try of the corresponding second-order Dirac Hamiltonian [8]. This aspect of the phenomenon appears to be a quite exceptional situation and indicates that the dynamical generation of mass is not so universal as one would expect by extrapolating the results obtained for homogeneous or unidirectional [9] magnetic field profiles. For instance, it has been demonstrated by Ragazzon [10] that the Nambu-Jona-Lasinio (NJL) model [11] minimally coupled to a background magnetic field with variable direction does not possess a massive phase until the coupling constant exceeds some critical value. Obviously, in this case one faces the conventional scenario of dynamical chiral symmetry breaking, where the magnetic field does not play an essential role. Conversely, having in mind that homogeneous magnetic fields can act as strong catalysts of chiral symmetry breaking, one might ask what is the effect caused by the strong interaction, when higher order multi-fermion interactions are present. These extensions of the NJL model are well known [12][13][14], for instance, the four-quark U(3) L × U(3) R chiral symmetric Lagr...