1 We consider long-range soliton interaction mediated by radiation in nonlinear 1 d system with frozen disorder. The problem is of a great importance for nonlinear fiber optics of the next generation (see, e.g., [1,2]), and it is also of general relevance for any of the traditional fields, like plasma physics, where propagation of solitary waves is possible. Our aim here is to answer the following sets of fundamental questions:What statistics describe the radiation emitted due to disorder by a single soliton or a pattern of solitons? How far do the radiation wings extend from the peak of the soliton(s)? What is the structure of the wings?How strong is the radiation mediating interaction between the solitons? How is the interaction modified if we vary the soliton positions and phases within a pattern of solitons?We focus on the dynamics of wave packets. The universal coarse-grained description of a wave packet envelope is given by the nonlinear Schrödinger equation (NLS) [3][4][5]. We consider the 1 d problem motivated mainly by applications to fiber optics [6] (1)The medium (fiber) is imperfect; i.e., various macroscopic characteristics of the fiber fluctuate in space. Fluctuations of the dispersion coefficient, d , are believed to be one of the major sources of disorder present in real fibers [7]. This disorder is frozen; i.e., d is a random function of z . We assume that d fluctuates on short spatial scales and that the fiber is homogeneous on larger scales. The averaged value of d is a constant, which can be rescaled to unity by changing the units 1 This article was submitted by the authors in English., where 〈ξ〉 = 0. According to the Central Limit Theorem [8], ξ at scales larger than the correlation length can be treated as a homogeneous Gaussian random process with zero mean and described by the quantity D = , which is the noise intensity. The pair correlation function of ξ is (2) We assume that the disorder is weak; i.e., D Ӷ 1.At z = 0, a sequence of well-separated solitons is launched. In an ideal medium ( ξ = 0), each of the solitons is preserved dynamically gaining, according to the exact single-soliton solution of Eq. (1), only a multiplicative phase factor. Because the medium is imperfect, the solitons, perturbed by impurities, shed radiation. The first problem is to describe the radiation. The soliton looses energy shedding radiation. Another problem is to describe the degradation of a single soliton. The tails of different solitons interfere with each other, forming a collective background. This fluctuating background affects all solitons. It results in the emergence of a long-range effective intersoliton interaction, which is the final (but not the least) problem to be addressed. The long-range interaction dominates the direct interaction due to overlapping of soliton tails, as this direct interaction decays exponentially with separation [9,10]. The emergence of the long-range interaction between the imperfect solitons in the pure ( ξ = 0) NLS, mediated by the emitted radiation, was noted in [11]....