We present a precise measurement of the CP violation parameter sin2phi_1 and the direct CP violation parameter A_f using the final data sample of 772x10^6 B\bar B pairs collected at the Upsilon(4S) resonance with the Belle detector at the KEKB asymmetric-energy e+e- collider. One neutral B meson is reconstructed in a J/psi K0S, psi(2S) K0S, chi_c1 K0S or J/psi K0L CP-eigenstate and its flavor is identified from the decay products of the accompanying B meson. From the distribution of proper time intervals between the two B decays, we obtain the following CP violation parameters: sin2phi_1=0.667+-0.023(stat)+-0.012(syst) and A_f=0.006+-0.016(stat)+-0.012(syst).

We present a numerical method that consistently implements thermal fluctuations and hydrodynamic interactions to the motion of Brownian particles dispersed in incompressible host fluids. In this method, the thermal fluctuations are introduced as random forces acting on the Brownian particles. The hydrodynamic interactions are introduced by directly resolving the fluid motions with the particle motion as a boundary condition to be satisfied. The validity of the method has been examined carefully by comparing the present numerical results with the fluctuation-dissipation theorem whose analytical form is known for dispersions of a single spherical particle. Simulations are then performed for more complicated systems, such as a dispersion composed of many spherical particles and a single polymeric chain in a solvent.