We consider a system formed by two different segments of particles, coupled to thermal baths, one at each end, modeled by Langevin thermostats. The particles in each segment interact harmonically and are subject to an on-site potential for which three different types are considered, namely, harmonic, ϕ^{4}, and Frenkel-Kontorova. The two segments are nonlinearly coupled, between interfacial particles, by means of a power-law potential with exponent μ, which we vary, scanning from subharmonic to superharmonic potentials, up to the infinite-square-well limit (μ→∞). Thermal rectification is investigated by integrating the equations of motion and computing the heat fluxes. As a measure of rectification, we use the difference of the currents, resulting from the interchange of the baths, divided by their average (all quantities taken in absolute value). We find that rectification can be optimized by a given value of μ that depends on the bath temperatures and details of the chains. But, regardless of the type of on-site potential considered, the interfacial potential that produces maximal rectification approaches the infinite square well (μ→∞) when reducing the average temperature of the baths. Our analysis of thermal rectification focuses on this regime, for which we complement numerical results with heuristic considerations.