Coupling functions for domino tilings of Aztec diamonds

The inverse Kasteleyn matrix of a bipartite graph holds much information about the perfect matchings of the system such as local statistics which can be used to compute local and global asymptotics. In this paper, we consider three different weightings of domino tilings of the Aztec diamond and show using recurrence relations, that we can compute the inverse Kasteleyn matrix. These weights are the one-periodic weighting where the horizontal edges have one weight and the vertical edges have another weight, the q^vol$q^{\mathrm{vol}}$ weighting which corresponds to multiplying the product of tile weights by q if we add a ‘box’ to the height function and the two-periodic weighting which exhibits a flat region with defects in the center.     

Laurent biorthogonal polynomials, q-Narayana polynomials and domino tilings of the Aztec diamonds

A Toeplitz determinant whose entries are described by a q-analogue of the Narayana polynomials is evaluated by means of Laurent biorthogonal polynomials which allow of a combinatorial interpretation in terms of Schröder paths. As an application, a new proof is given to the Aztec diamond theorem by Elkies, Kuperberg, Larsen and Propp concerning domino tilings of the Aztec diamonds. The proof is based on the correspondence with non-intersecting Schröder paths developed by Johansson.     

Domino tilings of the expanded Aztec diamond

The expanded Aztec diamond is a generalized version of the Aztec diamond, with an arbitrary number of long columns and long rows in the middle. In this paper, we count the number of domino tilings of the expanded Aztec diamond. The exact number of domino tilings is given by recurrence relations of state matrices by virtue of the state matrix recursion algorithm, recently developed by the author to solve various two-dimensional regular lattice model enumeration problems.     

A note on the structure of spaces of domino tilings

We study spaces of tilings, formed by tilings which are on a geodesic between two fixed tilings of the same domain (the distance is defined using local flips). We prove that each space of tilings is homeomorphic to an interval of tilings of a domain when flips are classically directed by height functions.     

A shuffling theorem for reflectively symmetric tilings

The author and Rohatgi recently proved a ‘shuffling theorem’ for doubly-dented hexagons. In particular, they showed that shuffling removed unit triangles along a horizontal axis in a hexagon changes the tiling number by only a simple multiplicative factor. In this paper, we consider a similar phenomenon for a symmetry class of tilings, namely, the reflectively symmetric tilings. We also prove several shuffling theorems for halved hexagons.     

Lozenge tilings of hexagons with arbitrary dents

Eisenkölbl gave a formula for the number of lozenge tilings of a hexagon on the triangular lattice with three unit triangles removed from along alternating sides. In earlier work, the first author extended this to the situation when an arbitrary set of unit triangles is removed from along alternating sides of the hexagon. In this paper we address the general case when an arbitrary set of unit triangles is removed from along the boundary of the hexagon.     

Symmetry classes of spanning trees of Aztec diamonds and perfect matchings of odd squares with a unit hole

We say that two graphs are similar if their adjacency matrices are similar matrices. We show that the square grid G _n of order n is similar to the disjoint union of two copies of the quartered Aztec diamond QAD _n−1 of order n−1 with the path P _n ⁽²⁾ on n vertices having edge weights equal to 2. Our proof is based on an explicit change of basis in the vector space on which the adjacency matrix acts. The arguments verifying that this change of basis works are combinatorial. It follows in particular that the characteristic polynomials of the above graphs satisfy the equality P(G _n )=P(P _n ⁽²⁾)[P(QAD _n−1)]². On the one hand, this provides a combinatorial explanation for the “squarishness” of the characteristic polynomial of the square grid—i.e., that it is a perfect square, up to a factor of relatively small degree. On the other hand, as formulas for the characteristic polynomials of the path and the square grid are well known, our equality determines the characteristic polynomial of the quartered Aztec diamond. In turn, the latter allows computing the number of spanning trees of quartered Aztec diamonds. We present and analyze three more families of graphs that share the above described “linear squarishness” property of square grids: odd Aztec diamonds, mixed Aztec diamonds, and Aztec pillowcases—graphs obtained from two copies of an Aztec diamond by identifying the corresponding vertices on their convex hulls. We apply the above results to enumerate all the symmetry classes of spanning trees of the even Aztec diamonds, and all the symmetry classes not involving rotations of the spanning trees of odd and mixed Aztec diamonds. We also enumerate all but the base case of the symmetry classes of perfect matchings of odd square grids with the central vertex removed. In addition, we obtain a product formula for the number of spanning trees of Aztec pillowcases. Research supported in part by NSF grant DMS-0500616.     

Mixed boundary value problems of diffraction by a half‐plane with an obstacle perpendicular to the boundary

The paper is devoted to the analysis of wave diffraction problems modeled by classes of mixed boundary conditions and the Helmholtz equation, within a half‐plane with a crack. Potential theory together with Fredholm theory, and explicit operator relations, are conveniently implemented to perform the analysis of the problems. In particular, an interplay between Wiener–Hopf plus/minus Hankel operators and Wiener–Hopf operators assumes a relevant preponderance in the final results. As main conclusions, this study reveals conditions for the well‐posedness of the corresponding boundary value problems in certain Sobolev spaces and equivalent reduction to systems of Wiener–Hopf equations. Copyright © 2013 John Wiley & Sons, Ltd.     

Solution for a problem of linear plane elasticity with mixed boundary conditions on an ellipse by the method of boundary integrals

A numerical boundary integral scheme is proposed for the solution of the system of field equations of plane, linear elasticity in stresses for homogeneous, isotropic media in the domain bounded by an ellipse under mixed boundary conditions. The stresses are prescribed on one half of the ellipse, while the displacements are given on the other half. The method relies on previous analytical work within the Boundary Integral Method [1], [2].The considered problem with mixed boundary conditions is replaced by two subproblems with homogeneous boundary conditions, one of each type, having a common solution. The equations are reduced to a system of boundary integral equations, which is then discretized in the usual way and the problem at this stage is reduced to the solution of a rectangular linear system of algebraic equations. The unknowns in this system of equations are the boundary values of four harmonic functions which define the full elastic solution inside the domain, and the unknown boundary values of stresses or displacements on proper parts of the boundary.On the basis of the obtained results, it is inferred that the tangential stress component on the fixed part of the boundary has a singularity at each of the two separation points, thought to be of logarithmic type. A tentative form for the singular solution is proposed to calculate the full solution in bulk directly from the given boundary conditions using the well-known Boundary Collocation Method. It is shown that this addition substantially decreases the error in satisfying the boundary conditions on some interval not containing the singular points.The obtained results are discussed and boundary curves for unknown functions are provided, as well as three-dimensional plots for quantities of practical interest. The efficiency of the used numerical schemes is discussed, in what concerns the number of boundary nodes needed to calculate the approximate solution.     

Optimal and approximate boundary controls of an elastic body with quasistatic evolution of damage

In this paper, we study an optimal control problem for the mixed boundary value problem for an elastic body with quasistatic evolution of an internal damage variable. We suppose that the evolution of microscopic cracks and cavities responsible for the damage is described by a nonlinear parabolic equation. A density of surface traction p acting on a part of boundary of an elastic body Ω is taken as a boundary control. Because the initial boundary value problem of this type can exhibit the Lavrentieff phenomenon and non‐uniqueness of weak solutions, we deal with the solvability of this problem in the class of weak variational solutions. Using the convergence concept in variable spaces and following the direct method in calculus of variations, we prove the existence of optimal and approximate solutions to the optimal control problem under rather general assumptions on the quasistatic evolution of damage. Copyright © 2014 John Wiley & Sons, Ltd.     

Modeling a swimming fish with an initial boundary value problem: Unsteady maneuvers of an elastic plate with internal force generation

In order to model unsteady maneuvers in swimming fish, we develop an initial-boundary value problem for a fourth-order hyperbolic partial differential equation in which the fish's body is treated as an inhomogeneous elastic plate. The model is derived from the three-dimensional equations of elastic dynamics, and is essentially a simple variant of the classical Kirchhoff model for a dynamic plate. The model incorporates body forces generating moment to simulate muscle force generation in fish. The initial-boundary value problem is reduced to a beam model in one spatial dimension and formulated computationally using finite differences. Interaction with the surrounding water is represented by nonlinear viscous damping. Two example applications using simple but physically reasonable physiological parameters are presented and interpreted. One models the acceleration from rest to steady swimming, the other a rapid turn from rest.     

An expansion theorem for an eigenvalue problem of an arbitrary multiply connected domain with an eigenparameter in a general type of boundary conditions

The object of this paper is to establish an expansion theorem for a regular right-definite eigenvalue problem with an eigenvalue parameter which is contained in the Schrödinger partial differential equation and in a general type of boundary conditions on the boundary of an arbitrary multiply connected bounded domain inR ⁿ (n2). We associate with this problem an essentially self-adjoint operator in a suitably defined Hilbert space and then we develop an associated eigenfunction expansion theorem.     

Self‐adjoint Sturm–Liouville problems with an infinite number of boundary conditions

We study Sturm–Liouville (SL) problems on an infinite number of intervals, adjacent endpoints are linked by means of boundary conditions, and characterize the conditions which determine self‐adjoint operators in a Hilbert space which is the direct sum of the spaces for each interval. These conditions can be regular or singular, separated or coupled. Furthermore, the inner products of the summand spaces may be multiples of the usual inner products with different spaces having different multiples. We also extend the GKN Theorem to cover the infinite number of intervals theory with modified inner products and discuss the connection between our characterization and the classical one with the usual inner products. Our results include the finite number of intervals case.     

Locating the peaks of the least energy solutions to an ellyptic system with Neumann boundary conditions

We consider a system of the form in Ω with Neumann boundary condition on ∂Ω, where Ω is a smooth bounded domain in and f,g are power-type nonlinearities having superlinear and subcritical growth at infinity. We prove that the least energy solutions to such a system concentrate, as ε goes to zero, at a point of the boundary which maximizes the mean curvature of the boundary of Ω.     

On the existence of nonnegative solutions to an integral equation with applications to boundary value problems

By the use of a Guo-Krasnoselskii theorem in cones, existence of positive eigenvalues yielding nonnegative or positive solutions to an integral equation is studied. The results are applied to a variety of boundary value problems concerning ordinary differential equations.     

具超前变元的二阶微分方程三点边值问题的正解

§ 1　 IntroductionRecently,certain three-point boundary value problems for nonlinear ordinarydifferential equations have been studied by many authors[1— 6] .However,few papers havebeen published on the same problems for nonlinear functional differential equations.In thispaper,we are concerned with the following second order differential equation with anadvanced argumentu″(t) +λa(t) f(u(h(t) ) ) =0 ,t∈ (0 ,1 ) (1 .1 )with the three-point boundary conditionsu(0 ) =0 ,αu(η) =u(1 ) ,(1 .2 )…     

Self-adjoint realization of a class of third-order differential operators with an eigenparameter contained in the boundary conditions

The present paper deals with a class of third-order differential operators with eigenparameter dependent boundary conditions. Using operator theoretic formulation, the self-adjointness of this operator is proved, the properties of spectrum are investigated, its Green function and the resolvent operator are also obtained.