Self-consistent simulations of beam-beam interaction in future <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:msup><mml:mi>e</mml:mi><mml:mo>+</mml:mo></mml:msup><mml:msup><mml:mi>e</mml:mi><mml:mo>−</mml:mo></mml:msup></mml:math> circular colliders including beamstrahlung and longitudinal coupling impedance

For the past generation ${e}^{+}{e}^{\ensuremath{-}}$ storage ring colliders, we usually used natural bunch length or its impedance lengthened value in beam-beam simulations instead of considering the impedance directly. In the future colliders, such as FCC-ee and CEPC, the beam-beam interaction becomes essentially three dimensional. In order to increase the luminosity, the future accelerators will collide very intense beams of high energy with low emittances and small beta functions at the collision points exploiting the crab waist collision scheme with a large Piwinski angle. For these extreme parameters several new effects become important for the collider performance such as beamstrahlung, coherent X-Z instability, 3D flip-flop so that the longitudinal beam dynamics should be also treated in a self-consistent manner. In this paper we describe the numerical code for the self-consistent 3D beam-beam simulations including beamstrahlung and the longitudinal beam coupling impedance and study interplay of different effects arising in beam-beam collisions of the future colliders.