Nearest-neighbor Interactions Research Articles

Highly controllable stabilization and switching of multiple colliding soliton sequences with generic Ginzburg–Landau gain–loss

We investigate propagation of J soliton sequences in a nonlinear optical waveguide array with generic weak Ginzburg–Landau (GL) gain–loss and nearest-neighbor (NN) interaction. The propagation is described by a system of J perturbed coupled nonlinear Schrödinger (NLS) equations. The NN interaction property leads to the elimination of collisional three-pulse interaction effects, which prevented the observation of stable multisequence soliton propagation with J>2 sequences in the presence of generic GL gain–loss in all previous studies. We show that the dynamics of soliton amplitudes can be described by a generalized J-dimensional Lotka–Volterra (LV) model. Stability and bifurcation analysis for the equilibrium points of the LV model, which is augmented by an application of the Lyapunov function method, is used to develop setups that lead to robust and scalable transmission stabilization and switching for a general J value. The predictions of the LV model are confirmed by extensive numerical simulations with the perturbed coupled-NLS model with J=3, 4, and 5 soliton sequences. Furthermore, soliton stability and the agreement between the LV model’s predictions and the simulations are not reduced with increasing value of J. Therefore, our study provides the first demonstration of robust control of multiple colliding sequences of NLS solitons in the presence of generic weak GL gain–loss with an arbitrary number of sequences. Due to the robustness and scalability of the results, they can have important applications in stabilization and switching of broadband soliton-based optical waveguide transmission.

Coupled Electronic and Magnonic Topological States in Two-Dimensional Ferromagnets.

Magnetic materials offer a fertile playground for fundamental physics discovery, with not only electronic but also magnonic topological states intensively explored. However, one natural material with both electronic and magnonic nontrivial topologies is still unknown. Here, we demonstrate the coexistence of first-order topological magnon insulators (TMIs) and electronic second-order topological insulators (SOTIs) in 2D honeycomb ferromagnets, giving rise to the nontrivial corner states being connected by the charge-free magnonic edge states. We show that, with symmetry, the phase factor ± ϕ caused by the next nearest-neighbor Dzyaloshinskii-Moriya interaction breaks the pseudo-spin time-reversal symmetry , which leads to the split of magnon bands, i.e., the emergence of TMIs with a nonzero Chern number of , in experimentally feasible candidates of MoI3, CrSiTe3, and CrGeTe3 monolayers. Moreover, protected by the symmetry, the electronic SOTIs characterized by nontrivial corner states are obtained, bridging the topological aspect of fermions and bosons with a high possibility of innovative applications in spintronics devices.

Magnetic-induced chiral dynamics in an extended two-leg bosonic ladder

The realization and detection of chiral physics with ultracold atomic gases provide a unique path for the exploration of topological phases. Here, we show that the interplay of magnetic field and interacting particles in an extended two-leg ladder leads to rich chiral Bloch dynamics. Considering both the on-site contact interaction and nearest-neighbor interactions, the ground state and Bloch dynamics of the system are studied analytically and numerically. When the system is in the ground state, the threshold and phase diagram for the transition between zero-momentum state and plane-wave state are analytically obtained, showing the nearest-neighbor interactions along the legs and rungs have opposite impact on the ground state transition, providing new opportunity to manipulate the ground state transition. When the ladder is perturbated under an external linear force, chiral dephasing of Bloch oscillations (BOs), i.e. symmetry breaking damped BOs (the damping rate of BOs on one leg is larger than the other), are observed. This chirality is absent for vanishing the magnetic field and atomic interaction. Particularly, the chirality of damped BOs is inversed when the magnetic field (chiral current) is inversed. In addition, the damping of BOs induced by different types of atomic interactions is different, and the strength and damping rate of BOs initialized in different ground states are distinct, offering dynamic ways to detect the different ground states. Furthermore, the persistent chiral Bloch oscillations observed in single particle case is predicted analytically, which is a crucial requirement for observation and application of chiral BOs in nonlinear regime. Our results provide an interesting path towards the exploration of topological atomic superfluids.

Classical spin-liquid behavior in interactionally distorted antiferromagnetic spin systems on the kagome lattice

The influence of an interaction distortion of the antiferromagnetic system on the kagome lattice on its magnetic and thermodynamic properties is investigated in the framework of the antiferromagnetic J1−J2−J3 spin-1/2 Ising system with three different nearest-neighbor antiferromagnetic interactions on the star kagome-like recursive lattice. The exact solution of the model is found and the phase diagram is determined and analyzed. It is shown that, depending on the model parameters, three antiferromagnetic phases can be identified, which are separated from the paramagnetic phase by the second order phase transitions. The magnetic and entropy properties of all ground states of the model are determined. Three different ground states that arise from the paramagnetic phase in the zero-temperature limit are identified. Due to the frustration, each of these ground states is highly macroscopically degenerated with different values of the residual entropy. In addition, since the paramagnetic phase in the vicinity of these ground states also exhibits high degeneracy with the presence of strong correlations, related to the proximity of the second-order phase transitions, these ground states together with their immediate paramagnetic surroundings can be considered as regions with the classical spin-liquid behavior. It is also shown that the typical anomalous (Schottky-type) behavior of the specific heat with the existence of a two-peak (in general multi-peak) structure in its temperature dependence in the vicinity of the highly macroscopically degenerated ground states can be suppressed by the presence of the second-order phase transitions.

Nonequilibrium generation of charge defects in kagome spin ice under slow cooling.

Kagome spin ice is one of the canonical examples of highly frustrated magnets. The effective magnetic degrees of freedom in kagome spin ice are Ising spins residing on a two-dimensional network of corner-sharing triangles. Due to strong geometrical frustration, nearest-neighbor antiferromagnetic interactions on the kagome lattice give rise to a macroscopic number of degenerate classical ground states characterized by ice rules. Elementary excitations at low temperatures are defect-triangles that violate the ice rules and carry an additional net magnetic charge relative to the background. We perform large-scale Glauber dynamics simulations to study the nonequilibrium dynamics of kagome ice under slow cooling. We show that the density of residual charge defects exhibits a power-law dependence on the quench rate for the class of algebraic cooling protocols. The numerical results are well captured by the rate equationfor the charge defects based on the reaction kinetics theory. As the relaxation time of the kagome ice phase remains finite, there is no dynamical freezing as in the Kibble-Zurek scenario. Instead, we show that the power-law behavior originates from a thermal excitation that decays algebraically with time at the late stage of the cooling schedule. Similarities and differences in quench dynamics of other spin ice systems are also discussed.

Multipolar Analysis in Symmetrical Meta-Atoms Sustaining Fano Resonances

We present an optical metasurface with symmetrical individual elements sustaining Fano resonances with high Q-factors. This study combines plane-wave illumination and modal analysis to investigate the resonant behavior that results in a suppression of the forward scattering, and we investigate the role of the lattice constant on the excited multipoles and on the spectral position and Q-factor of the Fano resonances, revealing the nonlocal nature of the resonances. The results show that the intrinsic losses play a crucial role in modulating the resonance amplitude in specific conditions and that the optical behavior of the device is extremely sensitive to the pitch of the metasurface. The findings highlight the importance of near-neighbor interactions to achieve high Q resonances and offer an important tool for the design of spectrally tunable metasurfaces using simple geometries.

Frustrated Extended Bose-Hubbard Model and Deconfined Quantum Critical Points with Optical Lattices at the Antimagic Wavelength.

The study of geometrically frustrated many-body quantum systems is of central importance to uncover novel quantum mechanical effects. We design a scheme where ultracold bosons trapped in a one-dimensional state-dependent optical lattice are modeled by a frustrated Bose-Hubbard Hamiltonian. A derivation of the Hamiltonian parameters based on Cesium atoms, further show large tunability of contact and nearest-neighbor interactions. For pure contact repulsion, we discover the presence of two phases peculiar to frustrated quantum magnets: the bond-order-wave insulator with broken inversion symmetry and a chiral superfluid. When the nearest-neighbor repulsion becomes sizable, a further density-wave insulator with broken translational symmetry can appear. We show that the phase transition between the two spontaneously symmetry-broken phases is continuous, thus representing a one-dimensional deconfined quantum critical point not captured by the Landau-Ginzburg-Wilson symmetry-breaking paradigm. Our results provide a solid ground to unveil the novel quantum physics induced by the interplay of nonlocal interactions, geometrical frustration, and quantum fluctuations.

Magnetic supersolid phases of two-dimensional extended Bose–Hubbard model with spin–orbit coupling

We investigate the quantum phases and phase transitions for spin–orbit coupled two-species bosons with nearest-neighbor (NN) interaction in a two-dimensional square lattice using inhomogeneous dynamical Guztwiller mean-field method. Under the effect of spin–orbit coupling and NN interaction, we uncover a rich variety of different magnetic supersolid (SS) phases. In the presence of intraspecies NN interaction, the phase diagram exhibits the phase-twisted double-checkerboard SS (PT-DCSS) and phase-striped double-checkerboard SS (PS-DCSS) phases. For both intra- and interspecies NN interactions, apart from the phase-twisted lattice SS (PT-LSS) and phase-striped lattice SS (PS-LSS) phases, some nontrivial SS phases with interesting properties occur. More importantly, we find that the emergences of these nontrivial SS phases are dependent of the interspecies on-site interaction. To further characterize the SS phases, we also discuss the spin-dependent momentum distributions and magnetic textures. The magnetic textures, such as antiferromagnetic, spiral and stripe orders are shown. Finally, we give the fully analytical insights into the numerical results.

Ordering kinetics of the two-dimensional voter model with long-range interactions.

We study analytically the ordering kinetics of the two-dimensional long-range voter model on a two-dimensional lattice, where agents on each vertex take the opinion of others at distance r with probability P(r)∝r^{-α}. The model is characterized by different regimes, as α is varied. For α>4, the behavior is similar to that of the nearest-neighbor model, with the formation of ordered domains of a typical size growing as L(t)∝sqrt[t], until consensus is reached in a time of the order of NlnN, with N being the number of agents. Dynamical scaling is violated due to an excess of interfacial sites whose density decays as slowly as ρ(t)∝1/lnt. Sizable finite-time corrections are also present, which are absent in the case of nearest-neighbor interactions. For 0<α≤4, standard scaling is reinstated and the correlation length increases algebraically as L(t)∝t^{1/z}, with 1/z=2/α for 3<α<4 and 1/z=2/3 for 0<α<3. In addition, for α≤3,L(t) depends on N at any time t>0. Such coarsening, however, only leads the system to a partially ordered metastable state where correlations decay algebraically with distance, and whose lifetime diverges in the N→∞ limit. In finite systems, consensus is reached in a time of the order of N for any α<4.

Double-well Bose-Hubbard model with nearest-neighbor and cavity-mediated long-range interactions

Double-well Bose-Hubbard model with nearest-neighbor and cavity-mediated long-range interactions

Topological Phase Diagram of an Interacting Kitaev Chain: Mean Field versus DMRG Study

In this work, we study the topological phase transitions of a Kitaev chain generalized by the addition of nearest-neighbor Coulomb interaction. We show the presence of a robust topological phase as a function of the interaction strength and of the on-site energy with associated non-zero energy Majorana states localized at the chain edges. We provide an effective mean-field model that allows for the self-consistent computation of the mean value of the local particle number operator, and we also perform Density Matrix Renormalization Group numerical simulations based on a tensor network approach. We find that the two methods show a good agreement in reporting the phase transition between trivial and topological superconductivity. Temperature robustness within a physically relevant threshold has also been demonstrated. These findings shed light on an entire class of topological interacting one-dimensional systems in which the effects of residual Coulomb interactions play a relevant role.

Entanglement-enabled symmetry-breaking orders

A spontaneous symmetry-breaking order is conventionally described by a tensor-product wavefunction of some few-body clusters; some standard examples include the simplest ferromagnets and valence bond solids. We discuss a type of symmetry-breaking orders, dubbed entanglement-enabled symmetry-breaking orders, which cannot be realized by any such tensor-product state. Given a symmetry breaking pattern, we propose a criterion to diagnose if the symmetry-breaking order is entanglement-enabled, by examining the compatibility between the symmetries and the tensor-product description. For concreteness, we present an infinite family of exactly solvable gapped models on one-dimensional lattices with nearest-neighbor interactions, whose ground states exhibit entanglement-enabled symmetry-breaking orders from a discrete symmetry breaking. In addition, these ground states have gapless edge modes protected by the unbroken symmetries. We also propose a construction to realize entanglement-enabled symmetry-breaking orders with spontaneously broken continuous symmetries. Under the unbroken symmetries, some of our examples can be viewed as symmetry-protected topological states that are beyond the conventional classifications.

A programmable hybrid digital chemical information processor based on the Belousov-Zhabotinsky reaction

The exponential growth of the power of modern digital computers is based upon the miniaturization of vast nanoscale arrays of electronic switches, but this will be eventually constrained by fabrication limits and power dissipation. Chemical processes have the potential to scale beyond these limits by performing computations through chemical reactions, yet the lack of well-defined programmability limits their scalability and performance. Here, we present a hybrid digitally programmable chemical array as a probabilistic computational machine that uses chemical oscillators using Belousov-Zhabotinsky reaction partitioned in interconnected cells as a computational substrate. This hybrid architecture performs efficient computation by distributing information between chemical and digital domains together with inbuilt error correction logic. The efficiency is gained by combining digital logic with probabilistic chemical logic based on nearest neighbour interactions and hysteresis effects. We demonstrated the computational capabilities of our hybrid processor by implementing one- and two-dimensional Chemical Cellular Automata demonstrating emergent dynamics of life-like entities called Chemits. Additionally, we demonstrate hybrid probabilistic logic as a viable logic for solving combinatorial optimization problems.

Investigation of Aluminum and Gold Flip-Chip Bonding for Quantum Device Integration

The majority of qubit chip integration is realized in a two-dimensional (2D) architecture. Whereas 3D architecture enables more advantages like efficient interconnect routing, allowing for more compact qubit coupling geometries, reducing form factor, and increased connectivity beyond nearest-neighbor interactions. Flipchip (FC) assembly has been demonstrated to enable 3D architecture connecting qubit chip, interposer, and readout in a sandwich-like structure. Moreover, 3D integration allows the fabrication of hosting chip circuitry without degrading the qubit performance [1–4]. Different material considerations have to be taken into account since qubit operational frequency is in the gigahertz range and operating temperature is in the mK range to avoid thermal excitation. Materials that possess superconducting characteristics like Indium (In), Titanium Nitride (TiN), Tantalum Nitride (TaN), and Niobium (Nb) have been discussed as potential interconnect between building blocks [3, 5]. The In bumping on Aluminum (Al) redistribution layers require under-bump metallization (UBM) layers thus introducing multiple fabrication steps before the bonding process. An alternative approach would be to use the existing Al surface to electrochemically grow Al bumps and bond the chip using thermosonic bonding (TSB) at below 150°C to form a homogeneous metal-metal interface [6, 7].

The collective dynamics of a stochastic Port-Hamiltonian self-driven agent model in one dimension

This paper studies the collective motion of self-driven agents in a one-dimensional space with periodic boundaries, using a stochastic Port-Hamiltonian system (PHS) with symmetric nearest-neighbor interactions and additive Brownian noise as an external input. In the case of a quadratic potential the PHS is an Ornstein-Uhlenbeck process for which we explicitly determine the distribution for any time t ≥ 0 and in the limit t → ∞. In particular, we characterize the collective motion by showing that the agents’ positions tend to build exactly one cluster. This is confirmed in simulations that show rapid and coordinated motion among agents, driven by noise, despite the absence of a preferred direction of motion in the model. Remarkably, the theoretical properties observed in the Ornstein-Uhlenbeck process also emerge in simulations of the nonlinear model incorporating a general interaction potential.

A quantization of interacting particle systems

Interacting particle systems studied in this paper are probabilistic cellular automata with nearest-neighbor interaction including the Domany-Kinzel model. A special case of the Domany-Kinzel model is directed percolation. We regard the interacting particle system as a Markov chain on a graph. Then we present a new quantization of the interacting particle system. After that, we introduce a zeta function of the quantized model and give its determinant expression. Moreover, we calculate the absolute zeta function of the quantized model for the Domany-Kinzel model.

The Effect of Polymer-Solvent Interaction on the Swelling of Polymer Matrix Tablets: A Magnetic Resonance Microscopy Study Complemented by Bond Fluctuation Model Simulations.

Polymer matrix tablets are an important drug-delivery system widely used for oral drug administration. Understanding the tablet hydration process, both experimentally and theoretically, is, thus, very important for the development of drug delivery systems that exhibit high drug loading capacity and controlled release potential. In this study, we used magnetic resonance microscopy (MRM) to nondestructively and dynamically analyze the water hydration process of xanthan-based tablets. The swelling process was characterized by well-resolved fronts of erosion, swelling, and penetration. The experimental results were complemented by numerical simulations of the polymer matrix hydration process. In the simulations, the polymer tablet matrix was modeled as an assembly of interacting chains with embedded drug particles, while its hydration process was mediated by interaction with solvent particles. The swelling dynamics were modeled within a Monte Carlo-based bond fluctuation model (BFM) that elegantly accounted for steric and nearest-neighbor interactions. This study provides an efficient experimental-theoretical approach for the study of polymer matrix swelling processes.

Persistence of weak ferromagnetism in antiferromagnetic systems on the body-centered octahedral lattice.

The possibility of existence of weak ferromagnetism in antiferromagnetic systems on the body-centered octahedral lattice is investigated in the framework of the corresponding spin-1/2J_{1}-J_{2} model in the recursive-lattice approach. The exact solution of the model is found and its phase diagram is determined. The magnetic and thermodynamic properties of all phases are studied, the nature of all phase transitions is established, and an equationthat determines the positions of all second-order phase transitions of the model is found. The magnetic and entropy properties of all ground states of the model are also determined. It is shown that the weak ferromagnetism predicted earlier in the pure spin-1/2 antiferromagnetic system on the octahedral lattice remains present even in the case of the model on the body-centered octahedral lattice with antiferromagnetic as well as ferromagnetic interactions between the central site of each elementary octahedron and each of its vertices. However, the presence of the interacting central site suppresses the weak ferromagnetism at very low temperatures, where the standard antiferromagnetic phase emerges. At the same time, the temperature region with the antiferromagnetic phase increases with simultaneous decreasing of the region, where the weak ferromagnetism can be observed, when the strength of the interaction between the central site of each elementary octahedron and each of its vertices increases. Moreover, the phenomenon of weak ferromagnetism disappears completely when this interaction becomes sufficiently stronger than the nearest-neighbor antiferromagnetic interaction between spin variables placed in vertices of each elementary octahedron of the lattice. Moreover, the possibility of the existence of the classical spin-liquid behavior in such magnetic systems is also discussed.

Digital quantum simulation of non-perturbative dynamics of open systems with orthogonal polynomials

Classical non-perturbative simulations of open quantum systems&amp;apos; dynamics face several scalability problems, namely, exponential scaling of the computational effort as a function of either the time length of the simulation or the size of the open system. In this work, we propose the use of the Time Evolving Density operator with Orthogonal Polynomials Algorithm (TEDOPA) on a quantum computer, which we term as Quantum TEDOPA (Q-TEDOPA), to simulate non-perturbative dynamics of open quantum systems linearly coupled to a bosonic environment (continuous phonon bath). By performing a change of basis of the Hamiltonian, the TEDOPA yields a chain of harmonic oscillators with only local nearest-neighbour interactions, making this algorithm suitable for implementation on quantum devices with limited qubit connectivity such as superconducting quantum processors. We analyse in detail the implementation of the TEDOPA on a quantum device and show that exponential scalings of computational resources can potentially be avoided for time-evolution simulations of the systems considered in this work. We applied the proposed method to the simulation of the exciton transport between two light-harvesting molecules in the regime of moderate coupling strength to a non-Markovian harmonic oscillator environment on an IBMQ device. Applications of the Q-TEDOPA span problems which can not be solved by perturbation techniques belonging to different areas, such as the dynamics of quantum biological systems and strongly correlated condensed matter systems.

Transition to period-3 synchronized state in coupled gauss maps.

We study coupled Gauss maps in one dimension with nearest-neighbor interactions. We observe transitions from spatiotemporal chaos to period-3 states in a coarse-grained sense and synchronized period-3 states. Synchronized fixed points are frequently observed in one dimension. However, synchronized periodic states are rare. The obvious reason is that it is very easy to create defects in one dimension. We characterize all transitions using the following order parameter. Let x∗ be the fixed point of the map. The values above (below) x∗ are classified as +1 (-1) spins. We expect all sites to return to the same band after three time steps for a coarse-grained periodic or three-period state. We define the flip rate F(t) as the fraction of sites i such that si(3t-3)≠si(t). It is zero in the coarse-grained periodic state. This state may or may not be synchronized. We observe three different transitions. (a) If different sites reach different bands, the transition is in the directed-percolation universality class. (b) If all sites reach the same band, we find an Ising-type transition. (c) A synchronized period-3 state where a new exponent is observed. We also study the finite-size scaling at critical points. The exponents obtained indicate that the synchronized period-3 transition is in a new universality class.