CONSTRUCTING QUANTUM OBSERVABLES AND SELF-ADJOINT EXTENSIONS OF SYMMETRIC OPERATORS. III. SELF-ADJOINT BOUNDARY CONDITIONS

This paper completes the review of the theory of self-adjoint extensions of symmetric operators for physicists as a basis for constructing quantum-mechanical observables. It contains a comparative presentation of the well-known methods and a newly proposed method for constructing ordinary self-adjoint differential operators associated with self-adjoint differential expressions in terms of self-adjoint boundary conditions. The new method has the advantage that it does not require explicitly evaluating deficient subspaces and deficiency indices (these latter are determined in passing) and that boundary conditions are of explicit character irrespective of the singularity of a differential expression. General assertions and constructions are illustrated by examples of well-known quantum-mechanical operators like momentum and Hamiltonian.

A Direct Numerov Sixth-order Numerical Scheme to Accurately Solve the Unidimensional Poisson Equation with Dirichlet Boundary Conditions

In this article, we present an analytical direct method, based on a Numerov three-point scheme, which is sixth order accurate and has a linear execution time on the grid dimension, to solve the discrete one-dimensional Poisson equation with Dirichlet boundary conditions. Our results should improve numerical codes used mainly in self-consistent calculations in solid state physics.

A complex network-based approach for boundary shape analysis

This paper introduces a novel methodology to shape boundary characterization, where a shape is modeled into a small-world complex network. It uses degree and joint degree measurements in a dynamic evolution network to compose a set of shape descriptors. The proposed shape characterization method has all efficient power of shape characterization, it is robust, noise tolerant, scale invariant and rotation invariant. A leaf plant classification experiment is presented on three image databases in order to evaluate the method and compare it with other descriptors in the literature (Fourier descriptors, Curvature, Zernike moments and multiscale fractal dimension). (C) 2008 Elsevier Ltd. All rights reserved.

The Sphene-Centered Ocellar Texture: An Effect of Grain-Supported Flow and Melt Migration in a Hyperdense Magma Mush

The sphene-centered ocellar texture consists of leucocratic ocelli with sphene (titanite) crystals at the center, enclosed in a biotite-rich matrix. This texture has been recognized worldwide in hybrid intermediate rocks. On the basis of structural, petrological, and geochronological data from selected outcrops of the Variscan Ribadelago pluton (NW Iberian Massif), we propose that the ocelli were formed by migration and accumulation of a residual melt through a plagioclase- and biotite-dominated crystalline framework. At the late stage of crystallization, the magma acted as a hyperdense suspension and reacted to the pressure gradient caused by the regional stress field, entering the domain of grain-supported flow. Microstructures reveal that aligned crystal domains arose in the crystal framework from the shearing and compaction of the crystal mush and behaved as magmatic microshears. Relative displacement of adjacent crystal clusters along these microshears corresponded to the onset of Reynolds dilatancy that generated an expansion of the crystal mush, involving melt migration and pore aperture. The mineralogy of the ocelli, dominated by andesine and sphene, represents the composition of the migrating melt. The chemistry of this late, Ti-rich melt stems from the incongruent melting of biotite. Magmatic sphene from the ocelli yields a U-Pb age of 317 +/- 1 Ma, which represents the final crystallization of the hybridized magmatic system. Moreover, this texture offers an opportunity to better understand the rheological behavior of highly crystallized magmas.

Controls of heavy minerals and grain size in a holocene regressive barrier (Ilha Comprida, southeastern Brazil)

Ilha Comprida is a regressive barrier island located in southeastern Brazil that was formed essentially by Quaternary sandy sediments. Ilha Comprida sediments were analyzed to assess heavy mineral indices and grain size variables. The spatial variation of heavy minerals and grain size was interpreted in terms of the present barrier dynamics and the barrier`s evolution since the Middle Holocene. These analyses allowed for the identification of the main factors and processes that control the variation of heavy minerals and grain size on the barrier. Rutile and zircon (RZi) and tourmaline and hornblende (THi) are significantly sensitive to provenance and exhibit the contributions of the Ribeira de Iguape River sediments, which reach the coast next to the northeastern end of Ilha Comprida. In addition to the influence of provenance, TZi responds mainly to hydraulic sorting processes. This agrees with a sediment transport pattern characterized by a divergence of two resultant net alongshore drifts southwest of the barrier. The sediments from the Ribeira de Iguape River reach the barrier directly through the river mouth and indirectly after temporary storage in the inner shelf. The combination of grain size and heavy mineral analyses is a reliable method for determining sediment transport patterns and provenance. (C) 2010 Elsevier Ltd. All rights reserved.

Grain size and heavy minerals of the Late Quaternary eolian sediments from the Imbituba-Jaguaruna coast, Southern Brazil: Depositional controls linked to relative sea-level changes

The stratigraphic subdivision and correlation of dune deposits is difficult, especially when age datings are not available. A better understanding of the controls on texture and composition of eolian sands is necessary to interpret ancient eolian sediments. The Imbituba-Jaguaruna coastal zone (Southern Brazil, 28 degrees-29 degrees S) stands out due to its four well-preserved Late Pleistocene (eolian generation 1) to Holocene eolian units (eolian generations 2, 3, and 4). In this study, we evaluate the grain-size and heavy-mineral characteristics of the Imbituba-Jaguartma eolian units through statistical analysis of hundreds of sediment samples. Grain-size parameters and heavy-mineral content allow us to distinguish the Pleistocene from the Holocene units. The grain size displays a pattern of fining and better sorting from generation 1 (older) to 4 (younger), whereas the content of mechanically stable (dense and hard) heavy minerals decreases from eolian generation 1 to 4. The variation in grain size and heavy-mineral content records shifts in the origin and balance (input versus output) of eolian sediment supply attributable mainly to relative sea-level changes. Dunefields submitted to relative sea-level lowstand conditions (eolian generation 1) are characterized by lower accumulation rates and intense post-depositional dissection by fluvial incision. Low accumulation rates favor deflation in the eolian system, which promotes concentration of denser and stable heavy minerals (increase of ZTR index) as well as coarsening of eolian sands. Dissection involves the selective removal of finer sediments and less dense heavy minerals to the coastal source area. Under a high rate of relative sea-level rise and transgression (eolian generation 2), coastal erosion prevents deflation through high input of sediments to the coastal eolian source. This condition favors dunefield growth. Coastal erosion feeds sand from local sources to the eolian system. including sands from previous dunefields (eolian generation 1) and from drowned incised valleys. Therefore, dunefields corresponding to transgressive phases inherit the grain-size and heavy-mineral characteristics of previous dunefields, leading to selective enrichment of finer sands and lighter minerals. Eolian generations 3 and 4 developed during a regressive-progradational phase (Holocene relative sea level highstand). The high rate of sediment supply during the highstand phase prevents deflation. The lack of coastal erosion favors sediment supply from distal sources (fluvial sediments rich in unstable heavy minerals). Thus, dunefields of transgressive and highstand systems tracts may be distinguished from dunefields of the lowstand systems tract through high rates of accumulation (low deflation) in the former. The sediment source of the transgressive dunefields (high input of previously deposited coastal sands) differs from that of the highstand dunefields (high input of fluvial distal sands). Based on this case study, we propose a general framework for the relation between relative sea level, sediment supply and the texture and mineralogy of eolian sediments deposited in siliciclastic wet coastal zones similar to the Imbituba-Jaguaruna coast (C) 2009 Elsevier B.V. All rights reserved.

Cascades of Hopf bifurcations from boundary delay

We consider a 1-dimensional reaction-diffusion equation with nonlinear boundary conditions of logistic type with delay. We deal with non-negative solutions and analyze the stability behavior of its unique positive equilibrium solution, which is given by the constant function u equivalent to 1. We show that if the delay is small, this equilibrium solution is asymptotically stable, similar as in the case without delay. We also show that, as the delay goes to infinity, this equilibrium becomes unstable and undergoes a cascade of Hopf bifurcations. The structure of this cascade will depend on the parameters appearing in the equation. This equation shows some dynamical behavior that differs from the case where the nonlinearity with delay is in the interior of the domain. (C) 2009 Elsevier Inc. All rights reserved.

Efficient solutions to robust, semi-implicit discretizations of the immersed boundary method

The immersed boundary method is a versatile tool for the investigation of flow-structure interaction. In a large number of applications, the immersed boundaries or structures are very stiff and strong tangential forces on these interfaces induce a well-known, severe time-step restriction for explicit discretizations. This excessive stability constraint can be removed with fully implicit or suitable semi-implicit schemes but at a seemingly prohibitive computational cost. While economical alternatives have been proposed recently for some special cases, there is a practical need for a computationally efficient approach that can be applied more broadly. In this context, we revisit a robust semi-implicit discretization introduced by Peskin in the late 1970s which has received renewed attention recently. This discretization, in which the spreading and interpolation operators are lagged. leads to a linear system of equations for the inter-face configuration at the future time, when the interfacial force is linear. However, this linear system is large and dense and thus it is challenging to streamline its solution. Moreover, while the same linear system or one of similar structure could potentially be used in Newton-type iterations, nonlinear and highly stiff immersed structures pose additional challenges to iterative methods. In this work, we address these problems and propose cost-effective computational strategies for solving Peskin`s lagged-operators type of discretization. We do this by first constructing a sufficiently accurate approximation to the system`s matrix and we obtain a rigorous estimate for this approximation. This matrix is expeditiously computed by using a combination of pre-calculated values and interpolation. The availability of a matrix allows for more efficient matrix-vector products and facilitates the design of effective iterative schemes. We propose efficient iterative approaches to deal with both linear and nonlinear interfacial forces and simple or complex immersed structures with tethered or untethered points. One of these iterative approaches employs a splitting in which we first solve a linear problem for the interfacial force and then we use a nonlinear iteration to find the interface configuration corresponding to this force. We demonstrate that the proposed approach is several orders of magnitude more efficient than the standard explicit method. In addition to considering the standard elliptical drop test case, we show both the robustness and efficacy of the proposed methodology with a 2D model of a heart valve. (C) 2009 Elsevier Inc. All rights reserved.

Boundary and internal conditions for adjoint fluid-flow problems

The ever-increasing robustness and reliability of flow-simulation methods have consolidated CFD as a major tool in virtually all branches of fluid mechanics. Traditionally, those methods have played a crucial role in the analysis of flow physics. In more recent years, though, the subject has broadened considerably, with the development of optimization and inverse design applications. Since then, the search for efficient ways to evaluate flow-sensitivity gradients has received the attention of numerous researchers. In this scenario, the adjoint method has emerged as, quite possibly, the most powerful tool for the job, which heightens the need for a clear understanding of its conceptual basis. Yet, some of its underlying aspects are still subject to debate in the literature, despite all the research that has been carried out on the method. Such is the case with the adjoint boundary and internal conditions, in particular. The present work aims to shed more light on that topic, with emphasis on the need for an internal shock condition. By following the path of previous authors, the quasi-1D Euler problem is used as a vehicle to explore those concepts. The results clearly indicate that the behavior of the adjoint solution through a shock wave ultimately depends upon the nature of the objective functional.

ON THE SUPERCRITICALITY OF THE FIRST HOPF BIFURCATION IN A DELAY BOUNDARY PROBLEM

The goal of this paper is to analyze the character of the first Hopf bifurcation (subcritical versus supercritical) that appears in a one-dimensional reaction-diffusion equation with nonlinear boundary conditions of logistic type with delay. We showed in the previous work [Arrieta et al., 2010] that if the delay is small, the unique non-negative equilibrium solution is asymptotically stable. We also showed that, as the delay increases and crosses certain critical value, this equilibrium becomes unstable and undergoes a Hopf bifurcation. This bifurcation is the first one of a cascade occurring as the delay goes to infinity. The structure of this cascade will depend on the parameters appearing in the equation. In this paper, we show that the first bifurcation that occurs is supercritical, that is, when the parameter is bigger than the delay bifurcation value, stable periodic orbits branch off from the constant equilibrium.

Reaction-diffusion systems coupled at the boundary and the Morse-Smale property

We study an one-dimensional nonlinear reaction-diffusion system coupled on the boundary. Such system comes from modeling problems of temperature distribution on two bars of same length, jointed together, with different diffusion coefficients. We prove the transversality property of unstable and stable manifolds assuming all equilibrium points are hyperbolic. To this end, we write the system as an equation with noncontinuous diffusion coefficient. We then study the nonincreasing property of the number of zeros of a linearized nonautonomous equation as well as the Sturm-Liouville properties of the solutions of a linear elliptic problem. (C) 2008 Elsevier Inc. All rights reserved.

Existence of orthogonal geodesic chords on Riemannian manifolds with concave boundary and homeomorphic to the N-dimensional disk

In this paper we give a proof of the existence of an orthogonal geodesic chord on a Riemannian manifold homeomorphic to a closed disk and with concave boundary. This kind of study is motivated by the link (proved in Giambo et al. (2005) [8]) of the multiplicity problem with the famous Seifert conjecture (formulated in Seifert (1948) [1]) about multiple brake orbits for a class of Hamiltonian systems at a fixed energy level. (C) 2010 Elsevier Ltd. All rights reserved.

A FREE BOUNDARY ISOPERIMETRIC PROBLEM IN HYPERBOLIC 3-SPACE BETWEEN PARALLEL HOROSPHERES

We investigate the isoperimetric problem of finding the regions of prescribed volume with minimal boundary area between two parallel horospheres in hyperbolic 3-space (the part of the boundary contained in the horospheres is not included). We reduce the problem to the study of rotationally invariant regions and obtain the possible isoperimetric solutions by studying the behavior of the profile curves of the rotational surfaces with constant mean curvature in hyperbolic 3-space. We also classify all the connected compact rotational surfaces M of constant mean curvature that are contained in the region between two horospheres, have boundary partial derivative M either empty or lying on the horospheres, and meet the horospheres perpendicularly along their boundary.

On the Multiplicity of Orthogonal Geodesics in Riemannian Manifold With Concave Boundary. Applications to Brake Orbits and Homoclinics

Let (M, g) be a complete Riemannian Manifold, Omega subset of M an open subset whose closure is diffeomorphic to an annulus. If partial derivative Omega is smooth and it satisfies a strong concavity assumption, then it is possible to prove that there are at least two geometrically distinct geodesics in (Omega) over bar = Omega boolean OR partial derivative Omega starting orthogonally to one connected component of partial derivative Omega and arriving orthogonally onto the other one. The results given in [6] allow to obtain a proof of the existence of two distinct homoclinic orbits for an autonomous Lagrangian system emanating from a nondegenerate maximum point of the potential energy, and a proof of the existence of two distinct brake orbits for a. class of Hamiltonian systems. Under a further symmetry assumption, it is possible to show the existence of at least dim(M) pairs of geometrically distinct geodesics as above, brake orbits and homoclinics.

GENERICITY OF NONDEGENERATE GEODESICS WITH GENERAL BOUNDARY CONDITIONS

Let M be a possibly noncompact manifold. We prove, generically in the C(k)-topology (2 <= k <= infinity), that semi-Riemannian metrics of a given index on M do not possess any degenerate geodesics satisfying suitable boundary conditions. This extends a result of L. Biliotti, M. A. Javaloyes and P. Piccione [6] for geodesics with fixed endpoints to the case where endpoints lie on a compact submanifold P subset of M x M that satisfies an admissibility condition. Such condition holds, for example, when P is transversal to the diagonal Delta subset of M x M. Further aspects of these boundary conditions are discussed and general conditions under which metrics without degenerate geodesics are C(k)-generic are given.