Source author record

Avadh Saxena

Avadh Saxena appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

82works
20topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

82 published item(s)

preprint2022arXiv

Kink Solutions With Power Law Tails

We present a comprehensive review about the various facets of kink solutions with a power law tail which have received considerable attention during the last few years. This area of research is in its early stages and while several aspects have become clear by now, there are a number of issues which have only been partially understood or not understood at all. We first discuss the aspects which are reasonably well known and then address in some detail the issues which are only partially or not understood at all. We present a wide class of higher (than sixth) order field theory models admitting implicit kink as well as mirror kink solutions where the two tails facing each other have a power law or a power-tower type fall off while the other two ends not facing each other could have either an exponential or a power law tail. The models admitting implicit kink solutions where the two ends facing each other have an exponential tail while the other two ends have a power law tail are also discussed. Moreover, we present several field theory models which admit explicit kink solutions with a power law fall off. We note that in all the polynomial models while the potential $V(ϕ)$ is continuous, its derivative is discontinuous. We also discuss one of the most important and only partially understood issues of the kink-kink and the kink-antikink forces in case the tails facing each other have a power law fall off. Finally, we briefly discuss the kink-antikink collisions at finite velocity and present some open questions.

preprint2022arXiv

Novel Superposed Kink and Pulse Solutions for $ϕ^4$, MKdV, NLS and Other Nonlinear Equations

We show that a number of nonlinear equations including symmetric as well as asymmetric $ϕ^4$, modified Korteweg de Vries (MKdV), mixed KdV-MKdV, nonlinear Schrödinger (NLS), quadratic-cubic NLS as well as higher order neutral scalar field theories, higher order KdV-MKdV and higher order quadratic-cubic NLS admit superposed periodic kink and pulse solutions and some of them also admit superposed hyperbolic kink solutions.

preprint2022arXiv

Spectral Filtering Induced by Non-Hermitian Evolution with Balanced Gain and Loss: Enhancing Quantum Chaos

The dynamical signatures of quantum chaos in an isolated system are captured by the spectral form factor, which exhibits as a function of time a dip, a ramp, and a plateau, with the ramp being governed by the correlations in the level spacing distribution. While decoherence generally suppresses these dynamical signatures, the nonlinear non-Hermitian evolution with balanced gain and loss (BGL) in an energy-dephasing scenario can enhance manifestations of quantum chaos. In the Sachdev-Ye-Kitaev model and random matrix Hamiltonians, BGL increases the span of the ramp, lowering the dip as well as the value of the plateau, providing an experimentally realizable physical mechanism for spectral filtering. The chaos enhancement due to BGL is optimal over a family of filter functions that can be engineered with fluctuating Hamiltonians.

preprint2022arXiv

Stability of exact solutions of the $(2+1)$-dimensional nonlinear Schrödinger equation with arbitrary nonlinearity parameter $κ$

In this work, we consider the nonlinear Schrödinger equation (NLSE) in $2+1$ dimensions with arbitrary nonlinearity exponent $κ$ in the presence of an external confining potential. Exact solutions to the system are constructed, and their stability over their "mass" (i.e., the $L^2$ norm) and the parameter $κ$ is explored. We observe both theoretically and numerically that the presence of the confining potential leads to wider domains of stability over the parameter space compared to the unconfined case. Our analysis suggests the existence of a stable regime of solutions for all $κ$ as long as their mass is less than a critical value $M^{\ast}(κ)$. Furthermore, we find that there are two different critical masses, one corresponding to width perturbations and the other one to translational perturbations. The results of Derrick's theorem are also obtained by studying the small amplitude regime of a four-parameter collective coordinate (4CC) approximation. A numerical stability analysis of the NLSE shows that the instability curve $M^{\ast}(κ)$ vs. $κ$ lies below the two curves found by Derrick's theorem and the 4CC approximation. In the absence of the external potential, $κ=1$ demarcates the separation between the blowup regime and the stable regime. In this 4CC approximation, for $κ<1$, when the mass is above the critical mass for the translational instability, quite complicated motions of the collective coordinates are possible. Energy conservation prevents the blowup of the solution as well as confines the center of the solution to a finite spatial domain. We call this regime the "frustrated" blowup regime and give some illustrations. In an appendix, we show how to extend these results to arbitrary initial ground state solution data and arbitrary spatial dimension $d$.

preprint2022arXiv

Superposed Hyperbolic Kink and Pulse Solutions of Coupled $ϕ^4$, NLS and MKdV Equations

We obtain novel solutions of a coupled $ϕ^4$, a coupled nonlinear Schrödinger (NLS) and a coupled modified Korteweg de Vries (MKdV) model which can be re-expressed as a linear superposition of either the sum or the difference of two hyperbolic kink or two hyperbolic pulse solutions. These results demonstrate that the notion of superposed solutions extends to coupled nonlinear equations as well.

preprint2021arXiv

Explicit Kink Solutions in Several One-Parameter Family of Higher Order Field Theory Models

We present several one-parameter family of higher order field theory models some of which admit explicit kink solutions with an exponential tail while others admit explicit kink solutions with a power-law tail. Various properties of these families of kink solutions are examined in detail. In particular for models harboring kink solutions with a power-law tail, we show that there is no gap between the zero mode and the beginning of the continuum in the kink stability potential. Further, by applying the recent Manton formalism, we provide estimates for the kink-kink and antikink-kink acceleration and hence the ratio of the corresponding antikink-kink and kink-kink forces.

preprint2021arXiv

Novel Deformation Function Creating or Destroying any Number of Even Kink Solutions

We present a one-parameter family of deformation functions $f(ϕ)$ which have novel properties. Firstly, the deformation function is its own inverse. We show that a class of potentials remains invariant under this deformation. Further, when applied to a certain class of kink bearing potentials, one obtains potentials which are not bounded from below. Besides, we show that there is a wide class of potentials having power law kink tails which are connected by this deformation but the corresponding kink tails of the two potentials have different asymptotic behavior. Finally, we show that when this deformation function is applied to an appropriate one-parameter family of potentials having two kink solutions, it creates new potentials with an arbitrary even number ($2n$) of kink solutions. Conversely, by starting from this potential with $2n$ kink solutions the deformation function annihilates $2n-2$ kinks and we get back the potential which we started with having only two kink solutions.

preprint2020arXiv

Alignment induced re-configurable walls for patterning and assembly of liquid crystal skyrmions

Skyrmions have attracted rapidly growing interest due to their topological properties and unique aspects for potential novel applications such as data storage and soft robotics. They can also serve as key elements for materials by design, self-assembly, and functional soft materials. While not real particles, these skyrmions behave like particles-they interact with each other and can be actuated by means of electric field, surface anchoring, and light. On the other hand, they are field configurations which have properties not possessed by real particles. Here, we show that, by means of alignment induced attractive and repulsive walls, skyrmions in chiral nematic liquid crystals can be precisely controlled and programmed to serve as suitable building blocks for the realization of the above goals. Our work may stimulate new experimental efforts and concomitant applications in this direction.

preprint2020arXiv

Anti-$\mathcal{PT}$-symmetric Qubit: Decoherence and Entanglement Entropy

We investigate a two-level spin system based anti-parity-time (anti-$\mathcal{PT}$)-symmetric qubit and study its decoherence as well as entanglement entropy properties. We compare our findings with that of the corresponding $\mathcal{PT}$-symmetric and Hermitian qubits. First we consider the time-dependent Dyson map to find the exact analytical dynamics for a general non-Hermitian qubit system coupled with a bath, then we specialize it to the case of the anti-$\mathcal{PT}$-symmetric qubit. We find that the decoherence function for the anti-$\mathcal{PT}$-symmetric qubit decays slower than the $\mathcal{PT}$-symmetric and Hermitian qubits. For the entanglement entropy we find that for the anti-$\mathcal{PT}$-symmetric qubit it grows more slowly compared to the $\mathcal{PT}$-symmetric and Hermitian qubits. Similarly, we find that the corresponding variance and area of Fisher information is much higher compared to the $\mathcal{PT}$-symmetric and Hermitian qubits. These results demonstrate that anti-$\mathcal{PT}$-symmetric qubits may be better suited for quantum computing and quantum information processing applications than conventional Hermitian or even $\mathcal{PT}$-symmetric qubits.

preprint2020arXiv

Domain wall diode based on functionally graded Dzyaloshinskii-Moriya interaction

We present a general approach for studying the dynamics of domain walls in biaxial ferromagnetic stripes with functionally graded Dzyaloshinskii-Moriya interaction (DMI). By engineering the spatial profile of the DMI parameter we propose the concept of a diode, which implements filtering of domain walls of certain topological charge and helicity. We base our study on phenomenological Landau-Lifshitz-Gilbert equations with additional Zhang-Li spin-transfer terms using a collective variable approach. In the effective equations of motion the gradients of DMI play the role of a driving force which competes with current driving. All analytical predictions are confirmed by numerical simulations.

preprint2020arXiv

Parity-time symmetric systems with memory

Classical open systems with balanced gain and loss, i.e. parity-time ($\mathcal{PT}$) symmetric systems, have attracted tremendous attention over the past decade. Their exotic properties arise from exceptional point (EP) degeneracies of non-Hermitian Hamiltonians that govern their dynamics. In recent years, increasingly sophisticated models of $\mathcal{PT}$-symmetric systems with time-periodic (Floquet) driving, time-periodic gain and loss, and time-delayed coupling have been investigated, and such systems have been realized across numerous platforms comprising optics, acoustics, mechanical oscillators, optomechanics, and electrical circuits. Here, we introduce a $\mathcal{PT}$-symmetric (balanced gain and loss) system with memory, and investigate its dynamics analytically and numerically. Our model consists of two coupled $LC$ oscillators with positive and negative resistance, respectively. We introduce memory by replacing either the resistor with a memristor, or the coupling inductor with a meminductor, and investigate the circuit energy dynamics as characterized by $\mathcal{PT}$-symmetric or $\mathcal{PT}$-symmetry broken phases. Due to the resulting nonlinearity, we find that energy dynamics depend on the sign and strength of initial voltages and currents, as well as the distribution of initial circuit energy across its different components. Surprisingly, at strong inputs, the system exhibits self-organized Floquet dynamics, including $\mathcal{PT}$-symmetry broken phase at vanishingly small dissipation strength. Our results indicate that $\mathcal{PT}$-symmetric systems with memory show a rich landscape.

preprint2019arXiv

Logarithmic Potential with Super-Super-Exponential Kink Profiles and Tails

We consider a novel one dimensional model of a logarithmic potential which has super-super-exponential kink profiles as well as kink tails. We provide analytic kink solutions of the model -- it has 3 kinks, 3 mirror kinks and the corresponding antikinks. While some of the kink tails are super-super-exponential, some others are super-exponential whereas the remaining ones are exponential. The linear stability analysis reveals that there is a gap between the zero mode and the onset of continuum. Finally, we compare this potential and its kink solutions with those of very high order field theories harboring seven degenerate minima and their attendant kink solutions, specifically $ϕ^{14}$, $ϕ^{16}$ and $ϕ^{18}$.

preprint2019arXiv

Noise amplification at spin-glass bottlenecks of quantum annealing: a solvable model

To gain better insight into the complexity theory of quantum annealing, we propose and solve a class of spin systems which contain bottlenecks of the kind expected to dominate the runtime of quantum annealing as it tries to solve difficult optimization problems. We uncover a noise amplification effect at these bottlenecks, whereby tunneling rates caused by flux-qubit noise scale in proportion to the number of qubits $N$ in the limit that $N\to \infty$. By solving the incoherent annealing dynamics exactly, we find a wide range of regimes where the probability that a quantum annealer remains in the ground-state upon exiting the bottleneck is close to one-half. We corroborate our analysis with detailed simulations of the performance of the D-Wave 2X quantum annealer on our class of computational problems.

preprint2019arXiv

Parametrically driven nonlinear Dirac equation with arbitrary nonlinearity

The damped and parametrically driven nonlinear Dirac equation with arbitrary nonlinearity parameter $κ$ is analyzed, when the external force is periodic in space and given by $f(x) =r\cos(K x)$, both numerically and in a variational approximation using five collective coordinates (time dependent shape parameters of the wave function). Our variational approximation satisfies exactly the low-order moment equations. Because of competition between the spatial period of the external force $λ=2 π/K$, and the soliton width $l_s$, which is a function of the nonlinearity $κ$ as well as the initial frequency $ω_0$ of the solitary wave, there is a transition (at fixed $ω_0$) from trapped to unbound behavior of the soliton, which depends on the parameters $r$ and $K$ of the external force and the nonlinearity parameter $κ$. We previously studied this phenomena when $κ=1$ (2019 J. Phys. A: Math. Theor. {\bf 52} 285201) where we showed that for $λ\gg l_s$ the soliton oscillates in an effective potential, while for $λ\ll l_s$ it moves uniformly as a free particle. In this paper we focus on the $κ$ dependence of the transition from oscillatory to particle behavior and explicitly compare the curves of the transition regime found in the collective coordinate approximation as a function of $r$ and $K$ when $κ=1/2,1,2$ at fixed value of the frequency $ω_0$. Since the solitary wave gets narrower for fixed $ω_0$ as a function of $κ$, we expect and indeed find that the regime where the solitary wave is trapped is extended as we increase $κ$.

preprint2019arXiv

Quantum Electron Transport in Degenerate Donor-Acceptor Systems

We develop a mathematically rigorous theory for the quantum transfer processes in degenerate donor-acceptor dimers in contact with a thermal environment. We calculate explicitly the transfer rates and the acceptor population efficiency. The latter depends critically on the initial donor state. We show that quantum coherence in the initial state enhances the transfer process. If the electron is initially shared coherently by the donor levels then the efficiency can reach values close to 100%, while an incoherent initial donor state will significantly suppress the efficiency. The results are useful for a better understanding of the quantum electron transport in many chemical, solid state, and biological systems with complex degenerate and quasi-degenerate energy landscapes.

preprint2019arXiv

Wide Class of Logarithmic Potentials with Power-Tower Kink Tails

We present a wide class of potentials which admit kinks and corresponding mirror kinks with either a power law or an exponential tail at the two extreme ends and a power-tower form of tails at the two neighbouring ends, i.e. of the forms $ette$ or $pttp$ where $e, p$ and $t$ denote exponential, power law and power-tower tail, respectively. We analyze kink stability equation in all these cases and show that there is no gap between the zero mode and the beginning of the continuum. Finally, we provide a recipe for obtaining logarithmic potentials with power-tower kink tails and estimate kink-kink interaction strength.

preprint2018arXiv

Connections Between Complex PT-Invariant Solutions and Complex Periodic Solutions of Several Nonlinear Equations

We point out novel connections between complex PT-invariant solutions of several nonlinear equations such as $ϕ^4$, $ϕ^6$, sine-Gordon, hyperbolic sine-Gordon, double sine-Gordon, double hyperbolic sine-Gordon, mKdV, etc. We then use these connections to obtain several new complex PT-invariant periodic solutions of some of these equations.

preprint2018arXiv

Noise-Assisted Quantum Exciton and Electron Transfer in Bio-Complexes with Finite Donor and Acceptor Bandwidths

We present an analytic and numerical study of noise-assisted quantum exciton (electron) transfer (ET) in a bio-complex, consisting of an electron donor and acceptor (a dimer), modeled by interacting continuous electron bands of finite widths. The interaction with the protein-solvent environment is modeled by a stationary stochastic process (noise) acting on all the donor and acceptor energy levels. We start with discrete energy levels for both bands. Then, by using a continuous {limit} for the electron spectra, we derive integro-differential equations for ET dynamics between two bands. Finally, we derive from these equations rate-type differential equations for ET dynamics. We formulate the conditions of validity of the rate-type equations. We consider different regions of parameters characterizing the widths of the donor and acceptor bands and the strength of the dimer-noise interaction. For a simplified model with a single energy level donor and a continuous acceptor band, we derive a generalized analytic expression and provide numerical simulations for the ET rate. They are consistent with Wigner-Weisskopf, Förster-type, and Marcus-type expressions, in their corresponding regime of parameters.Our analytic results are confirmed by numerical simulations. We demonstrate how our theoretical results are modified {when both the donor and the acceptor are described by finite bands}. We also show that, for a relatively wide acceptor band, the efficiency of the ET from donor to acceptor can be close to 100% for a broad range of noise amplitudes, for both "downhill" and "uphill" ET, for sharp and flat redox potentials, and for reasonably short times. We discuss possible experimental implementations of our approach with application to bio-complexes.

preprint2016arXiv

Dynamics of Dirac strings and monopole-like excitations in chiral magnets under a current drive

Skyrmion lines in metallic chiral magnets carry an emergent magnetic field experienced by the conduction electrons. The inflow and outflow of this field across a closed surface is not necessarily equal, thus it allows for the existence of emergent monopoles. One example is a segment of skyrmion line inside a crystal, where a monopole and antimonopole pair is connected by the emergent magnetic flux line. This is a realization of Dirac string-like excitations. Here we study the dynamics of monopoles in chiral magnets under an electric current. We show that in the process of creation of skyrmion lines, skyrmion line segments are first created via the proliferation of monopoles and antimonopoles. Then these line segments join and span the whole system through the annihilation of monopoles. The skyrmion lines are destroyed via the proliferation of monopoles and antimonopoles at high currents, resulting in a chiral liquid phase. We also propose to create the monopoles in a controlled way by applying an inhomogenous current to a crystal. Remarkably, an electric field component in the magnetic field direction proportional to the current squared in the low current region is induced by the motion of distorted skyrmion lines, in addition to the Hall and longitudinal voltage. The existence of monopoles can be inferred from transport or imaging measurements.

preprint2016arXiv

Helicoidal Graphene Nanoribbons: Chiraltronics

We present a calculation of the effective geometry-induced quantum potential for the carriers in graphene shaped as a helicoidal nanoribbon. In this geometry the twist of the nanoribbon plays the role of an effective transverse electric field in graphene and this is reminiscent of the Hall effect. However, this effective electric field has a different sign for the two iso-spin states and translates into a mechanism to separate the two chiral species on the opposing rims of the nanoribbon. Iso-spin transitions are expected with the emission or absorption of microwave radiation which could be adjusted to be in the THz region.

preprint2016arXiv

Non-hermitian quantum thermodynamics

Thermodynamics is the phenomenological theory of heat and work. Here we analyze to what extent quantum thermodynamic relations are immune to the underlying mathematical formulation of quantum mechanics. As a main result, we show that the Jarzynski equality holds true for all non-hermitian quantum systems with real spectrum. This equality expresses the second law of thermodynamics for isothermal processes arbitrarily far from equilibrium. In the quasistatic limit however, the second law leads to the Carnot bound which is fulfilled even if some eigenenergies are complex provided they appear in conjugate pairs. Furthermore, we propose two setups to test our predictions. Namely with strongly interacting excitons and photons in a semiconductor microcavity and in the non-hermitian tight-binding model.

preprint2016arXiv

Novel PT-invariant Kink and Pulse Solutions For a Large Number of Real Nonlinear Equations

For a large number of real nonlinear equations, either continuous or discrete, integrable or nonintegrable, uncoupled or coupled, we show that whenever a real nonlinear equation admits kink solutions in terms of $\tanh βx$, where $β$ is the inverse of the kink width, it also admits solutions in terms of the PT-invariant combinations $\tanh 2βx \pm i \sech 2 βx$, i.e. the kink width is reduced by half to that of the real kink solution. We show that both the kink and the PT-invariant kink are linearly stable and obtain expressions for the zero mode in the case of several PT-invariant kink solutions. Further, for a number of real nonlinear equations we show that whenever a nonlinear equation admits periodic kink solutions in terms of $\sn(x,m)$, it also admits periodic solutions in terms of the PT-invariant combinations $\sn(x,m) \pm i \cn(x,m)$ as well as $\sn(x,m)\pm i \dn(x,m)$. Finally, for coupled equations we show that one cannot only have complex PT-invariant solutions with PT eigenvalue $+1$ or $-1$ in both the fields but one can also have solutions with PT eigenvalue $+1$ in one field and $-1$ in the other field.

preprint2016arXiv

Repeatability of measurements: Equivalence of hermitian and non-hermitian observables

A non-commuting measurement transfers, via the apparatus, information encoded in a system's state to the external "observer". Classical measurements determine properties of physical objects. In the quantum realm, the very same notion restricts the recording process to orthogonal states as only those are distinguishable by measurements. Therefore, even a possibility to describe physical reality by means of non-hermitian operators should \emph{volens nolens} be excluded as their eigenstates are not orthogonal. Here, we show that non-hermitian operators with real spectrum can be treated within the standard framework of quantum mechanics. Furthermore, we propose a quantum canonical transformation that maps hermitian systems onto non-hermitian ones. Similar to classical inertial forces this transformation is accompanied by an energetic cost pinning the system on the unitary path.

preprint2016arXiv

Stability of exact solutions of the nonlinear Schroedinger equation in an external potential having supersymmetry and parity-time symmetry

We discuss the stability properties of the solutions of the general nonlinear Schroedinger equation (NLSE) in 1+1 dimensions in an external potential derivable from a parity-time (PT) symmetric superpotential $W(x)$ that we considered earlier [Kevrekedis et al Phys. Rev. E 92, 042901 (2015)]. In particular we consider the nonlinear partial differential equation $\{ i \partial_t + \partial_x^2 - V^{-}(x) +| ψ(x,t) |^{2κ} \} \, ψ(x,t) = 0$, for arbitrary nonlinearity parameter $κ$. We study the bound state solutions when $V^{-}(x) = (1/4- b^2)$ sech$^2(x)$, which can be derived from two different superpotentials $W(x)$, one of which is complex and $PT$ symmetric. Using Derrick's theorem, as well as a time dependent variational approximation, we derive exact analytic results for the domain of stability of the trapped solution as a function of the depth $b^2$ of the external potential. We compare the regime of stability found from these analytic approaches with a numerical linear stability analysis using a variant of the Vakhitov-Kolokolov (V-K) stability criterion. The numerical results of applying the V-K condition give the same answer for the domain of stability as the analytic result obtained from applying Derrick's theorem. Our main result is that for $κ>2$ a new regime of stability for the exact solutions appears as long as $b > b_{crit}$, where $b_{crit}$ is a function of the nonlinearity parameter $κ$. In the absence of the potential the related solitary wave solutions of the NLSE are unstable for $κ>2$.

preprint2016arXiv

Variational Approach to studying solitary waves in the nonlinear Schrodinger equation with Complex Potentials

We discuss the behavior of solitary wave solutions of the nonlinear Schr{ö}dinger equation (NLSE) as they interact with complex potentials, using a four parameter variational approximation based on a dissipation functional formulation of the dynamics. We concentrate on spatially periodic potentials with the periods of the real and imaginary part being either the same or different. Our results for the time evolution of the collective coordinates of our variational ansatz are in good agreement with direct numerical simulation of the NLSE. We compare our method with a collective coordinate approach of Kominis and give examples where the two methods give qualitatively different answers. In our variational approach, we are able to give analytic results for the small oscillation frequency of the solitary wave oscillating parameters which agree with the numerical solution of the collective coordinate equations. We also verify that instabilities set in when the slope of $dp(t)/dv(t)$ becomes negative when plotted parametrically as a function of time, where $p(t)$ is the momentum of the solitary wave and $v(t)$ the velocity.

preprint2015arXiv

Excitonic effects in two-dimensional semiconductors: Path integral Monte Carlo approach

One of the most striking features of novel 2D semiconductors (e.g., transition metal dichalcogenide monolayers or phosphorene) is a strong Coulomb interaction between charge carriers resulting in large excitonic effects. In particular, this leads to the formation of multi-carrier bound states upon photoexcitation (e.g., excitons, trions and biexcitons), which could remain stable at near-room temperatures and contribute significantly to optical properties of such materials. In the present work we have used the Path Integral Monte Carlo methodology to numerically study properties of multi-carrier bound states in 2D semiconductors. Specifically, we have accurately investigated and tabulated the dependence of single exciton, trion and biexciton binding energies on the strength of dielectric screening, including the limiting cases of very strong and very weak screening. The results of this work are potentially useful in the analysis of experimental data and benchmarking of theoretical and computational models.

preprint2015arXiv

Integrable and Superintegrable Klein-Gordon and Schrödinger Type Dimers

A $PT$-symmetric dimer is a two-site nonlinear oscillator or a nonlinear Schrödinger dimer where one site loses and the other site gains energy at the same rate. We present a wide class of integrable oscillator type dimers whose Hamiltonian is of arbitrary even order. Further, we also present a wide class of integrable and superintegrable nonlinear Schrödinger type dimers where again the Hamiltonian is of arbitrary even order.

preprint2015arXiv

Jarzynski equality in $\mathcal{PT}$-symmetric quantum mechanics

We show that the quantum Jarzynski equality generalizes to $\mathcal{PT}$-symmetric quantum mechanics with unbroken $\mathcal{PT}$-symmetry. In the regime of broken $\mathcal{PT}$-symmetry the Jarzynski equality does not hold as also the $\mathcal{CPT}$-norm is not preserved during the dynamics. These findings are illustrated for an experimentally relevant system -- two coupled optical waveguides. It turns out that for these systems the phase transition between the regimes of unbroken and broken $\mathcal{PT}$-symmetry is thermodynamically inhibited as the irreversible work diverges at the critical point.

preprint2015arXiv

More Opportunities than Wealth: A Network of Power and Frustration

We introduce a minimal agent-based model to qualitatively conceptualize the allocation of limited wealth among more abundant opportunities. We study the interplay of power, satisfaction and frustration in distribution, concentration, and inequality of wealth. Our framework allows us to compare subjective measures of frustration and satisfaction to collective measures of fairness in wealth distribution, such as the Lorenz curve and the Gini index. We find that a completely libertarian, law-of-the-jungle setting, where every agent can acquire wealth from, or lose wealth to, anybody else invariably leads to a complete polarization of the distribution of wealth vs. opportunity. The picture is however dramatically modified when hard constraints are imposed over agents, and they are limited to share wealth with neighbors on a network. We then propose an out of equilibrium dynamics {\it of} the networks, based on a competition between power and frustration in the decision-making of agents that leads to network coevolution. We show that the ratio of power and frustration controls different dynamical regimes separated by kinetic transitions and characterized by drastically different values of the indices of equality. The interplay of power and frustration leads to the emergence of three self-organized social classes, lower, middle, and upper class, whose interactions drive a cyclical regime.

preprint2015arXiv

Non-circular skyrmion and its anisotropic response in thin films of chiral magnets under tilted magnetic field

We study the equilibrium and dynamical properties of skyrmions in thin films of chiral magnets with oblique magnetic field. The shape of an individual skyrmion is non-circular and the skyrmion density decreases with the tilt angle from the normal of films. As a result, the interaction between two skyrmions depends on the relative angle between them in addition to their separation. The triangular lattice of skyrmions under a perpendicular magnetic field is distorted into a centered rectangular lattice for a tilted magnetic field. For a low skyrmion density, skyrmions form a chain like structure. The dynamical response of the non-circular skyrmions depends on the direction of external currents.

preprint2015arXiv

Nonlinear Dirac equation solitary waves in the presence of external driving forces

We consider the nonlinear Dirac (NLD) equation in 1+1 dimension with scalar-scalar self-interaction in the presence of external forces as well as damping of the form $ f(x,t) - i μγ^0 Ψ$, where both $f$ and $Ψ$ are two-component spinors. We develop an approximate variational approach using collective coordinates (CC) for studying the time dependent response of the solitary waves to these external forces. This approach predicts intrinsic oscillations of the solitary waves, i.e. the amplitude, width and phase all oscillate with the same frequency. The translational motion is also affected, because the soliton position oscillates around a mean trajectory. We then compare the results of the variational approximation with numerical simulations of the NLD equation, and find a good agreement, if we take into account a certain linear excitation with specific wavenumber that is excited together with the intrinsic oscillations such that the momentum in a transformed NLD equation is conserved. We also solve explicitly the CC equations of the variational approximation in the non-relativistic regime for a homogeneous external force and obtain excellent agreement with the numerical solution of the CC equations.

preprint2015arXiv

Novel PT-invariant Solutions For a Large Number of Real Nonlinear Equations

For a large number of real nonlinear equations, either continuous or discrete, integrable or nonintegrable, we show that whenever a real nonlinear equation admits a solution in terms of $\sech x$, it also admits solutions in terms of the PT-invariant combinations $\sech x \pm i \tanh x$. Further, for a number of real nonlinear equations we show that whenever a nonlinear equation admits a solution in terms $\sech^2 x$, it also admits solutions in terms of the PT-invariant combinations $\sech^2 x \pm i \sech x \tanh x$. Besides, we show that similar results are also true in the periodic case involving Jacobi elliptic functions.

preprint2015arXiv

PT-symmetric phase in kagome photonic lattices

Kagome lattice is a two-dimensional network of corner-sharing triangles and is often associated with geometrical frustration. In particular, the frustrated coupling between waveguide modes in a kagome array leads to a dispersionless flat band consisting of spatially localized modes. Here we propose a complex photonic lattice by placing $\mathcal{PT}$-symmetric dimers at the kagome lattice points. Each dimer corresponds to a pair of strongly coupled waveguides. With balanced arrangement of gain and loss on individual dimers, the system exhibits a $\mathcal{PT}$-symmetric phase for finite gain/loss parameter up to a critical value. The beam evolution in this complex kagome waveguide array exhibits a novel oscillatory rotation of optical power along the propagation distance. Long-lived local chiral structures originating from the nearly flat bands of the kagome structure are observed when the lattice is subject to a narrow beam excitation.

preprint2015arXiv

Quantum work statistics of charged Dirac particles in time-dependent fields

The quantum Jarzynski equality is an important theorem of modern quantum thermodynamics. We show that the Jarzynski equality readily generalizes to relativistic quantum mechanics described by the Dirac equation. After establishing the conceptual framework we solve a pedagogical, yet experimentally relevant, system analytically. As a main result we obtain the exact quantum work distributions for charged particles traveling through a time-dependent vector potential evolving under Schrödinger as well as under Dirac dynamics, and for which the Jarzynski equality is verified. Special emphasis is put on the conceptual and technical subtleties arising from relativistic quantum mechanics.

preprint2015arXiv

Recent Progress in the thermodynamics of ferrotoroidic materials

Recent theoretical and experimental progress on the study of ferrotoroidic materials is reviewed. The basic field equations are first described and then the expressions for magnetic toroidal moment and toroidization are derived. Relevant materials and experimental observation of magnetic toroidal moment and toroidal domains are summarized next. The thermodynamics of such magnetic materials is discussed in detail with examples of ferrorotoidic phase transition studied using Landau modelling. Specifically, an example of application of Landau modelling to the study of toroidocaloric effect is also provided. Recent results of polar nanostructures with electrical toroidal moment are finally reviewed.

preprint2015arXiv

Reply to the Comment by Zhang and Li [J. Math. Phys. 56, 084101, 2015] on the Paper [J. Math. Phys. {\bf 56}, 032104, 2015]}

In a recent paper Zhang and Li have doubted our claim that whenever a nonlinear equation has solutions in terms of the Jacobi elliptic functions $\cn(x,m)$ and $\dn(x,m)$, then the same nonlinear equation will necessarily also have solutions in terms of $\dn(x,m) \pm \sqrt{m} \cn(x,m)$. We point out the flaw in their argument and show why our assertion is indeed valid.

preprint2015arXiv

Skyrmion fractionalization and merons in chiral magnets with easy-plane anisotropy

We study the equilibrium phase diagram of ultrathin chiral magnets with easy-plane anisotropy $A$. The vast triangular skyrmion lattice phase that is stabilized by an external magnetic field evolves continuously as a function of increasing $A$ into a regime in which nearest-neighbor skyrmions start overlapping with each other. This overlap leads to a continuous reduction of the skyrmion number from its quantized value $Q=1$ and to the emergence of antivortices at the center of the triangles formed by nearest-neighbor skyrmions. The antivortices also carry a small "skyrmion number" $Q_A \ll 1$ that grows as a function of increasing $A$. The system undergoes a first order phase transition into a square vortex-antivortex lattice at a critical value of $A$. Finally, a canted ferromagnetic state becomes stable through another first order transition for a large enough anisotropy $A$. Interestingly enough, this first order transition is accompanied by {\it metastable} meron solutions.

preprint2015arXiv

Solitary waves of a PT-symmetric Nonlinear Dirac equation

In the present work, we consider a prototypical example of a PT-symmetric Dirac model. We discuss the underlying linear limit of the model and identify the threshold of the PT-phase transition in an analytical form. We then focus on the examination of the nonlinear model. We consider the continuation in the PT-symmetric model of the solutions of the corresponding Hamiltonian model and find that the solutions can be continued robustly as stable ones all the way up to the PT-transition threshold. In the latter, they degenerate into linear waves. We also examine the dynamics of the model. Given the stability of the waveforms in the PT-exact phase we consider them as initial conditions for parameters outside of that phase. We find that both oscillatory dynamics and exponential growth may arise, depending on the size of the corresponding "quench". The former can be characterized by an interesting form of bi-frequency solutions that have been predicted on the basis of the SU(1,1) symmetry. Finally, we explore some special, analytically tractable, but not PT-symmetric solutions in the massless limit of the model.

preprint2014arXiv

AC current generation in chiral magnetic insulators and skyrmion motion induced by the spin Seebeck effect

We show that a temperature gradient induces an ac electric current in multiferroic insulators when the sample is embedded in a circuit. We also show that a thermal gradient can be used to move magnetic skyrmions in insulating chiral magnets: the induced magnon flow from the hot to the cold region drives the skyrmions in the opposite direction via a magnonic spin transfer torque. Both results are combined to compute the effect of skyrmion motion on the ac current generation and demonstrate that skyrmions in insulators are a promising route for spin caloritronics applications.

preprint2014arXiv

Bent waveguides for matter-waves: supersymmetric potentials and reflectionless geometries

Non-zero curvature in a waveguide leads to the appearance of an attractive quantum potential which crucially affects the dynamics in matter-wave circuits. Using methods of supersymmetric quantum mechanics, pairs of bent waveguides are found whose geometry-induced potentials share the same scattering properties. As a result, reflectionless waveguides, dual to the straight waveguide, are identified. Strictly isospectral waveguides are also found by modulating the depth of the trapping potential. Numerical simulations are used to demonstrate the efficiency of these approaches in tailoring and controlling curvature-induced quantum-mechanical effects.

preprint2014arXiv

Dynamics of skyrmions in chiral magnets: dynamic phase transitions and equation of motion

We study the dynamics of skyrmions in a metallic chiral magnet. First we show that skyrmions can be created dynamically by destabilizing the ferromagnetic background state through a spin polarized current. We then treat skyrmions as rigid particles and derive the corresponding equation of motion. The dynamics of skyrmions is dominated by the Magnus force, which accounts for the weak pinning of skyrmions observed in experiments. Finally we discuss the quantum motion of skyrmions.

preprint2014arXiv

Internal modes of a skyrmion in the ferromagnetic background in chiral magnets

A spin texture called skyrmion has been recently observed in certain chiral magnets without inversion symmetry. The observed skyrmions are extended objects with typical linear sizes of 10 nm to 100 nm that contain $10^3$ to $10^5$ spins and can be deformed in response to external perturbations. Weak deformations are characterized by the internal modes which are localized around the skyrmion center. Knowledge of internal modes is crucial to assess the stability and rigidity of these topological textures. Here we compute the internal modes of a skyrmion in the ferromagnetic background state by numerical diagonalization of the dynamical matrix. We find several internal modes below the magnon continuum, such as the mode corresponding to the translational motion, and different kinds of breathing modes. The number of internal modes is larger for lower magnetic fields. Indeed, several modes become gapless in the low field region indicating that the single skyrmion solution becomes unstable, although a skyrmion lattice is thermodynamically stable. On the other hand, at high fields only three internal modes exist and the skyrmion is stable even at higher fields when the ferromagnetic state is thermodynamically stable. We also show that the presence of out-of-plane easy-axis anisotropy stabilizes the single skyrmion solution. Finally, we discuss the effects of damping and possible experimental observations of these internal modes.

preprint2014arXiv

Intrinsic Localized Modes of a Classical Discrete Anisotropic Heisenberg Ferromagnetic Spin Chain

We report several exact intrinsic localized mode solutions of the classical spin evolution equation of a one-dimensional anisotropic Heisenberg ferromagnetic spin chain in terms of Jacobian elliptic functions. These include one, two and three spin excitations. All these solutions have smooth anticontinuum limits. Their linear stability and semiclassical quantization are also discussed briefly.

preprint2014arXiv

Periodic and Hyperbolic Soliton Solutions of a Number of Nonlocal PT-Symmetric Nonlinear Equations

For a number of nonlocal nonlinear equations such as nonlocal, nonlinear Schrödinger equation (NLSE), nonlocal Ablowitz-Ladik (AL), nonlocal, saturable discrete NLSE (DNLSE), coupled nonlocal NLSE, coupled nonlocal AL and coupled nonlocal, saturable DNLSE, we obtain periodic solutions in terms of Jacobi elliptic functions as well as the corresponding hyperbolic soliton solutions. Remarkably, in all the six cases, we find that unlike the corresponding local cases, all the nonlocal models simultaneously admit both the bright and the dark soliton solutions. Further, in all the six cases, not only $\dn(x,m)$ and $\cn(x,m)$ but even their linear superposition is shown to be an exact solution. Finally, we show that the coupled nonlocal NLSE not only admits solutions in terms of Lamé polynomials of order 1, but it also admits solutions in terms of Lamé polynomials of order 2, even though they are not the solutions of the uncoupled nonlocal problem. We also remark on the possible integrability in certain cases.

preprint2014arXiv

PT-symmetry breaking with divergent potentials: lattice and continuum cases

We investigate the parity- and time-reversal ($\mathcal{PT}$)-symmetry breaking in lattice models in the presence of long-ranged, non-hermitian, $\mathcal{PT}$-symmetric potentials that remain finite or become divergent in the continuum limit. By scaling analysis of the fragile $\mathcal{PT}$ threshold for an open finite lattice, we show that continuum loss-gain potentials $V_α(x)\propto i |x|^α\mathrm{sign}(x)$ have a positive $\mathcal{PT}$-breaking threshold for $α>-2$, and a zero threshold for $α\leq -2$. When $α<0$ localized states with complex (conjugate) energies in the continuum energy-band occur at higher loss-gain strengths. We investigate the signatures of $\mathcal{PT}$-symmetry breaking in coupled waveguides, and show that the emergence of localized states dramatically shortens the relevant time-scale in the $\mathcal{PT}$-symmetry broken region.

preprint2014arXiv

Quasi-One-Dimensional Thermal Breakage

Breakage is generally understood in mechanical terms, yet nano-structures can rupture not only under external loads but also via thermal activation. Here we treat in a general framework the statistical mechanics of thermally induced breakage at the nano-scale for one-dimensional systems. We test it on a simple approximation and find that the probability of breakage controls distinct regimes, characterized by sharp crossovers and narrow peaks in the thermal fluctuations and specific heat. Our work provides predictions on clustering of new phases, of relevance in nano-fabrication.

preprint2014arXiv

Stability of solitary waves in the nonlinear Dirac equation with arbitrary nonlinearity

We consider the nonlinear Dirac equation in 1+1 dimension with scalar-scalar self interaction $ \frac{g^2}{κ+1} ({\bar Ψ} Ψ)^{κ+1}$ and with mass $m$. Using the exact analytic form for rest frame solitary waves of the form $Ψ(x,t) = ψ(x) e^{-i ωt}$ for arbitrary $ κ$, we discuss the validity of various approaches to understanding stability that were successful for the nonlinear Schrödinger equation. In particular we study the validity of a version of Derrick's theorem, the criterion of Bogolubsky as well as the Vakhitov-Kolokolov criterion, and find that these criteria yield inconsistent results. Therefore, we study the stability by numerical simulations using a recently developed 4th-order operator splitting integration method. For different ranges of $κ$ we map out the stability regimes in $ω$. We find that all stable nonlinear Dirac solitary waves have a one-hump profile, but not all one-hump waves are stable, while all waves with two humps are unstable. We also find that the time $t_c$, it takes for the instability to set in, is an exponentially increasing function of $ω$ and $t_c$ decreases monotonically with increasing $κ$.

preprint2014arXiv

Successive phase transitions and kink solutions in $ϕ^{8}$, $ϕ^{10}$, and $ϕ^{12}$ field theories

We obtain exact solutions for kinks in $ϕ^{8}$, $ϕ^{10}$ and $ϕ^{12}$ field theories with degenerate minima, which can describe a second-order phase transition followed by a first-order one, a succession of two first-order phase transitions and a second-order phase transition followed by two first-order phase transitions, respectively. Such phase transitions are known to occur in ferroelastic and ferroelectric crystals and in meson physics. In particular, we find that the higher-order field theories have kink solutions with algebraically-decaying tails and also asymmetric cases with mixed exponential-algebraic tail decay, unlike the lower-order $ϕ^4$ and $ϕ^6$ theories. Additionally, we construct distinct kinks with equal energies in all three field theories considered, and we show the co-existence of up to three distinct kinks (for a $ϕ^{12}$ potential with six degenerate minima). We also summarize phonon dispersion relations for these systems, showing that the higher-order field theories have specific cases in which only nonlinear phonons are allowed. For the $ϕ^{10}$ field theory, which is a quasi-exactly solvable (QES) model akin to $ϕ^6$, we are also able to obtain three analytical solutions for the classical free energy as well as the probability distribution function in the thermodynamic limit.

preprint2014arXiv

Superposition of Elliptic Functions as Solutions For a Large Number of Nonlinear Equations

For a large number of nonlinear equations, both discrete and continuum, we demonstrate a kind of linear superposition. We show that whenever a nonlinear equation admits solutions in terms of both Jacobi elliptic functions $\cn(x,m)$ and $\dn(x,m)$ with modulus $m$, then it also admits solutions in terms of their sum as well as difference. We have checked this in the case of several nonlinear equations such as the nonlinear Schrödinger equation, MKdV, a mixed KdV-MKdV system, a mixed quadratic-cubic nonlinear Schrödinger equation, the Ablowitz-Ladik equation, the saturable nonlinear Schrödinger equation, $λϕ^4$, the discrete MKdV as well as for several coupled field equations. Further, for a large number of nonlinear equations, we show that whenever a nonlinear equation admits a periodic solution in terms of $\dn^2(x,m)$, it also admits solutions in terms of $\dn^2(x,m) \pm \sqrt{m} \cn(x,m) \dn(x,m)$, even though $\cn(x,m) \dn(x,m)$ is not a solution of these nonlinear equations. Finally, we also obtain superposed solutions of various forms for several coupled nonlinear equations.

preprint2013arXiv

Driven skyrmions and dynamical transitions in chiral magnets

We study the dynamics of skyrmions in chiral magnets in the presence of a spin polarized current. The motion of skyrmions in the ferromagnetic background excites spin waves and contributes to additional damping. At a large current, the spin wave spectrum becomes gapless and skyrmions are created dynamically from the ferromagnetic state. At an even higher current, these skyrmions are strongly deformed due to the damping and become unstable at a threshold current, leading to a chiral liquid. We show how skyrmions can be created by increasing the current in the magnetic spiral state. We then construct a dynamic phase diagram for a chiral magnet with a current. The instability transitions between different states can be observed as experimentally clear signatures in the transport measurements, such as jumps and hysteresis.

preprint2013arXiv

Higher Order Spin Noise Statistics

The optical spin noise spectroscopy (SNS) is a minimally invasive route towards obtaining dynamical information about electrons and atomic gases by measuring mesoscopic time-dependent spin fluctuations. Recent improvements of the sensitivity of SNS should make it possible to observe higher order spin correlators at thermodynamic equilibrium. We develop theoretical methods to explore higher order (3rd and 4th) cumulants of the spin noise in the frequency domain. We make predictions for the possible functional form of these correlators in single quantum dot experiments and then apply the method of the stochastic path integral to estimate effects of many-body interactions.

preprint2013arXiv

Linear Superposition for a Large Number of Nonlinear Equations

We demonstrate a kind of linear superposition for a large number of nonlinear equations, both continuum and discrete. In particular, we show that whenever a nonlinear equation admits solutions in terms of Jacobi elliptic functions $\cn(x,m)$ and $\dn(x,m)$, then it also admits solutions in terms of their sum as well as difference, i.e. $\dn(x,m) \pm \sqrt{m}\, \cn(x,m)$. Further, we also show that whenever a nonlinear equation admits a solution in terms of $\dn^2(x,m)$, it also has solutions in terms of $\dn^2(x,m) \pm \sqrt{m}\, \cn(x,m)\, \dn(x,m)$ even though $\cn(x,m)\, \dn(x,m)$ is not a solution of that nonlinear equation. Finally, we obtain similar superposed solutions in coupled theories.

preprint2013arXiv

Manipulation of skyrmions in nanodisks with a current pulse and skyrmion rectifier

A skyrmion in a nanosized disk of a chiral magnet can be used as a bit of information. To this end, it is desirable to control the creation and removal of a skyrmion only by currents without using external magnetic fields. Here we propose to create a skyrmion by applying a current pulse to a nanodisk. The skyrmion can be removed from the disk by applying a dc current. We show that the dynamics of the created skyrmion can lead to a rectification effect, in which a dc voltage is generated by the motion of skyrmion in the presence of an ac current.

preprint2013arXiv

Nonlinear Dirac equation solitary waves in external fields

We consider the nonlinear Dirac equations (NLDE's) in 1+1 dimension with scalar-scalar self interaction $\frac{g^2}{κ+1} ({\bPsi} Ψ)^{κ+1}$ in the presence of various external electromagnetic fields. Starting from the exact solutions for the unforced problem we study the behavior of solitary wave solutions to the NLDE in the presence of a wide variety of fields in a variational approximation depending on collective coordinates which allows the position, width and phase of these waves to vary in time. We find that in this approximation the position $q(t)$ of the center of the solitary wave obeys the usual behavior of a relativistic point particle in an external field. For time independent external fields we find that the energy of the solitary wave is conserved but not the momentum which becomes a function of time. We postulate that similar to the nonlinear Schr{ö}dinger equation (NLSE) that a sufficient dynamical condition for instability to arise is that $ dP(t)/d \dq(t) < 0$. Here $P(t)$ is the momentum of the solitary wave, and $\dq$ is the velocity of the center of the wave in the collective coordinate approximation. We found for our choices of external potentials we always have $ dP(t)/d \dq(t) > 0$ so when instabilities do occur they are due to a different source. We investigate the accuracy of our variational approximation using numerical simulations of the NLDE and find that when the forcing term is small and we are in a regime where the solitary wave is stable, that the behavior of the solutions of the collective coordinate equations agrees very well with the numerical simulations.

preprint2013arXiv

Particle model for skyrmions in metallic chiral magnets: dynamics, pinning and creep

Recently spin textures called skyrmions have been discovered in certain chiral magnetic materials without spatial inversion symmetry, and have attracted enormous attention due to their promising application in spintronics since only a low applied current is necessary to drive their motion. When a conduction electron moves around the skyrmion, its spin is fully polarized by the spin texture and acquires a quantized phase; thus, the skyrmion yields an emergent electrodynamics that in turn determines skyrmion motion and gives rise to a finite Hall angle. While Skyrmions behave as particles, no particle level description of their dynamics exists to date. Such a model would have tremendous impact on understanding skyrmion dynamics by theoretical analysis and computational modeling. Here we derive the equation of motion from a microscopic continuum model and obtain the short-range interaction between skyrmions, and the interaction between skyrmions and defects. Skyrmions also experience a Magnus force perpendicular to their velocity due to the underlying emergent electromagnetic field. We validate the equation of motion by studying the depinning transition using both the particle and the continuum models. By using the particle description, we explain the recent experimental observations of the rotation of a skyrmion lattice in the presence of a temperature gradient. We also predict quantum and thermal creep motion of skyrmions in the pinning potential.

preprint2013arXiv

Spin noise spectroscopy of quantum dot molecules

We discuss advantages and limitations of the spin noise spectroscopy for characterization of interacting quantum dot systems on specific examples of individual singly and doubly charged quantum dot molecules (QDMs). It is shown that all the relevant parameters of the QDMs including tunneling amplitudes with spin-conserving and spin-non-conserving interactions, decoherence rates, Coulomb repulsions, anisotropic g-factors and the distance between the dots can be determined by measuring properties of the spin noise power spectrum.

preprint2013arXiv

Statistical Mechanics of a Discrete Schrödinger Equation with Saturable Nonlinearity

We study the statistical mechanics of the one-dimensional discrete nonlinear Schrödinger (DNLS) equation with saturable nonlinearity. Our study represents an extension of earlier work [Phys. Rev. Lett. {\bf 84}, 3740 (2000)] regarding the statistical mechanics of the one-dimensional DNLS equation with a cubic nonlinearity. As in this earlier study we identify the spontaneous creation of localized excitations with a discontinuity in the partition function. The fact that this phenomenon is retained in the saturable DNLS is non-trivial, since in contrast to the cubic DNLS whose nonlinear character is enhanced as the excitation amplitude increases, the saturable DNLS in fact becomes increasingly linear as the excitation amplitude increases. We explore the nonlinear dynamics of this phenomenon by direct numerical simulations.

preprint2012arXiv

Direct observation of magnetic phase coexistence and magnetization reversal in a Gd$_{0.67}$Ca$_{0.33}$MnO$_{3}$ thin film

We have investigated the ferrimagnetic domain structure in a Gd$_{0.67}$Ca$_{0.33}$MnO$_{3}$ thin film using magnetic force microscopy. We observe clear signs of phase separation, with magnetic islands embedded in a non-magnetic matrix. We also directly visualize the reversal of magnetization of ferrimagnetic domains as a function of temperature and attribute it to a change in the balance of magnetization of anti-aligned Mn and Gd sublattices.

preprint2012arXiv

Thermodynamics of Ferrotoroidic Materials: Toroidocaloric Effect

The three primary ferroics, namely ferromagnets, ferroelectrics and ferroelastics exhibit corresponding large (or even giant) magnetocaloric,electrocaloric and elastocaloric effects when a phase transition is induced by the application of an appropriate external field. Recently the suite of primary ferroics has been extended to include ferrotoroidic materials in which there is an ordering of toroidic moments in the form of magnetic vortex-like structures, examples being LiCo(PO_4)_3 and Ba_2CoGe_2O_7. In the present work we formulate the thermodynamics of ferrotoroidic materials. Within a Landau free energy framework we calculate the toroidocaloric effect by quantifying isothermal entropy change (or adiabatic temperature change) in the presence of an applied toroidic field when usual magnetization and polarization may also be present simultaneously. We also obtain a nonlinear Clausius-Clapeyron relation for phase coexistence.

preprint2011arXiv

Electronic properties of corrugated graphene, the Heisenberg principle and wormhole geometry in solid state

Adopting a purely two dimensional relativistic equation for graphene's carriers contradicts the Heisenberg uncertainty principle since it requires setting off-the-surface coordinate of a three-dimensional wavefunction to zero. Here we present a theoretical framework for describing graphene's massless relativistic carriers in accordance with this most fundamental of all quantum principles. A gradual confining procedure is used to restrict the dynamics onto a surface and normal to the surface parts and in the process the embedding of this surface into the three dimensional world is accounted for. As a result an invariant geometric potential arises in the surface part which scales linearly with the Mean curvature and shifts the Fermi energy of the material proportional to bending. Strain induced modification of the electronic properties or "straintronics" is clearly an important field of study in graphene. This opens a venue to producing electronic devices, MEMS and NEMS where the electronic properties are controlled by geometric means and no additional alteration of graphene is necessary. The appearance of this geometric potential also provides us with clues as to how quantum dynamics looks like in the curved space-time of general relativity. In this context, we explore a two-dimensional cross-section of the wormhole geometry realized with graphene as a solid state thought experiment.

preprint2011arXiv

Forced Nonlinear Schroedinger Equation with Arbitrary Nonlinearity

We consider the nonlinear Schr{ö}dinger equation (NLSE) in 1+1 dimension with scalar-scalar self interaction $\frac{g^2}{κ+1} (ψ^\star ψ)^{κ+1}$ in the presence of the external forcing terms of the form $r e^{-i(kx + θ)} -δψ$. We find new exact solutions for this problem and show that the solitary wave momentum is conserved in a moving frame where $v_k=2 k$. These new exact solutions reduce to the constant phase solutions of the unforced problem when $r \rightarrow 0.$ In particular we study the behavior of solitary wave solutions in the presence of these external forces in a variational approximation which allows the position, momentum, width and phase of these waves to vary in time. We show that the stationary solutions of the variational equations include a solution close to the exact one and we study small oscillations around all the stationary solutions. We postulate that the dynamical condition for instability is that $ dp(t)/d \dot{q} (t) < 0$, where $p(t)$ is the normalized canonical momentum $p(t) = \frac{1}{M(t)} \frac {\partial L}{\partial {\dot q}}$, and $\dot{q}(t)$ is the solitary wave velocity. Here $M(t) = \int dx ψ^\star(x,t) ψ(x,t)$. Stability is also studied using a "phase portrait" of the soliton, where its dynamics is represented by two-dimensional projections of its trajectory in the four-dimensional space of collective coordinates. The criterion for stability of a soliton is that its trajectory is a closed single curve with a positive sense of rotation around a fixed point. We investigate the accuracy of our variational approximation and these criteria using numerical simulations of the NLSE.

preprint2011arXiv

Mapping Between Generalized Nonlinear Schroedinger Equations and Neutral Scalar Field Theories and New Solutions of the Cubic-Quintic NLS Equation

We highlight an interesting mapping between the moving breather solutions of the generalized Nonlinear Schrodinger (NLS) equations and the static solutions of neutral scalar field theories. Using this connection, we then obtain several new moving breather solutions of the cubic-quintic NLS equation both with and without uniform phase in space. The stability of some stationary solutions is investigated numerically and the results confirmed via dynamical evolution.

preprint2011arXiv

On the properties of compacton-anticompacton collisions

We study the properties of compacton-anticompacton collision processes. We compare and con- trast results for the case of compacton-anticompacton solutions of the K(l, p) Rosenau-Hyman (RH) equation for l = p = 2, with compacton-anticompacton solutions of the L(l,p) Cooper-Shepard- Sodano (CSS) equation for p = 1 and l = 3. This study is performed using a Padé discretization of the RH and CSS equations. We find a significant difference in the behavior of compacton- anticompacton scattering. For the CSS equation, the scattering can be interpreted as "annihila- tion" as the wake left behind dissolves over time. In the RH equation, the numerical evidence is that multiple shocks form after the collision which eventually lead to "blowup" of the resulting waveform.

preprint2011arXiv

Solutions of Several Coupled Discrete Models in terms of Lame Polynomials of Arbitrary Order

Coupled discrete models abound in several areas of physics. Here we provide an extensive set of exact quasiperiodic solutions of a number of coupled discrete models in terms of Lamé polynomials of arbitrary order. The models discussed are (i) coupled Salerno model, (ii) coupled Ablowitz-Ladik model, (iii) coupled $ϕ^4$ model, and (iv) coupled $ϕ^6$ model. In all these cases we show that the coefficients of the Lamé polynomials are such that the Lamé polynomials can be reexpressed in terms of Chebyshev polynomials of the relevant Jacobi elliptic function.

preprint2011arXiv

Solutions of Several Coupled Discrete Models in terms of Lame Polynomials of Order One and Two

Coupled discrete models abound in several areas of physics. Here we provide an extensive set of exact quasiperiodic solutions of a number of coupled discrete models in terms of Lame polynomials of order one and two. Some of the models discussed are (i) coupled Salerno model, (ii) coupled Ablowitz-Ladik model, (iii) coupled saturated discrete nonlinear Schrodinger equation, (iv) coupled phi4 model, and (v) coupled phi6 model. Furthermore, we show that most of these coupled models in fact also possess an even broader class of exact solutions.

preprint2011arXiv

The robust $\mP\mT$-symmetric chain and properties of its Hermitian counterpart

We study the properties of a parity- and time-reversal- (PT) symmetric tight-binding chain of size N with position-dependent hopping amplitude. In contrast to the fragile PT-symmetric phase of a chain with constant hopping and imaginary impurity potentials, we show that, under very general conditions, our model is {\it always} in the PT-symmetric phase. We numerically obtain the energy spectrum and the density of states of such a chain, and show that they are widely tunable. By studying the size-dependence of inverse participation ratios, we show that although the chain is not translationally invariant, most of its eigenstates are extended. Our results indicate that tight-binding models with non-Hermitian PT-symmetric hopping have a robust PT-symmetric phase and rich dynamics which may be explored in coupled waveguides.

preprint2010arXiv

Exact Solutions of the Two-Dimensional Discrete Nonlinear Schrödinger Equation with Saturable Nonlinearity

We show that the two-dimensional, nonlinear Schrödinger lattice with a saturable nonlinearity admits periodic and pulse-like exact solutions. We establish the general formalism for the stability considerations of these solutions and give examples of stability diagrams. Finally, we show that the effective Peierls-Nabarro barrier for the pulse-like soliton solution is zero.

preprint2010arXiv

Robust and fragile PT-symmetric phases in a tight-binding chain

We study the phase-diagram of a parity and time-reversal (PT) symmetric tight-binding chain with $N$ sites and hopping energy $J$, in the presence of two impurities with imaginary potentials $\pm iγ$ located at arbitrary (P-symmetric) positions $(m, \bar{m}=N+1-m)$ on the chain where $m\leq N/2$. We find that except in the two special cases where impurities are either the farthest or the closest, the PT-symmetric region - defined as the region in which all energy eigenvalues are real - is algebraically fragile. We analytically and numerically obtain the critical impurity potential $γ_{PT}$ and show that $γ_{PT}\propto 1/N\rightarrow 0$ as $N\rightarrow\infty$ except in the two special cases. When the PT symmetry is spontaneously broken, we find that the maximum number of complex eigenvalues is given by $2m$. When the two impurities are the closest, we show that the critical impurity strength $γ_{PT}$ in the limit $N\rightarrow\infty$ approaches $J$ ($J/2$) provided that $N$ is even (odd). For an even $N$ the PT symmetry is maximally broken whereas for an odd $N$, it is sequentially broken. Our results show that the phase-diagram of a PT-symmetric tight-binding chain is extremely rich and that, in the continuum limit, this model may give rise to new PT-symmetric Hamiltonians.

preprint2010arXiv

Solitary wave in the Nonlinear Dirac Equation with arbitrary nonlinearity

We consider the nonlinear Dirac equations (NLDE's) in 1+1 dimension with scalar-scalar self interaction $\frac{g^2}{k+1} ({\bar Ψ} Ψ)^{k+1}$, as well as a vector-vector self interaction $\frac{g^2}{k+1} ({\bar Ψ} γ_μΨ\bPsi γ^μΨ)^{\frac{1}{2}(k+1)}$. We find the exact analytic form for solitary waves for arbitrary $k$ and find that they are a generalization of the exact solutions for the nonlinear Schrödinger equation (NLSE) and reduce to these solutions in a well defined nonrelativistic limit. We perform the nonrelativistic reduction and find the $1/2m$ correction to the NLSE, valid when $|ω-m |\ll 2m$, where $ω$ is the frequency of the solitary wave in the rest frame. We discuss the stability and blowup of solitary waves assuming the modified NLSE is valid and find that they should be stable for $k < 2$.

preprint2010arXiv

Stability and dynamical properties of Cooper-Shepard-Sodano compactons

Extending a Pade approximant method used for studying compactons in the Rosenau-Hyman (RH) equation, we study the numerical stability of single compactons of the Cooper-Shepard-Sodano (CSS) equation and their pairwise interactions. The CSS equation has a conserved Hamiltonian which has allowed various approaches for studying analytically the nonlinear stability of the solutions. We study three different compacton solutions and find they are numerically stable. Similar to the collisions between RH compactons, the CSS compactons reemerge with same coherent shape when scattered. The time evolution of the small-amplitude ripple resulting after scattering depends on the values of the parameters $l$ and $p$ characterizing the corresponding CSS equation. The simulation of the CSS compacton scattering requires a much smaller artificial viscosity to obtain numerical stability, than in the case of RH compactons propagation.

preprint2010arXiv

Stability and dynamical properties of Rosenau-Hyman compactons using Pade approximants

We present a systematic approach for calculating higher-order derivatives of smooth functions on a uniform grid using Padé approximants. We illustrate our findings by deriving higher-order approximations using traditional second-order finite-differences formulas as our starting point. We employ these schemes to study the stability and dynamical properties of K(2,2) Rosenau-Hyman (RH) compactons including the collision of two compactons and resultant shock formation. Our approach uses a differencing scheme involving only nearest and next-to-nearest neighbors on a uniform spatial grid. The partial differential equation for the compactons involves first, second and third partial derivatives in the spatial coordinate and we concentrate on four different fourth-order methods which differ in the possibility of increasing the degree of accuracy (or not) of one of the spatial derivatives to sixth order. A method designed to reduce roundoff errors was found to be the most accurate approximation in stability studies of single solitary waves, even though all derivates are accurate only to fourth order. Simulating compacton scattering requires the addition of fourth derivatives related to artificial viscosity. For those problems the different choices lead to different amounts of "spurious" radiation and we compare the virtues of the different choices.

preprint2010arXiv

Stable structures with high topological charge in nonlinear photonic quasicrystals

Stable vortices with topological charge of 3 and 4 are examined numerically and analytically in photonic quasicrystals created by interference of 5 as well as 8 beams, in the cases of cubic as well as saturable nonlinearities. These structures are experimentally realizable, including a prototypical example of a stable charge 4 vortex. Direct numerical simulations corroborate the analytical and numerical linear stability analysis predictions.

preprint2010arXiv

Staggered and short period solutions of the Saturable Discrete Nonlinear Schrödinger Equation

We point out that the nonlinear Schr{ö}dinger lattice with a saturable nonlinearity also admits staggered periodic as well as localized pulse-like solutions. Further, the same model also admits solutions with a short period. We examine the stability of these solutions and find that the staggered as well as the short period solutions are stable in most cases. We also show that the effective Peierls-Nabarro barrier for the pulse-like soliton solutions is zero.

preprint2010arXiv

Tuning electronic properties of corrugated graphene: confinement, curvature and band gap opening

It is shown that for monolayer graphene electrons are confined on a perfect two dimensional surface. The implications for the electronic properties of corrugated graphene are discussed in view of a derivation of the constrained relativistic dynamics for the massless carriers in two dimensions. Surface curvature is related to a series of phenomena with practical applications such as curvature induced p-n junctions, band gap opening and decoherence. We also establish a bending free energy by treating graphene as a soft electronic membrane.

preprint2009arXiv

Geometry induced potential on a 2D-section of a wormhole: catenoid

We show that a two dimensional wormhole geometry is equivalent to a catenoid, a minimal surface. We then obtain the curvature induced geometric potential and show that the ground state with zero energy corresponds to a reflectionless potential. By introducing an appropriate coordinate system we also obtain bound states for different angular momentum channels. Our findings can be realized in suitably bent bilayer graphene sheets with a neck or in a honeycomb lattice with an array of dislocations or in nanoscale waveguides in the shape of a catenoid.

preprint2009arXiv

Localized structures in Kagome lattices

We investigate the existence and stability of gap vortices and multi-pole gap solitons in a Kagome lattice with a defocusing nonlinearity both in a discrete case and in a continuum one with periodic external modulation. In particular, predictions are made based on expansion around a simple and analytically tractable anti-continuum (zero coupling) limit. These predictions are then confirmed for a continuum model of an optically-induced Kagome lattice in a photorefractive crystal obtained by a continuous transformation of a honeycomb lattice.

preprint2009arXiv

Thermal Stability of Strained Nanowires

Stranski-Krastanov strained islands undergo a shape anisotropy transition as they grow in size, finally evolving toward nanowires. This effect has been explained until now via simple energetic models that neglect thermodynamics. We investigate theoretically the stability of strained nanowires under thermal fluctuations of the long side. We find phase transitions from nanowires back to nanoislands as the temperature is increased and as the height of the nanostructure is raised or lowered and we predict regions of phase coexistence. Our results are general, but explain recent data on the growth of erbium silicide on a vicinal Si surface.

preprint2009arXiv

Thermally Induced Local Failures in Quasi-One-Dimensional Systems: Collapse in Carbon Nanotubes, Necking in Nanowires and Opening of Bubbles in DNA

We present a general framework to explore thermally activated failures in quasi one dimensional systems. We apply it to the collapse of carbon nanotubes, the formation of bottlenecks in nanowires, both of which limit conductance, and the opening of local regions or "bubbles" of base pairs in strands of DNA that are relevant for transcription and danaturation. We predict an exponential behavior for the probability of the opening of bubbles in DNA, the average distance between flattened regions of a nanotube or necking in a nanowire as a monotonically decreasing function of temperature, and compute a temperature below which these events become extremely rare. These findings are difficult to obtain numerically, however, they could be accessible experimentally.

preprint2008arXiv

Compactons in PT-symmetric generalized Korteweg-de Vries Equations

In an earlier paper Cooper, Shepard, and Sodano introduced a generalized KdV equation that can exhibit the kinds of compacton solitary waves that were first seen in equations studied by Rosenau and Hyman. This paper considers the PT-symmetric extensions of the equations examined by Cooper, Shepard, and Sodano. From the scaling properties of the PT-symmetric equations a general theorem relating the energy, momentum, and velocity of any solitary-wave solution of the generalized KdV equation is derived, and it is shown that the velocity of the solitons is determined by their amplitude, width, and momentum.

preprint2005arXiv

Exact Elliptic Compactons in Generalized Korteweg-DeVries Equations

We derive a general theorem relating the energy, momentum and velocity of any solitary wave solution of the generalized KdV equation which enables us to relate the amplitude, width, and momentum to the velocity of these solutions. We obtain the general condition for linear and Lyapunov stability. We then obtain a two parameter family of exact solutions to these equations which include elliptic and hyper-elliptic compacton solutions. For this general family we explicitly verify both the theorem and the stability criteria.

preprint2002arXiv

Glassy behavior in systems with Kac-type step-function interaction

We study a system with a weak, long-range repulsive Kac-type step-function interaction within the framework of a replicated effective $ϕ^4$-theory. The occurrence of extensive configurational entropy, or an exponentially large number of metastable minima in the free energy (characteristic of a glassy state), is demonstrated. The underlying mechanism of mesoscopic patterning and defect organizations is discussed.

preprint1995arXiv

Statistical Mechanics of Nonlinear Coherent Structures: Kinks in the Phi^6 Model

We study the thermodynamics of kinks in the Phi^6 model using a Langevin code implemented on a massively parallel computer. This code can be used to study first order dynamical phase transitions which exhibit multiple length and time scales. The classical statistical mechanics of a 1+1-dimensional field theory reduces to a time-independent quantum problem in one dimension via the transfer integral method. Exact solutions of the Schrodinger equation exist for the Phi^6 potential (unlike the case for Phi^4) and can be used to check results from the simulations. The Phi^6 model is also much richer than the Phi^4 model in terms of the variety of coherent structures and possible phases accompanying a phase transition. Specifically, we have calculated (in a one dimensional model) such quantities as the probability density function (PDF) and field-field correlation functions. These quantities help us understand the contribution to the specific heat from coherent structures such as domain walls (kinks) and other transformation structures as opposed to the contribution from lattice vibrations. We have calibrated our results against known exact solutions for limiting cases with very high accuracy. Having understood this problem, we are now extending our Langevin code to higher dimensions.