Many liquids have curves (isomorphs) in their phase diagrams along which structure, dynamics, and some thermodynamic quantities are invariant in reduced units. A substantial part of their phase diagrams is thus effectively one dimensional. The shape of these isomorphs is described by a material-dependent function of density h(ρ), which for real liquids is well approximated by a power law ρ γ . However, in simulations, a power law is not adequate when density changes are large; typical models such as the Lennard-Jones liquid show that γ(ρ) ≡ d ln h(ρ)/ d ln ρ is a decreasing function of density. This paper presents results from computer simulations using a new pair potential that diverges at a nonzero distance and can be tuned to give a more realistic shape of γ(ρ). Our results indicate that the finite size of molecules is an important factor to take into account when modeling liquids over a large density range.It is still an open question what controls the dynamics of viscous, glass-forming liquids [1][2][3][4][5][6]. Although the dynamics in general depends on both temperature T and density ρ, the dynamics of many organic supercooled liquids can be collapsed onto a single curve when plotted against a combined, material-specific variable h(ρ)/T [7][8][9]. It was found in many experiments that the scaling function h(ρ) is generally well approximated by a power law as h(ρ) = ρ γ , with γ being the materialspecific density-scaling exponent [10,11]; we refer to this as power-law density scaling. Another important development was the discovery that the dynamics of liquids are a function of the excess entropy [12,13].The isomorph theory [14] explains why both density scaling and excess-entropy scaling work for some liquids. Liquids that obey the isomorph theory have curves in their phase diagram, so-called isomorphs, along which not only the dynamics, but also the structure, excess entropy and other thermodynamic quantities are invariant. The development of the isomorph theory was initiated by the observation that in computer simulations some liquids have strongly correlated fluctuations in their energy and pressure. More specifically, if the energy E and pressure p are split in a kinetic part and a configurational part that only depends on the particle positions R ≡ (r 1 , . . . , r N ), as followsthe strong correlations are found between the thermal equilibrium fluctuations of the potential energy U and the virial W in the N V T ensemble [15], although strong correlations have also been found at high pressures in * a.a.veldhorst@gmail.com † tbs@ruc.dk ‡ dyre@ruc.dk the N pT ensemble [16]. Indeed, the standard correlation coefficientindicates whether a liquid obeys the isomorph theory: this is the case whenever R > 0.9 (although this value is of course somewhat arbitrary). The standard linear regression "slope" of the fluctuationsis the density-scaling exponent [14], and the theory thus provides a convenient way to determine the densityscaling exponent in computer simulations. Another empirical observation that can b...