This work was supported by the National Research Foundation (KRF) grant funded by the Korea gov- ernment (MEST) (2013R1A1A2057784) and Korea Ministry of Environment as “GAIA Project” (2012000550026). The main calculations were carried out at the Supercomputing Center/Korea Institute of Science and Technology Informa- tion with the supercomputing resources including technical support (KSC-2014-C1-028)