By using large-scale molecular dynamics simulations, the dynamics of two-dimensional (2D) supercooled liquids turns out to be dependent on the system size, while the size dependence is not pronounced in three dimensional (3D) systems. It is demonstrated that the strong system-size effect in 2D amorphous systems originates from the enhanced fluctuations at long wavelengths, which are similar to those of 2D crystal phonons. This observation is further supported by the frequency dependence of the vibrational density of states, consisting of the Debye approximation in the low-wavenumber-limit. However, the system-size effect in the intermediate scattering function becomes negligible when the length scale is larger than the vibrational amplitude. This suggests that the finite-size effect in a 2D system is transient and also that the structural relaxation itself is not fundamentally different from that in a 3D system. In fact, the dynamic correlation lengths estimated from the bond-breakage function, which do not suffer from those enhanced fluctuations, are not size dependent in either 2D or 3D systems. PACS numbers: 62.60.+v, 61.20.Lc, Dimensionality plays a key role in the physics of solids and liquids -from high to low dimensions -and fluctuation shows up differently, as typically observed in phase transitions [1,2]. Indeed, two-dimensional (2D) systems often exhibit enhanced fluctuations, leading to various anomalies that are not experienced in three-dimensional (3D) systems. The melting of a 2D solid is a marked example [3][4][5][6][7][8][9], where the long-wavelength structural correlation is induced by thermal fluctuations sthat span an infinite length. For the glass transition from supercooled liquids to amorphous solids, the dimensionality dependence of the fluctuation has become an issue only recently. Gigantic fluctuation in 2D supercooled liquids has been observed that is far stronger than that in their 3D counterparts [10][11][12]. The aim of this Letter is to elucidate the similarity of this fluctuation to that in crystals [13], and also to investigate the heterogeneous dynamics in both 2D and 3D systems.For a crystalline solid of monodisperse particle assemblies, the mean-squared thermal displacement (MSTD) is given by using the vibrational density of state (VDOS) g(ω) as a function of angular frequency ω aswhere m the particle mass, d the spatial dimension, and (k B T ) −1 the inverse temperature. Under the Debye approximation for the VDOS of acoustic plane waves, g(ω) becomes proportional to ω d−1 [14]. It leads to divergence of the integral in 2D systems owing to the low-frequency acoustic waves, while it converges in 3D systems. As a result, the long-range translation order is prohibited in 2D systems [15,16]. Integration of Eq. (1) over ω ≥ 2πc/L provides us with its dependence on the linear system size L aswhere µ and K are shear and bulk moduli, σ 0 is the particle radius, and c is the velocity of sound. Such fluctuation is the source of the size-dependent behavior of 2D solids undergoing melting...