Figure 2: Computational box for the experiment, shown in Fig. 1 (b) and (c), relays on the sinusoidal laser energy distribution on the target surface (a). The  total computational cell therefore is schematically represented for multiprocessing mode in (b). In each particular processor core an atomistic-continuum MD-TTM model for metal part and ordinary MD model for water part are solved in 3D space (internal mesh is shown for one processor core). Beneath the non-reflective boundaries (NRB) an ordinary TTM model is solved in continuum to account on the heat flux outward the metal surface region describing the cooling process [10].