The early stages of the bacteriorhodopsin photocycle, including the J625, K~w, and L550 intermediates and the role of water molecules within the protein interior, are studied by means of molecular dynamics simulations. Our calculations examine two models for the excited state potential surface governing the observed all-trans - 13-cis photoisomerization: one surface hindering a Cl4-CI5 single-bond corotation and the other surface allowing such corotation. The investigations use as a starting structure a model of bacteriorhodopsin based on electron-microscopy studies and subsequent molecular dynamics refinement. The following questions are addressed: How does the binding site guide retinal’s photoisomerization? How does the photoisomerization depend on features of the excited state potential surface? Can one recognize a J625 intermediate? How does water participate in the early part of the pump cycle? How is the initial photoreaction affected by a lowering of temperature? To model the quantum yield, i.e., the dependence of the dynamics on initial conditions, 50 separate isomerization teals are completed for each potential surface, at both 300 and 77 K, the trials being distinguished by different initial, random velocity distributions. From these trials emerge, besides all-transretinal, three different photoproducts as candidates for the K5w intermediate: (1) 13-cis-retinal, with the Schiff base proton oriented toward Asp-96; (2) 13-cis-retinal, highly twisted about the c6-c7 bond, with the Schiff base proton oriented perpendicular to the membrane normal; (3) 13,14-dicis-retinal, with the Schiff base proton oriented toward the extracellular side. Two candidates for the K~w intermediate, case 2 and case 3 above, were subjected to simulated annealing to determine corresponding L550 structures. We suggest that photoproduct 2 above most likely represents the true K~w intermediate. Water molecules near the Schiff base binding site are found to play a crucial role in stabilizing the K590 state and in establishing a pathway for proton transfer to Asp-85.