Thermostating by Deterministic Scattering: The Periodic Lorentz Gas

We present a novel mechanism for thermalizing a system of particles in equilibrium and nonequilibrium situations, based on specifically modeling energy transfer at the boundaries via a microscopic collision process. We apply our method to the periodic Lorentz gas, where a point particle moves diffusively through an ensemble of hard disks arranged on a triangular lattice. First, collision rules are defined for this system in thermal equilibrium. They determine the velocity of the moving particle such that the system is deterministic, time-reversible, and microcanonical. These collision rules can systematically be adapted to the case where one associates arbitrarily many degrees of freedom to the disk, which here acts as a boundary. Subsequently, the system is investigated in nonequilibrium situations by applying an external field. We show that in the limit where the disk is endowed by infinitely many degrees of freedom it acts as a thermal reservoir yielding a well-defined nonequilibrium steady state. The characteristic properties of this state, as obtained from computer simulations, are finally compared to those of the so-called Gaussian thermostated driven Lorentz gas.     

Smooth Dynamics and New Theoretical Ideas in Nonequilibrium Statistical Mechanics

This paper reviews various applications of the theory of smooth dynamical systems to conceptual problems of nonequilibrium statistical mecanics. We adopt a new point of view which has emerged progressively in recent years, and which takes seriously into account the chaotic character of the microscopic time evolution. The emphasis is on nonequilibrium steady states rather than the traditional approach to equilibrium point of view of Boltzmann. The nonequilibrium steady states, in presence of a Gaussian thermostat, are described by SRB measures. In terms of these one can prove the Gallavotti–Cohen fluctuation theorem. One can also prove a general linear response formula and study its consequences, which are not restricted to near-equilibrium situations. At equilibrium one recovers in particular the Onsager reciprocity relations. Under suitable conditions the nonequilibrium steady states satisfy the pairing theorem of Dettmann and Morriss. The results just mentioned hold so far only for classical systems; they do not involve large size, i.e., they hold without a thermodynamic limit.     

Surface-driven instability and enhanced relaxation in the dynamics of a nonequilibrium interface

We determine the stability of a nonequilibrium interface between two coexisting solid phases in the presence of a weak external field. Starting at the coarsegrained (Cahn-Hilliard) level, we use the method of matched asymptotics to derive the macroscopic interfacial dynamics. We then show that the external field leads to an instability due to flux along the interface, in contrast with the more common Mullins-Sekerka type instability, which involves fluxes normal to the interface. We also find that the external field produces an important modification of the Gibbs-Thomson relation. With these results, we perform the linear stability analysis for an approximately flat interface. If the field is tangent to the interface, the modification of the Gibbs-Thomson relation is important and the interface is stabilized. If the field is normal to the interface, the surface flux is important, and the effect can be stabilizing or destabilizing, but the orientational dependence is opposite what would be obtained if the Mullins-Sekerka instability dominates. Numerical simulations are performed to study the effect of the surface current and are in agreement with our analytical results.     

Phase transitions in a driven lattice gas at two temperatures

By suitably combining the uniformly driven lattice gas and the two-temperature kinetic Ising model, we obtain a generalized model that allows us to probe a variety of nonequilibrium phase transitions, including a type not previously observed. This new type of transition involves longitudinally ordered steady states, which are phase-segragated states with interface normalsparallel to the drive. Using computer simulations on a two-dimensional lattice gas, we map out the structure of the phase diagram, and the nature of the transitions, in the three-dimensional space of the drive and the two temperatures. While recovering anticipated results in most cases, we find one surprise, namely, that the transition from disorder to longitudinal order is continuous. Unless it turns out to be very weakly first order, this result is inconsistent with the expectation of field-theoretic renormalization group calculations.     

A Novel Measure Inspired by Lyapunov Exponents for the Characterization of Dynamics in State-Transition Networks

The combination of network sciences, nonlinear dynamics and time series analysis provides novel insights and analogies between the different approaches to complex systems. By combining the considerations behind the Lyapunov exponent of dynamical systems and the average entropy of transition probabilities for Markov chains, we introduce a network measure for characterizing the dynamics on state-transition networks with special focus on differentiating between chaotic and cyclic modes. One important property of this Lyapunov measure consists of its non-monotonous dependence on the cylicity of the dynamics. Motivated by providing proper use cases for studying the new measure, we also lay out a method for mapping time series to state transition networks by phase space coarse graining. Using both discrete time and continuous time dynamical systems the Lyapunov measure extracted from the corresponding state-transition networks exhibits similar behavior to that of the Lyapunov exponent. In addition, it demonstrates a strong sensitivity to boundary crisis suggesting applicability in predicting the collapse of chaos.     

Maximal Lyapunov exponent and rotation numbers for two coupled oscillators driven by real noise

Asymptotic expansions for the exponential growth rate, known as the Lyapunov exponent, and rotation numbers for two coupled oscillators driven by real noise are constructed. Such systems arise naturally in the investigation of the stability of steady-state motions of nonlinear dynamical systems and in parametrically excited linear mechanical systems. Almost-sure stability or instability of dynamical systems depends on the sign of the maximal Lyapunov exponent. Stability conditions are obtained under various assumptions on the infinitesimal generator associated with real noise provided that the natural frequencies are noncommensurable. The results presented here for the case of the infinitesimal generator having a simple zero eigenvalue agree with recent results obtained by stochastic averaging, where approximate ItÔ equations in amplitudes and phases are obtained in the sense of weak convergence.Dedicated to Thomas K. Caughey on the occasion of his 65th birthday.     

Computing the pressure for Axiom-A attractors by time series and large deviations for the Lyapunov exponent

For the Axiom-A attractors a relation is given between the topological pressure and the spectrum of the generalized Lyapunov exponents. As a consequence, a simple formula is found to compute the topological entropy of the attractor by means of a time series. The results are used to compute the large deviations for positive Lyapunov exponents.     

Upper Quantum Lyapunov Exponent and Anosov Relations for Quantum Systems Driven by a Classical Flow

We generalize the definition of quantum Anosov properties and the related Lyapunov exponents to the case of quantum systems driven by a classical flow, i.e. skew-product systems. We show that the skew Anosov properties can be interpreted as regular Anosov properties in an enlarged Hilbert space, in the framework of a generalized Floquet theory. This extension allows us to describe the hyperbolicity properties of almost-periodic quantum parametric oscillators and we show that their upper Lyapunov exponents are positive and equal to the Lyapunov exponent of the corresponding classical parametric oscillators. As second example, we show that the configurational quantum cat system satisfies quantum Anosov properties.     

Mixed Convection Heat and Mass Transfer in a Micropolar Fluid with Soret and Dufour Effects

A mathematical model for the steady, mixed convection heat and mass transfer along a semi-infinite vertical plate embedded in a micropolar fluid in the presence of Soret and Dufour effects is presented. The non-linear governing equations and their associated boundary conditions are initially cast into dimensionless forms using local similarity transformations. The resulting system of equations is then solved numerically using the Keller-box method. The numerical results are compared and found to be in good agreement with previously published results as special cases of the present investigation. The non-dimensional velocity, microrotation, temperature and concentration profiles are displayed graphically for different values of coupling number, Soret and Dufour numbers. In addition, the skin-friction coefficient, the Nusselt number and Sherwood number are shown in a tabular form.     

格子Boltzmann模型平衡分布边界条件和多棱柱排列渗流流场的数值模拟

介绍了适用于多种流场数值模拟的无滑动格子Boltzmann平衡分布边界条件，这一边界条件是以Bounce-Back方法为基础并满足质量、动量守恒的准则.数值计算结果表明平衡分布边界条件克服了Bounce-Back方法在边界上所产生的滑动速度误差效应.利用平衡分布边界条件数值模拟了由棱柱形充填粒子构成的微尺度渗流流场中的Darcy-Forchheimer方程，通过与Lee 和Yang的数值结果比较，该预测结果是足够可靠的. 关键词： 平衡分布边界条件 渗流介质 Darcy-Forchheimer阻力     

矩量法研究平面波导点缺陷引起的散射损耗

对平面波导工艺过程里出现的点缺陷,如气泡、灰尘颗粒等对平面波导内光场造成的散射提供了一种简单而又实用的计算方法。计算基于格林函数方法,并利用矩量法求解积分方程,采用分块矩阵优化算法,使得算法的计算时间直接与点缺陷大小相关,提高了效率。证明了一定大小的点缺陷在特定的波长处会产生很大的散射损耗,证明点缺陷的散射损耗依赖于缺陷大小和光学性质而与其在波导中的位置无直接关系。以二氧化硅平面波导为例,分析了对中心波长1．55μm的入射光场,分别存在气泡和灰尘颗粒两种点缺陷时,散射损耗随点缺陷大小的关系特性,指出工艺过程特别应该避免某些孔径的点缺陷。     

A Numerical Investigation of Nanocomposite of Copper and Titanium Dioxide in Water Based Fluid Influenced by Instigated Magnetic Region

Presence of external electrical field plays a vital role in heat transfer and fluid flow phenomena. Keeping this in view present article is a numerical investigation of stagnation point flow of water based nanoparticles suspended fluid under the influence of induced magnetic field. A detailed comparative analysis has been performed by considering Copper and Titanium dioxide nanoparticles. Utilization of similarity analysis leads to a simplified system of coupled nonlinear differential equations, which has been tackled numerically by means of shooting technique followed by Runge-Kutta of order 5. The solutions are computed correct up to 6 decimal places. Influence of pertinent parameters is examined for fluid flow, induced magnetic field, and temperature profile. One of the key findings includes that magnetic parameter plays a vital role in directing fluid flow and lowering temperature profile. Moreover, it is concluded that Cu-water based nanofluid high thermal conductivity contributes in enhancing heat transfer efficiently.     

Out‐of‐Equilibrium Fluctuation‐Dissipation Relations Verified by the Electrical and Thermoelectrical AC‐Conductances in a Quantum Dot

The electrical and heat currents flowing through a quantum dot are calculated in the presence of a time‐modulated gate voltage with the help of the out‐of‐equilibrium Green function technique. From the first harmonics of the currents, we extract the electrical and thermoelectrical trans‐admittances and ac‐conductances. Next, by a careful comparison of the ac‐conductances with the finite‐frequency electrical and mixed electrical‐heat noises, we establish the fluctuation‐dissipation relations linking these quantities, which are thus generalized out‐of‐equilibrium for a quantum system. It is shown that the electrical ac‐conductance associated to the displacement current is directly linked to the electrical noise summed over reservoirs, whereas the relation between the thermoelectrical ac‐conductance and the mixed noise contains an additional term proportional to the energy step that the electrons must overcome when traveling through the junction. A numerical study reveals however that a fluctuation‐dissipation relation involving a single reservoir applies for both electrical and thermoelectrical ac‐conductances when the frequency dominates over the other characteristic energies.     

Two-Body Scattering in(1+1) Dimensions by a Semi-relativistic Formalism and a Hulthén Interaction Potential

Scattering solutions of two-body Spinless Salpeter Equation(SSE) are investigated in the center of mass frame with a repulsive, symmetric Hulth′en potential in one spatial dimension. Transmission and reflection coefficients are calculated and discussed.     

多孔介质中流体流动及扩散的耦合格子Boltzmann模型

多孔介质中高Péclet数和大黏性比下混溶流体的流动和扩散广泛存在于二氧化碳驱油、化工生产等工业过程中.用数值方法对该问题进行研究时,关键在于如何正确描述高Péclet数和大黏性比下多孔介质内流体的行为.为此,提出了一种基于多松弛模型和格子动理模型的耦合格子Boltzmann模型.通过Chapman-Enskog分析,证明该模型能有效求解不可压Navier-Stokes方程和对流扩散方程.数值结果表明,该模型不仅具有二阶精度和良好的稳健性,而且对于高Péclet数和大黏性比的问题具有良好的数值稳定性,为模拟此类问题提供了有效工具.     

Rotation of magnetization in unison and langevin equations for a large spin

A single-domain ferromagnetic particle is represented as a large spin (model of rotation in unison) whose stochastic dynamics is derived from a spin-boson Hamiltonian. It is shown in the Markovian limit that thermal equilibrium exists provided that the fluctuation-dissipation theorem is supplemented by a symmetry constraint which for bilinear anisotropic and nonlinear (magnetoelastic) spin-bath coupling can only be satisfied in the underdamped limit. Only for bilinear isotropic coupling (Gilbert's theory) is it satisfied identically for arbitrary damping strength. Uniaxial and cubic symmetries are considered. For a model uniaxial crystal the thermal decay rate of M and the thermal enhancement of the macroscopic quantum tunneling rate are calculated for Gilbert and magnetoelastic dissipative couplings and compared. The effects of memory are discussed.     

Floquet study of ionization and high-order harmonic generation for a hydrogen atom in high-frequency laser fields

In this paper, we investigate the dynamics of a hydrogen atom in high-frequency (several atomic units) super strong (up to several tens of atomic units) laser fields within the high frequency Floquet theory framework. The ionization rate, ionization spectrum, angular distribution and high-order harmonic generation are all investigated. Our studies reveal the universal behavior of the total ionization rate, excess-photon ionization spectrum and angular distribution of the ionization rate in the stabilization regime, and achieve a deep insight into the dynamics of high-order harmonic generation in the stabilization regime. Received 6 June 2001 and Received in final form 31 August 2001     

Superfluidity and Anderson localisation for a weakly interacting Bose gas in a quasiperiodic potential

Using exact diagonalisation and Density Matrix Renormalisation group (DMRG) approach we analyse the transition to a localised state of a weakly interacting quasi-1D Bose gas subjected to a quasiperiodic potential. The analysis is performed by calculating the superfluid fraction, density profile, momentum distribution and visibility for different periodicities of the second lattice and in the presence (or not) of a weak repulsive interaction. It is shown that the transition is sharper towards the maximally incommensurate ratio between the two lattice periodicities, and shifted to higher values of the second lattice strength by weak repulsive interactions. We also relate our results to recent experiments.     

Generation of 3rd and 5th harmonics in a thin superconducting film by temperature oscillations and isothermal nonlinear current response

The generation of harmonics of the voltage response is considered when an AC current is applied through a superconducting film above T_c. It is shown that almost at all temperatures the mechanism of the temperature oscillations created by the AC current and the temperature dependence of the resistance dominates over the isothermal nonlinear electric conductivity. Only in a narrow critical region close to T_c the latter is essential for the generation of the harmonics. A detailed investigation of harmonics generation provides an accurate method for measuring the thermal boundary conductance between the film and the insulating substrate. The critical behaviour of the third harmonic will give a new method for the determination of the lifetime of metastable Cooper pairs above T_c. The comparison of the calculated fifth harmonics of the voltage with the experiment is proposed as an important test for the applicability of the employed theoretical models. Received 8 September 2001