Abstract. -The nonlinear dynamics of an elastic filament that is forced to rotate at its base is studied by hydrodynamic simulation techniques; coupling between stretch, bend, twist elasticity and thermal fluctuations is included. The twirling-overwhirling transition is located and found to be strongly discontinuous. For finite bend and twist persistence length, thermal fluctuations lower the threshold rotational frequency, for infinite persistence length the threshold agrees with previous analytical predictions.The dynamics and morphology of elastic filaments are of interest in various fields encompassing systems on many different length scales [1]. Relevant examples are provided by biopolymes and bio-assemblies, including DNA, actin filaments, microtubules, and multicellular organisms like Bacillus subtilis [2]. Modern micro-manipulation techniques allow to observe and analyze filament dynamics on the single-molecule level [3], prompting a detailed understanding of the combined effects of fluctuations, hydrodynamics and elasticity. Static properties of elastic rods under external force loads have a long history of study, dating back to Euler [4]. Of current interest is the time-dependent behavior of such driven filaments and the possibility of shape bifurcations, especially on the biologically relevant nano-to-micro length scales, where viscous hydrodynamic dissipation dominates over inertia.In this Letter, we describe a hydrodynamic simulation technique for elastic filaments with arbitrary shape and rigidity and subject to external forces or boundary conditions, including full coupling between thermal, elastic and hydrodynamic forces in low-Reynolds-number flow. We apply our method to the whirling dynamics of a slender filament that is axially rotated at one end at frequency ω (while the other end is free). This model system has first been studied analytically by Wolgemuth et al. [5], who showed that a critical frequency ω c separates whirling (steady-state crankshafting motion with axial spinning) from twirling (diffusion-dominated simple axial rotation) by a supercritical Hopf bifurcation (i.e., a continuous shape transition). Subsequently, a zero-temperature simulation has been performed using the immersed boundary method [6], where the microscopic structure of the filament was modelled by interconnected springs. In contrast to analytic predictions, the filament was shown to undergo a subcritical (i.e. discontinuous) shape transition from twirling to a strongly bent state where the filament almost folds back on itself (termed overwhirling). The transition frequency ω c was found to be smaller by a factor of 3.5 compared to the analytic estimate. Using our simulation technique, we first map out the stability diagram in the absence c EDP Sciences