The effect of bond randomness on the spin-gapped ground state of the spin-1 bond-alternating antiferromagnetic Heisenberg chain is discussed. By using the loop cluster quantum Monte Carlo method, we investigate the stability of topological order in terms of the recently proposed twist order parameter [M. Nakamura and S. Todo: Phys. Rev. Lett. 89 (2002) 077204]. It is observed that the dimer phases as well as the Haldane phase of the spin-1 Heisenberg chain are robust against a weak randomness, though the valence-bond-solid-like topological order in the latter phase is destroyed by introducing a disorder stronger than the critical value.KEYWORDS: Haldane chain, bond randomness, quantum phase transitions, random-singlet phase, quantum Monte Carlo, loop algorithm, topological order, twist order parameterDisorder effects on low-dimensional quantum magnets have been investigated extensively in recent theoretical studies. In particular, impurity effects on spin-gapped Heisenberg antiferromagnets 1 have aroused much interest in relation to the impurity-induced antiferromagnetic (AF) long-range order observed experimentally in real materials.2 It has been established by recent numerical simulations 3, 4 that in two dimensions or higher, there are two classes of disorder, that affect spin-gapped states in essentially different ways. Site dilution and bond dilution are representatives of each class. The former induces localized moments around impurity sites. There exist strong correlations between such effective spins retaining the staggeredness with respect to the original lattice, and therefore the AF long-range order emerges by an infinitesimal concentration of dilution. In the bonddilution case, on the other hand, localized moments are always induced in pairs and they form a singlet again by AF interactions through the two-or three-dimensional shortest paths as long as the concentration of bond dilution is smaller than a finite critical value.In one dimension, since quantum fluctuations are much stronger than those in higher-dimensional systems, novel quantum critical phenomena are observed under disorder at the magnitude of coupling constants (bond randomness). Theoretically, the decimation renormalization group (DRG) approaches have achieved great success in predicting rich physics, such as the random-singlet (RS) phase for spin-1 2 chains.5-7 Recently, this technique has been extended to higher-spin cases, [8][9][10][11] where two of the main debates are on the robustness of the Haldane gap 12 against disorder and on the presence of the spin-1 RS phase. A number of numerical studies have also been carried out [13][14][15][16] to establish a quantitative phase diagram. However, this problem has not been sufficiently clarified yet. One of the main difficulties in simulating random quantum systems is the extremely wide energy scale that has to be taken into account. Another difficulty * Present address: Speech Interface Technology Group, NEC Corporation, Kawasaki 211-8666, Japan. † E-mail address: wistaria@ap.t.u-to...