We demonstrate a systematic implementation of coupling between a scalar field and the geometry of the space (curve, surface, etc.) which carries the field. This naturally gives rise to a feedback mechanism between the field and the geometry. We develop a systematic model for the feedback in a general form, inspired by a specific implementation in the context of molecular dynamics (the so-called Rahman-Parrinello molecular dynamics, or RP-MD). We use a generalized Lagrangian that allows for the coupling of the space's metric tensor (the first fundamental form) to the scalar field, and add terms motivated by RP-MD. We present two implementations of the scheme: one in which the metric is only time-dependent [which gives rise to ordinary differential equation (ODE) for its temporal evolution], and one with spatio-temporal dependence [wherein the metric's evolution is governed by a partial differential equation (PDE)]. Numerical results are reported for the (1+1)-dimensional model with a nonlinearity of the sine-Gordon type.Recently, much attention has been focused on softcondensed-matter objects, such as vesicles, microtubules, and membranes [1][2][3][4]. Many nanoscale physical systems, including nanotubes and electronic and photonic waveguide structures [5,6], have nontrivial geometry and are influenced by substrate effects. These classes of systems, many of which are inherently nonlinear, raise the question of the interplay between nonlinearity and a substrate with variable curvature. Of particular interest is a possibility of developing curvature in the substrate due to forces generated by the nonlinear field. The resulting curvature can in turn affect the field.There is an increasing body of literature dealing with the interplay of nonlinearity and a curved substrate. Usually, however, the substrate geometry is assumed to be fixed, see, e.g., [7]. Nevertheless, for many applications, ranging from condensed matter to optics to biophysics, it is relevant to introduce models that admit a flexible substrate, which is affected by the field(s) that it carries, as well as feeding back into the field dynamics. In this situation, equations for the fields in a nonlinear system abutting on the flexible substrate should include both the field dynamics proper and the feedback coupling to the substrate. Equations for the evolution of the substrate should in turn be affected by the evolution of the field. A prototypical physical example of this type is Euler buckling [8], where the evolution of a thermal profile causes the underlying surface to buckle (and hence locally modify its curvature).In a discrete setting, a model of this type has recently been presented in [9]. However, it was limited to a system of masses coupled by nonlinear springs. Some studies have also been performed in a special case of the continuum limit of classical spin systems (such as the Heisenberg chain) coupled to the curvature; geometric frustration was found to arise in such settings [10].About twenty years ago, a problem similar to the theme of our...