abstract: The quest for large and low frequency band gaps is one of the principal objectives pursued in a number of engineering applications, ranging from noise absorption to vibration control, to seismic wave abatement. For this purpose, a plethora of complex architectures (including multi-phase materials) and multi-physics approaches have been proposed in the past, often involving difficulties in their practical realization.
  .
  To address this issue, in this work we propose an easy-to-manufacture design able to open large, low frequency complete Lamb band gaps exploiting a suitable arrangement of masses and stiffnesses produced by cavities in a monolithic material. The performance of the designed structure is evaluated by numerical simulations and confirmed by Scanning Laser Doppler Vibrometer (SLDV) measurements on an isotropic polyvinyl chloride plate in which a square ring region of cross-like cavities is fabricated. The full wave field reconstruction clearly confirms the ability of even a limited number of unit cell rows of the proposed design to efficiently attenuate Lamb waves. In addition, numerical simulations show that the structure allows to shift of the central frequency of the BG through geometrical modifications. The design may be of interest for applications in which large BGs at low frequencies are required.
address: $a$) University of Le Havre, Laboratoire Ondes et Milieux Complexes, UMR CNRS 6294, 75 Rue Bellot, 76600 Le Havre, France; $b$) Civil, Architectural & Environmental Engineering (CAEE) Department, Drexel University, 3141 Chestnut St., Philadelphia, PA 19104, USA; $c$) Institute of Fluid-Flow Machinery, Polish Academy of Science, Fiszera 14 st. 80-231 Gdańsk, Poland; $d$) Department of Physics and Nanostructured Interfaces and Surfaces Centre, University of Torino, Via Pietro Giuria 1, 10125 Torino, Italy; $e$) Laboratory of Bio-Inspired and Graphene Nanomechanics, Department of Civil, Environmental and Mechanical Engineering, University of Trento, Via Mesiano 77, 38123 Trento, Italy; $f$) School of Engineering and Materials Science, Queen Mary University of London, Mile End Road, London E1 4NS, United Kingdom; $g$) Ket Lab, Edoardo Amaldi Foudation, Italian Space Agency, Via del Politecnico snc, 00133 Rome, Italy; $*$ Correspondence to: email@example.com, firstname.lastname@example.com
author: M. Miniaci$^{a,*}$, M. Mazzotti$^{b}$, M. Radzieński$^c$, N. Kherraz$^a$, P. Kudela$^c$, W. Ostachowicz$^c$, B. Morvan$^a$, F. Bosia$^{d,*}$, N. M. Pugno$^{e,f,g}$
title: Experimental observation of a large, low frequency Lamb band gap induced by periodic cross-like cavities in a polymer waveguide

# Introduction

One of the main problems facing physicist and engineers working in the field of metamaterials is to achieve vibration damping and control over large, low frequency ranges. This is true in fields ranging from noise absorption to seismic wave abatement . Phononic crystals (PCs) and acoustic metamaterials (AMMs), generally made of periodically distributed inclusions in a matrix (or hosting material), are of particular interest because of their ability to act as stop-band filters, i.e. attenuate mechanical waves over entire frequency bands, commonly known as band gaps (BGs) .

Bragg scattering and local resonance are the mechanisms mainly exploited for BG nucleation. BGs due to Bragg scattering arise from the wave diffraction by periodic inclusions and thus occur at wavelengths of the order of the unit cell size, whilst local resonance is related to the vibration of individual elements within the medium, and is thus independent of the spatial periodicity of the lattice . In some cases, Bragg BGs can be coupled with hybridization BGs , due to resonating "inhomogeneities" (Fano interference effects , electrical resonances , etc.).

The concept of a BG naturally lends itself to applications involving vibration damping, acoustic filtering and wave attenuation, as well as waveguiding . Applications at high frequencies include wave filters, couplers, sensing devices, wave splitters, demultiplexers, etc. . A nonexhaustive list of low-frequency applications includes acoustic absorption, vibration shielding, subwavelength imaging, cloaking, etc .

With these applications in mind, a plethora of different architectures (including multi-phase materials) and multi-physics approaches have been proposed in recent years , but many involve considerable practical difficulties in their realization due to their inherent complexity.

Recently, particular attention has been focused on phononic plates because of their potential technological applications, ranging from microelectromechanical systems to nondestructive evaluation . Among others, phononic plates made of periodic distributions of studs or gratings on the surface as well as those realized by a periodic distribution of empty cavities perpendicular or parallel to the propagation plane have been studied.

In addition to theoretical predictions, experimental measurements of BGs in phononic plates have been performed. In this context, Brunet et al. reported BGs for Lamb waves propagating in rectangular and square arrays of holes drilled in a silicon phononic plate, whilst Bonello et al. proved the existence of BGs for Lamb waves in silicon plates coated by a very thin phononic film. Publications dealing with BGs for Lamb waves generated in a plate with a periodic grating on the surface have also appeared with the aim of quantitatively verifying the relation between the width of the BG and the depth of the grooves. Finally, several studies have considered and shown that filtering and waveguiding properties are achievable in plates with stubs, including bi-phase materials .

In these approaches, geometrical/physical complexity is often inevitable, leading to considerable complications in their practical realization, which requires unit cells capable of providing the desired dynamic performance, with realistic structural features, such as being compact, lightweight and easy to manufacture. In this work, we propose phononic plate that is simple to fabricate and is capable of opening large, low frequency complete Lamb BGs, exploiting a suitable arrangement of masses and stiffnesses produced by cavities in a monolithic material. This simple but efficient design capable of producing large and low frequency Lamb band gaps when cross-like holes are considered, showing that a sort of "mass centrifugation" process is beneficial. The idea is supported by numerical calculations and experimental Scanning Laser Doppler Vibrometer (SLDV) measurements on an isotropic PolyVinyl Chloride (PVC) plate with a square ring region of cross-like holes. The full wave field reconstruction existing below, within and above the BG frequencies clearly confirms the ability of the proposed design to stop Lamb waves in the desired frequency range.

The paper is organized as follows: first we present extensive numerical simulations aimed at finding suitable plate configurations capable of nucleating large BGs within the optimal operative frequency range of the SLDV using Finite Element simulations (FEM).

Through-the-thickness holes with different geometries are considered as scatterers (inclusions). Next, finite element time-transient analysis is performed to study the influence of the phononic region length (number of unit cells) on wave propagation. Finally, SLDV measurements on ad-hoc machined phononic polyvinyl chloride (PVC) plate are presented and compared to the numerically predicted results. The calculated BGs are verified experimentally and the direct observation of scattering phenomena below, within and above the BGs is described.

# Design and numerical analysis

In this section, the description of the phononic plate design is provided. Different geometries of the scatterers, made of through-the-thickness cavities, are considered to find Bragg-type BGs in the frequency range $0 - 70$ kHz.

In view of the experimental phase, PVC has been chosen as the material matrix since it can be easily machined. The PVC mechanical properties are the density $\rho_{PVC} = 1430$ kg$/$m$^3$, the Young's modulus $E_{PVC} = 3$ GPa and the Poisson's ratio $\nu_{PVC} = 0.4$ .

A lattice parameter of $a = 20$ mm is chosen for the unit cell and plate thickness $H = 0.4 a$, so that the lowest Bragg BGs relative to the $A0$ mode are expected at around $14$ and $36$ kHz, corresponding to the half and to the full wavelength of the $A0$ propagating mode, respectively. The $A0$ mode is chosen as it is easy to excite via piezoelectric (PZT) sensors and to measure via SLDV. Once the lattice parameter has been set, band structures for different types of structures are computed. The unit cells are shown in Fig. : relative to (a) circular, (b) square and (c) rounded cross-like cavities. Geometrical parameters (see Table ) are chosen so as to guarantee a constant filling fraction $f_f = 41\%$ for all the three models.

The corresponding band structures are shown in Figs. d-f along with the mode shapes for some specific branches (g-j). The dispersion diagram is computed along the three high symmetry directions of the first irreducible Brillouin zone (Fig. ) exploiting the Bloch-Floquet theorem by using COMSOL MultyPhysics. Band structures are reported in terms of reduced wavevector $k^*$ with wavenumber components $k^* = [k_x a/\pi; k_y a/\pi]$. A detailed description of the finite element procedure can be found in Ref. .

Phononic plates with circular and square holes do not nucleate any BG, whereas the unit cell with a rounded cross-like hole exhibits two complete BGs (highlighted as light grey regions). In order to shed light on this different behaviour we highlight the in-plane or out-of-plane polarization of the modes in the dispersion diagrams by defining a polarization factor:

$$p=\frac{\int_V (|u_z |)^2 dV}{\int_V (|u_x |^2+|u_y |^2+ |u_z |^2)dV},$$

where $V$ is the volume of the unit cell, $u_x$, $u_y$ and $u_z$ are the displacement components along $x$, $y$ and $z$ axes, respectively. The points of the dispersion curves are shaded accordingly, with colours varying from $p = 0$ (blue) to $p = 1$ (red). Thus, colours close to red indicate vibration modes that are predominantly polarized out-of-plane, while colours close to blue are predominantly polarized in-plane.

This allows us to infer that the different arrangement of the cavities is responsible for the shift of some dispersion bands while leaving others unaltered. For instance, Figs. d,e highlight the shift occurring for two modes: one with a mainly in-plane deformation mechanism (blue rectangle) and another one with a mainly out-of-plane deformation mechanism (red rectangle). The shift is clearly due to the equivalent stiffness variations occurring due to the cavity geometry (see Fig. g,h). Introducing a cross-like cavity can accentuate this process and lead to the opening of two BGs (Fig. f) with a smaller deformation of the unit cell for the same branches (Fig. j).

Results hold true for varying plate thickness with the same type of mechanical effects. Fig. shows a parametric analysis in which band structures for the unit cell with a cross-like cavity have been computed increasing the thicknesses from $H = 0.2 a$ to $H = 1.2 a$ in steps of $0.2 a$. Results show that the plate thickness plays a critical role for (i) the nucleation of BGs, (ii) the frequency shift of the curves, and thus of the upper and/or lower BG boundaries, as well as (iii) polarization of the unit cell modes (as for the previous case, the color scale denotes here the in-plane$/$out-of-plane contribution to the total modal displacement). If the thickness of the plate is relatively small (see Fig. a for $H \le 0.4 a$) the modes are characterized by mainly in-plane or -out-of-plane deformation (curves are either red or blue) whereas for thicknesses $H \ge 0.8 a$) (see Fig. b) deformation modes with non-negligible in-plane and out-of plane components are present.

The dependence of the BGs with respect to the plate thickness are summarized in Fig. where the lower and upper bounds and width of the BGs are reported as a function of $H/a$ (in steps of $0.1 a$). No BG appears in the system when $H/a = 0.2$ whereas a single large BG from approximately $20.70$ kHz to $28.00$ kHz is present when $H/a = 0.7$. In all the other cases two total BGs are present. The plate thickness giving rise to the widest total BG is for $H/a = 0.6$, corresponding to a relative width of $36 \%$($f_{BG_L} \in [20.60 - 27.95]$ kHz and $f_{BG_U} \in [28.42 - 32.52]$ kHz, (highlighted in red in Fig. ). In this case, the two BGs are connected, except for an almost flat band between them.

In conclusion, results show that as the plate thickness increases, the frequency range covered by BGs increases up to the largest value at $H/a = 0.6$, which is therefore the chosen thickness for the manufacturing of experimental samples. To determine the necessary number of repeated unit cell to obtain efficient signal damping in experiments, finite element time transient analysis are performed (using ABAQUS). Numerical models with unit cells such as the one presented in Fig. c, arranged in concentric square rings with $1, 2, 3$ and $4$ rows in a specific portion of the PVC plate, as shown in Fig. , are considered.

First, a wave propagation analysis is performed to evaluate the screening power of the phononic region as a function of the number of rows of unit cells over a large frequency range. Guided waves are excited by means of an imposed orthogonally displacement ($1 \times 10^{-6}$ mm in the $z$-direction) at point $A$ (see Fig. a). A pulse of $2$ sinusoidal cycles centred at $50$ kHz and modulated by a Hanning window (Fig. a) is chosen since it generates Lamb waves over a large frequency spectrum.

The out-of-plane displacement detected at point $R2$ for the four considered cases along with their Fourier transforms are shown in Fig . The system is modelled as linear elastic, so that the signal attenuation is totally due to the Bragg scattering mechanism. The results clearly show that waves with a frequency content falling outside the BG (gray rectangle) propagate through the phononic region reaching the monitoring point $R2$ without substantial attenuation as the number of cross-like cavities increases, whilst when the frequency content of the propagating waves falls inside the BG, the phononic crystal ring region inhibits the wave propagation so that the displacement at point $R2$ is much smaller as the number of cross-like cavities increases. Results show that the destructive interference mechanism starts taking place after two unit cells and that as the number of unit cells increases the edges of the attenuation frequency regions become sharper. The results are presented in semi-logarithmic scale and prove that each additional row of unit cells introduces an attenuation factor of about $4$, allowing to infer that a plate with four rows of unit cells can reduce incident wave amplitudes by more than one order of magnitude.

The dynamics of the plate is further investigated assigning as excitation pulse $21$ sinusoidal cycles centred at $27.5$ kHz (i.e., inside the BG) and modulated by a Hanning window (Fig. b). The full wavefield reconstruction of the Von Mises stress after $780 \mu s$ from the excitation is reported in Fig. . Green and red colours correspond to zero and maximum stresses induced in the plate. The area surrounded by four rows of cross-like cavities is practically free from stresses. For this reason, the phononic plate with four rows of cross-like cavities is chosen for the experimental phase of the study. In particular, the cavities are distributed over a square frame of width $4 \times a = 80$ mm. An homogeneous area of $120 \times 120$ mm$^2$ is enclosed by the annular phononic crystal region.

# Experimental measurements and full wavefield reconstruction

The sample used for experimental analysis consists of a standard $1000 \times 500 \times 12$ mm$^3$ (PVC) plate subjected to a machining process whereby $160$ hollow rounded cross-cylinder inclusions were drilled, as shown in Fig. a. The cavities cover a total area of $64$ cm$^2$ and are arranged in $4$ concentric rings which leave a homogeneous $120 \times 120$ mm$^2$ PVC region at the centre. The manufacturing process requireds a tolerance of $0.1$ mm, which slightly modifies the band structure of the system .

The adopted experimental set-up consists of a scanning laser Doppler vibrometer (SLDV) PSV 400 3D by Polytec (see Fig. b) which allows waveform digitization up to $70$ kHz, which is sufficient for detecting and visualizing ultrasonic waves in the frequency range of interest since the probe beam size is small enough to detect the motion at a single point on the surface of the phononic plate . Elastic guided waves are induced by means of a ceramic piezoelectric disk of diameter $10$ mm made of Sonox$^\text{{\textregistered}}$ by CeramTec$^\text{{\textregistered}}$ bonded to the surface of the investigated sample (point $A$ in Fig. a) using commercial super-glue. The pulse is excited by a TGA1241 generator by Thurlby Thandar Instruments through an EPA-104 amplifier by Piezo Systems$^\text{{\textregistered}}$ Inc, inducing a $200$ V$_{\text{pp}}$ signal. In order to improve measurement accuracy the investigated specimen is covered with self-adhesive retro-reflective film by ORALITE$^\text{{\textregistered}}$ enhancing the laser vibrometer signal level in each measurement point regardless of the incidence angle of the measurement beam.

The investigated spatial grid covers the entire phononic surface (see Fig. a) and consists of $500 \times 500$ equally spaced grid points. The laser vibrometer is perpendicularly positioned $1$ m away from the investigated surface and only one SLDV head is used to perform the out-of-plane measurements of the velocity amplitudes over the target area (Fig. b). Multiple ($512$) measurements are performed and averaged for each node, to filter out part of the noise.

A set of verification experiments is conducted to evaluate the elastic wavefield of the phononic plate at some frequencies of interest highlighted by the numerical study. Based on this analysis, sine functions with $21$ Hanning modulated sine cycles with central frequencies of $15$, $27.5$ and $42.5$ kHz are used as input signals (Fig. b-c) to reconstruct the wavefield for excitation below, within and above the BGs, respectively.

The excitation signal is used to trigger the data acquisition process. Out-of-plane velocity time-histories are acquired starting from laser beam reflections back to the scanning head and their Doppler frequency shift detected by interferometers. Upon completion of measurements at all grid points, the recorded responses are post-processed to obtain full images of the propagating wavefield aggregating single-point measurements within the region of inspection. The knowledge of the velocity time histories at all grid points (which can be integrated to obtain displacement time histories) allows the reconstruction of the time-evolving wavefields established in the scanning domain .

Figure shows the comparison of the experimental wavefield reconstruction and the corresponding calculated ones for the three aforementioned excitations (Figs. a-c). At frequencies below the BGs (Figs. a,b), Lamb waves travel through the periodic lattice relatively undisturbed. Limited reflection, scattering, or other losses can be observed, indicating that the phononic region does not interfere significantly with the wave motion.

When operating at a frequency within the BGs, instead, the destructive interferences due to the Bragg scattering occurring within this region becomes large. Figs. c,d clearly show that waves are reflected in the plate portion around the phononic region, mainly between the transmitting PZT and the lower edge of the unit cell ring. This behaviour is accompanied by an extremely low transmission due to the absence of detectable wave amplitudes inside the phononic region.

Above the BG, the SLDV again registers transmission inside the phononic region, allowing the wavefield reconstruction at the same intensity scale. However, in this case unit cells scatter the wave field, resulting in an observable delay in the wave propagation (Fig. e,f). In this case, despite the scattering, the phononic region does not cause significant attenuation of the wave field.

Thus, numerical and experimental results show excellent agreement. In particular, the clear attenuation of the wave field in the phononic region due to Bragg scattering can be observed in the case of a pulse centred within the BG, i.e. at $27.5$ kHz, confirming the presence of the numerically calculated BG. It is worth highlighting that as the excitation is provided at a single point, the waves impinge on the phononic region from a wide range of angles proving the "complete" nature of the BG.

# Conclusions

In this paper, the existence of a complete BG for Lamb waves propagating in a phononic plate is provided by means of full wave field measurements. Extensive numerical simulations are presented to optimize the attenuating performance of a $12$-mm thick PVC plate with $4$ rows of unit cells of lattice parameter $a = 20$ mm, characterized by through-the-thickness cross-like cavities arranged in a square ring. Ultrasonic measurements on an ad-hoc machined phononic PVC plate are performed via SLDV measurements and compared to the numerical results. Excellent agreement between the numerical prediction and experimental measurements is found both in terms of the BG frequency range as well as the decreasing trend of the transmission power spectrum. Besides confirming the existence of the numerically predicted complete BGs, measurements provide a direct observation of the scattering phenomena observed for waves characterized by frequency content below, within and above the BG.

This design shows promise due to its simplicity and effectiveness, and could be exploited at different size scales, e.g. at large scales for seismic shielding applications. In we have shown that for sandy type soils, similar structures in the range of about 10 m are sufficient for effective shielding at frequencies below 5 Hz. Experimental verification remains to be performed on scaled models of these shields. In future, it would be interesting to extend the experiments to other geometries of practical interest, and especially to weakly disordered phononic materials or stubbed phononic plates, in order to further extend the study of metamaterial geometries providing efficient attenuation over wide low-frequency ranges.

**Acknowledgements**

M.M. acknowledges funding from the European Union's Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement no. 658483. N.M.P. is supported by the European Research Council PoC 2015 "Silkene" No. 693670, by the European Commission H2020 under the Graphene Flagship Core 1 No. 696656 (WP14 "Polymer composites") and FET Proactive "Neurofibres" grant No. 732344. F.B. is supported by H2020 FET Proactive "Neurofibres" grant No. 732344.

**References**

|  Geometrical  | Length \[mm\] |             |              |
|:-------------:|:-------------:|:-----------:|:------------:|
| 2-4 parameter |  $case (i)$   | $case (ii)$ | $case (iii)$ |
|       a       |      20       |     20      |      20      |
|       h       |      12       |     12      |      12      |
|       b       |     14.48     |    12.83    |      18      |
|       c       |       /       |      /      |      3       |


abstract: A material which is an insulator entirely because of interaction effects is called a correlated insulator. Examples are trans-polyacetylene and the cuprate high temperature superconductors. Whereas doping of a band insulator results in a shift of the chemical potential into the conduction or valence band, doping of a correlated insulator produces fundamental changes in the electronic density of states itself. We have found that a general feature of doping a correlated insulator is the generation of topological defects; solitons in one-dimension and anti-phase domain walls in higher dimensions. We review the well known features of this process in polyacetylene, and describe the experimental evidence that the analogous features are seen in the cuprate superconductors. We also distinguish the case in which the doping-induced features can be viewed as a Fermi surface instability, as in polyacetylene, and the more usual case in which they are a consequence of a Coulomb frustrated electronic tendency to phase separation.
address: Dept. of Physics, UCLA, Los Angeles, CA 90024;  Dept. of Physics, Brookhaven National Laboratory, Upton, NY 11973
author: S. A.  Kivelson; V. J.  Emery
title:  Topological Doping of Correlated Insulators

= 10000

# Introduction

Perhaps the most intensively studied problem in condensed matter physics in recent years is the doping of a Mott insulator, a system with one electron per unit cell, which is insulating rather than metallic by virtue of the strong electron-electron repulsion. This long-standing unsolved problem has been revitalized by the discovery, almost a decade ago, of the cuprate high temperature superconductors, which are themselves doped Mott insulators. The scale of temperatures below which a system is insulating is set by the gap in the charge excitation spectrum; the interactions are sufficiently strong that this scale is much larger than the temperature of any transition to an electronically-ordered state. Nonetheless, at least on a bipartite lattice, such systems typically exhibit ordered ground states in which the size of the unit cell doubles. Of course, the most-studied example has a ground state with Néel, or antiferromagnetic order, as in the high temperature superconductors.

This problem is a subclass of the somewhat broader problem of doping a correlated insulator, which is relevant for a host of synthetic metals that have been studied in the past three decades. By a correlated insulator, we simply mean a material that would be a metal in a single-electron theory, but is rendered insulating by interaction effects. As with the Mott insulator, the commensurability between the density of conduction electrons and the underlying lattice is important. Frequently the unit cell contains one conduction electron or hole in the high-temperature phase, and is doubled in the broken symmetry state.

It is our contention that a certain more or less universal motif, which we have named "topological doping", characterizes the evolution of the electronic structure of such materials. The basic idea is that it is energetically favorable for the system to maintain the local commensurate insulating structure even when lightly doped, and for the additional charges to be incorporated as local defects. Moveover we find, quite generally, that the preferred defects are topological, in that the commensurate structure experiences a phase shift across the defect. When there is one conduction electron or hole per unit cell, the phase shift is equal to $\pi$. In one dimension such defects are $\pi$ solitons, while in higher dimensions they are the various forms of discommensuration, familiar from studies of the commensurate-incommensurate transition. It is a natural consequence of this view that, as the doping is increased, there is a crossover from a "doped insulator" regime, where the defects are far separated, to an "overdoped" or more or less conventionally metallic regime, when the spacing between defects is comparable to their width.

Many materials display the characteristics of topological doping, and it is useful to take a broad perspective in order to understand what they have in common. In this paper, we use the language of topological doping to describe the well known polyacetylene story, whose understanding was largely based on the work of A. J. Heeger and collaborators. We then show that analagous considerations govern the evolution of the electronic structure of the high temperature superconductors.

# Trans-polyacetylene

*The Undoped State:* Undoped trans-polyacetylene has a half-filled $\pi$ band on the order of 10 eV wide, and it would be a rather good metal were it not for the Peierls instability that causes it to dimerize; thus it is a prime example of a correlated insulator. The dimerized state is insulating and characterized by a broken symmetry, *i.e.* the dimerization itself, which can be observed as a Bragg peak in X-ray scattering. (The existence and magnitude of the static dimerization also has been inferred from fancy NMR techniques.) The undoped polymer exhibits a gap to both charge and spin excitations. The charge gap can be seen, for example, in the optical absorption spectrum and also has been measured by electrochemical means. The existence of a spin gap can be deduced from the vanishingly small value of the spin susceptibility, although its magnitude has been obtained only by indirect means, such as third-harmonic generation spectroscopy. The gap in the excitation spectrum of a single electron, which of course combines charge and spin, also has been measured by photoemission spectroscopy. All of these measurements are more or less consistent with a semiconducting model of the insulating state with a gap of order 1.8 eV. There are, however, two features of the insulating state that suggest the inadequacy of the semiconductor model: **i)** Long sub-gap tails, observed by for example optical absorption spectroscopy, suggest that there are excitations which have a substantially lower energy and are weakly accessible in the insulating state. These tails are now thought to be due to the direct photoproduction of soliton pairs, and are weak since they rely on the existence of virtual solitons as quantum fluctuations in the ground state. **ii)** There are remarkable photo-induced changes in the electronic structure as seen, for example, in photo-absorption experiments. These changes are known to be directly analagous to the dopant-induced changes in the electronic structure, and hence these experiments can be viewed as exploring the effects of photodoping the polymer.

*The Doped Insulator:* Upon doping, the conductivity of trans-polyacetylene rises by many orders of magnitude, signifying that there is no longer a gap (or at most a very small gap) in the charge excitation spectrum. This can be seen in more detail from the dopant-induced changes in the optical absorption spectrum. Here, it is found that, upon doping, the spectral weight of the features associated with the insulating gap is transferred continuously (and rapidly) to a broad "mid-gap" absorption feature, with substantial weight extending to much lower frequencies. At the same time, the spin susceptibility remains very small, suggesting that the gap in the spin excitation spectrum is not destroyed. As far as we know, the symmetry of the ground state is restored even at the lowest doping levels.

The basic explanation of these features was set forward in the seminal paper by W-P. Su, J. R. Schrieffer, and A. J. Heeger in 1979. They realized that doping does not proceed uniformly, as in a semiconductor, but rather it is energetically favorable to destroy the dimerization, and hence the gap, in local "defect" regions, associated with antiphase domain walls or discommensurations in the dimerization order. These are the famous solitons in polyacetylene. The local nature of the dopant-induced changes in the dimerization is responsible for most of the salient features of the doping of polyacetylene: **1)** The spin gap remains intact because the dimerization (hence the gap) persists between solitons. **2)** In the core of the soliton, the dimerization vanishes, hence there exists a zero energy (midgap) electronic bound state associated with each soliton. This mid-gap state leads to the shift of density of states into the middle of the insulating gap and hence to the associated peak in the optical absoption, with weight proportional to the dopant concentration. (The fact that all bound states are either fully occupied or unoccupied implies that there are no low-energy spin excitations associated with the soliton cores.) **3)** Low energy collective charge excitations are associated with the translational motion of the domain walls and appear in the far infrared absorption spectrum with oscillator strength proportional to the dopant concentration.

It is a more subtle matter to deduce the topological character of the dopant-induced defects directly from experiments. Soliton doping (as opposed to, for example, bipolaron doping) produces different distributions of states in the gap, although this is complicated by interaction effects and disorder. But, most importantly, the rapid restoration of the symmetry of the ground state upon doping is a consequence of the topological character of the defects; dopant-induced non-topological defects would produce only a modest reduction in the magnitude of the order parameter in proportion to the dopant concentration. One additional consequence of the topological character of the solitons is that they are anomalously one-dimensional, since a single soliton cannot tunnel between two polymer chains. This is a form of "confinement", and is a general feature of the one-dimensional electron gas in the spin gap regime. Confinement has not been directly observed experimentally, but it is fair to say that no evidence of three-dimensional quantum mechanical coherence has ever been observed in these materials, despite the fairly substantial magnitude of the wave function overlap between adjacent polyacetylene chains.

Of course, as a consequence of the dynamics of the solitons, the above picture of an inhomogeneous dopant-induced state is valid only at intermediate time scales and length scales. In the absence of disorder, and for weak enough interchain coupling, thermal and quantum flucutations of the soliton positions produce a "soliton liquid" state. Here, dimerization is completely absent in the static structure factor, but would be evident at intermediate times (or frequencies). In the more realistic case in which disorder plays a dominant role, the solitons are pinned at random positions, producing a "soliton glass" state. In either case, all configuration-averaged static correlation functions ultimately are translationally invariant or, equivalently, the dimerization order parameter is zero when averaged over distances long compared to the mean spacing between solitons.

Finally, we observe that the conductivity of lightly-doped polyacetylene, although enormous compared to the undoped polymer, is still considerably too small to be thought of as metallic in the usual sense, since it violates the so-called Ioffe-Regel criterion. Basically, the idea is that the smallest value of the conductivity consistent with Boltzmann transport theory occurs when the mean-free-path, $\ell$, is equal to the Fermi wavelength, $\lambda_F=2\pi/k_F$, and this constitutes a sort of minimum conductivity, $\sigma_{m}$, for a conventional metal. Any material with a smaller conductivity cannot be understood in terms of the occasional scattering of coherently propagating quasiparticles. For a quasi one-dimensional material, such as polyacetylene, $\sigma_{m}=4(e^2/\hbar)
(n_{\perp}/k_F)\approx 1.5 m\Omega-cm$, where $n_{\perp}$ is the areal density of polyacetylene chains. In the entire doped-insulator regime, the conductivity at room temperature is substantially smaller than $\sigma_{m}$, and it decreases with decreasing temperature. This poor conductivity is a consequence of soliton doping. The fact that the solitons are confined to one dimension implies that they are easily pinned by disorder. This is exacerbated by their local character, in that the majority of the electrons remain more or less condensed in the insulating state, and hence cannot participate in screening the impurity potential. If, by contrast, doping supressed all dimerization entirely, the room temperature conductivity would be large, as it is in the overdoped regime.

*The Overdoped Metal:* All of the above properties are features of the lightly-doped polymer which one can think of as a "doped insulator". Clearly, there is a characteristic dopant concentration $y_c$, which is roughly equal to the reciprocal of the soliton width, at which the defects begin to overlap. For higher dopant concentrations, one expects a crossover or transition to new physics. This crossover is clearly evident in the abrupt destruction of the spin gap above a critical dopant concentration which is of order $y_c \approx 5\%$, although it depends on the dopant species. The overdoped polymer behaves grossly like a three-dimensional Fermi liquid: **1)** It has a Pauli susceptibility and a T-linear term in the specific heat of roughly the magnitude that would be expected from band theory. **2)** Its D.C. conductivity can be as large as that of typical metals, and certainly much larger than $\sigma_{m}$. **3)** It exhibits three dimensional coherent transport, as seen for instance in the measured negative magnetoresistance .

The essential known features of the topological doping of trans polyacetylene are summarized in Table I.

# Coulomb frustrated electronic phase separation

In a series of previous publications, we have worked out the theory of topological doping of correlated insulators in higher dimensions, with special emphasis on the high temperature superconductors . While the ideas are quite general, we shall use language specific to the high temperature superconductors in summarizing the basic ideas, *i.e.* we will take the commensurate insulating state to be an ordered spin 1/2 Heisenberg antiferromagnet, and the doping to be p-type, *i.e.* hole doping:

**1)** Given a short-range model of a doped antiferromagnet, such as the $t-J$ model, we have shown that the most common behavior at low doping concentrations is phase separation. That is to say, the system separates into two macroscopically distinct regions: a commensurate ordered insulating region and a hole-rich metallic region. It is important to recognise that the driving force for this electronic phase separation is kinetic; it derives from the fact that the hole motion is frustrated in the antiferromagnetic state. Thus, phase separation is the best compromise between the hole kinetic energy and antiferromagnetic ordering. Since superexchange is itself a kinetic process, the strongly-correlated system has chosen phase separation as a means of minimizing the total zero-point energy. Elsewhere we have reviewed the set of model systems which can be proven to exhibit phase separation, and the experimental evidence that this physics is operative in the high temperature superconductors. Among other things, high temperature superconductors exhibit macroscopic phase separation whenever the dopants are mobile enough to counterbalance the effect of the long-range Coulomb interaction, as in the case of photodoping or oxygen doping. In other words "turning off" the long-range Coulomb interaction reveals the microscopic tendency to phase separation in these materials.

**2)** In the presence of long-range Coulomb interactions, the short-range tendency to phase separation is frustrated. So long as the Coulomb interaction is not too strong, this always results in a state which is inhomogeneous on intermediate time scales or length scales, *i.e.* the material forms a state consisting of a mixture of local undoped antiferromganetic regions and heavily-doped metallic regions. While many mesoscopic motifs can result from this physics, we have shown that the typical consequence is an ordered "stripe crystal" or fluctuating "stripe liquid" phase. Within the stripe liquid phase, there is also the possibility of a "nematic stripe liquid" phase in which the stripes are positionally disordered, but are, on average, oriented along a particular crystallographic direction, thus spontaneous breaking the four-fold rotational symmetry of the crystal. The stripes are antiphase domain walls in the antiferromagnetic order, and the doped holes are concentrated in the domain walls, where the antiferromagnetic order vanishes.

# Concerning the microscopic origin of the stripes

It is necessary to distinguish the various microscopic routes to stripe phases, since they have rather different consequences. In one mechanism, stripes (*i.e.* long period charge-density waves) arise from Fermi-surface nesting in a weakly incommensurate system . So far as is known, in the Hartree-Fock approximation which captures the essence of the Fermi surface instability, the phase separated state and stripes give the best variational solutions of the two-dimensional Hubbard model at small values of doping.

The Hartree-Fock stripes are, in fact, extremely close analogues of the solitons in polyacetylene, in that they are a lattice of antiphase domain walls whose mean separation is determined self-consistently from the hole concentration so that the mid-gap bands are completely empty. It is, of course, a general feature of a Fermi-surface instability that the order develops in such a way as to open gaps on the Fermi surface. Moreover in the Hartree-Fock solution, the instability is spin-driven, in the sense that there is a single transition temperature below which the broken symmetry solution of the Hartree-Fock equations is stable, and that the magnitude of the spin-order parameter $m$ turns on with the usual mean-field exponent, $m \sim (T_c-T)^{1/2}$, while the charge modulation turns on more weakly, (**i.e.** $\sim m^2$). In short, there are three clear features of stripes that arise from a Fermi surface instability: **i)** The transition is spin driven. **ii)** In the low temperature phase, there are gaps (or pseudogaps) on the Fermi surface. **iii)** The spacing between domain walls is equal to $1/x$, where $x$ is the hole concentration. (We have checked, by explicit calculation, that inclusion of long-range Coulomb interactions in the Hartree-Fock caclulation has little effect on this result .) **iv)** Although it is not easy to quantify, the notion of a Fermi-surface instability requires that the high temperature phase should be a Fermi liquid, with quasiparticle scattering rates no greater than the ultimate low-temperature energy gap.

The situation is quite different if the stripes arise from Coulomb-frustrated phase separation: **i)** One typically expects the transition to be charge-driven . That is to say, local spin order between the anti-phase domain walls can only develop *after* the holes are expelled from these regions. In this case, general Landau-Ginzburg considerations, lead one to expect either a first order transition, in which spin and charge order turn on simultaneously, or a sequence of transitions, in which first the charge order and then the spin order appears as the temperature is lowered. **ii)** Since the stripe concentration is determined primarily by the competition between the Coulomb interaction and the local tendency to phase separation, we do not necessarily expect that the spacing between stripes should be a simple function of $x$ and, as a consequence, there is no reason to expect the Fermi energy to lie in a gap or psuedogap; in fact, the Fermi energy typically will lie in a region of large density of states. **iii)** A high temperature Fermi liquid phase is not a prerequisite, since it is not an essential part of the physics of frustrated phase separation.

Although the formation of stripes is driven by frustrated phase separation, their conformation is influenced by other more subtle effects. **i)** As with any other charge-density wave system, commensurability energies are likely to be important. For instance, in the 214 family of superconductors (discussed below) either the stripe separation is nearly equal to 4 lattice spacings or the concentration of holes along a stripe is near to one per two sites. These two commensurabilities are compatible at $x=1/8$; otherwise the system must choose which one (if either) to satisfy. Commensurability effects have been explored in some detail in doped La$_2$NiO$_4$ by Tranquada and co-workers . Since systems of this sort can easily support a disorder line, commensurability effects can be important even in a stripe liquid phase. **ii)** The underlying crystalline lattice possesses, at best, fourfold rotational symmetry, not full rotational symmetry. This gives rise to a lattice anisotropy energy which strongly affects the ordering of stripe phases, and can also stabilize a nematic stripe-liquid phase. Moreover, the strength of the commensurability and anisotropy energies may be strongly influenced by structural phase transformations .

# The La$_2$CuO$_4$ (214) family

The "high temperature superconductors" are a group of materials which contain essentially equivalent CuO$_2$ planes separated by a variety of interstitial charge reservoir regions which are the source of the doping. It is generally believed that the basic physics is the same for all these materials, although subtle differences sometimes make it difficult to distinguish the essential and universal from the inessential or material dependent.

We choose to focus on the so-called 214 sub-family of these materials. Not only are they among the best studied but neutron scattering structure factors provide direct "photographic" evidence of the topological character of the doping process. In particular, we will contrast two slightly different materials of the family: La$_{2-x}$Sr$_x$CuO$_4$ and La$_{1.6-x}$Nd$_{0.4}$Sr$_x$CuO$_4$. In both cases, the undoped system ($x=0$) is an antiferromagnetic insulator and the concentration of doped holes in the CuO$_2$ planes is equal to the Sr concentration, $x$. However, at the temperatures and dopant concentrations of interest to us here, La$_{2-x}$Sr$_x$CuO$_4$ has an orthorhombic structure, whereas La$_{1.6-x}$Nd$_{0.4}$Sr$_x$CuO$_4$ undergoes a structural phase transition to a low temperature tetragonal (LTT) phase below about 70K. Above 70K, the electronic properties of the two materials are essentially indistinguishable for the same dopant concentration $x$ but the LTT phase supresses superconductivity at dopant concentrations $x < 20\%$ and, as has been recently discovered, stabilizes an ordered "striped" state, in which dopant-induced antiphase domain walls crystallize.

*The Undoped State:* The undoped state of all the high temperature superconductors is a quasi two-dimensional antiferromagnetic insulator. At all temperatures below 1000K, there are substantial antiferromagnetic correlations, although the system does not finally order into a Néel state until a temperature of the order of 400K. Since the low-energy charge excitations involve the movement of charge from copper to oxygen, the material is known as a charge-transfer insulator. However, in common with a Mott insulator, the gap in the charge excitation spectrum (roughly 2eV) is much greater than the scale of magnetic ordering .

*The Doped Insulator:* A breakthrough in understanding the properties of the doped insulating state was achieved recently through neutron scattering studies of a crystal of La$_{1.6-x}$Nd$_{0.4}$Sr$_x$CuO$_4$ with $x=0.12$. As the temperature is lowered, there is a succession of transitions to 1) the LTT phase, 2) a charged-ordered state (charge-stripes) and 3) at a slightly lower temperature, a period-doubling magnetically-ordered state. This is a photograph of an ordered stripe phase of the sort predicted above. The factor of two between the periods of the charge and spin orders implies that the charge is localized in magnetic antiphase domain walls, and the fact that the charge order appears first implies that the physics is driven by charge (*i.e.* frustrated phase separation), rather than by Fermi surface or Hartree-Fock physics. Moreover the charge density along the stripes corresponds to a half-filled band (one doped hole per unit cell containing two Cu ions), which is consistent with the mobility required by kinetic phase separation and commensurability with the underlying lattice, but half as large as the prediction of Hartree-Fock theory. Therefore, the origin of the stripe crystal phase of this material may be attributed uniquely to frustrated phase separation.

La$_{1.6-x}$Nd$_{0.4}$Sr$_x$CuO$_4$ with $x=0.12$ is thought not to exhibit any bulk superconductivity. In a manner that we believe we understand well , the LTT phase stabilizes an ordered stripe phase and, the same time, supresses superconductivity. That the ordered stripe phase is in direct competition with superconductivity already suggests that it is an extremely important piece of the physics. However, there is more direct evidence. La$_{1.6-x}$Nd$_{0.4}$Sr$_x$CuO$_4$ (with $x=0.12$) has static magnetic order below about 50K, but La$_{2-x}$Sr$_x$CuO$_4$ (with $x=0.15$) is magnetically disordered at all temperatures. Nevertheless the dynamical neutron scattering structure factors of both materials, measured at a frequency $\omega =$ 3 meV and T=40K (which effectively take a snapshot of the magnetic structures of the systems) are somewhat broadened versions of the elastic (Bragg) scattering observed in La$_{1.6-x}$Nd$_{0.4}$Sr$_x$CuO$_4$. Thus, there can be no question that there exist substantial dynamic stripe fluctuations in the optimally superconducting samples of La$_{2-x}$Sr$_x$CuO$_4$. Moreover, the fact that these fluctuations appear at such low frequencies, implies that they are very nearly frozen, even in the superconducting material!

This evidence is merely the latest and most graphic demonstration of topological doping in the high temperature superconductors. It is a dramatic confirmation of our interpretation of other experiments, as discussed in our earlier reviews, which includes analogues of most of the properties of polyacetylene, discussed above. A noteworthy property of the lightly-doped superconducting materials (known in the field as "underdoped") is the existence of a spin gap in the normal state, which we have associated with local pairing correlations without global phase coherence. Recent photoemission experiments appear to confirm this interpretation. The only known ways of realizing this situation are in highly anisotropic systems (*e.g.* quasi one-dimensional where the superconducting phase coherence is suppressed by thermal and quantum fluctuations) or in granular, spatially inhomogeneous, systems in which global superconductivity is suppressed by the small Josephson coupling between grains. To some extent, both of these features are present in the high temperature superconductors.

*The Overdoped Metal:* It is conventional to classify the high temperature superconductors as "optimally doped" (*i.e.* having the highest transition temperature in a given class of materials), "underdoped" (*i.e.* having a lower dopant concentration than optimal), and "overdoped" (*i.e.* having a larger than optimal dopant concentration). By now, the characteristics of the members of these three classes are sufficiently well known that it is possible to make an assignment to a given material, even if systematic studies of the transition temperature as a function of dopant concentration are unavailable. We have identified "underdoped" materials as doped insulators, in the sense of this article, and have suggested that "optimally doped" is simply a crossover between underdoped and overdoped. As in polyacetylene, the overdoped material appears to be more conventionally metallic. Certainly the conductivity at all temperatures is much greater than $\sigma_{m}$, and it is much more isotropic (*i.e.* three dimensional) than in under and optimally doped materials, although it still has a strange temperature dependence. Moreover, the magnetic properties appear to be those of a conventional metal. Thus, while it is not certain that the overdoped state can be classified as a Fermi liquid, it surely is a more conventional metallic state than at lower doping.

The essential features of topological doping in the 214 family of high temperature superconductors are summarized in Table II, which includes some additional properties discussed in our published papers.

# Other Cuprate Superconductors

The existence of charge-density wave fluctuations adds a new dimension to discussions of the physics of the high temperature superconductors which, almost from the outset, have been dominated by the interplay of antiferromagnetism and superconductivity. It is clear that the three phenomena are simply different manifestations of a common underlying physical theme, and that the best way to investigate any one of them is to find the material in which it shows up the most clearly. Nevertheless, although charge-density fluctuations are most clearly manifested in the 214 family, it is natural to ask how their effects might be revealed in other materials. Evidently it is desirable to have a space-sensitive probe and, in the absence of suitable neutron or X-ray scattering data, the obvious technique is angle-resolved photoemission spectroscopy (ARPES), which, until recently, was most reliable for Bi$_2$Sr$_2$CaCu$_2$O$_{8+x}$. We have undertaken a detailed study of the implications of ordered and fluctuating charge stripes for the single-particle properties of high temperature superconductors, and have found that it is possible to understand the most striking features of the ARPES data on Bi$_2$Sr$_2$CaCu$_2$O$_{8+x}$, especially the shape of the Fermi surface and the existence of nearly dispersionless states in its neighborhood. The essential ingredient in this interpretation of the data is a background of slowly-fluctuating stripes whose dynamics is determined by collective effects (the competition between phase separation and the long-range Coulomb force), rather than the single-particle behavior.

Recently a detailed analysis of the ARPES experiments has led to an understanding of the magnitude and temperature-dependence of the width $\Gamma$ of the single-particle spectral functions, with two important consequences. First of all, in the normal state, $\Gamma$ is too large to be consistent with any kind of quasiparticle picture; this observation strongly supports our analysis of the resistivity of the normal state , which was based on the violation of the Ioffe-Regel condition . Secondly, it was found that $\Gamma$ decreases rapidly as the temperature falls below T$_c$, indicating that quasiparticles are resurrected in the superconducting state. In our picture, quasiparticle propagation is inhibited in the normal state by the charge inhomogeneity, and especially by the frustration of hole motion in the direction perpendicular to the stripes, but phase order restores coherent particle motion, even though the stripes persist in the superconducting state.

# What does this have to do with high temperature superconductivity?

Our discussion has scarcely mentioned superconductivity which, after all, provoked the worldwide interest in the oxide superconductors. But, as might be expected, an appreciation of the electronic structure and the physics that gives rise to it also leads to a natural mechanism of high temperature superconductivity. The fundamental driving force is the frustration of the motion of holes in an antiferromagnet. This implies that when the holes on the stripes tunnel into the intervening region, which has a low charge density and strong antiferromagnetic correlations, they experience the very same effective attraction that would have caused them to coalesce (phase separate) in the absence of the longer range Coulomb interaction. Moreover, the Coulomb force, which so effectively inhibits pairing with a short coherence length and yet is widely ignored in discussions of the mechanism of high temperature superconductivity, has been incorporated into the picture from the outset. In fact, the holes moving on a stripe form a quasi one-dimensional electron gas in an "active" environment with low-energy excitations in the spin (and possibly charge) degrees of freedom, which renormalize the kinetic energy of the holes and mediate effective interactions between them . We have investigated this problem as a generalization of the standard theory of the one-dimensional electron gas, and have shown that it provides a possible microscopic explanation of the phenomenon of high temperature superconductivity . Whatever the mechanism, it must survive (and even thrive) on significant charge-density wave (stripe) fluctuations, and account for their intimate competition with superconductivity. In this regard, we find the mechanism described above to be an attractive possibility.

Empirically and on general theoretical grounds , it is clear that the optimal place for high temperature superconductivity is in the crossover regime between the doped insulator and overdoped regimes. This suggests a general strategy for looking for superconductivity in other doped correlated insulators.

Acknowledgements: SK was supported in part by NSF grant \#DMR -93-12606. This work also was supported by the Division of Materials Sciences, U.S. Department of Energy, under contract DE-AC02-76CH00016. We have enjoyed a continuing interaction with Dr. John Tranquada on the physics of oxide materials, and especially on charge ordering.

| **General Property** | **Particular Manifestation** |
|:--:|:--:|
| **Correlated insulator** | Commensurate (dimerized) CDW |
|  | Spontaneously broken reflection |
|  | and glide plane symmetry |
|  | Charge gap |
|  | Spin gap |
|  | Single-particle gap |
| **Topological Doping** |  Soliton liquid or glass |
|  | Rapid restoration of symmetry |
| Doped insulator for $y < y_c$ |  $y_c \sim$ 1/ (soliton width) |
|  | Spin gap survives doping |
|  | Low-energy optical oscillator strength proportional to $x$ |
|  | Charge gap feature persists with reduced oscillator strength |
| Inhomogeneous on intermediate | Transfer of oscillator strength to "mid-gap" states |
|  scales | Pinning of chemical potential in the gap |
|  | Conductivity $\ll \sigma_{m}$ |
|  |  |
|  |  |
|  |  |
| **"Overdoped" metallic state** | 3D metallic state |
| **for $y > y_c$** | Pauli spin susceptibility |
|  | 3D Coherent quasiparticle propagation |
|  | Metallic magnitude of the conductivity |
|  | Density of states at E$_F$ consistent with band theory |

Topological Doping in Polyacetylene (CH X$_y$)$_x$

| **General Property** | **Particular Manifestation** |
|:--:|:--:|
| **Correlated insulator** | Commensurate (Néel) SDW |
|  | Spontaneously broken spin rotational |
|  |      and translational symmetry |
|  | Charge gap |
|  | Spin waves |
|  | Single particle gap |
| **Topological doping** |  Stripe liquid, glass, or crystal |
|  | Rapid restoration of spin-rotational symmetry |
| Doped insulator for $x < x_c$ |  $x_c \sim$ 1/ (stripe width) |
|  | Spin gap appears at intermediate doping |
|  | Low-energy optical oscillator strength proportional to $x$ |
|  | Charge gap feature persists with reduced oscillator strength |
| Inhomogeneous on intermediate | Transfer of oscillator strength to "mid-gap" states |
|  scales | Pinning of chemical potential in the gap |
|  | Conductivity at high T $< \sigma_{m}$ |
|  | Ordered striped phase in LTT structure |
|  | Macroscopic phase separation when dopants are mobile |
|  | Fluctuating stripes observed in neutron scattering |
| **"Overdoped" metallic state** | 3D metallic state |
| **for $x > x_c$** | Destruction of spin gap |
|  | 3D Coherent quasiparticle propagation |
|  | Metallic magnitude of the conductivity |

Topological Doping in La$_{2-x}$Sr$_x$CuO$_4$


abstract: A Prolog-based framework for fully automated verification currently under development for heap-based object-oriented data is introduced. Dynamically allocated issues are discussed, recent approaches and criteria are analysed. The architecture and its components are introduced by example. Finally, propositions to further and related work are given.
author:
title: Dynamically Allocated Memory Verification in Object-Oriented Programs using Prolog

# Introduction

The main interest of this work is dedicated to the correctness of a program according to its memory consumption behaviour. It may, however, also be extended to performance considerations based on results coming from the dynamic memory verification, particularly but not only during the alias analysis or during the garbage collection phase, for instance. A Prolog-based verifier is suggested.

The structure of this paper is to give a short overview on current approaches in section 2. Section 3 presents the languages being used as programming and specification languages, and introduces the overall architecture. The architecture is designed to be open. Currently the project is under progress, the final section gives an outlook on related and future work.
discusses why the aliasing issue has still not been solved yet. provides a more detailed mark on the issues related to aliasing issues in a commercial Unix-environment. Despite its age and its partial closeness because of the commercial background, it is still often cited in very recent publications and numerous technical reports on the same topic, and in open-source projects, and surprisingly enough most of the issues found initially are still there almost with no changes. One key aspect that makes aliasing such a hard issue is that its local changes in a program listing may effect other regions unexpectedly –– but at the same time this is its strength, since no additional copying is required. had particularly analysed previously fixed bugs for a long-term period over last releases and found out - according to the bug distribution over time - that undesired memory behaviour is one of the most expensive bug reasons in terms of time and efforts to locate and fix.
Object-orientation , based on concepts such as encapsulation, polymorphism and inheritance, has been one of the most successful and widely adapted programming paradigms by now for a long time in industry, hence its combination with pointer structures remains a relevant research task up to date .
Prolog is considered for program verification for several reasons. First, it is a logical and declarative programming language which offers a high abstraction in writing Horn-clauses as they correspond with defined axioms and rules. A proof tree occasionally insists on a back-tracking strategy, which Prolog supports for free as one of its core-language features. The hope is Prolog's *generate-and-test* goal strategy may be found useful in simplifying and abstracting a proof significantly. Second, programs and internal states can be represented as terms. Terms can be easily processed in Prolog. The hope is, abduction and general symbolic term evaluation will allow generating lemmas more efficiently and make the reasoning terminate and terminate earlier. As a previous successful verification attempt, shall be noted. The authors solved a fairly hard problem from mathematical numerics using Prolog elegantly and straight. Since proofs are rule-centric, a proof contradiction will eventually help generating counter-examples easily by simply matching terms from the memory and from a rule or axiom.
**Example1** –– memory leak

    MyClass object1=new MyClass();
    ...
    object1=new MyClass();

**Example2** — unachievable memory

    // object1 has been created
    MyClass object2=new MyClass();
    object2.ref=object1;

**Example3** –– invalid memory access

    // object1.ref==null
    value = (object1.ref).attribute1;

**Example4** –– data structure with cycle

    object1.next=object1;
    ...
    root=object1;
    while(root.next!=null){
      printf("%d", object.data);
      root=root.next;
    }

Example 1 demonstrates the case where a fresh memory region is allocated, and without freeing it, allocates it again, which may cause the previously created region becoming unachievable. The second example demonstrates where `object2` is linked to an occupied `object1`, but `object2` itself remains unused. A very common problem in practice might be the third example, when an object reference is not set, but later referenced causing either an abnormal runtime failure or continues execution, which might be even worse in realistic scenarios because the further program execution becomes totally unpredictable with invalid value settings. The forth example might not immediately be seen as a problem, but if, for whatever reason, there is a cycle in `root` the program will not terminate. Apart from direct consequences like crashes or non-termination, one more side effect is there are spontaneous allocations/deallocations taking place on runtime which may eventually become a performance bottleneck.

# Current Approaches

One approach to verify correctness of dynamic memory is to get the program run and to record all memory cells that will be referenced and allocated/deallocated. This is what the *valgrind* tool does. This open-source tool requires on compilation guarding memory-checking code is injected to the assembly code. Not only that the enhanced program runs with huge delays, the general problem underneath this approach as well as SAFECode , which on runtime checks whether programmer-inserted assertions are fulfilled, is that only a small subset of all possible execution paths can be tested and that it requires additional code is inserted. For this reason, only static approaches are considered further that analyse the incoming program listing prior to running it.
In order to address the problems mentioned in the introduction part, several approaches exist: (i) Shape-based Analysis , (ii) Separating heap. For sake of completeness a heap-free alternative proposed by (iii) Tufte and Talpin and Meyer shall be mentioned, who both appeal to a stack-based approach, if any possible, to avoid expensive heap allocation and deallocation operations. Automatic handling of stack-based locations is essential for both, where the stack sizes are determined during compilation. Thus, control passing which happens during a call will allow to allocate objects almost for *free*, since a stack frame needs to be created in any case, and no more expensive heap operations are needed in fact, for instance garbage collection. The disadvantage of (iii) is, however, there is often a platform-dependent restriction on stack sizes and number of entries, so in practise there are tough frame restrictions, e.g. a maximum offset, which should not be exceeded without getting a severe performance penalty on concrete target architectures at the same time. Meyer asks to turn garbage collection steadily on during program execution. This implies for efficient execution runtime critical parts will not trigger dynamic memory operations and those operations are opted out as efficiently as it can possibly be done.

Approaches (i) and (ii) are similar, both describe the memory state, although (i) describes the entire dynamic heap as an entire graph, where edges are region dependencies and vertices are locations. The problem with (i) is *locality*, because if a particular function is called, the entire graph has to be specified before, after and during the call, where approach (ii) allows to hide all non-affected heaps (*framed heaps*) – this is what is meant by *locality principle* in terms of *Separation Logic* . The most important concept behind Separation Logic is the specification of two non-interleaving memory regions. Heaps might be composed, and programs may change heaps. If two heaps are connected, then the dependency has to be added explicitly to a current heap's specification. If a heap depends on some other heap data, then this is called *aliasing* (or "*big brother property*" as found in ).

The first implementation which makes use of Separation Logic is Smallfoot . In order to extend deductive reasoning capabilities, an abductive approach was proposed, called bi-abduction , for Separation Logic, which is a constructive guess of unchanged heaps by a greedy symbolic table-construction algorithm that chooses bigger rules first. The extension of Smallfoot is called *SpaceInvader*.

Hurlin extends in the classic Separation Logic proposed in by classes for a Java-like language. He suggests a *heap factorization*, an attempt to normalise heaps in order to remove redundant heap specification fragments which are considered as *noise*, even if the main goal of his thesis is focused on multi-threaded applications. He re-uses the same concept of *abstract predicates* as it was introduced by , and generates unchanged parts during deduction with a parallel algorithm.

Parkinson introduces a Java-like language with object-orientated features. Nevertheless, many problems are not being addressed yet: abstraction mismatch on encapsulation and inheritance, particularly, the problem of expanding specifications in subclasses seems to be a real hinder in simple and elegant specifications. The most essential contribution of is the introduction of *abstract predicates*, although there are currently tough restrictions concerning expressibility. Super calls, static fields, reflection, inner classes and quantified predicates, for example, are currently missing language features.

Verifast is another forward verifier based on Separation Logic. In comparison to all previously introduced verifiers which do very similar operations on the heap, all introduced conventions per tool differ strongly, and it does not automatically process loop invariants nor predicates -– here it depends entirely on user-interaction or requires explicit injections within specification annotations inside the program which are used as internal reorganisation commands.
In objects as class-instances are treated as records, typing and verification rules are introduced and a soundness proof is provided. Problems which neither in nor are addressed are that objects may have references to other objects and that a lack in abstraction causes a dramatic increase in specification length which makes it in practice impossible to read and understand specification to a full extend. The memory state is specified by temporal predicates and a result-register for the previous computation step's result. There is no general recursive definition allowed, although attempts to relax this hard restriction by an algebraic ideal-construction. Still there are hard restrictions, such as no aliasing nor late binding at all, and object-records only which even may become unsound for eager type evaluation.

# Architecture and Design

Before going on with more details on the architecture, the prerequisites on the architecture shall be summarised. The proposed architecture may be considered for teaching purposes in the future:

1.  Automatic proof. The program and its annotations shall be sufficient in order to get the verification run. If there is an endless cycle in the proof however, there shall be no mandatory recognition, since termination is beyond the main focus of this work.

2.  Openness. The provided architecture shall be open for extensions and variability, and the attached models shall be exportable, so it might eventually be passed through to another arbitrary model transformer, if needed.

3.  Extensibility. The target language shall be fixed but interchangeable with an imperative programming language in the front-end. The rules and user-defined predicates shall be designed amendable, so the user may want to add rules directly to the rule set.

4.  Plausibility: There shall be configurable visualisation facilities, so the incoming annotated program may be retrospected on each stage of the verification process. If, for instance, a proof fails or stops abruptly the user would perhaps like to see the proof tree and a counter-example.

In figure the architecture of the Prolog-based verification system is shown. The input is a C-program with object-orientated extension that is annotated with assertions specifying the dynamic memory. The shortened syntax can be described by the Extended Backus-Naur form in figure . Not mentioned definitions as actual parameters, blocks, class methods, variable declarations have been skipped here for the sake of readability and follow mostly ANSI C. For readability purposes the expression sub-grammar has not been expanded according to its precedence hierarchy nor for optional assertions. `new` and `delete` reserve/free new chunks in the heap associated with previously defined locations. The access to heap memory is performed by `[<location>]`, where `<location>` denotes either a field variable, another object's field or a either of those with an offset in order specify non-aligned memory regions, for instance. The rule `<funcall>` denotes the syntax for a method call, which may have a object specifier optionally and a method name which is required to exist with the matching total number and types of parameters being passed as expressions.

C-programs are annotated by assertions which are injected as usual Prolog terms into blocks. Blocks are encoded as lists of statement-terms. Assertions are inductively defined and can be found in figure . Keep in mind the expression might request object references and $\alpha(\vec{p})$" assumes predicate named $\alpha$ was defined prior to using it, and $\vec{p}$ contains as many actual parameters as the arity of predicate $\alpha$ require there are.

For example,

    int f(int a, int b) @ a<10 @ {
      id=2; a=1; b=6;
    } @ a->5 * b->c * c->object(myClass1,15) @

is transformed into this Prolog-term:

    function(f, int,
      [param(a,int), param(b,int)],
      [assert(le(a,10)),
       assign(id,2), assign(a,1), assign(b,6),
       assert(a->5 * b->c *
              c->object(myClass1,15))])

An important specification fragment of the intermediate Prolog-term syntax can be found in figure , where the remaining part is close to the syntax of figure . Apart from *'ite'* which represents a *if-then-else*-construct with at least one block for the if-case and one more optional block for the else-block – as long as it was provided, there are also while-loops and further constructs, like class definitions, which will not be mentioned here for simplicity purposes. All expressions, particularly with binary operators, are encoded as terms where the literal operator becomes the functor, for example `add(i,7)`.

Since the architecture is designed flexible, it allows the user to interchange the compiler front-end for a different language, so the user has the possibility to write own Prolog-terms directly without even having an ordinary C-program. In this case syntax and semantic constraints remain on full responsibility to the user. Prolog-terms are internally checked and may also be directed to a graphical output, e.g. for proof tree visualisation. Antlr 4 is currently used as compilation front-end.
Once the Prolog term is constructed, it can be passed to the verification. Hereby, the term is now processed while the internal environment, which has to keep the states of the memory, needs to be updated after every statement. All locals are residing in stack, where dynamically allocated memory locations may remain in memory –– even if a stack-based variable stores a dynamic address it would be freed at the end of a block.

If we decide to specify a list concatenation of two lists, we have several opportunities to describe the heap. If `list(s,e)` denotes a heap predicate where `s` is the location of the beginning root element of a list, and `e` denotes the last element in that list, then having two lists `x` and `y` with `x->a,b,c` and `y->d,e,f` will concatenate for instance to either (i) `x->a,b,c,d,e,f * y->f` or to (ii) `x->a,b,c * y->d,e,f * z->a,b,c,d,e,f`. Remark: The ','-operator is defined as a list constructor with variable input amount for all consecutive objects currently in memory linked together to a simply-linked list . The main difference between (i) and (ii) is that (i) requires only a single assignment if the end of `x` is known, therefore `x` and `y` are no more as they used to be before concatenation. (ii) creates an entirely new copy of all element from `x` and `y` and does not touch neither `x` nor `y`. (ii) is safer from a general reuse perspective, but it is considerably slower and consumes more memory due to additional copies to be generated.

Finally, the SMT-solver is required whenever taking out trivial calculations, for instance in basic arithmetics. For instance, if there is an expression that might be reduced to a value, then this should in general be tried first before triggering a certain rule. Beside finding solutions to basic arithmetic and other theories, re-arrangement needs to be taken into consideration. Formal rules will usually also not deal too much about heap permutation, although a strategy must be found and is crucial in fact for the overall performance.

# Conclusion

So far the open architecture was presented providing several suggestions for further research activities. Prolog was proposed as specification and proof platform for memory-specific research, e.g. on extending the expressibility of abstract predicates or abduction. We believe, questions related to abduction in Separation Logic with objects still have not been profoundly investigated yet, as well as some object-oriented features like polymorphism in Separation Logic.

The platform might be used to incorporate with existing compiler packages in order to research improvement on code optimization during the alias analysis phase, but also garbage collection, based on knowledge obtained during the dynamic memory verification.

Further rules of normalisation and re-arrangement will be applied to cover more real world scenarios, particularly in order to resolve arithmetic equivalency by the integration of a SMT-solver ().

Related work includes Jacobs who suggests to investigate Banerjee's Regional Logic approach as substitute for the Symbolic Execution approach .