function animate(sol::Solution, args, name_part::String; framerate::Int64=10)
    println("Generating animation...")

    set_theme!(theme_black())

    fig = Figure()
    ax = Axis(
        fig[1, 1];
        limits=(-args.l, args.l, -args.l, args.l),
        aspect=AxisAspect(1),
        xlabel=L"x",
        ylabel=L"y",
    )

    Colorbar(fig[1, 2]; limits=(0, 2), colormap=ColorSchemes.rainbow, label=L"\frac{\varphi}{\pi}")

    animation_path = "exports/$name_part.mkv"

    record(fig, animation_path; framerate=framerate) do io
        segments_x = zeros(2 * args.N)
        segments_y = zeros(2 * args.N)

        circles = Vector{Circle}(undef, args.N)
        interaction_circles = Vector{Circle}(undef, args.N)
        skin_circles = Vector{Circle}(undef, args.N)

        colors = Vector{RGBAf}(undef, args.N)
        interaction_colors = Vector{RGBAf}(undef, args.N)
        skin_colors = Vector{RGBAf}(undef, args.N)

        @showprogress 0.5 for frame in 1:args.n_frames
            if frame > 1
                empty!(ax)

                @simd for i in 1:(args.N)
                    segments_x[2 * i - 1] = sol.center[i, frame - 1][1]
                    segments_x[2 * i] = sol.center[i, frame][1]

                    segments_y[2 * i - 1] = sol.center[i, frame - 1][2]
                    segments_y[2 * i] = sol.center[i, frame][2]
                end

                linesegments!(ax, segments_x, segments_y; color=1:(args.N))
            else
                println("Started recording!")
            end

            @simd for i in 1:args.N
                circles[i] = Circle(Point2(sol.center[i, frame]), args.particle_diameter / 2)
                interaction_circles[i] = Circle(Point2(sol.center[i, frame]), args.interaction_r)
                skin_circles[i] = Circle(Point2(sol.center[i, frame]), args.skin_r)

                color = get(ColorSchemes.rainbow, rem2pi(sol.φ[i, frame] / (2 * π), RoundDown))

                colors[i] = RGBAf(color)
                interaction_colors[i] = RGBAf(color, 0.2)
                skin_colors[i] = RGBAf(color, 0.05)
            end

            poly!(ax, circles; color=colors)
            poly!(ax, interaction_circles; color=interaction_colors)
            poly!(ax, skin_circles; color=skin_colors)

            ax.title = "t = $(round(sol.t[frame], digits=3))"

            recordframe!(io)
        end
    end

    println("Animation done and saved to $animation_path.")

    return nothing
end