Full-dimensional quantum calculations of ground-state tunneling splitting of malonaldehyde using an accurate ab initio potential energy surface

Quantum calculations of the ground vibrational state tunneling splitting of H-atom and D-atom transfer in malonaldehyde are performed on a full-dimensional ab initio potential energy surface (PES). The PES is a fit to 11 147 near basis-set-limit frozen-core CCSD(T) electronic energies. This surface properly describes the invariance of the potential with respect to all permutations of identical atoms. The saddle-point barrier for the H-atom transfer on the PES is 4.1 kcal/mol, in excellent agreement with the reported ab initio value. Model one-dimensional and "exact" full-dimensional calculations of the splitting for H- and D-atom transfer are done using this PES. The tunneling splittings in full dimensionality are calculated using the unbiased "fixed-node" diffusion Monte Carlo (DMC) method in Cartesian and saddle-point normal coordinates. The ground-state tunneling splitting is found to be 21.6 cm(-1) in Cartesian coordinates and 22.6 cm(-1) in normal coordinates, with an uncertainty of 2-3 cm(-1). This splitting is also calculated based on a model which makes use of the exact single-well zero-point energy (ZPE) obtained with the MULTIMODE code and DMC ZPE and this calculation gives a tunneling splitting of 21-22 cm(-1). The corresponding computed splittings for the D-atom transfer are 3.0, 3.1, and 2-3 cm(-1). These calculated tunneling splittings agree with each other to within less than the standard uncertainties obtained with the DMC method used, which are between 2 and 3 cm(-1), and agree well with the experimental values of 21.6 and 2.9 cm(-1) for the H and D transfer, respectively. (C) 2008 American Institute of Physics.

A uniformly valid far field asymptotic expansion of the green function for two-dimensional propagation above a homogeneous impedance plane

A generalized asymptotic expansion in the far field for the problem of cylindrical wave reflection at a homogeneous impedance plane is derived. The expansion is shown to be uniformly valid over all angles of incidence and values of surface impedance, including the limiting cases of zero and infinite impedance. The technique used is a rigorous application of the modified steepest descent method of Ot

On the existence of two-dimensional Euler flows satisfying energy-Casimir stability criteria

The energy-Casimir stability method, also known as the Arnold stability method, has been widely used in fluid dynamical applications to derive sufficient conditions for nonlinear stability. The most commonly studied system is two-dimensional Euler flow. It is shown that the set of two-dimensional Euler flows satisfying the energy-Casimir stability criteria is empty for two important cases: (i) domains having the topology of the sphere, and (ii) simply-connected bounded domains with zero net vorticity. The results apply to both the first and the second of Arnold’s stability theorems. In the spirit of Andrews’ theorem, this puts a further limitation on the applicability of the method. © 2000 American Institute of Physics.

On Arnol'd's second nonlinear stability theorem for two-dimensional quasi-geostrophic flow

Arnol'd's second hydrodynamical stability theorem, proven originally for the two-dimensional Euler equations, can establish nonlinear stability of steady flows that are maxima of a suitably chosen energy-Casimir invariant. The usual derivations of this theorem require an assumption of zero disturbance circulation. In the present work an analogue of Arnol'd's second theorem is developed in the more general case of two-dimensional quasi-geostrophic flow, with the important feature that the disturbances are allowed to have non-zero circulation. New nonlinear stability criteria are derived, and explicit bounds are obtained on both the disturbance energy and potential enstrophy which are expressed in terms of the initial disturbance fields. While Arnol'd's stability method relies on the second variation of the energy-Casimir invariant being sign-definite, the new criteria can be applied to cases where the second variation is sign-indefinite because of the disturbance circulations. A version of Andrews' theorem is also established for this problem.

Eigenvalues of a one-dimensional Dirac operator pencil

We study the spectrum of a one-dimensional Dirac operator pencil, with a coupling constant in front of the potential considered as the spectral parameter. Motivated by recent investigations of graphene waveguides, we focus on the values of the coupling constant for which the kernel of the Dirac operator contains a square integrable function. In physics literature such a function is called a confined zero mode. Several results on the asymptotic distribution of coupling constants giving rise to zero modes are obtained. In particular, we show that this distribution depends in a subtle way on the sign variation and the presence of gaps in the potential. Surprisingly, it also depends on the arithmetic properties of certain quantities determined by the potential. We further observe that variable sign potentials may produce complex eigenvalues of the operator pencil. Some examples and numerical calculations illustrating these phenomena are presented.

Finite-dimensional non-associative algebras and codimension growth

Let A be a (non-necessarily associative) finite-dimensional algebra over a field of characteristic zero. A quantitative estimate of the polynomial identities satisfied by A is achieved through the study of the asymptotics of the sequence of codimensions of A. It is well known that for such an algebra this sequence is exponentially bounded. Here we capture the exponential rate of growth of the sequence of codimensions for several classes of algebras including simple algebras with a special non-degenerate form, finite-dimensional Jordan or alternative algebras and many more. In all cases such rate of growth is integer and is explicitly related to the dimension of a subalgebra of A. One of the main tools of independent interest is the construction in the free non-associative algebra of multialternating polynomials satisfying special properties. (C) 2010 Elsevier Inc. All rights reserved.

The Gauss-Kronecker curvature of minimal hypersurfaces in four-dimensional space forms

LetQ(4)( c) be a four-dimensional space form of constant curvature c. In this paper we show that the infimum of the absolute value of the Gauss-Kronecker curvature of a complete minimal hypersurface in Q(4)(c), c <= 0, whose Ricci curvature is bounded from below, is equal to zero. Further, we study the connected minimal hypersurfaces M(3) of a space form Q(4)( c) with constant Gauss-Kronecker curvature K. For the case c <= 0, we prove, by a local argument, that if K is constant, then K must be equal to zero. We also present a classification of complete minimal hypersurfaces of Q(4)( c) with K constant.

Modeling of interfaces in two-dimensional problems using solid finite elements with high aspect ratio

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Numerical solution of the two-dimensional Gross-Pitaevskii equation for trapped interacting atoms

We present a numerical scheme for solving the time-independent nonlinear Gross-Pitaevskii equation in two dimensions describing the Bose-Einstein condensate of trapped interacting neutral atoms at zero temperature. The trap potential is taken to be of the harmonic-oscillator type and the interaction both attractive and repulsive. The Gross-Pitaevskii equation is numerically integrated consistent with the correct boundary conditions at the origin and in the asymptotic region. Rapid convergence is obtained in all cases studied. In the attractive case there is a limit Co the maximum number of atoms in the condensate. (C) 2000 Published by Elsevier B.V. B.V. All rights reserved.

First study of the radiation-amplitude zero in W gamma production and limits on anomalous WW gamma couplings at root s=1.96 TeV

We present results from a study of p (p) over bar -> W gamma+X events utilizing data corresponding to 0.7 fb(-1) of integrated luminosity at root s = 1.96 TeV collected by the D0 detector at the Fermilab Tevatron Collider. We set limits on anomalous WW gamma couplings at the 95% C.L. The one-dimensional 95% C.L. limits are 0.49

SEMICLASSICAL THEORY OF MAGNETIZATION FOR A 2-DIMENSIONAL NONINTERACTING ELECTRON-GAS

We compute the semiclassical magnetization and susceptibility of non-interacting electrons, confined by a smooth two-dimensional potential and subjected to a uniform perpendicular magnetic field, in the general case when their classical motion is chaotic. It is demonstrated that the magnetization per particle m(B) is directly related to the staircase function N(E), which counts the single-particle levels up to energy E. Using Gutzwiller's trace formula for N, we derive a semiclassical expression for m. Our results show that the magnetization has a non-zero average, which arises from quantum corrections to the leading-order Weyl approximation to the mean staircase and which is independent of whether the classical motion is chaotic or not. Fluctuations about the average are due to classical periodic orbits and do represent a signature of chaos. This behaviour is confirmed by numerical computations for a specific system.

Spin, charge, and orbital correlations in the one-dimensional t(2g)-orbital Hubbard model

We present the zero-temperature phase diagram of the one-dimensional t(2g)-orbital Hubbard model, obtained using the density-matrix renormalization group and Lanczos techniques. Emphasis is given to the case of the electron density n=5 corresponding to five electrons per site, while several other cases for electron densities between n=3 and 6 are also studied. At n=5, our results indicate a first-order transition between a paramagnetic (PM) insulator phase, with power-law slowly decaying correlations, and a fully polarized ferromagnetic (FM) state by tuning the Hund's coupling. The results also suggest a transition from the n=5 PM insulator phase to a metallic regime by changing the electron density, either via hole or electron doping. The behavior of the spin, charge, and orbital correlation functions in the FM and PM states are also described in the text and discussed. The robustness of these two states against varying parameters suggests that they may be of relevance in quasi-one-dimensional Co-oxide materials, or even in higher dimensional cobaltite systems as well.

Testing the resolving power of two-dimensional K+K+ interferometry

We show results from an analysis performed to test the resolving power of a two-dimensional χ2 method proposed previously when applied to the case of kaon interferometry, where no significant contribution from long-lived resonances is expected. For that purpose, use is made of the preliminary E859 K+K+ interferometry data from Si+Au collisions at 14.6/4 GeV/c. Although less sensitivity is achieved in the present case, this analysis seems to favor scenarios with no resonance formation at the AGS energy range. The possible compatibility of data with zero decoupling proper time interval, conjectured by the three-dimensional experimental analysis, is also investigated and is ruled out when considering more realistic dynamical models with expanding sources. Furthermore, these results strongly emphasize that the static Gaussian parametrization cannot be trusted under more realistic conditions, leading to a distorted or even wrong interpretation of the source parameters.

First study of the radiation-amplitude zero in W gamma production and limits on anomalous WW gamma couplings at root s=1.96 TeV

We present results from a study of pp̄→Wγ+X events utilizing data corresponding to 0.7fb-1 of integrated luminosity at s=1.96TeV collected by the D0 detector at the Fermilab Tevatron Collider. We set limits on anomalous WWγ couplings at the 95% C.L. The one-dimensional 95% C.L. limits are 0.49<κγ<1.51 and -0.12<λγ<0.13. We make the first study of the charge-signed rapidity difference between the lepton and the photon and find it to be indicative of the standard model radiation-amplitude zero in the Wγ system. © 2008 The American Physical Society.

Critical currents and melting temperature of a two-dimensional vortex lattice with periodic pinning

The critical current and melting temperature of a vortex system are analyzed. Calculations are made for a two-dimensional film at finite temperature with two kinds of periodic pinning: hexagonal and Kagomé. A transport current parallel and perpendicular to the main axis of the pinning arrays is applied and molecular dynamics simulations are used to calculate the vortex velocities to obtain the critical currents. The structure factor and displacements of vortices at zero transport current are used to obtain the melting temperature for both pinning arrays. The critical currents are higher for the hexagonal pinning lattice and anisotropic for both pinning arrays. This anisotropy is stronger with temperature for the hexagonal array. For the Kagomé pinning lattice, our analysis shows a multi stage phase melting; that is, as we increase the temperature, each different dynamic phase melts before reaching the melting temperature. Both the melting temperature and critical currents are larger for the hexagonal lattice, indicating the role for the interstitial vortices in decreasing the pinning strength. © 2012 Springer Science+Business Media New York.