  We use hydrodynamic simulations to predict correlations between Lya forest
absorption and galaxies at redshift z~3. The probability distribution function
(PDF) of Lya flux decrements shifts systematically towards higher values in the
vicinity of galaxies, reflecting the overdense environments in which these
galaxies reside. The predicted signal remains strong in spectra smoothed over
50-200 km/s, allowing tests with moderate resolution quasar spectra. The strong
bias of high redshift galaxies towards high density regions imprints a clear
signature on the flux PDF, but the predictions are not sensitive to galaxy
baryon mass or star formation rate, and they are similar for galaxies and for
dark matter halos. The dependence of the flux PDF on galaxy proximity is
sensitive to redshift determination errors, with rms errors of 150-300 km/s
substantially weakening the predicted trends. On larger scales, the mean galaxy
overdensity in a cube of 5 or 10 Mpc/h (comoving) is strongly correlated with
the mean Lya flux decrement on a line of sight through the cube center. The
slope of the correlation is ~3 times steeper for galaxies than for dark matter
as a result of galaxy bias. The predicted large scale correlation is in
qualitative agreement with recently reported observational results. However,
observations also show a drop in absorption in the immediate vicinity of
galaxies, which our models do not predict even if we allow the galaxies or AGNs
within them to be ionizing sources. This decreased absorption could be a
signature of galaxy feedback on the surrounding IGM, perhaps via galactic
winds. Peculiar velocities often allow gas at comoving distances ~1.5 Mpc/h to
produce saturated absorption at the galaxy redshift, so any feedback mechanism
must suppress neutral hydrogen out to these radii to match the data. (Abridged)
