On solvability of linear differential equations in ${\Bbb R}^{\Bbb N}$

We construct $x^0$ in ${\Bbb R}^{\Bbb N}$ and a row-finite matrix $T=\{T_{i,j}(t)\}_{i,j\in\N}$ of polynomials of one real variable $t$ such that the Cauchy problem $\dot x(t)=T_tx(t)$, $x(0)=x^0$ in the Fr\'echet space $\R^\N$ has no solutions. We also construct a row-finite matrix $A=\{A_{i,j}(t)\}_{i,j\in\N}$ of $C^\infty(\R)$ functions such that the Cauchy problem $\dot x(t)=A_tx(t)$, $x(0)=x^0$ in ${\Bbb R}^{\Bbb N}$ has no solutions for any $x^0\in{\Bbb R}^{\Bbb N}\setminus\{0\}$. We provide some sufficient condition of solvability and of unique solvability for linear ordinary differential equations $\dot x(t)=T_tx(t)$ with matrix elements $T_{i,j}(t)$ analytically dependent on $t$.

Some results on solvability of ordinary linear differential equations in locally convex spaces

Let $\Gamma$ be the class of sequentially complete locally convex spaces such that an existence theorem holds for the linear Cauchy problem $\dot x = Ax$, $x(0) = x_0$ with respect to functions $x: R\to E$. It is proved that if $E\in \Gamma$, then $E\times R^A$ is-an-element-of $\Gamma$ for an arbitrary set $A$. It is also proved that a topological product of infinitely many infinite-dimensional Frechet spaces, each not isomorphic to $\omega$, does not belong to $\Gamma$.

R-matrix-Floquet theory of multiphoton processes:VIII. A linear equations method

A new linear equations method for calculating the R-matrix, which arises in the R-matrix-Floquet theory of multiphoton processes, is introduced. This method replaces the diagonalization of the Floquet Hamiltonian matrix by the solution of a set of linear simultaneous equations which are solved, in the present work, by the conjugate gradient method. This approach uses considerably less computer memory and can be readily ported onto parallel computers. It will thus enable much larger problems of current interest to be treated. This new method is tested by applying it to three-photon ionization of helium at frequencies where double resonances with a bound state and autoionizing states are important. Finally, an alternative linear equations method, which avoids the explicit calculation of the R-matrix by incorporating the boundary conditions directly, is described in an appendix.

Blind Characterisation of Sensors with Second-Order Dynamic Response

Compensation for the dynamic response of a temperature sensor usually involves the estimation of its input on the basis of the measured output and model parameters. In the case of temperature measurement, the sensor dynamic response is strongly dependent on the measurement environment and fluid velocity. Estimation of time-varying sensor model parameters therefore requires continuous textit{in situ} identification. This can be achieved by employing two sensors with different dynamic properties, and exploiting structural redundancy to deduce the sensor models from the resulting data streams. Most existing approaches to this problem assume first-order sensor dynamics. In practice, however second-order models are more reflective of the dynamics of real temperature sensors, particularly when they are encased in a protective sheath. As such, this paper presents a novel difference equation approach to solving the blind identification problem for sensors with second-order models. The approach is based on estimating an auxiliary ARX model whose parameters are related to the desired sensor model parameters through a set of coupled non-linear algebraic equations. The ARX model can be estimated using conventional system identification techniques and the non-linear equations can be solved analytically to yield estimates of the sensor models. Simulation results are presented to demonstrate the efficiency of the proposed approach under various input and parameter conditions.

An approach to determine the local boundary of small disturbance voltage stability region in a cut-set power space:基于微扰的割集电压稳定域局部边界求解方法

A new approach to determine the local boundary of voltage stability region in a cut-set power space (CVSR) is presented. Power flow tracing is first used to determine the generator-load pair most sensitive to each branch in the interface. The generator-load pairs are then used to realize accurate small disturbances by controlling the branch power flow in increasing and decreasing directions to obtain new equilibrium points around the initial equilibrium point. And, continuous power flow is used starting from such new points to get the corresponding critical points around the initial critical point on the CVSR boundary. Then a hyperplane cross the initial critical point can be calculated by solving a set of linear algebraic equations. Finally, the presented method is validated by some systems, including New England 39-bus system, IEEE 118-bus system, and EPRI-1000 bus system. It can be revealed that the method is computationally more efficient and has less approximation error. It provides a useful approach for power system online voltage stability monitoring and assessment. This work is supported by National Natural Science Foundation of China (No. 50707019), Special Fund of the National Basic Research Program of China (No. 2009CB219701), Foundation for the Author of National Excellent Doctoral Dissertation of PR China (No. 200439), Tianjin Municipal Science and Technology Development Program (No. 09JCZDJC25000), National Major Project of Scientific and Technical Supporting Programs of China During the 11th Five-year Plan Period (No. 2006BAJ03A06). ©2009 State Grid Electric Power Research Institute Press.

Sequentially continuous non-linear fundamental systems of solutions of affine equations in locally convex spaces

According to the Mickael's selection theorem any surjective continuous linear operator from one Fr\'echet space onto another has a continuous (not necessarily linear) right inverse. Using this theorem Herzog and Lemmert proved that if $E$ is a Fr\'echet space and $T:E\to E$ is a continuous linear operator such that the Cauchy problem $\dot x=Tx$, $x(0)=x_0$ is solvable in $[0,1]$ for any $x_0\in E$, then for any $f\in C([0,1],E)$, there exists a continuos map $S:[0,1]\times E\to E$, $(t,x)\mapsto S_tx$ such that for any $x_0\in E$, the function $x(t)=S_tx_0$ is a solution of the Cauchy problem $\dot x(t)=Tx(t)+f(t)$, $x(0)=x_0$ (they call $S$ a fundamental system of solutions of the equation $\dot x=Tx+f$). We prove the same theorem, replacing "continuous" by "sequentially continuous" for locally convex spaces from a class which contains strict inductive limits of Fr\'echet spaces and strong duals of Fr\'echet--Schwarz spaces and is closed with respect to finite products and sequentially closed subspaces. The key-point of the proof is an extension of the theorem on existence of a sequentially continuous right inverse of any surjective sequentially continuous linear operator to some class of non-metrizable locally convex spaces.

Boundary conditions in the simplest model of linear and second harmonic magneto-optical effects

This paper is concerned with linear and nonlinear magneto- optical effects in multilayered magnetic systems when treated by the simplest phenomenological model that allows their response to be represented in terms of electric polarization, The problem is addressed by formulating a set of boundary conditions at infinitely thin interfaces, taking into account the existence of surface polarizations. Essential details are given that describe how the formalism of distributions (generalized functions) allows these conditions to be derived directly from the differential form of Maxwell's equations. Using the same formalism we show the origin of alternative boundary conditions that exist in the literature. The boundary value problem for the wave equation is formulated, with an emphasis on the analysis of second harmonic magneto-optical effects in ferromagnetically ordered multilayers. An associated problem of conventions in setting up relationships between the nonlinear surface polarization and the fundamental electric field at the interfaces separating anisotropic layers through surface susceptibility tensors is discussed. A problem of self- consistency of the model is highlighted, relating to the existence of resealing procedures connecting the different conventions. The linear approximation with respect to magnetization is pursued, allowing rotational anisotropy of magneto-optical effects to be easily analyzed owing to the invariance of the corresponding polar and axial tensors under ordinary point groups. Required representations of the tensors are given for the groups infinitym, 4mm, mm2, and 3m, With regard to centrosymmetric multilayers, nonlinear volume polarization is also considered. A concise expression is given for its magnetic part, governed by an axial fifth-rank susceptibility tensor being invariant under the Curie group infinityinfinitym.

Ion-acoustic waves in a plasma consisting of adiabatic warm ions, non-isothermal electrons and a weakly relativistic electron beam: linear and higher-order nonlinear effects

The nonlinear propagation of finite amplitude ion acoustic solitary waves in a plasma consisting of adiabatic warm ions, nonisothermal electrons, and a weakly relativistic electron beam is studied via a two-fluid model. A multiple scales technique is employed to investigate the nonlinear regime. The existence of the electron beam gives rise to four linear ion acoustic modes, which propagate at different phase speeds. The numerical analysis shows that the propagation speed of two of these modes may become complex-valued (i.e., waves cannot occur) under conditions which depend on values of the beam-to-background-electron density ratio , the ion-to-free-electron temperature ratio , and the electron beam velocity v0; the remaining two modes remain real in all cases. The basic set of fluid equations are reduced to a Schamel-type equation and a linear inhomogeneous equation for the first and second-order potential perturbations, respectively. Stationary solutions of the coupled equations are derived using a renormalization method. Higher-order nonlinearity is thus shown to modify the solitary wave amplitude and may also deform its shape, even possibly transforming a simple pulse into a W-type curve for one of the modes. The dependence of the excitation amplitude and of the higher-order nonlinearity potential correction on the parameters , , and v0 is numerically investigated.

Linear and nonlinear dynamics of a dust bicrystal consisting of positive and negative dust particles

A dusty plasma crystalline configuration consisting of charged dust grains of alternating charge sign (.../+/-/+/-/+/...) and mass is considered. Both charge and mass of each dust species are taken to be constant. Considering the equations of longitudinal motion, a dispersion relation for linear longitudinal vibrations is derived from first principles and then analyzed. Two harmonic modes are obtained, namely, an acoustic mode and an inverse-dispersive optic-like one. The nonlinear aspects of acoustic longitudinal dust grain motion are addressed via a generalized Boussinesq (and, alternatively, a generalized Korteweg-de Vries) description. (C) 2005 American Institute of Physics.

Relativistic laser pulse compression in plasmas with a linear axial density gradient

The self-compression of a relativistic Gaussian laser pulse propagating in a non-uniform plasma is investigated. A linear density inhomogeneity (density ramp) is assumed in the axial direction. The nonlinear Schrodinger equation is first solved within a one-dimensional geometry by using the paraxial formalism to demonstrate the occurrence of longitudinal pulse compression and the associated increase in intensity. Both longitudinal and transverse self-compression in plasma is examined for a finite extent Gaussian laser pulse. A pair of appropriate trial functions, for the beam width parameter (in space) and the pulse width parameter (in time) are defined and the corresponding equations of space and time evolution are derived. A numerical investigation shows that inhomogeneity in the plasma can further boost the compression mechanism and localize the pulse intensity, in comparison with a homogeneous plasma. A 100 fs pulse is compressed in an inhomogeneous plasma medium by more than ten times. Our findings indicate the possibility for the generation of particularly intense and short pulses, with relevance to the future development of tabletop high-power ultrashort laser pulse based particle acceleration devices and associated high harmonic generation. An extension of the model is proposed to investigate relativistic laser pulse compression in magnetized plasmas.

Discretization and solution of the inhomogeneous Bogoliubov-de Gennes equations

In the presence of inhomogeneities, defects and currents, the equations describing a Bose-condensed ensemble of alkali atoms have to be solved numerically. By combining both linear and nonlinear equations within a Discrete Variable Representation framework, we describe a computational scheme for the solution of the coupled Bogoliubov-de Gennes (BdG) and nonlinear Schrodinger (NLS) equations for fields in a 3D spheroidal potential. We use the method to calculate the collective excitation spectrum and quasiparticle mode densities for excitations of a Bose condensed gas in a spheroidal trap. The method is compared against finite-difference and spectral methods, and we find the DVR computational scheme to be superior in accuracy and efficiency for the cases we consider. (C) 2004 Elsevier B.V. All rights reserved.

Performance of polarisation functionals for linear and nonlinear optical properties of bulk zinc chalcogenides ZnX (X=S,Se,Te).

We calculated the frequency dependent macroscopic dielectric function and second-harmonic generation of cubic ZnS, ZnSe and ZnTe within time-dependent density-polarisation functional theory. The macroscopic dielectric function is calculated in a linear response framework, and second-harmonic generation in a real-time framework. The macroscopic exchange–correlation electric field that enters the time-dependent Kohn–Sham equations and accounts for long range correlation is approximated as a simple polarisation functional αP, where P is the macroscopic polarisation. Expressions for α are taken from the recent literature. The performance of the resulting approximations for the exchange–correlation electric field is analysed by comparing the theoretical spectra with experimental results and results obtained at the levels of the independent particle approximation and the random-phase approximation. For the dielectric function we also compare with state-of-the art calculations at the level of the Bethe–Salpeter equation.