Si lo he entendido, no necesitas complicarte demasiado, simplemente debes desplazar el cálculo de la población:

    pop = 100000
    t = 5 # tiempo
    t_m = 0.008
    t_n = 0.01
    
    
    población = []
    nac = []
    mort = []
    
    for i in range(t + 1):
        población.append(pop)
        nacimientos = int(round(pop * t_n))
        muertes = int(round(pop * t_m))
        pop = pop + nacimientos - muertes
        mort.append(muertes)
        nac.append(nacimientos)
        
    print(población)
    print(nac)
    print(mort)

>     [100000, 100200, 100400, 100601, 100802, 101004]
>     [1000, 1002, 1004, 1006, 1008, 1010]
>     [800, 802, 803, 805, 806, 808]

El cálculo en este caso es:
    
        población[0] = pop
        poblacion[1] = poblacion[0] + nac[0] - mort[0]
        poblacion[2] = poblacion[1] + nac[1] - mort[1]
        ...
    
Otra opción sería que fuera:

        población[0] = pop
        poblacion[1] = poblacion[0] + nac[1] - mort[1]
        poblacion[2] = poblacion[1] + nac[2] - mort[2]
        ...

pero en este caso lo apropiado es que `mort[0]` y `nac[0]` sean ambos 0:

    
    pop = 100000
    t = 5 # tiempo
    t_m = 0.008
    t_n = 0.01
    
    
    población = [pop]
    nac = [0]
    mot = [0]
    
    for i in range(t):
        nacimientos = int(round(pop * t_n))
        muertes = int(round(pop * t_m))
        pop = pop + nacimientos - muertes
        población.append(pop)
        mot.append(muertes)
        nac.append(nacimientos)
        
    print(población)
    print(nac)
    print(mot)

>     [100000, 100200, 100400, 100601, 100802, 101004]
>     [0, 1000, 1002, 1004, 1006, 1008]
>     [0, 800, 802, 803, 805, 806]

Asumo que estás trabajando con una población, tasa de natalidad y tasa de mortalidad, por eso he redondeado al entero más cercano dado que en principio no existen individuos medio nacidos o medio muertos ni poblaciones con una fracción de un individuo :). Además resto `mort` en vez de sumarlo a la población.