Molecular dynamics simulation of entropy and surface tension for grain boundary of α-Fe

The grain boundary is an interface and the surface tension is one of its important thermodynamic properties. In this paper, the surface tension of the ∑9 grain boundary for α-Fe at various temperatures and pressures is calculated by means of Computer Molecular Dynamics (CMD). The results agree satisfactorily with the experimental data. It is shown that the contribution of entropy to surface tension of grain boundary can be ignored. The project supported by the National Natural Science Foundation of China and the Science Foundation of Chinese Academy of Sciences.     

Clustering of spectra and fractals of regular graphs

We exhibit a characteristic structure of the class of all regular graphs of degree d that stems from the spectra of their adjacency matrices. The structure has a fractal threadlike appearance. Points with coordinates given by the mean and variance of the exponentials of graph eigenvalues cluster around a line segment that we call a filar. Zooming-in reveals that this cluster splits into smaller segments (filars) labeled by the number of triangles in graphs. Further zooming-in shows that the smaller filars split into subfilars labeled by the number of quadrangles in graphs, etc. We call this fractal structure, discovered in a numerical experiment, a multifilar structure. We also provide a mathematical explanation of this phenomenon based on the Ihara-Selberg trace formula, and compute the coordinates and slopes of all filars in terms of Bessel functions of the first kind.     

HIGH ACCURACY NUMERICAL METHOD OF THIN-FILM PROBLEMS IN MICROMAGNETICS

In this paper, a new high accuracy numerical method for the thin-film problems of micron and submicron size ferromagnetic elements is proposed. For the computation of stray field, we use the finite element method(FEM) by introducing a semi-discrete artificial boundary condition [1, 2]. In our numerical experiments about the domain patterns and their movement, we can see that the results are accordant to that of experiments and other numerical methods. Our method are very convenient to deal with arbitrary shape of thin films such as a polygon with high accuracy.     

On the properties of solutions to the Goursat-Darboux problem with boundary and distributed controls

We consider a control system described by the Goursat-Darboux equation. The system is controlled by distributed and boundary controls. The controls are subject to the constraints given as multivalued mappings with closed, possibly nonconvex, values depending on the phase variable. Alongside the initial constraints, we consider the convexified constraints and the constraints whose values are the extreme points of the convexified constraints. We study the questions of existence of solutions and establish connections between the solutions under various constraints.     

Front roughening in three-dimensional imbibition

We investigate the structure and dynamics of the interface between two immiscible liquids in a three-dimensional disordered porous medium. We apply a phase-field model that includes explicitly disorder and discuss both spontaneous and forced imbibition. The structure of the interface is dominated by a length scale ξ_× which arises from liquid conservation. We further show that disorder in the capillary and permeability act on different length scales and give rise to different scalings and structures of the interface properties. We conclude with a range of applications.     

Fixed-sign solutions for a system of singular focal boundary value problems

We consider a system of focal boundary value problems where the nonlinearities may be singular in the independent variable and may also be singular in the dependent arguments. Using Schauder fixed point theorem, we establish criteria such that the system of boundary value problems has at least one fixed-sign solution.     

Generalized 3G theorem and application to relativistic stable process on non-smooth open sets

Let G(x,y) and GD(x,y) be the Green functions of rotationally invariant symmetric α-stable process in Rd and in an open set D, respectively, where 0<α<2. The inequality GD(x,y)GD(y,z)/GD(x,z)?c(G(x,y)+G(y,z)) is a very useful tool in studying (local) Schrödinger operators. When the above inequality is true with c=c(D)∈(0,∞), then we say that the 3G theorem holds in D. In this paper, we establish a generalized version of 3G theorem when D is a bounded κ-fat open set, which includes a bounded John domain. The 3G we consider is of the form GD(x,y)GD(z,w)/GD(x,w), where y may be different from z. When y=z, we recover the usual 3G. The 3G form GD(x,y)GD(z,w)/GD(x,w) appears in non-local Schrödinger operator theory. Using our generalized 3G theorem, we give a concrete class of functions belonging to the non-local Kato class, introduced by Chen and Song, on κ-fat open sets. As an application, we discuss relativistic α-stable processes (relativistic Hamiltonian when α=1) in κ-fat open sets. We identify the Martin boundary and the minimal Martin boundary with the Euclidean boundary for relativistic α-stable processes in κ-fat open sets. Furthermore, we show that relative Fatou type theorem is true for relativistic stable processes in κ-fat open sets. The main results of this paper hold for a large class of symmetric Markov processes, as are illustrated in the last section of this paper. We also discuss the generalized 3G theorem for a large class of symmetric stable Lévy processes.     

Non-self-adjoint sturm-liouville operators with matrix potentials

We obtain asymptotic formulas for non-self-adjoint operators generated by the Sturm-Liouville system and quasiperiodic boundary conditions. Using these asymptotic formulas, we obtain conditions on the potential for which the system of root vectors of the operator under consideration forms a Riesz basis.     

On the convergence rate and optimization of a numerical method with splitting of boundary conditions for the stokes system in a spherical layer in the axisymmetric case: Modification for thick layers

The convergence rate of a fast-converging second-order accurate iterative method with splitting of boundary conditions constructed by the authors for solving an axisymmetric Dirichlet boundary value problem for the Stokes system in a spherical gap is studied numerically. For R/r exceeding about 30, where r and R are the radii of the inner and outer boundary spheres, it is established that the convergence rate of the method is lower (and considerably lower for large R/r) than the convergence rate of its differential version. For this reason, a really simpler, more slowly converging modification of the original method is constructed on the differential level and a finite-element implementation of this modification is built. Numerical experiments have revealed that this modification has the same convergence rate as its differential counterpart for R/r of up to 5 × 10³. When the multigrid method is used to solve the split and auxiliary boundary value problems arising at iterations, the modification is more efficient than the original method starting from R/r ～ 30 and is considerably more efficient for large values of R/r. It is also established that the convergence rates of both methods depend little on the stretching coefficient η of circularly rectangular mesh cells in a range of η that is well sufficient for effective use of the multigrid method for arbitrary values of R/r smaller than ～ 5 × 10³.     

四阶微分方程奇异边值问题的正解

本文利用锥上的不动点指数理论,在更一般的边界条件下讨论了一类四阶奇异边值问题正解的存在性并给出了应用．