Nonlinear breathing modes at a defect

Recent molecular dynamics (MD) simulations of Cubero et al (1999) of a DNA duplex containing the 'rogue' base difluorotoluene (F) in place of a thymine (T) base show that breathing events can occur on the nanosecond timescale, whereas breathing events in a normal DNA duplex take place on the microsecond timescale. The main aim of this paper is to analyse a nonlinear Klein-Gordon lattice model of the DNA duplex including both nonlinear interactions between opposing bases and a defect in the interaction at one lattice site; each of which can cause localisation of energy. Solutions for a breather mode either side of the defect are derived using multiple-scales asymptotics and are pieced together across the defect to form a solution which includes the effects of the nonlinearity and the defect. We consider defects in the inter-chain interactions and in the along chain interactions. In most cases we find in-phase breather modes and/or out-of-phase breather modes, with one case displaying a shifted mode.

Energy norm a-posteriori error estimation for hp-adaptive discontinuous Galerkin methods for elliptic problems in three dimensions

We develop the energy norm a-posteriori error estimation for hp-version discontinuous Galerkin (DG) discretizations of elliptic boundary-value problems on 1-irregularly, isotropically refined affine hexahedral meshes in three dimensions. We derive a reliable and efficient indicator for the errors measured in terms of the natural energy norm. The ratio of the efficiency and reliability constants is independent of the local mesh sizes and weakly depending on the polynomial degrees. In our analysis we make use of an hp-version averaging operator in three dimensions, which we explicitly construct and analyze. We use our error indicator in an hp-adaptive refinement algorithm and illustrate its practical performance in a series of numerical examples. Our numerical results indicate that exponential rates of convergence are achieved for problems with smooth solutions, as well as for problems with isotropic corner singularities.