We describe and implement a novel methodology for Monte Carlo simulation of one-dimensional killed diffusions. The proposed estimators represent an unbiased and efficient alternative to current Monte Carlo estimators based on discretization methods for the cases when the finite-dimensional distributions of the process are unknown. For barrier option pricing in finance, we design a suitable Monte Carlo algorithm both for the single barrier case and the double barrier case. Results from numerical investigations are in excellent agreement with the theoretical predictions.