We present a review of recent works on clustering phenomena in unstable nuclei studied by antisymmetrized molecular dynamics (AMD). The AMD studies in these decades have uncovered novel types of clustering phenomena brought about by the excess neutrons. Among them, this review focuses on the molecule-like structure of unstable nuclei.One of the earliest discussions on the clustering in unstable nuclei was made for neutron-rich Be and B isotopes. AMD calculations predicted that the ground state clustering is enhanced or reduced depending on the number of excess neutrons. Today, the experiments are confirming this prediction as the change of the proton radii. Behind this enhancement and reduction of the clustering, there are underlying shell effects called molecular-and atomic-orbits. These orbits form covalent and ionic bonding of the clusters analogous to the atomic molecules. It was found that this "molecular-orbit picture" reasonably explains the low-lying spectra of Be isotopes. The molecular-orbit picture is extended to other systems having parity asymmetric cluster cores and to the three cluster systems. O and Ne isotopes are the candidates of the former, while the 3α linear chains in C isotopes are the latter. For both subjects, many intensive studies are now in progress.We also pay a special attention to the observables which are the fingerprint of the clustering. In particular, we focus on the monopole and dipole transitions which are recently regarded as good probe for the clustering. We discuss how they have and will reveal the exotic clustering.where X denotes Z i , ν or χ i . λ and µ are arbitrary numbers, but µ must have negative sign. The energy of the system H is defined as (8) As time t being evolved by the Eq. (7), the energy of the system is decreased and we obtain the set of parameters which minimizes the energy of the system. The calculation which employs the eigenstate of the parity Φ π to evaluate the energy is called variation before angular momentum projection (VBP), while the calculation employing the eigenstate of parity and angular momentum Φ Jπ M K is called variation after angular momentum projection (VAP) [42]. In the case of VBP, the angular momentum projection is performed after the energy minimization.To obtain the wave functions which have the configurations different from the energy minimum, the constrained energy minimization is often performed in VBP calculations. For example, the constraint on the quadrupole deformation parameter β is imposed by adding the constraint potential to the energy asHere, the magnitude of constraint potential v β is chosen large enough so that the deformation parameter β of the wave function approximately equals to β 0 after the energy minimization.Another way to obtain the excited configurations is orthogonalization method [42] that is used in VAP calculatons in which the variational wave function is orthogonalized to the energy minimum state By applying this procedure successively, the second, third and more excited configurations are obtained. He...