Let λ be a positive number, and let ( x j : j ∈ Z ) ⊂ R be a fixed Riesz-basis sequence, namely, ( x j ) is strictly increasing, and the set of functions { R ∋ t ↦ e i x j t : j ∈ Z } is a Riesz basis ( i.e., unconditional basis) for L 2 [ − π , π ] . Given a function f ∈ L 2 ( R ) whose Fourier transform is zero almost everywhere outside the interval [ − π , π ] , there is a unique sequence ( a j : j ∈ Z ) in ℓ 2 ( Z ) , depending on λ and f , such that the function I λ ( f ) ( x ) ≔ ∑ j ∈ Z a j e − λ ( x − x j ) 2 , x ∈ R , is continuous and square integrable on ( − ∞ , ∞ ) , and satisfies the interpolatory conditions I λ ( f ) ( x j ) = f ( x j ) , j ∈ Z . It is shown that I λ ( f ) converges to f in L 2 ( R ) , and also uniformly on R , as λ → 0 + . In addition, the fundamental functions for the univariate interpolation process are defined, and some of their basic properties, including their exponential decay for large argument, are established. It is further shown that the associated interpolation operators are bounded on ℓ p ( Z ) for every p ∈ [ 1 , ∞ ] .