Spectral Simplicity and Asymptotic Separation of Variables

We describe a method for comparing the spectra of two real-analytic families, (a _t ) and (q _t ), of quadratic forms that both degenerate as a positive parameter t tends to zero. We suppose that the family (a _t ) is amenable to ‘separation of variables’ and that each eigenspace of a _t is 1-dimensional for some t. We show that if (q _t ) is asymptotic to (a _t ) at first order as t → 0, then the eigenspaces of (q _t ) are also 1-dimensional for all but countably many t. As an application, we prove that for the generic triangle (simplex) in Euclidean space (constant curvature space form) each eigenspace of the Laplacian acting on Dirichlet functions is 1-dimensional.     

The modified Korteweg-de Vries equation on the quarter plane with t-periodic data

We study the modified Korteweg-de Vries equation posed on the quarter plane with asymptotically t-periodic Dirichlet boundary datum u(0,t) in the sense that u(0,t) tends to a periodic function g?₀ (t) with period τ as t → ∞. We consider the perturbative expansion of the solution in a small ε > 0. Here we show that if the unknown boundary data u_x(0,t) and u_xx(0,t) are asymptotically t-periodic with period τ which tend to the functions g?₁ (t) and g?₂ (t) as t → ∞, respectively, then the periodic functions g?₁ (t) and g?₂ (t) can be uniquely determined in terms of the function g?₀ (t). Furthermore, we characterize the Fourier coefficients of g?₁ (t) and g?₂ (t) to all orders in the perturbative expansion by solving an infinite system of algebraic equations. As an illustrative example, we consider the case of a sine-wave as Dirichlet datum and we explicitly determine the coefficients for large t up to the third order in the perturbative expansion.     

Possible Loss and Recovery of Gibbsianness¶During the Stochastic Evolution of Gibbs Measures

We consider Ising-spin systems starting from an initial Gibbs measure ν and evolving under a spin-flip dynamics towards a reversible Gibbs measure μ≠ν. Both ν and μ are assumed to have a translation-invariant finite-range interaction. We study the Gibbsian character of the measure νS(t) at time t and show the following: (1) For all ν and μ, νS(t) is Gibbs for small t. (2) If both ν and μ have a high or infinite temperature, then νS(t) is Gibbs for all t > 0. (3) If ν has a low non-zero temperature and a zero magnetic field and μ has a high or infinite temperature, then νS(t) is Gibbs for small t and non-Gibbs for large t. (4) If ν has a low non-zero temperature and a non-zero magnetic field and μ has a high or infinite temperature, then νS(t) is Gibbs for small t, non-Gibbs for intermediate t, and Gibbs for large t. The regime where μ has a low or zero temperature and t is not small remains open. This regime presumably allows for many different scenarios. Received: 26 April 2001 / Accepted: 10 October 2001     

Solitonic properties for a forced generalized variable-coefficient Korteweg-de Vries equation for the atmospheric blocking phenomenon

In this paper, investigation is given to a forced generalized variable-coefficient Korteweg-de Vries equation for the atmospheric blocking phenomenon. Based on the Lax pair, under certain variable-coefficient-dependent constraints, we present an infinite sequence of the conservation laws. Through the Riccati equations obtained from the Lax pair, a Wahlquist-Estabrook-type Bäcklund transformation (BT) is derived, based on which the nonlinear superposition formula as well as one- and two-soliton-like solutions are obtained. Via the truncated Painlevé expansion, we give a Painlevé BT, along with the one-soliton-like solutions. With the Painlevé BT, bilinear forms are constructed, and we get a bilinear BT as well as the corresponding one-soliton-like solutions. Bell-type bright and dark soliton-like waves and kink-type soliton-like waves are observed, respectively. Graphic analysis shows that (1) the velocities of the soliton-like waves are related to h(t), d(t), f(t) and R(t), while the soliton-like wave amplitudes just depend on f(t), and (2) with the nonzero f(t) and R(t), soliton-like waves propagate on the varying backgrounds, where h(t), d(t) and f(t) are the dispersive, dissipative and line-damping coefficients, respectively, R(t) is the external-force term, and t is the scaled time coordinate.     

Growth and Roughness of the Interface for Ballistic Deposition

In ballistic deposition (BD), (d+1)-dimensional particles fall sequentially at random towards an initially flat, large but bounded d-dimensional surface, and each particle sticks to the first point of contact. For both lattice and continuum BD, a law of large numbers in the thermodynamic limit establishes convergence of the mean height and surface width (sample standard deviation of the height) of the interface to constants h(t) and w(t), respectively, depending on time t. We show that h(t) is asymptotically linear in t, while (w(t))² grows at least logarithmically in t when d=1. We use duality results showing that w(t) can be interpreted as the standard deviation of the height for deposition onto a surface growing from a single point.     

Laser thermal response of a finite slab as a function of the laser pulse parameters

This paper deals with the problem of heating a finite slab using laser radiation in relation to the parameters characterizing the laser pulse, namely: q_max(W/m²), the maximum laser power density, t₀ the time interval required to reach q_max and t_d, the pulse time duration. The pulse shape q(t) is suggested in the form: q(t)=βq_max(t/t_d)(1-(t/t_d))exp-B(t-t₀/t_d), where β and B are parameters. Fitting with published experimental pulse [Ready JF. Effects due to absorption of laser radiation. J Appl Phys 1965;36:462–68] is made. Fourier series expansion technique is considered to solve the problem. The critical time required to initiate melting t_m is estimated for four metallic elements and five semiconductors, namely: Al, Cu, Ag, Au (aluminum, copper, silver, and gold), cadmium sulfide, germanium, silicon, alpha beryllium oxide, and silicon carbide. Five pulses with different characteristic parameters are considered.Computations revealed that the thermal response of the targets is highly affected by q_max and t_o, while the pulse time duration is less effective in determining the value of t_m. Moreover, it is revealed that the relation between t_m and the melting temperature for the same laser pulse is nonlinear for the considered targets under the indicated conditions.     

Reacting polymers with highly correlated initial conditions

We propose and theoretically study an experiment designed to measure short-time polymer reaction kinetics in melts or dilute solutions. The photolysis of groups centrally located along chain backbones, one group per chain, creates pairs of spatially highly correlated macroradicals. We calculate time-dependent rate coefficients κ(t) governing their first-order recombination kinetics, which are novel on account of the far-from-equilibrium initial conditions. In dilute solutions (good solvents) reaction kinetics are intrinsically weak, despite the highly reactive radical groups involved. This leads to a generalised mean-field kinetics in which the rate of radical density decay - ∼S(t), where S(t) ∼t ^{- (1 + g/3)} is the equilibrium return probability for 2 reactive groups, given initial contact. Here g≈ 0.27 is the correlation hole exponent for self-avoiding chain ends. For times beyond the longest coil relaxation time τ, - ∼S(t) remains true, but center of gravity coil diffusion takes over with rms displacement of reactive groups x(t) ∼t ^1/2 and S(t) ∼ 1/x ³(t). At the shortest times ( t 10^-6s), recombination is inhibited due to spin selection rules and we find ∼tS(t). In melts, kinetics are intrinsically diffusion-controlled, leading to entirely different rate laws. During the regime limited by spin selection rules, the density of radicals decays linearly, n(0) - n(t) ∼t. At longer times the standard result - ∼d ³(t)/d (for randomly distributed ends) is replaced by ∼d2x ³(t)/d ² for these correlated initial conditions. The long-time behavior, t > τ, has the same scaling form in time as for dilute solutions. Received 18 May 2000     

The Modified Korteweg-de Vries Equation on the Half-Line with a Sine-Wave as Dirichlet Datum

Boundary value problems for integrable nonlinear evolution PDEs, like the modified KdV equation, formulated on the half-line can be analyzed by the so-called unified transform method. For the modified KdV equation, this method yields the solution in terms of the solution of a matrix Riemann-Hilbert problem uniquely determined in terms of the initial datum q(x,0), as well as of the boundary values {q(0, t),qx(0, t),qxx(0, t)}. For the Dirichlet problem, it is necessary to characterize the unknown boundary values qx(0, t) and qxx(0, t) in terms of the given data q(x, 0) and q(0, t). It is shown here that in the particular case of a vanishing initial datum and of a sine wave as Dirichlet datum, qx(0, t) and qxx(0, t) can be computed explicitly at least up to third order in a perturbative expansion and that at least up to this order, these functions are asymptotically periodic for large t.     

Ballistic diffusion induced by non-Gaussian noise

In this letter,we have analyzed the diffusive behavior of a Brownian particle subject to both internal Gaussian thermal and external non-Gaussian noise sources.We discuss two time correlation functions C(t) of the non-Gaussian stochastic process,and find that they depend on the parameter q,indicating the departure of the non-Gaussian noise from Gaussian behavior:for q ≤ 1,C(t) is fitted very well by the first-order exponentially decaying curve and approaches zero in the longtime limit,whereas for q 1,C(t) can be approximated by a second-order exponentially decaying function and converges to a non-zero constant.Due to the properties of C(t),the particle exhibits a normal diffusion for q ≤ 1,while for q 1 the non-Gaussian noise induces a ballistic diffusion,i.e.,the long-time mean square displacement of the free particle reads [x(t)-x(t)]2∝t2.     

On integrability of generalized Veronese curves of distributions

Given a 1-parameter family of 1-forms γ(t) = γ₀+tγ₁+ ···+tⁿψ_n, consider the condition dγ(t)γ(t) = 0 (of integrability for the annihilated by γ(t) distribution w(t)). We prove that in order that this condition is satisfied for any t it is sufficient that it is satisfied for N = n + 3 different values of t (the corresponding implication for N = 2n + 1 is obvious). In fact we give a stronger result dealing with distributions of higher codimension. This result is related to the so-called Veronese webs and can be applied in the theory of bihamiltonian structures.     

Transport Equation Evaluation of Coupled Continuous Time Random Walks

The transport behavior of a migrating particle in a disordered medium is exhibited in the solution of a transport equation derived from a coupled continuous time random walk (CTRW). A core aspect of CTRW is the spectrum of transitions in displacement s and time t, ψ(s,t), that characterizes the disordered system, which determine the transport. In many applications the CTRW approach has successfully accounted for the anomalous or non-Fickian nature of the particle plume propagation based on a power-law dependence ψ(t) in a decoupled p(s)ψ(t) approximation to ψ(s,t). For example, this power-law dependence in t derives from the complex Darcy flow fields in geological formations. Recently, the fully coupled CTRW was analyzed using a particle tracking approach, demonstrating that the decoupled approximation is valid only for a compact distribution of s. In this paper we solve the nonlocal-in-time transport equation with a ψ(s,t) containing a power-law dependence in both s (a Lévy-like distribution) and t, which necessitates the strong s,t coupling. We show enhanced transport behavior (relative to the plume propagation behavior reported in the literature) that derives from the rare large displacements in s (limited by the transition t). The interplay between the two coupled power laws is clearly shown in the changes in the breakthrough curves in the arrival times, dispersion and dependence on the velocity (v=s/t) distribution. Similar enhancements are exhibited in the particle tracking results.     

Smoothening of a random surface as it results from the edwards-wilkinson kinetic equation

The time evolution of a random surfacez=h(r, t) (r=x, y) formed by a deposition process of the Edwards-Wilkinson type is discussed. The discussion is based on the author’s former derivation of the autocorrelation functionA _h(|r − r′|,t, t′)=〈h(r,t)h(r′,t′)〉 of the height functionh(r,t) under the assumption of a stochastic initial condition [V. Bezák: Acta Physica Univ. Comenianae39 (1998) 135]. Under the assumption of a steady (non-zero) deposition rate, the varianceσ _h ² (t)=〈[h(r,t)]²〉 increases logarithmically in time whilst the correlation lengthl _h(t) of the height functionh(r,t) increases as ∼t ^1/2. Therefore, the ratioσ _h(t)/l _h (t) tends to zero and the surfacez=h(r,t) does always tend towards a smoothened appearance. This work has been supported by the Slovak Grant Agency VEGA under contract No. 1/4319/97.     

Muon charge exchange and muonium spin exchange in gases

The spin dynamics of the positive muon undergoing cyclic charge exchange (μ⁺⇌Mu) while it is slowing down is compared with that in the case of repeated muonium spin exchange after thermalization. The expectation value of the spin polarization at timet aftern spin exchange collisions (att ₁,t ₂, ...,t _n ) are calculated explicitly from time dependent wave functions, and the quantity is averaged over the statistical distribution of the times of collisions and over all possible numbers of collisions betweent=0 and the observation timet. This result is complementary and equivalent to the conventional density matrix formalism, but offers an insight into the roles of spin flip and spin non-flip processes. The neutral fraction during slowing down is also discussed.     

Density-matrix renormalization study of the Hubbard model [4]on a Bethe lattice

The half-filled Hubbard model on the Bethe lattice with coordination number z=3 is studied using the density-matrix renormalization group (DMRG) method. Ground-state properties such as the energy per site E, average local magnetization , its fluctuations and various spin correlation functions are determined as a function of the Coulomb interaction strength U/t. The local magnetic moments increase monotonically with increasing Coulomb repulsion U/t showing antiferromagnetic order between nearest neighbors []. At large U/t, is strongly reduced with respect to the saturation value 1/2 due to exchange fluctuations between nearest neighbors (NN) spins [ for ]. shows a maximum for U/t=2.4-2.9 that results from the interplay between the usual increase of with increasing U/t and the formation of important permanent moments at large U/t. While NN sites show antiferromagnetic spin correlations that increase with increasing Coulomb repulsion, the next NN sites are very weakly correlated over the whole range of U/t. The DMRG results are discussed and compared with tight-binding calculations for U=0, independent DMRG studies for the Heisenberg model and simple first-order perturbation estimates. Received 8 February 1999 and Received in final form 14 June 1999     

Non-linear scission/recombination kinetics of living polymerization

Living polymers are formed by reversible association of primary units (unimers). Generally the chain statistical weight involves a factor σ < 1 suppressing short chains in comparison with free unimers. Living polymerization is a sharp thermodynamic transition for σ ≪ 1 which is typically the case. We show that this sharpness has an important effect on the kinetics of living polymerization (one-dimensional association). The kinetic model involves i) the unimer activation step (a transition to an assembly-competent state); ii) the scission/recombination processes providing growth of polymer chains and relaxation of their length distribution. Analyzing the polymerization with no chains but unimers at t = 0 , with initial concentration of unimers M ≳ M ^* (M^* is the critical polymerization concentration), we determine the time evolution of the chain length distribution and find that: 1) for M ^* ≪ M ≪ M ^*/σ the kinetics is characterized by 5 distinct time stages demarcated by 4 characteristic times t₁, t₂, t₃ and t^*; 2) there are transient regimes (t ₁ ≲ t ≲ t ₃) when the molecular-weight distribution is strongly non-exponential; 3) the chain scissions are negligible at times shorter than t₂. The chain growth is auto-accelerated for t ₁ ≲ t ≲ t ₂ : the cut-off chain length (= polymerization degree 〈n〉_w N ₁ ∝ t ² in this regime. 4) For t ₂ < t < t ₃ the length distribution is characterized by essentially 2 non-linear modes; the shorter cut-off length N₁ is decreasing with time in this regime, while the length scale N₂ of the second mode is increasing. (5) The terminal relaxation time of the polymer length distribution, t^*, shows a sharp maximum in the vicinity of M^*; the effective exponent is as high as ∼ σ^-1/3 just above M^*.     

Berechnung der Amplitudenfluktuation eines mit Laserlicht gepumpten optischen parametrischen Oszillators

The intensity fluctuations of the radiation of an optical parametric oscillator have been calculated by taking into account the corresponding fluctuations of the pump light. The latter are shown to have a strong influence on the fluctuations of signal and idler waves. In contrary to the laser case, they are essentially reproduced in the signal and idler radiation; more quantitatively, the correlation function 〈?(t₁) ?(t₂)〉 where ? denotes the deviation of the amplitude of the signal (or idler) wave from its mean value is, for t₁ = t₂, equal to the corresponding quantity for the pump wave apart from a reduction factor, whose value depends on the resonator losses and the relative excitation of both the pumping laser and the optical parametric oscillator and is found to be of the order of one for the normal case. It is interesting to note that for sufficiently strong pumping, the correlation function 〈?(t₁) ?(t₂)〉 approaches the value zero for |t₁ – t₂| → ∞ not monotonically but in form of damped oscillations.     

利用杨辉三角形对称性推导高阶运动微分方程

利用Mathematica数学软件计算函数r=r(q(t),t)各变量之间偏导和高阶导数的关系，发现具有杨辉三角形对称性.结合杨辉三角形的对称性规律和牛顿第二定律推导出了高阶运动微分方程，并讨论了理想约束系统下的高阶运动微分方程. 关键词： 杨辉三角形 牛顿第二定律 高阶运动微分方程 高阶力变率 高阶速度能量 理想约束     

Vicious random walkers in the limit of a large number of walkers

The vicious random walker problem on a line is studied in the limit of a large number of walkers. The multidimensional integral representing the probability that thep walkers will survive a timet (denotedP _t ^(p) ) is shown to be analogous to the partition function of a particular one-component Coulomb gas. By assuming the existence of the thermodynamic limit for the Coulomb gas, one can deduce asymptotic formulas forP _t ^(p) in the large-p, large-t limit. A straightforward analysis gives rigorous asymptotic formulas for the probability that after a timet the walkers are in their initial configuration (this event is termed a reunion). Consequently, asymptotic formulas for the conditional probability of a reunion, given that all walkers survive, are derived. Also, an asymptotic formula for the conditional probability density that any walker will arrive at a particular point in timet, given that allp walkers survive, is calculated in the limittp.     

The second torsional state of acetic acid

We have fitted to within experimental accuracy a data set for acetic acid CH₃COOH consisting of 2103 v_t=0, 1111 v_t=1 and 634 v_t=2 microwave lines, using the so-called “rho axis method” and a model extended to include perturbation terms through eighth order. The previously published 2109 v_t=0 and 409 v_t=1 microwave lines have been extended by new measurements from Kharkov. The final fit requires only 62 parameters to achieve a unitless weighted standard deviation for the whole fit of 0.85 for a total of 3848 lines and includes transitions with J27 (15) and |K_a|13 (8) for the v_t=2 A (E)-species, respectively. This result represents a significant improvement over our previous fit which was dominated by v_t=0 transitions. A calculation of the linestrengths of torsion–rotation transitions is also provided.