Shock waves propagate in falling coupled harmonic oscillators. The bottom end of coupled harmonic oscillators does not fall downwards until a shock wave reaches the bottom end. The exact solution can be expressed by the Fourier series expansion, and an approximate solution can be expressed by the integral of the Airy function. The width of the shock wave increases slowly in accordance with a power law.Shock waves are generated in compressive fluids. Typical shock waves appear in air compressed by supersonic planes or meteorites. There is a jump in the fields of pressure, temperature, and fluid velocity. The Rankine-Hugoniot relation is satisfied for the jump under normal shock [1]. There have been numerous investigations of shock waves [2]. It is considered that a shock wave is a typical nonlinear wave. The simplest model of a shock wave is the Burgers equation [3]. Nonlinearity and dissipation are essential for shock waves.We consider a linear chain of coupled harmonic oscillators under gravity. It is a linear system and there is no dissipation. It is a typical system of particles considered in a basic course of mechanics. However, there is a nontrivial phenomenon similar to a shock wave in this simple system. A similar phenomenon was discussed in falling elastic bars, using a partial differential equation [4]. We will focus on the effect of the discreteness in this paper.The model equation is written aswhere N is the total number of particles of mass m, k is the spring constant, a is the natural length of the spring, x i is the height of the ith particle, and g denotes the acceleration of gravity. The heights of the bottom and top particles are expressed respectively as x 1 and x N . If x N is fixed to a constant value x N 0 by holding the top particle, the stationary positions of the other particles are determined from the relation x i+1 − x i = a + (mgi/k) aswhere the position of the bottom particle, x 1 , is expressed asWe study the free-fall motion of this system by releasing the top particle with an initial velocity of 0 from the stationary state. There are five parameters, i.e., m, k, a, g, and N , in this system, but N is the only essential parameter. The other parameters can be set to a unit value by changing the scales of x and t. Figure 1 shows the positions of ten particles at t = 0, 1, · · · , 10 for N = 10, k = 1, g = 0.1, m = 1, and a = 1. The top particle falls with a nearly constant velocity. The bottom particle does not move until t ∼ 8 in the gravity field.