Search

Achieving parametric excitation in an oscillating physical system involves periodically adjusting one of its parameters to modulate the oscillator's natural frequency. This phenomenon has been observed in numerous systems within physics and engineering, profoundly transforming modern science and technology. Despite rapid progress, the parametric control of collective excitations, such as phonons, remains a challenge while promising to generate novel and intriguing effects in a largely unexplored field. Here, we investigate the terahertz (THz) field-induced dynamics of Raman-active phonons in LaAlO$_3$ (LAO), an archetypal quantum material for its electronic and structural properties. Utilizing intense THz pulses, we demonstrate a novel mechanism of parametric phonon excitation marked by substantial subharmonic components. Theoretical analysis can successfully capture the hallmarks of the observed phenomena in a physical scenario with the THz field inducing a parametric coupling between the Raman mode and pairs of acoustic phonon excitations., Comment: Main text: 12 pages, 3 figures. Supplementary Information: 3 figures

Self-oscillating systems, described in classical dynamics as limit cycles, are emerging as canonical models for driven dissipative nonequilibrium open quantum systems, and as key elements in quantum technology. We consider a family of models that interpolates between the classical textbook examples of the Rayleigh and the van der Pol oscillators, and follow their transition from the classical to the quantum domain, while properly formulating their corresponding quantum descriptions. We derive an exact analytical solution for the steady-state quantum dynamics of the simplest of these models, applicable to any bosonic system---whether mechanical, optical, or otherwise---that is coupled to its environment via single-boson and double-boson emission and absorption. Our solution is a generalization to arbitrary temperature of existing solutions for very-low, or zero, temperature, often misattributed to the quantum van der Pol oscillator. We closely explore the classical to quantum transition of the bifurcation to self-oscillations of this oscillator, while noting changes in the dynamics and identifying features that are uniquely quantum.

We present a method for analyzing the phase noise of oscillators based on feedback driven high quality factor resonators. Our approach is to derive the phase drift of the oscillator by projecting the stochastic oscillator dynamics onto a slow time scale corresponding physically to the long relaxation time of the resonator. We derive general expressions for the phase drift generated by noise sources in the electronic feedback loop of the oscillator. These are mixed with the signal through the nonlinear amplifier, which makes them {cyclostationary}. We also consider noise sources acting directly on the resonator. The expressions allow us to investigate reducing the oscillator phase noise thereby improving the frequency precision using resonator nonlinearity by tuning to special operating points. We illustrate the approach giving explicit results for a phenomenological amplifier model. We also propose a scheme for measuring the slow feedback noise generated by the feedback components in an open-loop driven configuration in experiment or using circuit simulators, which enables the calculation of the closed-loop oscillator phase noise in practical systems.

We study the role of amplifier saturation in eliminating feedback noise in self-sustained oscillators. We extend previous works that use a saturated amplifier to quench fluctuations in the feedback magnitude, while simultaneously tuning the oscillator to an operational point at which the resonator nonlinearity cancels fluctuations in the feedback phase. We consider a generalized model which features an amplitude-dependent amplifier gain function. This allows us to determine the total oscillator phase noise in realistic configurations due to noise in both quadratures of the feedback, and to show that it is not necessary to drive the resonator to large oscillation amplitudes in order to eliminate noise in the phase of the feedback.

We investigate the synchronization of oscillators based on anharmonic nanoelectromechanical resonators. Our experimental implementation allows unprecedented observation and control of parameters governing the dynamics of synchronization. We find close quantitative agreement between experimental data and theory describing reactively coupled Duffing resonators with fully saturated feedback gain. In the synchronized state we demonstrate a significant reduction in the phase noise of the oscillators, which is key for sensor and clock applications. Our work establishes that oscillator networks constructed from nanomechanical resonators form an ideal laboratory to study synchronization given their high-quality factors, small footprint, and ease of co-integration with modern electronic signal processing technologies., Comment: 34 pages, 8 figures

Two-dimensional lattices of N synchronized oscillators with reactive coupling are considered as high-precision frequency sources in the case where a spiral pattern is formed. The improvement of the frequency precision is shown to be independent of N for large N, unlike the case of purely dissipative coupling where the improvement is proportional to N, but instead depends on just those oscillators in the core of the spiral that acts as the source region of the waves. Our conclusions are based on numerical simulations of up to N=29929 oscillators, and analytic results for a continuum approximation to the lattice in an infinite system. We derive an expression for the dependence of the frequency precision on the reactive component of the coupling constant, depending on a single parameter given by fitting the frequency of the spiral waves to the numerical simulations., Comment: 7 pages, 6 figures

We demonstrate an analytical method for calculating the phase sensitivity of a class of oscillators whose phase does not affect the time evolution of the other dynamic variables. We show that such oscillators possess the possibility for complete phase noise elimination. We apply the method to a feedback oscillator which employs a high Q weakly nonlinear resonator and provide explicit parameter values for which the feedback phase noise is completely eliminated and others for which there is no amplitude-phase noise conversion. We then establish an operational mode of the oscillator which optimizes its performance by diminishing the feedback noise in both quadratures, thermal noise, and quality factor fluctuations. We also study the spectrum of the oscillator and provide specific results for the case of 1/f noise sources.

Self-sustained oscillators are ubiquitous and essential for metrology, communications, time reference, and geolocation. In its most basic form an oscillator consists of a resonator driven on-resonance, through feedback, to create a periodic signal sustained by a static energy source. The generation of a stable frequency, the basic function of oscillators, is typically achieved by increasing the amplitude of motion of the resonator while remaining within its linear, harmonic, regime. Contrary to this conventional paradigm, in this Letter we show that by operating the oscillator at special points in the resonators anharmonic regime we can overcome fundamental limitations of oscillator performance due to thermodynamic noise as well as practical limitations due to noise from the sustaining circuit. We develop a comprehensive model that accounts for the major contributions to the phase noise of the nonlinear oscillator. Using a nanoelectromechanical system (NEMS)-based oscillator, we experimentally verify the existence of a special region in the operational parameter space that enables a significant reduction of the oscillators phase noise, as predicted by our model., Comment: 14 pages, 2 figures, 1 table

We study the difference between quantum and classical behavior in a pair of nonidentical cavities with second-harmonic generation. In the classical limit, each cavity has a limit-cycle solution, in which the photon number oscillates periodically in time. Coupling between the cavities leads to synchronization of the oscillations and classical correlations between the cavities. In the quantum limit, there are quantum correlations due to entanglement. The quantum correlations persist even when the cavities are far off resonance with each other, in stark contrast with the classical case. We also find that the quantum and classical limits are connected by an intermediate regime of almost no correlations. Our results can be extended to a wide variety of quantum models., Comment: 6 pages, 7 figures

We introduce a new method for reducing phase noise in oscillators, thereby improving their frequency precision. The noise reduction device consists of a pair of coupled nonlinear resonating elements that are driven parametrically by the output of a conventional oscillator at a frequency close to the sum of the linear mode frequencies. Above the threshold for parametric response, the coupled resonators exhibit self-oscillation at an inherent frequency. We find operating points of the device for which this periodic signal is immune to frequency noise in the driving oscillator, providing a way to clean its phase noise. We present results for the effect of thermal noise to advance a broader understanding of the overall noise sensitivity and the fundamental operating limits.

We study the nonequilibrium dynamics of quantum jumps in a one-dimensional chain of atoms. Each atom is driven on a strong transition to a short-lived state and on a weak transition to a metastable state. We choose the metastable state to be a Rydberg state so that when an atom jumps to the Rydberg state, it inhibits or enhances jumps in the neighboring atoms. This leads to rich spatiotemporal dynamics that are visible in the fluorescence of the strong transition., Comment: 10 pages

The study of collective nonlinear dynamics of coupled mechanical resonators is regaining attention in recent years thanks to rapid developments in the fields of microelectromechanical and nanoelectromechanical systems (MEMS and NEMS). We review a wide range of collective dynamical phenomena, while highlighting the common concepts and theoretical tools that we have developed for treating them. We provide detailed derivations of amplitude equations, which allow us to obtain reduced descriptions for the relevant dynamics of our complex systems. We apply these amplitude equations to study (a) resonant response to parametric excitation; (b) pattern selection, or the nonlinear competition between extended modes in situations of multistability; (c) formation and dynamics of intrinsically localized modes (ILM); and (d) spontaneous synchronization of oscillators with differing frequencies. All the predictions obtained from analyzing the different amplitude equations are in excellent agreement with numerical solutions of the underlying equations of motion, suggesting that the predicted effects can be observed in arrays of real micromechanical or nanomechanical resonators, thus motivating new experiments in these systems., Comment: To appear in "Fluctuating nonlinear oscillators", Ed. Mark I. Dykman (Oxford University Press, Oxford, 2012)

We study an open quantum system of atoms with long-range Rydberg interaction, laser driving, and spontaneous emission. Over time, the system occasionally jumps between a state of low Rydberg population and a state of high Rydberg population. The jumps are inherently collective and in fact exist only for a large number of atoms. We explain how entanglement and quantum measurement enable the jumps, which are otherwise classically forbidden., Comment: 4 pages

Improving the frequency precision by synchronizing a lattice of oscillators is studied in the phase reduction limit. For the most commonly studied case of purely dissipative phase coupling (the Kuramoto model) I confirm that the frequency precision of N oscillators perturbed by independent noise sources is improved by a factor N as expected from simple averaging arguments. In the presence of reactive coupling, such as will typically be the case for non-dissipatively coupled oscillators based on high-Q resonators, the synchronized state consists of target like waves radiating from a local source which is a region of higher frequency oscillators. In this state all the oscillators evolve with the same frequency, however I show that the improvement of the frequency precision is independent of N for large N, but instead depends on the disorder and reflects the dependence of the frequency of the synchronized state on just those oscillators in the source region of the waves., Comment: 18 pages, 8 figures; clarified wording of text and added references

We study a driven-dissipative system of atoms in the presence of laser excitation to a Rydberg state and spontaneous emission. The atoms interact via the blockade effect, whereby an atom in the Rydberg state shifts the Rydberg level of neighboring atoms. We use mean-field theory to study how the Rydberg population varies in space. As the laser frequency changes, there is a continuous transition between the uniform and antiferromagnetic phases. The nonequilibrium nature also leads to a novel oscillatory phase and bistability between the uniform and antiferromagnetic phases., Comment: 4 pages + appendix

Ion traps are a versatile tool to study nonequilibrium statistical physics, due to the tunability of dissipation and nonlinearity. We propose an experiment with a chain of trapped ions, where dissipation is provided by laser heating and cooling, while nonlinearity is provided by trap anharmonicity and beam shaping. The collective dynamics are governed by an equation similar to the complex Ginzburg-Landau equation, except that the reactive nature of the coupling leads to qualitatively different behavior. The system has the unusual feature of being both oscillatory and excitable at the same time. We account for noise from spontaneous emission and find that the patterns are observable for realistic experimental parameters. Our scheme also allows controllable experiments with noise and quenched disorder., Comment: 4 pages + appendix

We study synchronization in the two-dimensional lattice of coupled phase oscillators with random intrinsic frequencies. When the coupling $K$ is larger than a threshold $K_E$, there is a macroscopic cluster of frequency-synchronized oscillators. We explain why the macroscopic cluster disappears at $K_E$. We view the system in terms of vortices, since cluster boundaries are delineated by the motion of these topological defects. In the entrained phase ($K>K_E$), vortices move in fixed paths around clusters, while in the unentrained phase ($K

We present analytical calculations and numerical simulations for the synchronization of oscillators interacting via a long range power law interaction on a one dimensional lattice. We have identified the critical value of the power law exponent $\alpha_c$ across which a transition from a synchronized to an unsynchronized state takes place for a sufficiently strong but finite coupling strength in the large system limit. We find $\alpha_c=3/2$. Frequency entrainment and phase ordering are discussed as a function of $\alpha \geq 1$. The calculations are performed using an expansion about the aligned phase state (spin-wave approximation) and a coarse graining approach. We also generalize the spin-wave results to the {\it d}-dimensional problem., Comment: Final published version

We apply a recently developed renormalization group (RG) method to study synchronization in a one-dimensional chain of phase-coupled oscillators in the regime of weak randomness. The RG predicts how oscillators with randomly distributed frequencies and couplings form frequency-synchronized clusters. Although the RG was originally intended for strong randomness, i.e. for distributions with long tails, we find good agreement with numerical simulations even in the regime of weak randomness. We use the RG flow to derive how the correlation length scales with the width of the coupling distribution in the limit of large coupling. This leads to the identification of a universality class of distributions with the same critical exponent $\nu$. We also find universal scaling for small coupling. Finally, we show that the RG flow is characterized by a universal approach to the unsynchronized fixed point, which provides physical insight into low-frequency clusters., Comment: 14 pages, 10 figures

We study intrinsic localized modes (ILMs), or solitons, in arrays of parametrically-driven nonlinear resonators with application to microelectromechanical and nanoelectromechanical systems (MEMS and NEMS). The analysis is performed using an amplitude equation in the form of a nonlinear Schroedinger equation with a term corresponding to nonlinear damping (also known as a forced complex Ginzburg-Landau equation), which is derived directly from the underlying equations of motion of the coupled resonators, using the method of multiple scales. We investigate the creation, stability, and interaction of ILMs, show that they can form bound states, and that under certain conditions one ILM can split into two. Our findings are confirmed by simulations of the underlying equations of motion of the resonators, suggesting possible experimental tests of the theory., Comment: New version contains a new section discussing the dynamical formation of solitons and the stability of non-zero uniform states

Two elastically coupled nanomechanical resonators driven independently near their resonance frequencies show intricate nonlinear dynamics. The dynamics provide a scheme for realizing a nanomechanical system with tunable frequency and nonlinear properties. For large vibration amplitudes the system develops spontaneous oscillations of amplitude modulation that also show period doubling transitions and chaos. The complex nonlinear dynamics are quantitatively predicted by a simple theoretical model., Comment: Acknowledgment and references added. Slightly shortened to fit journal requirements

We develop a renormalization group method to investigate synchronization clusters in a one-dimensional chain of nearest-neighbor coupled phase oscillators. The method is best suited for chains with strong disorder in the intrinsic frequencies and coupling strengths. The results are compared with numerical simulations of the chain dynamics and good agreement in several characteristics is found. We apply the renormalization group and simulations to Lorentzian distributions of intrinsic frequencies and couplings and investigate the statistics of the resultant cluster sizes and frequencies, as well as the dependence of the characteristic cluster length upon parameters of these Lorentzian distributions., Comment: 13 pages, 7 figures: 8 figure files in *.eps format and one *.tex file

We study the problem of pattern selection in an array of parametrically-driven nonlinear resonators with application to microelectromechanical and nanoelectromechanical systems (MEMS & NEMS), using an amplitude equation recently derived by Bromberg, Cross, and Lifshitz [PRE 73, 016214 (2006)]. We describe the transitions between standing-wave patterns of different wave numbers as the drive amplitude is varied either quasistatically, abruptly, or as a linear ramp in time. We find novel hysteretic effects, which are confirmed by numerical integration of the original equations of motion of the interacting nonlinear resonators.

The stochastic dynamics of micron and nanoscale cantilevers immersed in a viscous fluid are quantified. Analytical results are presented for long slender cantilevers driven by Brownian noise. The spectral density of the noise force is not assumed to be white and the frequency dependence is determined from the fluctuation-dissipation theorem. The analytical results are shown to be useful for the micron scale cantilevers that are commonly used in atomic force microscopy. A general thermodynamic approach is developed that is valid for cantilevers of arbitrary geometry as well as for arrays of multiple cantilevers whose stochastic motion is coupled through the fluid. It is shown that the fluctuation-dissipation theorem permits the calculation of stochastic quantities via straightforward deterministic methods. The thermodynamic approach is used with deterministic finite element numerical simulations to quantify the autocorrelation and noise spectrum of cantilever fluctuations for a single micron scale cantilever and the cross-correlations and noise spectra of fluctuations for an array of two experimentally motivated nanoscale cantilevers as a function of cantilever separation. The results are used to quantify the noise reduction possible using correlated measurements with two closely spaced nanoscale cantilevers., Comment: Submitted to Nanotechnology April 26, 2006

Experiments and simulations from a variety of sample sizes indicated that the centrifugal force significantly affects rotating Rayleigh-B\'enard convection-patterns. In a large-aspect-ratio sample, we observed a hybrid state consisting of domain chaos close to the sample center, surrounded by an annulus of nearly-stationary nearly-radial rolls populated by occasional defects reminiscent of undulation chaos. Although the Coriolis force is responsible for domain chaos, by comparing experiment and simulation we show that the centrifugal force is responsible for the radial rolls. Furthermore, simulations of the Boussinesq equations for smaller aspect ratios neglecting the centrifugal force yielded a domain precession-frequency $f\sim\epsilon^\mu$ with $\mu\simeq1$ as predicted by the amplitude-equation model for domain chaos, but contradicted by previous experiment. Additionally the simulations gave a domain size that was larger than in the experiment. When the centrifugal force was included in the simulation, $\mu$ and the domain size closely agreed with experiment., Comment: 8 pages, 11 figures

We present a detailed analysis of a model for the synchronization of nonlinear oscillators due to reactive coupling and nonlinear frequency pulling. We study the model for the mean field case of all-to-all coupling, deriving results for the initial onset of synchronization as the coupling or nonlinearity increase, and conditions for the existence of the completely synchronized state when all the oscillators evolve with the same frequency. Explicit results are derived for Lorentzian, triangular, and top-hat distributions of oscillator frequencies. Numerical simulations are used to construct complete phase diagrams for these distributions.

Magnetostatic modes of Yttrium Iron Garnet (YIG) films are investigated by ferromagnetic resonance force microscopy (FMRFM). A thin film ``probe'' magnet at the tip of a compliant cantilever introduces a local inhomogeneity in the internal field of the YIG sample. This influences the shape of the sample's magnetostatic modes, thereby measurably perturbing the strength of the force coupled to the cantilever. We present a theoretical model that explains these observations; it shows that tip-induced variation of the internal field creates either a local ``potential barrier'' or ``potential well'' for the magnetostatic waves. The data and model together indicate that local magnetic imaging of ferromagnets is possible, even in the presence of long-range spin coupling, through the induction of localized magnetostatic modes predicted to arise from sufficiently strong tip fields., Comment: 5 pages, 4 figures, 1 table, 15 references

We study the response of a large array of coupled nonlinear oscillators to parametric excitation, motivated by the growing interest in the nonlinear dynamics of microelectromechanical and nanoelectromechanical systems (MEMS and NEMS). Using a multiscale analysis, we derive an amplitude equation that captures the slow dynamics of the coupled oscillators just above the onset of parametric oscillations. The amplitude equation that we derive here from first principles exhibits a wavenumber dependent bifurcation similar in character to the behavior known to exist in fluids undergoing the Faraday wave instability. We confirm this behavior numerically and make suggestions for testing it experimentally with MEMS and NEMS resonators., Comment: Version 2 is an expanded version of the article, containing detailed steps of the derivation that were left out in version 1, but no additional results

To understand how spatiotemporal chaos may modify material transport, we use direct numerical simulations of the three-dimensional Boussinesq equations and of an advection-diffusion equation to study the transport of a passive tracer by the spiral defect chaos state of a convecting fluid. The simulations show that the transport is diffusive and is enhanced by the spatiotemporal chaos. The enhancement in tracer diffusivity follows two regimes. For large Peclet numbers (that is, small molecular diffusivities of the tracer), we find that the enhancement is proportional to the Peclet number. For small Peclet numbers, the enhancement is proportional to the square root of the Peclet number. We explain the presence of these two regimes in terms of how the local transport depends on the local wave numbers of the convection rolls. For large Peclet numbers, we further find that defects cause the tracer diffusivity to be enhanced locally in the direction orthogonal to the local wave vector but suppressed in the direction of the local wave vector., Comment: 11 pages, 12 figures

We analyze a model for the synchronization of nonlinear oscillators due to reactive coupling and nonlinear frequency pulling motivated by the physics of arrays of nanoscale oscillators. We study the model for the mean field case of all-to-all coupling, deriving results for the onset of synchronization as the coupling or nonlinearity increase, and the fully locked state when all the oscillators evolve with the same frequency.

The coarsening and wavenumber selection of striped states growing from random initial conditions are studied in a non-relaxational, spatially extended, and far-from-equilibrium system by performing large-scale numerical simulations of Rayleigh-B\'{e}nard convection in a large-aspect-ratio cylindrical domain with experimentally realistic boundaries. We find evidence that various measures of the coarsening dynamics scale in time with different power-law exponents, indicating that multiple length scales are required in describing the time dependent pattern evolution. The translational correlation length scales with time as $t^{0.12}$, the orientational correlation length scales as $t^{0.54}$, and the density of defects scale as $t^{-0.45}$. The final pattern evolves toward the wavenumber where isolated dislocations become motionless, suggesting a possible wavenumber selection mechanism for large-aspect-ratio convection., Comment: 5 pages, 6 figures

The stochastic response of nanoscale oscillators of arbitrary geometry immersed in a viscous fluid is studied. Using the fluctuation-dissipation theorem it is shown that deterministic calculations of the governing fluid and solid equations can be used in a straightforward manner to directly calculate the stochastic response that would be measured in experiment. We use this approach to investigate the fluid coupled motion of single and multiple cantilevers with experimentally motivated geometries., Comment: 5 pages, 5 figures

We investigate a scheme that makes a quantum non-demolition measurement of the excitation level of a mesoscopic mechanical oscillator by utilizing the anharmonic coupling between two elastic beam bending modes. The non-linear coupling between the two modes shifts the resonant frequency of the readout oscillator proportionate to the excitation of the system oscillator. This frequency shift may be detected as a phase shift of the readout oscillation when driven on resonance. We show that in an appropriate regime this measurement approaches a quantum non-demolition measurement of the phonon number of the system oscillator. As phonon energies in micromechanical oscillators become comparable to or greater than the thermal energy, the individual phonon dynamics within each mode can be resolved. As a result it should be possible to monitor jumps between Fock states caused by the coupling of the system to the thermal reservoirs., Comment: revised, 21 pages, 9 figures

We describe a numerical procedure to construct a modified velocity field that does not have any mean flow. Using this procedure, we present two results. Firstly, we show that, in the absence of mean flow, spiral defect chaos collapses to a stationary pattern comprising textures of stripes with angular bends. The quenched patterns are characterized by mean wavenumbers that approach those uniquely selected by focus-type singularities, which, in the absence of mean flow, lie at the zig-zag instability boundary. The quenched patterns also have larger correlation lengths and are comprised of rolls with less curvature. Secondly, we describe how mean flow can contribute to the commonly observed phenomenon of rolls terminating perpendicularly into lateral walls. We show that, in the absence of mean flow, rolls begin to terminate into lateral walls at an oblique angle. This obliqueness increases with Rayleigh number., Comment: 14 pages, 19 figures

Rayleigh-B\'{e}nard convection is studied and quantitative comparisons are made, where possible, between theory and experiment by performing numerical simulations of the Boussinesq equations for a variety of experimentally realistic situations. Rectangular and cylindrical geometries of varying aspect ratios for experimental boundary conditions, including fins and spatial ramps in plate separation, are examined with particular attention paid to the role of the mean flow. A small cylindrical convection layer bounded laterally either by a rigid wall, fin, or a ramp is investigated and our results suggest that the mean flow plays an important role in the observed wavenumber. Analytical results are developed quantifying the mean flow sources, generated by amplitude gradients, and its effect on the pattern wavenumber for a large-aspect-ratio cylinder with a ramped boundary. Numerical results are found to agree well with these analytical predictions. We gain further insight into the role of mean flow in pattern dynamics by employing a novel method of quenching the mean flow numerically. Simulations of a spiral defect chaos state where the mean flow is suddenly quenched is found to remove the time dependence, increase the wavenumber and make the pattern more angular in nature., Comment: 9 pages, 10 figures