We have constructed a class of perturbative dynamical black hole solutions in presence of cosmological constant. We have done our calculation in large number of dimensions. The inverse power of dimension has been used as the perturbation parameter and our calculation is valid upto the first subleading order. The solutions are in one to one correspondence with a dynamical membrane and a velocity field embedded in the asymptotic geometry. Our method is manifestly covariant with respect to the asymptotic geometry. One single calculation and the same universal result works for both dS and AdS geometry or in case of AdS for both global AdS and Poincare patch. We have checked our final answer with various known exact solutions and the known spectrum of Quasi Normal modes in AdS/dS.

We have extended the results of [1] upto second subleading order in an expansion around large dimension D. Unlike the previous case, there are non-trivial metric corrections at this order. Due to our 'background-covariant' formalism, the dependence on Ricci and the Riemann curvature tensor of the background is manifest here. The gravity system is dual to a dynamical membrane coupled with a velocity field. The dual membrane is embedded in some smooth background geometry that also satisfies the Einstein equation in presence of cosmological constant. We explicitly computed the corrections to the equation governing the membrane-dynamics. Our results match with earlier derivations in appropriate limits. We calculated the spectrum of QNM from our membrane equations and matched them against similar results derived from gravity.