Search

The computational investigation of Fermi, Pasta, Ulam, and Tsingou of arrays of nonlinearly coupled oscillators has led to a wealth of studies in nonlinear dynamics. Most studies of oscillator arrays have considered homogeneous oscillators, even though there are inherent heterogeneities between {individual} oscillators in real-world arrays. Well-known FPUT phenomena, such as energy recurrence, can break down in such heterogeneous systems. In this paper, we present an approach -- the use of structured heterogeneities -- to recover recurrence in FPUT systems in the presence of oscillator heterogeneities. We examine oscillator variabilities in FPUT systems with cubic nonlinearities, and we demonstrate that centrosymmetry in oscillator arrays may be an important source of recurrence., Comment: revised version

We consider generalizations of nonlinear Schr\"odinger equations, which we call "Karpman equations", that include additional linear higher-order derivatives. Singularly-perturbed Karpman equations produce generalized solitary waves (GSWs) in the form of solitary waves with exponentially small oscillatory tails. Nanoptera are a special case of GSWs in which these oscillatory tails do not decay. Previous research on continuous third-order and fourth-order Karpman equations has shown that nanoptera occur in specific settings. We use exponential asymptotic techniques to identify traveling nanoptera in singularly-perturbed continuous Karpman equations. We then study the effect of discretization on nanoptera by applying a finite-difference discretization to continuous Karpman equations and studying traveling-wave solutions. The finite-difference discretization turns a continuous Karpman equation into an advance--delay equation, which we study using exponential asymptotic analysis. By comparing nanoptera in these discrete Karpman equations with nanoptera in their continuous counterparts, we show that the oscillation amplitudes and periods in the nanoptera tails differ in the continuous and discretized equations. We also show that the parameter values at which there is a bifurcation between nanopteron and decaying oscillatory solutions depends on the choice of discretization. Finally, by comparing different higher-order discretizations of the fourth-order Karpman equation, we show that the bifurcation value tends to a nonzero constant for large orders, rather than to $0$ as in the associated continuous Karpman equation., Comment: 14 figures, 42 pages, 4 tables

We study ``nanoptera'', which are non-localized solitary waves with exponentially small but non-decaying oscillations, in two singularly-perturbed Hertzian chains with precompression. These two systems are woodpile chains (which we model as systems of Hertzian particles and springs) and diatomic Hertzian chains with alternating masses. We demonstrate that nanoptera arise from Stokes phenomena and appear as special curves, called Stokes curves, are crossed in the complex plane. We use techniques from exponential asymptotics to obtain approximations of the oscillation amplitudes. Our analysis demonstrates that traveling waves in a singularly perturbed woodpile chain have a single Stokes curve, across which oscillations appear. Comparing these asymptotic predictions with numerical simulations reveals that this accurately describes the non-decaying oscillatory behavior in a woodpile chain. We perform a similar analysis of a diatomic Hertzian chain, that the nanpteron solution has two distinct exponentially small oscillatory contributions. We demonstrate that there exists a set of mass ratios for which these two contributions cancel to produce localized solitary waves. This result builds on prior experimental and numerical observations that there exist mass ratios that support localized solitary waves in diatomic Hertzian chains without precompression. Comparing asymptotic and numerical results in a diatomic Hertzian chain with precompression reveals that our exponential asymptotic approach accurately predicts the oscillation amplitude for a wide range of system parameters, but it fails to identify several values of the mass ratio that correspond to localized solitary-wave solutions., Comment: 29 pages, 9 figures

We conduct an extensive study of nonlinear localized modes (NLMs), which are temporally periodic and spatially localized structures, in a two-dimensional array of repelling magnets. In our experiments, we arrange a lattice in a hexagonal configuration with a light-mass defect, and we harmonically drive the center of the chain with a tunable excitation frequency, amplitude, and angle. We use a damped, driven variant of a vector Fermi- Pasta-Ulam-Tsingou lattice to model our experimental setup. Despite the idealized nature of the model, we obtain good qualitative agreement between theory and experiments for a variety of dynamical behaviors. We find that the spatial decay is direction-dependent and that drive amplitudes along fundamental displacement axes lead to nonlinear resonant peaks in frequency continuations that are similar to those that occur in one-dimensional damped, driven lattices. However, driving along other directions leads to the creation of asymmetric NLMs that bifurcate from the main solution branch, which consists of symmetric NLMs. When we vary the drive amplitude, we observe such behavior both in our experiments and in our simulations. We also demonstrate that solutions that appear to be time-quasi-periodic bifurcate from the branch of symmetric time-periodic NLMs.

In 1955, Fermi, Pasta, Ulam, and Tsingou reported recurrence over time of energy between modes in a one-dimensional array of nonlinear oscillators. Subsequently, there have been myriad numerical experiments using homogenous FPUT arrays, which consist of chains of ideal, nonlinearly-coupled oscillators. However, inherent variations --- e.g., due to manufacturing tolerance --- introduce heterogeneity into the parameters of any physical system. We demonstrate that such tolerances degrade the observance of recurrence, often leading to complete loss in moderately sized arrays. We numerically simulate heterogeneous FPUT systems to investigate the effects of tolerances on dynamics. Our results illustrate that tolerances in real nonlinear oscillator arrays may limit the applicability of results from numerical experiments on them to physical systems, unless appropriate heterogeneities are taken into account., Comment: 11 pages in main text + short supplement (2 extra figures and terse accompanying text)

We study quasiperiodicity-induced localization of waves in strongly precompressed granular chains. We propose three different setups, inspired by the Aubry--Andr\'e (AA) model, of quasiperiodic chains; and we use these models to compare the effects of on-site and off-site quasiperiodicity in nonlinear lattices. When there is purely on-site quasiperiodicity, which we implement in two different ways, we show for a chain of spherical particles that there is a localization transition (as in the original AA model). However, we observe no localization transition in a chain of cylindrical particles in which we incorporate quasiperiodicity in the distribution of contact angles between adjacent cylinders by making the angle periodicity incommensurate with that of the chain. For each of our three models, we compute the Hofstadter spectrum and the associated Minkowski--Bouligand fractal dimension, and we demonstrate that the fractal dimension decreases as one approaches the localization transition (when it exists). Finally, in a suite of numerical computations, we demonstrate both localization and that there exist regimes of ballistic, superdiffusive, diffusive, and subdiffusive transport. Our models provide a flexible set of systems to study quasiperiodicity-induced analogs of Anderson phenomena in granular chains that one can tune controllably from weakly to strongly nonlinear regimes., Comment: 16 pages

We study - experimentally, theoretically, and numerically - nonlinear excitations in lattices of magnets with long-range interactions. We examine breather solutions, which are spatially localized and periodic in time, in a chain with algebraically-decaying interactions. It was established two decades ago [S. Flach, Phys. Rev. E 58, R4116 (1998)] that lattices with long-range interactions can have breather solutions in which the spatial decay of the tails has a crossover from exponential to algebraic decay. In this Letter, we revisit this problem in the setting of a chain of repelling magnets with a mass defect and verify, both numerically and experimentally, the existence of breathers with such a crossover.

The study of energy transport properties in heterogeneous materials has attracted scientific interest for more than a century, and it continues to offer fundamental and rich questions. One of the unanswered challenges is to extend Anderson theory for uncorrelated and fully disordered lattices in condensed-matter systems to physical settings in which additional effects compete with disorder. Specifically, the effect of strong nonlinearity has been largely unexplored experimentally, partly due to the paucity of testbeds that can combine the effect of disorder and nonlinearity in a controllable manner. Here we present the first systematic experimental study of energy transport and localization properties in simultaneously disordered and nonlinear granular crystals. We demonstrate experimentally that disorder and nonlinearity --- which are known from decades of studies to individually favor energy localization --- can in some sense "cancel each other out", resulting in the destruction of wave localization. We also report that the combined effect of disorder and nonlinearity can enable the manipulation of energy transport speed in granular crystals from subdiffusive to superdiffusive ranges., Comment: main text + supplementary information

The study of granular crystals, metamaterials that consist of closely packed arrays of particles that interact elastically, is a vibrant area of research that combines ideas from disciplines such as materials science, nonlinear dynamics, and condensed-matter physics. Granular crystals, a type of nonlinear metamaterial, exploit geometrical nonlinearities in their constitutive microstructure to produce properties (such as tunability and energy localization) that are not conventional to engineering materials and linear devices. In this topical review, we focus on recent experimental, computational, and theoretical results on nonlinear coherent structures in granular crystals. Such structures --- which include traveling solitary waves, dispersive shock waves, and discrete breathers --- have fascinating dynamics, including a diversity of both transient features and robust, long-lived patterns that emerge from broad classes of initial data. In our review, we primarily discuss phenomena in one-dimensional crystals, as most research has focused on such scenarios, but we also present some extensions to two-dimensional settings. Throughout the review, we highlight open problems and discuss a variety of potential engineering applications that arise from the rich dynamic response of granular crystals., Comment: 41 pages, 11 figures, review article (comments welcome)

We study asymptotic solutions to a singularly-perturbed, period-2 Toda lattice and use exponential asymptotics to examine `nanoptera', which are nonlocal solitary waves with constant-amplitude, exponentially small wave trains. With this approach, we isolate the exponentially small, constant-amplitude waves, and we elucidate the dynamics of these waves in terms of the Stokes phenomenon. We find a simple asymptotic expression for the waves, and we study configurations in which these waves vanish, producing localized solitary-wave solutions. In the limit of small mass ratio, we derive a simple anti-resonance condition for the manifestation these wave-free solutions., Comment: 31 journal pages

We study scattering of waves by impurities in strongly precompressed granular chains. We explore the linear scattering of plane waves and identify a closed-form expression for the reflection and transmission coefficients for the scattering of the waves from both a single impurity and a double impurity. For single-impurity chains, we show that, within the transmission band of the host granular chain, high-frequency waves are strongly attenuated (such that the transmission coefficient vanishes as the wavenumber $k\rightarrow \pm\pi$), whereas low-frequency waves are well-transmitted through the impurity. For double-impurity chains, we identify a resonance --- enabling full transmission at a particular frequency --- in a manner that is analogous to the Ramsauer--Townsend (RT) resonance from quantum physics. We also demonstrate that one can tune the frequency of the RT resonance to any value in the pass band of the host chain. We corroborate our theoretical predictions both numerically and experimentally, and we directly observe complete transmission for frequencies close to the RT resonance frequency. Finally, we show how this RT resonance can lead to the existence of reflectionless modes even in granular chains (including disordered ones) with multiple double impurities., Comment: 12 pages, 9 figures

We study the spreading of initially localized excitations in 1D disordered granular crystals. We thereby investigate localization phenomena in strongly nonlinear systems, which we demonstrate to be fundamentally different from localization in linear and weakly nonlinear systems. We compare wave dynamics in chains with 3 different types of disorder: an uncorrelated (Anderson-like) disorder and 2 types of correlated disorders (random dimer arrangements), and for 2 types of initial conditions: displacement excitations and velocity excitations. For strongly precompressed chains, the dynamics depend strongly on the initial condition. For displacement excitations, the long-time behavior of the second moment $\tilde{m}_2$ has oscillations that depend on the type of disorder, with a complex trend that differs markedly from a power law and which is particularly evident for an Anderson disorder. For velocity excitations, we find a scaling $\tilde{m}_2\sim t^{\gamma}$ (for a constant $\gamma$) for all 3 types of disorder. For weakly precompressed (strongly nonlinear) chains, $\tilde{m}_2$ and the inverse participation ratio $P^{-1}$ satisfy $\tilde{m}_2\sim t^{\gamma}$ and $P^{-1}\sim t^{-\eta}$, and the dynamics is superdiffusive for all examined cases. When precompression is strong, the IPR decreases slowly for all 3 types of disorder, and we observe a partial localization around the core and the leading edge of the wave. For an Anderson disorder, displacement perturbations lead to localization of energy primarily in the core, and velocity perturbations cause the energy to be divided between the core and the leading edge. This localization phenomenon does not occur in the sonic-vacuum regime, which yields the surprising result that the energy is no longer contained in strongly nonlinear waves but instead is spread across many sites. In this regime, the exponents are very similar in all cases., Comment: 23 pages, 11 figures; the abstract in this website field is a shortened version of the paper's full abstract because of spatial limits

Force chains form heterogeneous physical structures that can constrain the mechanical stability and acoustic transmission of granular media. However, despite their relevance for predicting bulk properties of materials, there is no agreement on a quantitative description of force chains. Consequently, it is difficult to compare the force-chain structures in different materials or experimental conditions. To address this challenge, we treat granular materials as spatially-embedded networks in which the nodes (particles) are connected by weighted edges that represent contact forces. We use techniques from community detection, which is a type of clustering, to find sets of closely connected particles. By using a geographical null model that is constrained by the particles' contact network, we extract chain-like structures that are reminiscent of force chains. We propose three diagnostics to measure these chain-like structures, and we demonstrate the utility of these diagnostics for identifying and characterizing classes of force-chain network architectures in various materials. To illustrate our methods, we describe how force-chain architecture depends on pressure for two very different types of packings: (1) ones derived from laboratory experiments and (2) ones derived from idealized, numerically-generated frictionless packings. By resolving individual force chains, we quantify statistical properties of force-chain shape and strength, which are potentially crucial diagnostics of bulk properties (including material stability). These methods facilitate quantitative comparisons between different particulate systems, regardless of whether they are measured experimentally or numerically., Comment: 15 pages, 13 figures, and 2 appendices

Coupled map lattices (CMLs) are often used to study emergent phenomena in nature. It is typically assumed (unrealistically) that each component is described by the same map, and it is important to relax this assumption. In this paper, we characterize periodic orbits and the laminar regime of type-I intermittency in heterogeneous weakly coupled map lattices (HWCMLs). We show that the period of a cycle in an HWCML is preserved for arbitrarily small coupling strengths even when an associated uncoupled oscillator would experience a period-doubling cascade. Our results characterize periodic orbits both near and far from saddle--node bifurcations, and we thereby provide a key step for examining the bifurcation structure of heterogeneous CMLs.

We study statically homogeneous Bose-Einstein condensates with spatially inhomogeneous interactions and outline an experimental realization of compensating linear and nonlinear potentials that can yield constant-density solutions. We illustrate how the presence of a step in the nonlinearity coefficient can only be revealed dynamically and consider, in particular, how to reveal it by exploiting the inhomogeneity of the sound speed with a defect-dragging experiment. We conduct computational experiments and observe the spontaneous emergence of dark solitary waves. We use effective-potential theory to perform a detailed analytical investigation of the existence and stability of solitary waves in this setting, and we corroborate these results computationally using a Bogoliubov-de Gennes linear stability analysis. We find that dark solitary waves are unstable for all step widths, whereas bright solitary waves can become stable through a symmetry-breaking bifurcation as one varies the step width. Using phase-plane analysis, we illustrate the scenarios that permit this bifurcation and explore the dynamical outcomes of the interaction between the solitary wave and the step., Comment: 7 pages (published version), 4 figures

We study the structure, stability, and dynamics of dark solitary waves in parabolically trapped, collisionally inhomogeneous Bose-Einstein condensates (BECs) with spatially periodic variations of the scattering length. This collisional inhomogeneity yields a nonlinear lattice, which we tune from a small-amplitude, approximately sinusoidal structure to a periodic sequence of densely spaced spikes. We start by investigating time-independent inhomogeneities, and we subsequently examine the dynamical response when one starts with a collisionally homogeneous BEC and then switches on an inhomogeneity either adiabatically or nonadiabatically. Using Bogoliubov-de Gennes linearization as well as direct numerical simulations of the Gross-Pitaevskii equation, we observe dark solitary waves, which can become unstable through oscillatory or exponential instabilities. We find a critical wavelength of the nonlinear lattice that is comparable to the healing length. Near this value, the fundamental eigenmode responsible for the stability of the dark solitary wave changes its direction of movement as a function of the strength of the nonlinearity. When it increases, it collides with other eigenmodes, leading to oscillatory instabilities; when it decreases, it collides with the origin and becomes imaginary, illustrating that the instability mechanism is fundamentally different in wide-well versus narrow-well lattices. When starting from a collisionally homogeneous setup and switching on inhomogeneities, we find that dark solitary waves are preserved generically for aligned lattices. We briefly examine the time scales for the onset of solitary-wave oscillations in this scenario., Comment: 9 pages, 10 figures (with many subfigures), published in PRA

We present a systematic study of the existence and stability of discrete breathers that are spatially localized in the bulk of a one-dimensional chain of compressed elastic beads that interact via Hertzian contact. The chain is diatomic, consisting of a periodic arrangement of heavy and light spherical particles. We examine two families of discrete gap breathers: (1) an unstable discrete gap breather that is centered on a heavy particle and characterized by a symmetric spatial energy profile and (2) a potentially stable discrete gap breather that is centered on a light particle and is characterized by an asymmetric spatial energy profile. We investigate their existence, structure, and stability throughout the band gap of the linear spectrum and classify them into four regimes: a regime near the lower optical band edge of the linear spectrum, a moderately discrete regime, a strongly discrete regime that lies deep within the band gap of the linearized version of the system, and a regime near the upper acoustic band edge. We contrast discrete breathers in anharmonic FPU-type diatomic chains with those in diatomic granular crystals, which have a tensionless interaction potential between adjacent particles, and highlight in that the asymmetric nature of the latter interaction potential may lead to a form of hybrid bulk-surface localized solutions.

In this chapter, we discuss experiments that realize the discrete nonlinear Schr\"odinger (DNLS) equations. The relevance of such descriptions arises from the competition of three common features: nonlinearity, dispersion, and a medium to large level of (periodic, quasiperiodic, or random) discreteness in space. DNLS equations have been especially prevalent in atomic and molecular physics in the study of Bose-Einstein condensates in optical lattices or superlattices; and in nonlinear optics in the description of pulse propagation in waveguide arrays and photorefractive crystals. New experiments in both nonlinear optics and Bose-Einstein condensation provide new challenges for DNLS models, and DNLS and related equations have also recently been used to make important predictions in novel physical settings such as the study of composite metamaterials and arrays of superconducting devices., Comment: 17 pages, review article; (document is same as before, but I forgot to include the reference to the published version in Journal-ref); in P. G. Kevrekidis, Discrete Nonlinear Schr\"odinger Equation: Mathematical Analysis, Numerical Computations, and Physics Perspectives, 175-189, Springer Tracts in Modern Physics (Springer-Verlag)

We investigate nonlinear localized modes at light-mass impurities in a one-dimensional, strongly-compressed chain of beads under Hertzian contacts. Focusing on the case of one or two such "defects", we analyze the problem's linear limit to identify the system eigenfrequencies and the linear defect modes. We then examine the bifurcation of nonlinear defect modes from their linear counterparts and study their linear stability in detail. We identify intriguing differences between the case of impurities in contact and ones that are not in contact. We find that the former bears similarities to the single defect case, whereas the latter features symmetry-breaking bifurcations with interesting static and dynamic implications.

We investigate the propagation and scattering of highly nonlinear waves in disordered granular chains composed of diatomic (two-mass) units of spheres that interact via Hertzian contact. Using ideas from statistical mechanics, we consider each diatomic unit to be a "spin", so that a granular chain can be viewed as a spin chain composed of units that are each oriented in one of two possible ways. Experiments and numerical simulations both reveal the existence of two different mechanisms of wave propagation: In low-disorder chains, we observe the propagation of a solitary pulse with exponentially decaying amplitude. Beyond a critical level of disorder, the wave amplitude instead decays as a power law, and the wave transmission becomes insensitive to the level of disorder. We characterize the spatio-temporal structure of the wave in both propagation regimes and propose a simple theoretical interpretation for such a transition. Our investigation suggests that an elastic spin chain can be used as a model system to investigate the role of heterogeneities in the propagation of highly nonlinear waves., Comment: 10 pages, 8 figures (some with multiple parts), to appear in Physical Review E; summary of changes: new title, one new figure, additional discussion of several points (including both background and results)

We provide a quantitative characterization of dissipative effects in one-dimensional granular crystals. We use the propagation of highly nonlinear solitary waves as a diagnostic tool and develop optimization schemes that allow one to compute the relevant exponents and prefactors of the dissipative terms in the equations of motion. We thereby propose a quantitatively-accurate extension of the Hertzian model that encompasses dissipative effects via a discrete Laplacian of the velocities. Experiments and computations with steel, brass, and polytetrafluoroethylene reveal a {\em common} dissipation exponent with a material-dependent prefactor., Comment: 4 pages, 4 figures, accepted in Phys. Rev. Lett

Motivated by recent experiments in optics and atomic physics, we derive an averaged nonlinear partial differential equation describing the dynamics of the complex field in a nonlinear Schroedinger model in the presence of a periodic nonlinearity and a periodically-varying dissipation coefficient. The incorporation of dissipation is motivated by experimental considerations. We test the numerical behavior of the derived averaged equation by comparing it to the original nonautonomous model in a prototypical case scenario and observe good agreement between the two.

We employ an evolutionary algorithm to investigate the optimal design of composite protectors using one-dimensional granular chains composed of beads of various sizes, masses, and stiffnesses. We define a fitness function using the maximum force transmitted from the protector to a "wall" that represents the body to be protected and accordingly optimize the {topology} (arrangement), {size}, and {material} of the chain. We obtain optimally randomized granular protectors characterized by high-energy equipartition and the transformation of incident waves into interacting solitary pulses. We consistently observe that the pulses traveling to the wall combine to form an extended (long-wavelength), small-amplitude pulse., Comment: 13 pages, 23 figures (many with multiple parts), to appear in Mechanics of Advanced Materials and Structures

We use experiments, numerical simulations, and theoretical analysis to investigate the propagation of highly nonlinear solitary waves in periodic arrangements of dimer (two-mass) and trimer (three-mass) cell structures in one-dimensional granular lattices. To vary the composition of the fundamental periodic units in the granular chains, we utilize beads of different materials (stainless steel, bronze, glass, nylon, polytetrafluoroethylene, and rubber). This selection allows us to tailor the response of the system based on the masses, Poisson ratios, and elastic moduli of the components. For example, we examine dimer configurations with two types of heavy particles, two types of light particles, and alternating light and heavy particles. Employing a model with Hertzian interactions between adjacent beads, we find very good agreement between experiments and numerical simulations. We find equally good agreement between these results and a theoretical analysis of the model in the long-wavelength regime that we derive for heterogeneous environments (dimer chains) and general bead interactions. Our analysis encompasses previously-studied examples as special cases and also provides key insights on the influence of heterogeneous lattices on the properties (width and propagation speed) of the nonlinear wave solutions of this system., Comment: 17 pages, 2 tables, 12 figures (several with multiple panels)

Motivated by recent proposals of ``collisionally inhomogeneous'' Bose-Einstein condensates (BECs), which have a spatially modulated scattering length, we study the existence and stability properties of bright and dark matter-wave solitons of a BEC characterized by a periodic, piecewise-constant scattering length. We use a ``stitching'' approach to analytically approximate the pertinent solutions of the underlying nonlinear Schr\"odinger equation by matching the wavefunction and its derivatives at the interfaces of the nonlinearity coefficient. To accurately quantify the stability of bright and dark solitons, we adapt general tools from the theory of perturbed Hamiltonian systems. We show that solitons can only exist at the centers of the constant regions of the piecewise-constant nonlinearity. We find both stable and unstable configurations for bright solitons and show that all dark solitons are unstable, with different instability mechanisms that depend on the soliton location. We corroborate our analytical results with numerical computations., Comment: 16 pages, 7 figures (some with multiple parts), to appear in Physical Review A

We report the propagation of highly nonlinear solitary waves in heterogeneous, periodic granular media using experiments, numerical simulations, and theoretical analysis. We examine periodic arrangements of particles in experiments in which stiffer/heavier beads (stainless steel) are alternated with softer/lighter ones (PTFE beads). We find excellent agreement between experiments and numerics in a model with Hertzian interactions between adjacent beads, which in turn agrees very well with a theoretical analysis of the model in the long-wavelength regime that we derive for heterogeneous environments and general bead interactions. Our analysis encompasses previously-studied examples as special cases and also provides key insights on the influence of the dimer lattice on the properties (width and propagation speed) of the obtained highly nonlinear wave solutions., Comment: 4.2 pages, 3 figures (each with multiple panels); adjustments in title, text, and figures; typos corrected; to appear in Physical Review E (rapid communication)

The vortex density of a rotating superfluid, divided by its particle mass, dictates the superfluid's angular velocity through the Feynman relation. To find how the Feynman relation applies to superfluid mixtures, we investigate a rotating two-component Bose-Einstein condensate, composed of bosons with different masses. We find that in the case of sufficiently strong interspecies attraction, the vortex lattices of the two condensates lock and rotate at the drive frequency, while the superfluids themselves rotate at two different velocities, whose ratio is the ratio between the particle mass of the two species. In this paper, we characterize the vortex-locked state, establish its regime of stability, and find that it surives within a disk smaller than a critical radius, beyond which vortices become unbound, and the two Bose-gas rings rotate together at the frequency of the external drive., Comment: 6 pages, 2 figures

We investigate analytically, numerically, and experimentally the modulational instability in a layered, cubically-nonlinear (Kerr) optical medium that consists of alternating layers of glass and air. We model this setting using a nonlinear Schr\"odinger (NLS) equation with a piecewise constant nonlinearity coefficient and conduct a theoretical analysis of its linear stability, obtaining a Kronig-Penney equation whose forbidden bands correspond to the modulationally unstable regimes. We find very good {\it quantitative} agreement between the theoretical analysis of the Kronig-Penney equation, numerical simulations of the NLS equation, and the experimental results for the modulational instability. Because of the periodicity in the evolution variable arising from the layered medium, we find multiple instability regions rather than just the one that would occur in uniform media., Comment: 13 pages, 12 figures (several with multiple parts); some important changes from the page proof stage implemented in this preprint version

We present the first experimental investigation of modulational instability in a layered Kerr medium. The particularly interesting and appealing feature of our configuration, consisting of alternating glass-air layers, is the piecewise-constant nature of the material properties, which allows a theoretical linear stability analysis leading to a Kronig-Penney equation whose forbidden bands correspond to the modulationally unstable regimes. We find very good {\it quantitative} agreement between theoretical, numerical, and experimental diagnostics of the modulational instability. Because of the periodicity in the evolution variable arising from the layered medium, there are multiple instability regions rather than just one as in the uniform medium., Comment: 4 pages, 4 figures, contains experimental + computational + theoretical results, to appear in Physical Review Letters

We investigate the dynamics of an effectively one-dimensional Bose-Einstein condensate (BEC) with scattering length $a$ subjected to a spatially periodic modulation, $a=a(x)=a(x+L)$. This "collisionally inhomogeneous" BEC is described by a Gross-Pitaevskii (GP) equation whose nonlinearity coefficient is a periodic function of $x$. We transform this equation into a GP equation with constant coefficient $a$ and an additional effective potential and study a class of extended wave solutions of the transformed equation. For weak underlying inhomogeneity, the effective potential takes a form resembling a superlattice, and the amplitude dynamics of the solutions of the constant-coefficient GP equation obey a nonlinear generalization of the Ince equation. In the small-amplitude limit, we use averaging to construct analytical solutions for modulated amplitude waves (MAWs), whose stability we subsequently examine using both numerical simulations of the original GP equation and fixed-point computations with the MAWs as numerically exact solutions. We show that "on-site" solutions, whose maxima correspond to maxima of $a(x)$, are significantly more stable than their "off-site" counterparts., Comment: 25 pages, 10 figures (many with several parts), to appear in Physica D; higher resolution versions of some figures are available at http://www.its.caltech.edu/~mason/papers

We conduct an experimental investigation of nonlinearity management in optics using femtosecond pulses and layered Kerr media consisting of glass and air. By examining the propagation properties over several diffraction lengths, we show that wave collapse can be prevented. We corroborate these experimental results with numerical simulations of the (2 + 1)-dimensional focusing cubic nonlinear Schr\"odinger equation with piecewise constant coefficients and a theoretical analysis of this setting using a moment method., Comment: 4 pages, 4 figures (some with multiple parts; color online), to appear in Physical Review Letters; changes from previous version: revised figures, discussion of time-dependent simulations, other touch-ups

We investigate the generation of fractional-period states in continuum periodic systems. As an example, we consider a Bose-Einstein condensate confined in an optical-lattice potential. We show that when the potential is turned on non-adiabatically, the system explores a number of transient states whose periodicity is a fraction of that of the lattice. We illustrate the origin of fractional-period states analytically by treating them as resonant states of a parametrically forced Duffing oscillator and discuss their transient nature and potential observability., Comment: 10 pages, 6 figures (some with multiple parts); revised version: minor clarifications of a couple points, to appear in Physical Review A

We analyze the existence and stability of bright, dark, and gap matter-wave solitons in optical superlattices. Then, using these properties, we show that (time-dependent) ``dynamical superlattices'' can be used to controllably place, guide, and manipulate these solitons. In particular, we use numerical experiments to displace solitons by turning on a secondary lattice structure, transfer solitons from one location to another by shifting one superlattice substructure relative to the other, and implement solitonic ``path-following'', in which a matter wave follows the time-dependent lattice substructure into oscillatory motion., Comment: 6 pages, revtex, 6 figures, to appear in Physics Letters A; minor modifications from last version

We analyze gap solitons in trapped Bose-Einstein condensates (BECs) in optical lattice potentials under Feshbach resonance management. Starting with an averaged Gross-Pitaevsky (GP) equation with a periodic potential, we employ an envelope wave approximation to derive coupled-mode equations describing the slow BEC dynamics in the first spectral gap of the optical lattice. We construct exact analytical formulas describing gap soliton solutions and examine their spectral stability using the Chebyshev interpolation method. We show that these gap solitons are unstable far from the threshold of local bifurcation and that the instability results in the distortion of their shape. We also predict the threshold of the power of gap solitons near the local bifurcation limit., Comment: 8 pages, 4 figures (1 with six parts, 3 with two parts)

The Fermi-Pasta-Ulam (FPU) model, which was proposed 50 years ago to examine thermalization in non-metallic solids and develop ``experimental'' techniques for studying nonlinear problems, continues to yield a wealth of results in the theory and applications of nonlinear Hamiltonian systems with many degrees of freedom. Inspired by the studies of this seminal model, solitary-wave dynamics in lattice dynamical systems have proven vitally important in a diverse range of physical problems--including energy relaxation in solids, denaturation of the DNA double strand, self-trapping of light in arrays of optical waveguides, and Bose-Einstein condensates (BECs) in optical lattices. BECS, in particular, due to their widely ranging and easily manipulated dynamical apparatuses--with one to three spatial dimensions, positive-to-negative tuning of the nonlinearity, one to multiple components, and numerous experimentally accessible external trapping potentials--provide one of the most fertile grounds for the analysis of solitary waves and their interactions. In this paper, we review recent research on BECs in the presence of deep periodic potentials, which can be reduced to nonlinear chains in appropriate circumstances. These reductions, in turn, exhibit many of the remarkable nonlinear structures (including solitons, intrinsic localized modes, and vortices) that lie at the heart of the nonlinear science research seeded by the FPU paradigm., Comment: 10 pages, revtex, two-columns, 3 figs, accepted fpr publication in Chaos's focus issue on the 50th anniversary of the publication of the Fermi-Pasta-Ulam problem; minor clarifications (and a couple corrected typos) from previous version

We consider the Gross-Pitaevskii (GP) equation in the presence of periodic and quasiperiodic superlattices to study cigar-shaped Bose-Einstein condensates (BECs) in such potentials. We examine spatially extended wavefunctions in the form of modulated amplitude waves (MAWs). With a coherent structure ansatz, we derive amplitude equations describing the evolution of spatially modulated states of the BEC. We then apply second-order multiple scale perturbation theory to study harmonic resonances with respect to a single lattice wavenumber as well as ultrasubharmonic resonances that result from interactions of both wavenumbers of the superlattice. In each case, we determine the resulting system's equilibria, which represent spatially periodic solutions, and subsequently examine the stability of the corresponding solutions by direct simulations of the GP equation, identifying them as typically stable solutions of the model. We then study subharmonic resonances using Hamiltonian perturbation theory, tracing robust, spatio-temporally periodic patterns., Comment: 30 pages, 5 color figures (some with multiple parts), to appear in SIAM Journal on Applied Dynamical Systems; expanded work on spatially extended solutions; added introductory material and discussions of relevant experimental work; section on localized solutions moved to another paper

We consider a system of two Gross-Pitaevskii (GP) equations, in the presence of an optical-lattice (OL) potential, coupled by both nonlinear and linear terms. This system describes a Bose-Einstein condensate (BEC) composed of two different spin states of the same atomic species, which interact linearly through a resonant electromagnetic field. In the absence of the OL, we find plane-wave solutions and examine their stability. In the presence of the OL, we derive a system of amplitude equations for spatially modulated states which are coupled to the periodic potential through the lowest-order subharmonic resonance. We determine this averaged system's equilibria, which represent spatially periodic solutions, and subsequently examine the stability of the corresponding solutions with direct simulations of the coupled GP equations. We find that symmetric (equal-amplitude) and asymmetric (unequal-amplitude) dual-mode resonant states are, respectively, stable and unstable. The unstable states generate periodic oscillations between the two condensate components, which is possible only because of the linear coupling between them. We also find four-mode states, but they are always unstable. Finally, we briefly consider ternary (three-component) condensates., Comment: 16 pages, 4 figures (some of which have multiple parts), to appear in Physica D; streamlined paper; added some references and discussion concerning experimental realizations of this work; higher-resolution copies of a couple figures are available on the version of the document downloadable from http://www.math.gatech.edu/~mason/

We apply Lindstedt's method and multiple scale perturbation theory to analyze spatio-temporal structures in nonlinear Schr\"odinger equations and thereby study the dynamics of quasi-one-dimensional Bose-Einstein condensates with mean-field interactions. We determine the dependence of the intensity of modulated amplitude waves on their wave number. We also explore the band structure of Bose-Einstein condensates in detail using Hamiltonian perturbation theory and supporting numerical simulations., Comment: 24 pages, 20 figs (numbered to 9), 6 tables, to appear in Chaos

We analyze spatio-temporal structures in the Gross-Pitaevskii equation to study the dynamics of quasi-one-dimensional Bose-Einstein condensates (BECs) with mean-field interactions. A coherent structure ansatz yields a parametrically forced nonlinear oscillator, to which we apply Lindstedt's method and multiple-scale perturbation theory to determine the dependence of the intensity of periodic orbits (``modulated amplitude waves'') on their wave number. We explore BEC band structure in detail using Hamiltonian perturbation theory and supporting numerical simulations., Comment: 5 pages, 4 figs, revtex, final form of paper, to appear in PRE (forgot to include \bibliography command in last update, so this is a correction of that; the bibliography is hence present again)