We present the first calculation of an intrinsically relativistic quantity in fully non-linear cosmological large-scale structure studies. Traditionally, non-linear structure formation in standard ΛCDM cosmology is studied using N-body simulations, based on Newtonian gravitational dynamics on an expanding background. When one derives the Newtonian regime in a way that is a consistent approximation to the Einstein equations, a gravito-magnetic vector potential -giving rise to frame dragging -is present in the metric in addition to the usual Newtonian scalar potential. At leading order, this vector potential does not affect the matter dynamics, thus it can be computed from Newtonian N-body simulations. We explain how we compute the vector potential from simulations in ΛCDM and examine its magnitude relative to the scalar potential. We also discuss some possible observable effects.Introduction. Modern Cosmology is usually studied in two limits. On the largest scales, a perturbative approach is used in a general-relativistic framework. On small scales, non-linearities are treated in a Newtonian fashion, often with the use of N-body simulations.Few attempts have been made to go beyond the Newtonian approximation on non-linear scales by including post-Newtonian type corrections [1][2][3][4][5][6]. However, no attempt has been made to include post-Newtonian corrections in N-body simulations of cosmological large scale structure. Investigations have been carried out into the interpretation of N-body simulations on large scales, of the order of the Hubble length [7,8]. In [8], they go further and examine the dictionary between Newtonian and relativistic cosmologies on all scales and how accurately Newtonian cosmology satisfy the Einstein equations. Of course, no matter how well the Newtonian dynamics capture the full GR dynamics, there are GR quantities on all scales that have no counterpart in Newtonian theory.Recently, a new approximation scheme has been developed, dubbed the post-Friedmann approach [9], with the aim of providing a unified framework for all scales, from the fully non-linear Newtonian regime to the largest scales where relativistic effects become important [10]. It is based on an expansion in inverse powers of the speed of light, c, in a post-Newtonian [11] fashion, adapted to cosmology. When linearised, this approach correctly reproduces the linear general-relativistic perturbation theory. When one derives the Newtonian regime in this approach, in a way that is a consistent approximation to the Einstein equations, a vector potential must be present in the metric in addition to the usual Newtonian scalar gravitational potential.This vector potential is non-dynamical at leading order, therefore it does not affect the matter dynamics.