Here P ij (t) represent the off-diagonal components of the stress-tensor (i.e. P xy , P xz , and P yz ), . . . denotes ensemble averaging, and V represents the volume. The Dbased method was shown to compare favourably with EMD approaches, since due to the slow relaxation in highly viscous fluids, and due to large fluctuations in the components of the stress tensor, very long MD simulations are required to provide sufficient sampling. Alternatively, nonequilibrium MD (NEMD) approaches employ the response of the system to an external shear field to determine the viscosity . NEMD approaches have been shown to be quite efficient, when using large external fields. There is, however, always the inherent problem, that the determined viscosities might depend on the shear rate. In this contribution, we are testing the "OrthoBoXY"approach for a variety of systems, such as molecular fluids like liquid water, dimethyl ether, methanol, and triglyme, as well as methanol/water and triglyme/water liquid mixtures, and imidazolium-based ionic liquids, thus probing viscosities varying over almost three orders of magnitude. Our aim is to assess the efficiency of the method for computing selfdiffusion and viscosity data, and to provide simple formulae for estimating the MD simulation lengths and sizes, and to offer guidelines for delivering reliable data with targeted uncertainties for both, the system size independent self-diffusion coefficient D 0 , and the shear viscosity η.