Warm cascades and anomalous scaling in a diffusion model of turbulence

A phenomenological turbulence model in which the energy spectrum obeys a nonlinear diffusion equation is analyzed. The general steady state contains a nonlinear mixture of the constant-flux Kolmogorov and fluxless thermodynamic components. Such "warm cascade" solutions describe a bottleneck phenomenon of spectrum stagnation near the dissipative scale. Transient self-similar solutions describing a finite-time formation of steady cascades are analyzed and found to exhibit nontrivial scaling behavior.     

用连续法计算五维对流模型的定常解和周期解

利用连续算法(Continuation algorithm)对五维对流非线性动力系统的定常解和周期解进行了数值计算。在参数平面Ri-Re上计算出实分岔点曲线、极限点曲线、Hopf分岔点曲线,绘出了分岔图。在分岔图上的不同区域,存在性质不同的稳定解如定常吸引子、周期吸引子等。分析了定常解、周期解的分岔过程。计算结果很好地说明大气中由基本态到对流态再到波动态最后到湍流态的物理转换过程。 连续算法对研究非线性动力系统的分岔以及耗散结构是很有效的计算方法。     

Influence of a transport current on a domain wall in an antiferromagnetic metal

We consider the influence of an electric current on the position of a domain wall in an antiferromagnetic metal. We first microscopically derive an equation of motion for the Néel vector in the presence of current by performing, in the transport steady state, a linear-response calculation in the deviation from collinearity of the antiferromagnet. This equation of motion is then solved variationally for an antiferromagnetic domain wall. We find that, in the absence of dissipative or non-adiabatic coupling between magnetization and current, the current displaces the domain wall by a finite amount and that the domain wall is then intrinsically pinned by the exchange interactions. In the presence of dissipative or non-adiabatic current-to-domain-wall coupling, the domain wall velocity is proportional to the current and is no longer pinned.     

Variational Quantum Algorithms for the Steady States of Open Quantum Systems

The solutions of the problems related to open quantum systems have attracted considerable interest.We propose a variational quantum algorithm to find the steady state of open quantum systems.In this algorithm,we employ parameterized quantum circuits to prepare the purification of the steady state and define the cost function based on the Lindblad master equation,which can be efficiently evaluated with quantum circuits.We then optimize the parameters of the quantum circuit to find the steady state.Numerical simulations are performed on the one-dimensional transverse field Ising model with dissipative channels.The result shows that the fidelity between the optimal mixed state and the true steady state is over 99%.This algorithm is derived from the natural idea of expressing mixed states with purification and it provides a reference for the study of open quantum systems.     

Quantization of Hénon's map with dissipation II — numerical results

The master equation for a quantized version of Hénon's map with dissipation derived in a preceding paper is here solved numerically for the Wigner quasi-probability density, under conditions of period doubling and classical chaos both in the transient regime and in the dissipative steady state. Approximations of the quantum map by a classical stochastic process are also considered and compared with solutions incorporating non-classical quantum fluctuations.     

Stable vortex tori in the three-dimensional cubic-quintic Ginzburg-Landau equation

We demonstrate the existence of stable toroidal dissipative solitons with the inner phase field in the form of rotating spirals, corresponding to vorticity S=0, 1, and 2, in the complex Ginzburg-Landau equation with the cubic-quintic nonlinearity. The stable solitons easily self-trap from pulses with embedded vorticity. The stability is corroborated by accurate computation of growth rates for perturbation eigenmodes. The results provide the first example of stable vortex tori in a 3D dissipative medium, as well as the first example of higher-order tori (with S=2) in any nonlinear medium. It is found that all stable vortical solitons coexist in a large domain of the parameter space; in smaller regions, there coexist stable solitons with either S=0 and S=1, or S=1 and S=2.     

Cellular instabilities,sublimit structures and edge-flames in premixed counterflows

We examine twin premixed flames in a plane counterflow and uncover, in the parameter space, a hitherto unknown domain of cellular instability. This leads us to hypothesize that for small Lewis numbers a two-dimensional (2D) steady solution branch bifurcates from the one-dimensional (1D) solution branch at a neutral stability point located near the strain-induced quenching point. Solutions on this 2D branch are constructed indirectly by solving an initial-value problem in the edge-flame context defined by the multiple-valued bistable 1D solution. Three kinds of solution are found: a periodic array of flame-strings, a single isolated flame-string and a pair of interacting flame-strings. These structures can exist for values of strain greater than the 1D quenching value, corresponding to sublimit solutions.     

Continuous symmetry reduction and return maps for high-dimensional flows

We present two continuous symmetry reduction methods for reducing high-dimensional dissipative flows to local return maps. In the Hilbert polynomial basis approach, the equivariant dynamics is rewritten in terms of invariant coordinates. In the method of moving frames (or method of slices) the state space is sliced locally in such a way that each group orbit of symmetry-equivalent points is represented by a single point. In either approach, numerical computations can be performed in the original state space representation, and the solutions are then projected onto the symmetry-reduced state space. The two methods are illustrated by reduction of the complex Lorenz system, a five-dimensional dissipative flow with rotational symmetry. While the Hilbert polynomial basis approach appears unfeasible for high-dimensional flows, symmetry reduction by the method of moving frames offers hope.     

Existence of dual solutions and three-dimensional instability in helical pipe flow

Three-dimensional (3D) direct numerical simulations (DNS) of the viscous incompressible fluid flow through a helical pipe with circular cross section were performed. The flow is governed by three parameters: the Dean number (or the Reynolds number), curvature, and torsion. First, we obtained steady solutions by steady 3D calculations, where dual solutions were found, one was uniform in the pipe axial direction and the other varied very slowly, if torsion exceeded a critical value. Then, the instability of the steady solutions obtained was studied by unsteady 3D calculations. We obtained critical Reynolds numbers of steady to unsteady transition by observing the behaviors of the unsteady solutions. The present results of the critical Reynolds number nearly agreed with those by the 2D linear stability analysis (Yamamoto et al. [9]) except for the lowest critical Reynolds number region, where the present study gave the critical Reynolds number much less than that obtained by the 2D linear stability analysis. We found the vortical structures in the form of a standing wave slightly above the marginal instability state, which is a trigger of explosive 3D instability.     

A new type of global bifurcation in Hénon map

A simple branch of solution on a bifurcation diagram, which begins at static bifurcation and ends at boundary crisis (or interior crisis in a periodic window), is generally a period-doubling cascade. A domain of solution in parameter space, enclosed by curves of static bifurcation and that of boundary crisis (or the interior of a periodic window), is the trace of branches of solution. A P-n branch of solution refers to the one starting from a period-n (n≥1) solution, and the corresponding domain in parameter space is named the P-n domain of solution. Because of the co-existence of attractors, there may be several branches within one interval on a bifurcation diagram, and different domains of solution may overlap each other in some areas of the parameter space. A complex phenomenon, concerned both with the co-existence of attractors and the crises of chaotic attractors, was observed in the course of constructing domains of steady state solutions of the Hénon map in parameter space by numerical methods. A narrow domain of period-m solutions firstly co-exists with (lies on) a big period-n (m≥3n) domain. Then it enters the chaotic area of the big domain and becomes period-m windows. The co-existence of attractors disappears and is called the landing phenomenon. There is an interaction between the two domains in the course of landing: the chaotic area in the big domain is enlarged, and there is a crisis step near the landing area.     

Thermal Robustness of Entanglement in a Dissipative Two-Dimensional Spin System in an Inhomogeneous Magnetic Field

Recently new novel magnetic phases were shown to exist in the asymptotic steady states of spin systems coupled to dissipative environments at zero temperature. Tuning the different system parameters led to quantum phase transitions among those states. We study, here, a finite two-dimensional Heisenberg triangular spin lattice coupled to a dissipative Markovian Lindblad environment at finite temperature. We show how applying an inhomogeneous magnetic field to the system at different degrees of anisotropy may significantly affect the spin states, and the entanglement properties and distribution among the spins in the asymptotic steady state of the system. In particular, applying an inhomogeneous field with an inward (growing) gradient toward the central spin is found to considerably enhance the nearest neighbor entanglement and its robustness against the thermal dissipative decay effect in the completely anisotropic (Ising) system, whereas the beyond nearest neighbor ones vanish entirely. The spins of the system in this case reach different steady states depending on their positions in the lattice. However, the inhomogeneity of the field shows no effect on the entanglement in the completely isotropic (XXX) system, which vanishes asymptotically under any system configuration and the spins relax to a separable (disentangled) steady state with all the spins reaching a common spin state. Interestingly, applying the same field to a partially anisotropic (XYZ) system does not just enhance the nearest neighbor entanglements and their thermal robustness but all the long-range ones as well, while the spins relax asymptotically to very distinguished spin states, which is a sign of a critical behavior taking place at this combination of system anisotropy and field inhomogeneity.     

New,Spherical Solutions of Non-Relativistic,Dissipative Hydrodynamics

We present a new family of exact solutions of dissipative fireball hydrodynamics for arbitrary bulk and shear viscosities. The main property of these solutions is a spherically symmetric, Hubble flow field. The motivation of this paper is mostly academic: we apply non-relativistic kinematics for simplicity and clarity. In this limiting case, the theory is particularly clear: the non-relativistic Navier–Stokes equations describe the dissipation in a well-understood manner. From the asymptotic analysis of our new exact solutions of dissipative fireball hydrodynamics, we can draw a surprising conclusion: this new class of exact solutions of non-relativistic dissipative hydrodynamics is asymptotically perfect.     

Covariant Non-Equilibrium Transport Theory Solutions for RHIC

New numerical solutions of 3+1D covariant kinetic theory are reported for nuclear collisions in the energy domain E_cm200 AGeV. They were obtained using the MPC 0.1.2 parton transport code employing high parton subdivision to retain Lorentz covariance. The solutions are compared to those of relativistic hydrodynamics employing Cooper–Frye isotherm freeze-out. The transport solutions follow a different dynamical path than hydrodynamics due to large dissipative effects when pQCD scattering rates and HIJING initial conditions are assumed. The transport freeze-out four-volume is sensitive to the reaction rates. The final transverse momentum distributions are found to deviate by up to an order of magnitude from those of Cooper–Frye frozen hydrodynamics.     

Low-dimensional physics of ultracold gases with bound states and the sine-Gordon model

We present 2D steady concentration profiles of confined layers of off-critical polymer blends. The layer rests on a solid substrate and has a flat free surface due to very high surface tension. The profiles correspond to non-linear steady solutions of the Cahn-Hilliard equation in a rectangular domain. The free polymer-gas interface is considered to be sharp, while the internal interfaces are diffuse. We explore the rich solution structure (including laterally structured layers, stratified layers, checkerboard structures, oblique states and droplets) as a function of mean concentration.     

Dissipative Dynamics of Quantum Discord of Two Strongly Driven Qubits

The exact dynamics of quantum discord (QD) of two strongly driven qubits, which are initially prepared in the X-type quantum states and inserted in two independent dissipative cavities or in a common dissipative cavity, are studied. The results indicate that both in the two cases, the evolution of QD is independent of the initial cavity state. For the two independent dissipative cavities, it is found that the phenomenon of sudden transition between classical and quantum decoherence exists and the transition time can be greatly delayed by suitably choosing the initial state parameter of the two qubits, the cavity mode-driving field detunning and the decay rate of the cavity. For the common dissipative cavity, it is shown that for some initial states of the two qubits, the QD can increase for a finite time at first, and then it decreases to a steady value, while for some other initial states, the QD can increase monotonously or with oscillation till a stable value is reached. Moreover, the creation of QD for the two qubits in a common cavity is discussed.     

Slow cracklike dynamics at the onset of frictional sliding

We propose a friction model which incorporates interfacial elasticity and whose steady state sliding relation is characterized by a generic nonmonotonic behavior, including both velocity weakening and strengthening branches. In 1D and upon the application of sideway loading, we demonstrate the existence of transient cracklike fronts whose velocity is independent of sound speed, which we propose to be analogous to the recently discovered slow interfacial rupture fronts. Most importantly, the properties of these transient inhomogeneously loaded fronts are determined by steady state front solutions at the minimum of the sliding friction law, implying the existence of a new velocity scale and a "forbidden gap" of rupture velocities. We highlight the role played by interfacial elasticity and supplement our analysis with 2D scaling arguments.     

Pulsating, creeping, and erupting solitons in dissipative systems

We present three novel pulsating solutions of the cubic-quintic complex Ginzburg-Landau equation. They describe some complicated pulsating behavior of solitons in dissipative systems. We study their main features and the regions of parameter space where they exist.     

Entropy balance in distributed reversible Gray-Scott model

A practical way to calculate the entropy change in the distributed media composed of reversible Gray-Scott model is demonstrated. The entropy change is given as the sum of the entropy production and the divergence of entropy flow. The divergence of entropy is calculated based on the chemical potential of steady state. It becomes evident that: (i) the entropy change for the emergence of dissipative structures in the open system can be positive or negative, (ii) most of the entropy produced inside the system is thrown out to the environment when dissipative structures are developing, (iii) the entropy production and the divergence of entropy flow balance completely, when the system shows static steady states, (iv) the entropy change behaves as if it is the time derivative of the entropy production. Prior to these calculations of entropy balance, the features of emergent patterns in the two-dimensional system are examined in terms of entropy production solely. The results imply that the entropy production can be an index for us to discriminate spatial patterns, but is not a global thermodynamic potential for the evolution of dissipative structures.     

Rotating,charged, and uniformly accelerating mass in general relativity

A new general class of solutions of the Einstein-Maxwell equations is presented. It depends on seven arbitrary parameters that group in a natural way into three complex parameters m + in, a + ib, e + ig, and the cosmological constant λ. They correspond to mass, NUT parameter, angular momentum per unit mass, acceleration, and electric and magnetic charge. The metric is in general stationary and axially symmetric. These solutions are of type D and contain as special cases all known solutions of type D belonging to this class. The known solutions are recovered by performing limiting transitions. An appropriate limit of our solutions describes an electromagnetic field in flat spacetime. We investigate the properties of that field. Its singular region corresponds in general to two circles moving with uniform acceleration in the positive and negative directions along the axis of symmetry. One can easily extend our solutions to the complex domain. Then it turns out that the metric can be written in a double Kerr-Schild form.