We perform classical molecular dynamics simulations based on the Tersoff interatomic potential to investigate the size and surface orientation dependence of lattice constant and thermal expansion coefficient of one-dimensional silicon nanostructures. Three different surface orientations of silicon are considered, i.e., Si͑110͒, Si͑111͒, and Si͑100͒ with 2 ϫ 1 reconstruction. For each surface orientation, we investigate nanostructures with thicknesses ranging from 0.3 to 5.0 nm. We compute the vibrational amplitude of surface atoms, lattice constant, and thermal expansion coefficient as a function of size and temperature, and compare them for different surface orientations. An analytical expression is developed to compute the variation of the thermal expansion coefficient with size of the nanostructure.