Mergers of compact object binaries are one of the most powerful sources of gravitational waves (GWs) in the frequency range of second-generation ground-based gravitational wave detectors (Advanced LIGO and Virgo). Dynamical simulations of young dense star clusters (SCs) indicate that ~27 per cent of all double compact object binaries are members of hierarchical triple systems (HTs). In this paper, we consider 570 HTs composed of three compact objects (black holes or neutron stars) that formed dynamically in N-body simulations of young dense SCs. We simulate them for a Hubble time with a new code based on the Mikkola's algorithmic regularization scheme, including the 2.5 post-Newtonian term. We find that ~88 per cent of the simulated systems develop Kozai-Lidov (KL) oscillations. KL resonance triggers the merger of the inner binary in three systems (corresponding to 0.5 per cent of the simulated HTs), by increasing the eccentricity of the inner binary. Accounting for KL oscillations leads to an increase of the total expected merger rate by ~50 per cent. All binaries that merge because of KL oscillations were formed by dynamical exchanges (i.e. none is a primordial binary) and have chirp mass >20 Msun. This result might be crucial to interpret the formation channel of the first recently detected GW events.

In this paper, we study the formation and dynamical evolution of black hole-black hole (BH-BH) binaries in young star clusters (YSCs), by means of N-body simulations. The simulations include metallicity-dependent recipes for stellar evolution and stellar winds, and have been run for three different metallicities (Z = 0.01, 0.1 and 1 Zsun). Following recent theoretical models of wind mass-loss and core-collapse supernovae, we assume that the mass of the stellar remnants depends on the metallicity of the progenitor stars. We find that BH-BH binaries form efficiently because of dynamical exchanges: in our simulations, we find about 10 times more BH-BH binaries than double neutron star binaries. The simulated BH-BH binaries form earlier in metal-poor YSCs, which host more massive black holes (BHs) than in metal-rich YSCs. The simulated BH-BH binaries have very large chirp masses (up to 80 Msun), because the BH mass is assumed to depend on metallicity, and because BHs can grow in mass due to the merger with stars. The simulated BH-BH binaries span a wide range of orbital periods (10^-3-10^7 yr), and only a small fraction of them (0.3 per cent) is expected to merge within a Hubble time. We discuss the estimated merger rate from our simulations and the implications for Advanced VIRGO and LIGO.