In this Letter we investigate the rupture instability of thin liquid films by means of a bifurcation analysis in the vicinity of the short-scale instability threshold. The rupture time estimate obtained in closed form as a function of the relevant dimensionless groups is in striking agreement with the results of the numerical simulations of the original nonlinear evolution equations. This suggests that the weakly nonlinear theory captures the adequate physics of the instability. When antagonistic (attractive/repulsive) molecular forces are considered, nonlinear saturation of the instability becomes possible. We show that the stability boundaries are determined by the van der Waals potential alone.PACS numbers: 47.20.-k, 68.15+e It is well known that a liquid film on a planar solid surface may become unstable due to long-range molecular forces. The forces originating from van der Waals attractions [1] accelerate thinning in regions of film depression leading to film rupture and "spinodal dewetting" [2]. On the other hand, electrical double layers on the solid surface may give rise to intermolecular repulsions stabilizing thin films against rupture [3].In recent years, much effort has been put into theoretical modelling the dewetting phenomena [4,5,6,7,8,9,11,14,15]. A nonlinear theory of the film evolution based on the long-wave nature of the response was first posed in Ref. 4. This approach, which has already been considered for different situations [5], yields nonlinear partial differential equations that describes the evolution of the interface shape, surfactant concentration, etc. Linear stability analysis is routinely applied to predict the onset of the instability and the characteristic wavelength, but the rupture time estimate obtained from the linear theory turns out to be rather poor: it underestimates the rupture time due to highly nonlinear nature of response. The most common and straightforward approach is to solve the evolution equations numerically [6], [7] [11], [13,14,15,16]. The obvious disadvantage of the numerical simulation is that for a complex problem that involves many parameters, full parametric study of the rupture is quite elaborate.A bifurcation technique was first applied in [8] to arrive at the nonlinear estimate for the rupture time in the vicinity of a steady bifurcation point. It was demonstrated that nonlinear terms owing to van der Waals attractions contribute to rapid acceleration of the rupture beyond the linear regime. Analysis of the nonlinear evolution of small disturbances leads to a dynamic Landau equation for the perturbation amplitude. The closed form solution of the amplitude equation provides a time for "blowup" of the initial small-amplitude disturbance that was proposed to be a good estimate of the nonlinear rupture time. The approach has never been given enough attention perhaps because the analysis involves rather tedious algebra and can only be done "by hand" for some simple cases. It has been demonstrated in [9] that the derivation of the amplitude equation ca...