The application of the "inverse problem" method for constructing confining potentials that make N-soliton waveforms exact solutions in the Gross-Pitaevskii equation.

In this work, we discuss an application of the "inverse problem" method to find the external trapping potential, which has particular N trapped soliton-like solutions of the Gross-Pitaevskii equation (GPE) also known as the cubic nonlinear Schrödinger equation (NLSE). This inverse method assumes particular forms for the trapped soliton wave function, which then determines the (unique) external (confining) potential. The latter renders these assumed waveforms exact solutions of the GPE (NLSE) for both attractive (g<0) and repulsive (g>0) self-interactions. For both signs of g, we discuss the stability with respect to self-similar deformations and translations. For g<0, a critical mass Mc or equivalently the number of particles for instabilities to arise can often be found analytically. On the other hand, for the case with g>0 corresponding to repulsive self-interactions which is often discussed in the atomic physics realm of Bose-Einstein condensates, the bound solutions are found to be always stable. For g<0, we also determine the critical mass numerically by using linear stability or Bogoliubov-de Gennes analysis, and compare these results with our analytic estimates. Various analytic forms for the trapped N-soliton solutions in one, two, and three spatial dimensions are discussed, including sums of Gaussians or higher-order eigenfunctions of the harmonic oscillator Hamiltonian.

A nonlinear journey from structural phase transitions to quantum annealing.

Motivated by an exact mapping between equilibrium properties of a one-dimensional chain of quantum Ising spins in a transverse field (the transverse field Ising (TFI) model) and a two-dimensional classical array of particles in double-well potentials (the "ϕ4 model") with weak inter-chain coupling, we explore connections between the driven variants of the two systems. We argue that coupling between the fundamental topological solitary waves in the form of kinks between neighboring chains in the classical ϕ4 system is the analog of the competing effect of the transverse field on spin flips in the quantum TFI model. As an example application, we mimic simplified measurement protocols in a closed quantum model system by studying the classical ϕ4 model subjected to periodic perturbations. This reveals memory/loss of memory and coherence/decoherence regimes, whose quantum analogs are essential in annealing phenomena. In particular, we examine regimes where the topological excitations control the thermal equilibration following perturbations. This paves the way for further explorations of the analogy between lower-dimensional linear quantum and higher-dimensional classical nonlinear systems.

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.

Kink-Kink and Kink-Antikink Interactions with Long-Range Tails.

In this Letter, we address the long-range interaction between kinks and antikinks, as well as kinks and kinks, in φ^{2n+4} field theories for n>1. The kink-antikink interaction is generically attractive, while the kink-kink interaction is generically repulsive. We find that the force of interaction decays with the 2n/(n-1)th power of their separation, and we identify the general prefactor for arbitrary n. Importantly, we test the resulting mathematical prediction with detailed numerical simulations of the dynamic field equation, and obtain good agreement between theory and numerics for the cases of n=2 (φ^{8} model), n=3 (φ^{10} model), and n=4 (φ^{12} model).

Stability of Solitary Waves and Vortices in a 2D Nonlinear Dirac Model.

We explore a prototypical two-dimensional massive model of the nonlinear Dirac type and examine its solitary wave and vortex solutions. In addition to identifying the stationary states, we provide a systematic spectral stability analysis, illustrating the potential of spinor solutions to be neutrally stable in a wide parametric interval of frequencies. Solutions of higher vorticity are generically unstable and split into lower charge vortices in a way that preserves the total vorticity. These conclusions are found not to be restricted to the case of cubic two-dimensional nonlinearities but are found to be extended to the case of quintic nonlinearity, as well as to that of three spatial dimensions. Our results also reveal nontrivial differences with respect to the better understood nonrelativistic analogue of the model, namely the nonlinear Schrödinger equation.

PT-symmetric phase in kagome-based photonic lattices.

The 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 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 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.

Jarzynski Equality in PT-Symmetric Quantum Mechanics.

We show that the quantum Jarzynski equality generalizes to PT-symmetric quantum mechanics with unbroken PT symmetry. In the regime of broken PT symmetry, the Jarzynski equality does not hold as also the 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 PT symmetry is thermodynamically inhibited as the irreversible work diverges at the critical point.

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.

Quantum fluctuations drive nonmonotonic correlations in a qubit lattice.

Fluctuations may induce the degradation of order by overcoming ordering interactions, consequently leading to an increase of entropy. This is particularly evident in magnetic systems characterized by nontrivial, constrained disorder, where thermal or quantum fluctuations can yield counterintuitive forms of ordering. Using the proven efficiency of quantum annealers as programmable spin system simulators, we present a study based on entropy postulates and experiments on a platform of programmable superconducting qubits to show that a low level of uncertainty can promote ordering in a system impacted by both thermal and quantum fluctuations. A set of experiments is proposed on a lattice of interacting qubits arranged in a triangular geometry with precisely controlled disorder, effective temperature, and quantum fluctuations. Our results demonstrate the creation of ordered ferrimagnetic and layered anisotropic disordered phases, displaying characteristics akin to the elegant order-by-disorder phenomenon. Extensive experimental evidence is provided for the role of quantum fluctuations in lowering the total energy of the system by increasing entropy and defect clustering. Our thorough and comprehensive application of an intentionally introduced noise on a quantum platform provides insight into the dynamics of defects and fluctuations in quantum devices, which may help to reduce the cost associated with quantum processing.

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.

Uniform-density Bose-Einstein condensates of the Gross-Pitaevskii equation found by solving the inverse problem for the confining potential.

In this work, we study the existence and stability of constant density (flat-top) solutions to the Gross-Pitaevskii equation (GPE) in confining potentials. These are constructed by using the "inverse problem" approach which corresponds to the identification of confining potentials that make flat-top waveforms exact solutions to the GPE. In the one-dimensional case, the exact solution is the sum of stationary kink and antikink solutions, and in the overlapping region, the density is constant. In higher spatial dimensions, the exact solutions are generalizations of this wave function. In the absence of self-interactions, the confining potential is similar to a smoothed-out finite square well with minima also at the edges. When self-interactions are added, terms proportional to ±gψ^{*}ψ and ±gM with M representing the mass or number of particles in Bose-Einstein condensates get added to the confining potential and total energy, respectively. In the realm of stability analysis, we find (linearly) stable solutions in the case with repulsive self-interactions which also are stable to self-similar deformations. For attractive interactions, however, the minima at the edges of the potential get deeper and a barrier in the center forms as we increase the norm. This leads to instabilities at a critical value of M. Comparing the stability criteria from Derrick's theorem with Bogoliubov-de Gennes (BdG) analysis stability results, we find that both predict stability for repulsive self-interactions and instability at a critical mass M for attractive interactions. However, the numerical analysis gives a much lower critical mass. This is due to the emergence of symmetry-breaking instabilities that were detected by the BdG analysis and violate the symmetry x→-x assumed by Derrick's theorem.

Cross-polarized excitons in carbon nanotubes.

Polarization of low-lying excitonic bands in finite-size semiconducting single-walled carbon nanotubes (SWNTs) is studied by using quantum-chemical methodologies. Our calculations elucidate properties of cross-polarized excitons, which lead to the transverse optical absorption of nanotubes and presumably couple to intermediate-frequency modes recently observed in resonance Raman excitation spectroscopy. We identify up to 12 distinct excitonic transitions below the second fundamental band associated with the E(22) van Hove singularity. Calculations for several chiral SWNTs distinguish the optically active "bright" excitonic band polarized parallel to the tube axis and several optically "weak" cross-polarized excitons. The rest are optically (near) forbidden "dark" transitions. An analysis of the transition density matrices related to excitonic bands provides detailed information about delocalization of excitonic wavefunction along the tube. Utilization of the natural helical coordinate system accounting for the tube chirality allows one to disentangle longitudinal and circumferential components. The distribution of the transition density matrix along a tube axis is similar for all excitons. However, four parallel-polarized excitons associated with the E(11) transition are more localized along the circumference of a tube, compared with others related to the E(12) and E(21) cross-polarized transitions. Calculated splitting between optically active parallel- and cross-polarized transitions increases with tube diameter, which compares well with experimental spectroscopic data.

Thermalization in the one-dimensional Salerno model lattice.

The Salerno model constitutes an intriguing interpolation between the integrable Ablowitz-Ladik (AL) model and the more standard (nonintegrable) discrete nonlinear Schrödinger (DNLS) one. The competition of local on-site nonlinearity and nonlinear dispersion governs the thermalization of this model. Here, we investigate the statistical mechanics of the Salerno one-dimensional lattice model in the nonintegrable case and illustrate the thermalization in the Gibbs regime. As the parameter interpolating between the two limits (from DNLS toward AL) is varied, the region in the space of initial energy and norm densities leading to thermalization expands. The thermalization in the non-Gibbs regime heavily depends on the finite system size; we explore this feature via direct numerical computations for different parametric regimes.

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 affect conductance, and the opening of local regions or "bubbles" of base pairs in strands of DNA that are relevant for transcription and denaturation. 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.

Superconductivity in graphene induced by the rotated layer.

Recent discoveries in graphene bilayers have revealed that when one of the layers is rotated by a specific angle, superconductivity emerges. We provide an explanation for this phenomenon. We find that due to the layer rotations, the spinors are modified in such way that a repulsive interaction becomes attractive in certain directions. We also find that due to rotations the nodal points become angle dependent. The spinor in layer $ i=2 $ depends on the twisting angle in contrast to the spinor in layer $i=1$. As a result, the physics in the two layers depends on the twist and is identified with a twisted phase. In order to observe the twist we use an interaction term which changes sign. The change from a repulsive interaction to an attractive one gives rise to a one dimensional charge-density-wave. Due to tunneling between the two layers, the proximity of layer $i=1$ induces superconductivity in the charge-density-wave phase in layer $i=2$. This result is obtained by following a sequence of steps: when layer $2$ is rotated by an angle $\theta$, this rotation is equivalent to a rotation of an angle $-\theta$ of the linear momentum. Due to the discrete lattice, in layer $1$ the Fourier transform conserves the linear momentum $modulo$ the hexagonal reciprocal lattice vector. In layer $2$, due to the rotation, the linear momentum is conserved $modulo$ the {\it Moir\'e} reciprocal lattice vector. Periodicity is achieved at the $magic $ angles obtained from the condition of commensuration of the two lattices.

Spin Transport Based on Exchange Coupling in Doped Organic Polymers.

We develop a spin diffusion theory based on the exchange mechanism among polarons to understand the organic pure spin current. It is demonstrated that the exchange coupling is strong enough to induce spin transport within the organic layer with impurity concentrations higher than 1018 cm-3. By calculating the inverse spin Hall voltage in an organic spin device, we predict that the voltage depends nonmonotonically on the impurity concentration of the organic material. By tuning the doping concentration, one can achieve a maximum inverse spin Hall voltage. Our results not only explain some recent experimental data but also inspire further experimental investigation on pure spin current in organic devices with variable impurity doping.

Surface-wave-assisted nonreciprocity in spatio-temporally modulated metasurfaces.

Emerging photonic functionalities are mostly governed by the fundamental principle of Lorentz reciprocity. Lifting the constraints imposed by this principle could circumvent deleterious effects that limit the performance of photonic systems. Most efforts to date have been limited to waveguide platforms. Here, we propose and experimentally demonstrate a spatio-temporally modulated metasurface capable of complete violation of Lorentz reciprocity by reflecting an incident beam into far-field radiation in forward scattering, but into near-field surface waves in reverse scattering. These observations are shown both in nonreciprocal beam steering and nonreciprocal focusing. We also demonstrate nonreciprocal behavior of propagative-only waves in the frequency- and momentum-domains, and simultaneously in both. We develop a generalized Bloch-Floquet theory which offers physical insights into Lorentz nonreciprocity for arbitrary spatial phase gradients, and its predictions are in excellent agreement with experiments. Our work opens exciting opportunities in applications where free-space nonreciprocal wave propagation is desired.

Influence of the number of orientational domains on avalanche criticality in ferroelastic transitions.

We study the relationship between avalanche criticality and the number of orientational domains in ferroelastic transitions. To this end, we use a general Ginzburg-Landau model appropriate for displacive transitions of the square lattice. The model includes disorder as a quenched distribution of local transition temperatures. We focus on the square-to-rectangle and the square-to-oblique ferroelastic transitions, which have two and four orientational domains, respectively, which in turn determine the corresponding degeneracy of the ground state of the system. The phase transitions are driven by temperature under the assumption of a strict athermal behavior. That is, we assume that thermal fluctuations do not play any role. Numerical results are obtained using a purely relaxational dynamics, and it is shown that both the square-to-rectangle and the square-to-oblique transitions occur intermittently in the form of avalanches. Avalanche sizes and avalanche energies are found to display power-law distributions, which corroborates avalanche criticality. We compare and contrast the dependence of avalanche criticality on the number of orientational domains of the low-symmetry phase. It is found that the critical exponents depend on that number, in agreement with recent experimental results.

Speed-of-light pulses in a massless nonlinear Dirac equation.

In this work, we explore a massless nonlinear Dirac equation, i.e., a nonlinear Weyl equation. We study the dynamics of its pulse solutions and find that a localized one-hump initial condition splits into a localized two-hump pulse, while an associated phase structure emerges in suitable components of the spinor field. For times larger than a transient time t_{s} this pulse moves with the speed of light, effectively featuring linear wave dynamics and maintaining its shape (both in two and three dimensions). We show that for the considered nonlinearity, this pulse represents an exact solution of the nonlinear equation. Finally, we briefly comment on the generalization of the results to a broader class of nonlinearities.

Spin-current rectification in an organic magnetic/nonmagnetic device.

We propose a spin diode based on an organic magnetic co-oligomer or a magnetic/nonmagnetic heterojunction structure. The current and its spin polarization in the device are calculated with the spin-dependent Landauer-Büttiker formula. It is found that, by reversing the applied bias, the charge current and the spin current (SC) may be rectified at the same time or separately. A normal charge-current rectification usually takes place if the spatial electric structure is asymmetric. While a spin-current rectification may appear in two forms or their combination: one is that the spin-polarized orientation keeps unchanged but the magnitude of the SC is asymmetric with the bias; another is that only the spin orientation of the SC flips when the bias is reversed. By designing a suitable organic spin device, either of the two kinds of spin-current rectifications is obtained in our calculations. Finally, the effects of the properties of the organic interlayer and the structural asymmetry on the rectification are discussed.