A mixed spectral and finite difference model to study baroclinic annulus waves

A mixed spectral and finite difference model to study finite amplitude baroclinic waves in a differentially heated rotating annulus is presented. The model consists of the full Navier-Stokes equations and the heat equation. The field variables f = f(r, φ z; t) are decomposed into zonally averaged components f_o(r, z; t) and eddy components f(r, φ, z; t), the latter being periodic in f and represented in terms of Fourier series. The unknowns f_o(r, z; t) and f^{c, s}(r, z; t), which are Fourier amplitudes of f′(r, φ, z; t) are governed by two-dimensional primitive equations with the addition of source terms. These equations are solved semi-implicitly by the alternating direction implicit method on variable grids.A simplified model with two Fourier components which permits self-interaction of the chosen wave and the interaction of the wave and the mean fields had been used to repeat a computation done by G. P. Williams, who used a fully three-dimensional finite difference algorithm. We can reproduce almost all of Williams' results in 1/20 of the computing time with the present model. It only requires 1/30 the additional computer storage of Williams' finite difference model over the axisymmetric problem.The potential of the present model for investigation of multiwave interaction as well as the advantages and disadvantages of the two different approaches is discussed.     

Dynamical theory for photon and photoelectric pulse distributions in single molecule fluorescence

A single two-level molecule driven by CW-laser field and a photomultiplier tube (PMT) are considered as two parts of the united dynamical system connected with each other by photons of molecular fluorescence. Each PMT is characterized by a rate α of photo-effect and by a rate β of PMT recovery. A theory for the photon distribution function w_N(t) and for the photoelectric pulse distribution function f_n(t) for such a system is built up. If times 1/ α and 1/ β characterizing PMT are much shorter as compared to the average time interval 1/ k between two successively emitted photons of fluorescence, the photon and the photoelectron distribution functions coincide with each other, i.e. f_n(t) ≅ w_N(t). A relation between w_N(t) and f_n(t) is studied in detail for the case in which PMT works slower as compared to the rate k of photon emission, i.e. at 1/ α, 1/ β ≥ 1/ k.     

Kinetic Limit of N-Body Description of Wave–Particle Self-Consistent Interaction

A system of N particles $\xi ^N = x_1 ,\upsilon_1,...,x_N ,\upsilon _N )$ interacting self-consistently with one wave Z = A exp(iφ) is considered. Given initial data (Z ^(N)(0), ξ ^N (0)), it evolves according to Hamiltonian dynamics to (Z ^(N)(t), ξ ^N (t)). In the limit N → ∞, this generates a Vlasov-like kinetic equation for the distribution function f(x, v, t), abbreviated as f(t), coupled to the envelope equation for Z: initial data (Z ^(∞)(0), f(0)) evolve to (Z ^(∞)(t), f(t)). The solution (Z, f) exists and is unique for any initial data with finite energy. Moreover, for any time T>0, given a sequence of initial data with N particles distributed so that the particle distribution f ^N(0) → f(0) weakly and with Z ^(N)(0) → Z(0) as N → ∞, the states generated by the Hamiltonian dynamics at all times 0 ≤ t ≤ T are such that (Z ^(N)(t), f ^N(t)) converges weakly to (Z ^(∞)(t), f(t)).     

Finite and infinite systems of nonlinearly-coupled ordinary differential equations,the solutions of which feature remarkable Diophantine findings

We use previous results concerning the time evolution of the zeros x_n(t) of time-dependent polynomials p_N (z;t) or entire functions F(z;t) of the complex variable z, in order to identify lots of nonlinearly-coupled, finite or infinite, systems of Ordinary Differential Equations the solutions of which feature remarkable Diophantine properties.     

Asymptotic behavior of orbits of C 0-semigroups and solutions of linear and semilinear abstract differential equations

In this paper, we study the asymptotic behavior of solutions of semilinear abstract differential equations (*) u′(t) = Au(t) + t ⁿ f(t, u(t)), where A is the generator of a C ₀-semigroup (or group) T(·), f(·, x) ∈ A for each x ∈ X, A is the class of almost periodic, almost automorphic or Levitan almost periodic Banach space valued functions ϕ: ℝ → X and n ∈ {0, 1, 2, ...}. We investigate the linear case when T(·)x is almost periodic for each x ∈ X; and the semilinear case when T(·) is an asymptotically stable C ₀-semigroup, n = 0 and f(·, x) satisfies a Lipschitz condition. Also, in the linear case, we investigate (*) when ϕ belongs to a Stepanov class S ^p-A defined similarly to the case of S ^p-almost periodic functions. Under certain conditions, we show that the solutions of (*) belong to A _u:= A ∩ BUC(ℝ, X) if n = 0 and to t ⁿ A _u ⊕ w _n C ₀ (ℝ, X) if n ∈ ℕ, where w _n(t) = (1 + |t|)ⁿ. The results are new for the case n ∈ ℕ and extend many recent ones in the case n = 0. Dedicated to the memory of B. M. Levitan     

The use of unitarity bounds for a stable extrapolation of low energy data

The properties of the scattering amplitude allow to define a functionf(z) satisfying the following conditions:

1. f(z) is holomorphic in a simply connected domain , which can be mapped conformally onto the unit disk;
2. |Imf(z)| is bounded by some constantM in ;
3. |Ref(z)| is known not to exceed some constantm on a certain part Γ₁ of the boundary Γ of ;f(z) is continuously extensible onto Γ.

Using these properties, constraints are derived on the real part off(z) valid at any point . The result is used for performing a stable extrapolation of low energy pion-pion scattering data to any finite energy. We derive a bound on energy averaged values of the real part of the scattering amplitude. The bound depends onm, M, on the energy variabless and on the energy average intervals ₂?s ₁. Generalizations of the method are discussed.     

Correlation functions of the transverse Ising chain at the critical field for large temporal and spatial separations

We study the behavior of 〈σ₀^x(t)σ_n^x(0)〉 and 〈σ₀^y(t)σ_n^y(0)〉 for the transverse Ising chain at the critical magnetic field at T = 0. Explicit results are obtained for the three distinct regions where t → ∞ and n → ∞with $\text{0 ?}\text{n}\text{t}\text{<1}$, $\text{1 <}\text{n}\text{t}$, or $\text{t = n + n}{}^{\mathrm{<mtext>1</mtext><mtext>3</mtext>}}\text{(}\text{z}\text{2}\text{)}$ where z is fixed of order one. In this latter region the general Painlevé V solution is shown to reduce to a Painlevé II function. We use our results to discuss the general problem of long-time behavior of Toda equations with slowly decaying initial values.     

On the fragmentation functions of heavy quarks into hadrons

An explicit form for the charmed quark fragmentation function D^C_c(z) into hadrons has been obtained with the help of the “reciprocity relation” and the c-quark distribution function in charmed mesons (the function calculated in terms of the Kuti-Weisskopf model). D^C_c(z) turns out to peak mainly at z close to 1. The analysis of new data on muon pair production in neutrino reactions points to such a behaviour of the D^C_c(z) function. The obtained fragmentation function, contrary to those, used earlier, leads to a charmed particle (average) multiplicity in e⁺e^?-annihilation independent of energy.     

The Inviscid Burgers Equation with Initial Value of Poissonian Type

We study the discontinuities (shocks) of the solution to the Burgers equation in the limit of vanishing viscosity (the inviscid limit) when the initial value is the opposite of the standard Poisson process p. We show that this solution is only defined for t ε (0, 1). Let T ₀ = 0 and T _n , n≧1, be the successive jumps of p. We prove that for all M > 0 the inviscid limit is characterized on the region x ε (-∞, M], t ε (0, 1) by the increasing process $N(t) = \sup \{ n \in \mathbb{N} {\text{| }}M + nt > T_n \} $ and the random set I(x) = {n ε {0,..., N(t)}‖T _n -nt≦x<T _n+1 - nt}. The positions of shocks are given in a precise manner. We give the distribution of N(t) and also the distribution of its first jump. We also prove similar results when the initial value is u _μ(y, 0) = -μp(y/μ²) + μ^-1 max(y, 0), μ ε (0, 1).     

Spin autocorrelation of Ni in K2Mn0.975Ni0.025F4 studied by nuclear spin-lattice relaxation

From the nuclear spin-lattice relaxation of the out-of-layer ¹⁹F nuclei in magnetic fields perpendicular to the c-axis the low-frequency component of the autocorrelation function 〈S^z(t)S^z(O)〉 of Ni in ordered K₂Mn_0.975Ni_0.025F₄ is found to be substantially reduced relative to the Mn host. The experimental rates vs temperature are in accord with those for relaxation involving two spin excitations calculated with local Green's functions.     

Multiplicity distributions in the generalized Feynman gas model

Multiplicity distributions Ψ_n^(k) in the generalized Feynman gas model of order k (defined by saying that all integrated correlation functions f_n except f₁,…,f_k are zero) are derived and expressed in terms of Poisson distributions with different ”average multiplicities”, which are related to the integrated correlation functions. The relations between Ψ_n^(k) and Ψ_n^(j) for arbitrary positive integers k,j are found. An intuitive picture to gain a better feeling for these relations is developed.On the basis of our formulae we show that the experimentally observed multiplicity distributions (between 50 GeV/c and 303 GeV/c incoming momentum) can be well reproduced by those of a Feynman gas model of order two. Other applications of our formulae are suggested.     

The reactions N14(n, α)B11 and N14(n,t)C12 observed with a gridded ionization chamber

The cross sections for the reactions N¹⁴(n, α)B¹¹ and N¹⁴(n, t)C¹² have been measured in the neutron energy range 4.0 to 6.4 MeV and at 2.5 MeV. Mono-energetic neutrons were produced in the D(d, n) He³ reaction using a gas target. The (n, α) and (n, t) disintegrations were detected in a gridded ionization chamber filled with an argonnitrogen mixture. The response of the chamber under different operation conditions is described. The excitation functions, measured with a neutron energy resolution of 40 to 50 keV, are given for theα ₀ group from the N¹⁴(n,α)B¹¹ reaction over the entire neutron energy range and for theα ₁ group and the t₀ group from N¹⁴(n, t) C¹² for neutron energies above 4.3 and 5.6 MeV, respectively.     

Applications of exact solution for strongly interacting one-dimensional Bose-Fermi mixture: Low-temperature correlation functions, density profiles, and collective modes

We consider one-dimensional interacting Bose-Fermi mixture with equal masses of bosons and fermions, and with equal and repulsive interactions between Bose-Fermi and Bose-Bose particles. Such a system can be realized in current experiments with ultracold Bose-Fermi mixtures. We apply the Bethe ansatz technique to find the exact ground state energy at zero temperature for any value of interaction strength and density ratio between bosons and fermions. We use it to prove the absence of the demixing, contrary to prediction of a mean-field approximation. Combining exact solution with local density approximation in a harmonic trap, we calculate the density profiles and frequencies of collective modes in various limits. In the strongly interacting regime, we predict the appearance of low-lying collective oscillations which correspond to the counterflow of the two species. In the strongly interacting regime, we use exact wavefunction to calculate the single particle correlation functions for bosons and fermions at low temperatures under periodic boundary conditions. Fourier transform of the correlation function is a momentum distribution, which can be measured in time-of-flight experiments or using Bragg scattering. We derive an analytical formula, which allows to calculate correlation functions at all distances numerically for a polynomial time in the system size. We investigate numerically two strong singularities of the momentum distribution for fermions at k_f and k_f + 2k_b. We show, that in strongly interacting regime correlation functions change dramatically as temperature changes from 0 to a small temperature ∼E_f/γ ? E_f, where E_f = (π?n)²/(2m), n is the total density and γ = mg/(?²n) ? 1 is the Lieb-Liniger parameter. A strong change of the momentum distribution in a small range of temperatures can be used to perform a thermometry at very small temperatures.     

Reconstruction of Five-Dimensional Bounce Cosmological Models from Deceleration Factor

In this paper, we consider a class of five-dimensional Ricci-flat vacuum solutions, which contain two arbitrary functions μ(t) and ν(t). It is shown that μ(t) can be rewritten as a new arbitrary function f(z) in terms of redshift z and the f(z) can be determined by choosing particular deceleration parameters q(z) which gives early deceleration and late time acceleration. In this way, the 5D cosmological model can be reconstructed and the evolution of the universe can be determined. PACS: 04.50.+h, 98.80.-k     

The generalized exponential-integral V(x,y)=∫1∞ exp(−xt)ln(t+y)dtt and computer algorithms for y = 0, 1

The generalized exponential-integral function V(x, y) defined here includes as special cases the function E⁽²⁾₁(x) = V(x, 0) introduced by van de Hulst and functions M₀(x) = V(x, 1) and N₀(x) = V(x, -1) introduced by Kourganoff in connection with integrals of the form ∫ E_n)t)E_m(t±x), which play an important role in the theory of monochromatic radiative transfer. Series and asymptotic expressions are derived and, for the most important special cases, y = 0 and y = 1, Chebyshev expansions and rational approximations are obtained that permit the function to be evaluated to at least 10 sf on 0<x<∞ using 16 sf arithmetic.     

Long-range spin correlations in a honeycomb spin model with a magnetic field

We consider the spin-1/2 model on the honeycomb lattice [A. Kitaev, Ann. Phys. 321, 2 (2006)] in the presence of a weak magnetic field h _α ? J. Such a perturbation treated in the lowest nonvanishing order over h _α leads [K.S. Tikhonov, M.V. Feigel’man, and A.Yu. Kitaev, Phys. Rev. Lett. 106, 067203 (2011)] to a powerlaw decay of irreducible spin correlations 《s ^z (t, r)s ^z (0, 0)》 ∝ h _z ² f(t, r), where f(t, r) ∝ [max(t, Jr)]^–4. We have studied the effects of the next order of perturbation in h _z and found an additional term of the order h _z ⁴ in the correlation function 《s ^z (t, r)s ^z (0, 0)》 which scales as h _z ⁴ cosγ/r ³ at Jt? r, where γ is the polar angle in the 2D plane. We demonstrate that such a contribution can be understood as a result of a perturbation of the effective Majorana Hamiltonian by the weak imaginary vector potential A _x ∝ i h _z ² .     

On the derivation of the Boltzmann-Landau equation from the quantum mechanical hierarchy

The general solution of the equation of motion for the quantum mechanical distribu tion functionf ₂(r ₁ P ₁,r ₂ p ₂;t)in the two particle space is given by means of the Schrödinger scattering functions. A special initial condition leads to the usual Boltzmann equation plus density correction terms, which depend on the scattering matrixt(p′,p). In the long wavelength limit and in lowest order oft(p′,p) the Landau corrections to the simple Boltzmann streaming part are obtained.     

Computable upper bounds for the adiabatic approximation errors

For a given Hermitian Hamiltonian H(s)(s∈[0,1])with eigenvalues Ek(s)and the corresponding eigenstates|Ek(s)(1 k N),adiabatic evolution described by the dilated Hamiltonian HT(t):=H(t/T)(t∈[0,T])starting from any fixed eigenstate|En(0)is discussed in this paper.Under the gap-condition that|Ek(s)-En(s)|λ0 for all s∈[0,1]and all k n,computable upper bounds for the adiabatic approximation errors between the exact solution|ψT(t)and the adiabatic approximation solution|ψadi T(t)to the Schr¨odinger equation i|˙ψT(t)=HT(t)|ψT(t)with the initial condition|ψT(0)=|En(0)are given in terms of fidelity and distance,respectively.As an application,it is proved that when the total evolving time T goes to infinity,|ψT(t)-|ψadi T(t)converges uniformly to zero,which implies that|ψT(t)≈|ψadi T(t)for all t∈[0,T]provided that T is large enough.     

L ∞ estimates for the space-homogeneous Boltzmann equation

This paper studies the boundedness of solutionsf of the initial-value problem for the space-homogeneous Boltzmann equation for inverse kth power forces, whenk>5, and under angular cutoff. The main result is that if the initial value isf ₀ ? 0 with (1 + ¦υ¦²)φ₀ εL ¹ and (1 + ¦υ¦)^s f ₀ε L^∞ for somes > 2, then (1 + ¦υ¦^s'f _tεL ^∞ fort>0 and ess_υ,t sup(1 + ¦υ¦)^s'f(υ, t,) < ∞ for anys′ ? s whens ? 5, and anys′ ? s ifs > 5.