In this paper, we study the nonequilibrium dynamics of the Bose-Hubbard model with the nearestneighbor repulsion by using time-dependent Gutzwiller (GW) methods. In particular, we vary the hopping parameters in the Hamiltonian as a function of time, and investigate the dynamics of the system from the density wave (DW) to the superfluid (SF) crossing a first-order phase transition and vice versa. From the DW to SF, we find scaling laws for the correlation length and vortex density with respect to the quench time. This is a reminiscence of the Kibble-Zurek scaling for continuous phase transitions and contradicts the common expectation. We give a possible explanation for this observation. On the other hand from SF to DW, the system evolution depends on the initial SF state. When the initial state is the ground state obtained by the static GW methods, a coexisting state of the SF and DW domains forms after passing through the critical point. Coherence of the SF order parameter is lost as the system evolves. This is a phenomenon similar to the glass transition in classical systems. When the state starts from the SF with small local phase fluctuations, the system obtains a large size DW domain structure with thin domain walls. them, works in [27, 28] questioned the applicability of the KZ scaling theory to the QPT, whereas [29,30] concluded that the observed results were in good agreement with the KZ scaling law.In this paper, we focus on the two-dimensional (2D) Bose-Hubbard model (BHM) [33,34], which is a canonical model of the bosonic ultra-cold atomic gas systems in an optical lattice. In particular, we add nearestneighbor (NN) repulsions between atoms. Then, the resultant system is described by an extended Bose-Hubbard model (EBHM). As a result, a parameter region corresponding to the density wave (DW) appears in the ground-state phase diagram, in addition to the Mott insulator and SF. Near the half filling, there exists a firstorder phase transition between the SF and DW [35]. We shall study the quench dynamics of the EBHM on passing across the SF and DW phase boundary. There are only a few works for the dynamical properties of quantum systems at first-order phase transitions under a quench [36][37][38], and therefore detailed study on that problem is desired.This paper is organized as follows. In section 2, we introduce the EBHM and explain the Gutzwiller (GW) methods, which are used in the present work. In section 3, quench dynamics of the first-order phase transition from the DW to SF is studied. Behavior of SF and DW orders are investigated by solving the Schrödinger equation by means of time-dependent GW (tGW) methods. We focus on the order parameters, correlation length, vortex number, etc, in particular, scaling laws of these quantities with respect to the quench time τ Q . Contrary to the common expectation, we find that scaling laws hold for the correlation length and vortex density. In section 4, we give a possible explanation of the observed results from viewpoint of the SF-bubblenucleation process. W...