Compute minimum energy paths from arbitrary positions on the potential energy surface
There are some possibilities here:


The negative of the gradient of a function always points to the direction of optimal descent (for a discussion, see here).
A method that simply follows the gradient is called first order.
As permeakra pointed out, extrapolation techniques are commonly used for optimizations.
This is because we take finite steps in the PES and this introduces the risk of "jumping" over the minima.
Better approximations are either second order in nature (e.g., the Hessian matrix is calculated in each step) or quasi-second order (the Hessian matrix is approximated or approximately updated).

In order to simply follow the gradient you want the steepest descent-like, first order method with probably a very small maximum step and this can be accomplished by setting the Hessian or the Hessian approximation equal to the identity matrix.
In ORCA (I have never used Molcas, sorry) for example you may use (source):

! Opt
%geom
  inhess unit
  update noupdate
  step qn
end

permeakra also cited the Nudged Elastic Band, which is used for finding minimum energy paths connecting two structures.
This might be very useful and is available for instance in VASP and in NWChem (free).
You may also want to take a direction from one Hessian eigenvector.
This is commonly associated with a molecular vibrational degree of freedom.
This may be accomplished with NWChem.


If you just want to go down the PES from an arbitrary geometry without any jumping, solution 1. is probably the one you want.