  We investigate non-perturbative supersymmetry breaking in various models of
quantum mechanics, including an interesting class of $PT$-invariant models,
using lattice path integrals. These theories are discretized on a temporal
Euclidean lattice with anti-periodic boundary conditions. Hybrid Monte Carlo
algorithm is used to update the field configurations to their equilibrium
values. We used the Ward identities and expectation values of superpotentials
as tools for probing supersymmetry breaking.
