  We have further developed and extended a method for calculation of atomic
properties based on a combination of the configuration interaction and
coupled-cluster approach. We have applied this approach to the calculation of
different properties of atomic lead, including the energy levels, hyperfine
structure constants, electric-dipole transition amplitudes, and E1 parity
nonconserving (PNC) amplitude for the (6p^2) ^3P_0 - (6p^2) ^3P_1 transition.
The uncertainty of the E1 PNC amplitude was reduced by a factor of two in
comparison with the previous most accurate calculation [V.~A.~Dzuba et.al.,
Europhys. Lett. 7, 413 (1988)]. Our value for the weak charge Q_W=-117(5) is in
agreement with the standard model prediction.
