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.

An efficient numerical scheme for the two-dimensional shallow water equations using arithmetic averaging

A finite difference scheme based on flux difference splitting is presented for the solution of the two-dimensional shallow water equations of ideal fluid flow. A linearised problem, analogous to that of Riemann for gas dynamics is defined, and a scheme, based on numerical characteristic decomposition is presented for obtaining approximate solutions to the linearised problem, and incorporates the technique of operator splitting. An average of the flow variables across the interface between cells is required, and this average is chosen to be the arithmetic mean for computational efficiency leading to arithmetic averaging. This is in contrast to usual ‘square root’ averages found in this type of Riemann solver, where the computational expense can be prohibitive. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second order scheme which avoids nonphysical, spurious oscillations. An extension to the two-dimensional equations with source terms is included. The scheme is applied to the one-dimensional problems of a breaking dam and reflection of a bore, and in each case the approximate solution is compared to the exact solution of ideal fluid flow. The scheme is also applied to a problem of stationary bore generation in a channel of variable cross-section. Finally, the scheme is applied to two other dam-break problems, this time in two dimensions with one having cylindrical symmetry. Each approximate solution compares well with those given by other authors.

Four-dimensional variational data assimilation for high resolution nested models

Four-dimensional variational data assimilation (4D-Var) is used in environmental prediction to estimate the state of a system from measurements. When 4D-Var is applied in the context of high resolution nested models, problems may arise in the representation of spatial scales longer than the domain of the model. In this paper we study how well 4D-Var is able to estimate the whole range of spatial scales present in one-way nested models. Using a model of the one-dimensional advection–diffusion equation we show that small spatial scales that are observed can be captured by a 4D-Var assimilation, but that information in the larger scales may be degraded. We propose a modification to 4D-Var which allows a better representation of these larger scales.

A model of the three-dimensional evolution of Arctic melt ponds on first-year and multiyear sea ice

During winter the ocean surface in polar regions freezes over to form sea ice. In the summer the upper layers of sea ice and snow melts producing meltwater that accumulates in Arctic melt ponds on the surface of sea ice. An accurate estimate of the fraction of the sea ice surface covered in melt ponds is essential for a realistic estimate of the albedo for global climate models. We present a melt-pond–sea-ice model that simulates the three-dimensional evolution of melt ponds on an Arctic sea ice surface. The advancements of this model compared to previous models are the inclusion of snow topography; meltwater transport rates are calculated from hydraulic gradients and ice permeability; and the incorporation of a detailed one-dimensional, thermodynamic radiative balance. Results of model runs simulating first-year and multiyear sea ice are presented. Model results show good agreement with observations, with duration of pond coverage, pond area, and ice ablation comparing well for both the first-year ice and multiyear ice cases. We investigate the sensitivity of the melt pond cover to changes in ice topography, snow topography, and vertical ice permeability. Snow was found to have an important impact mainly at the start of the melt season, whereas initial ice topography strongly controlled pond size and pond fraction throughout the melt season. A reduction in ice permeability allowed surface flooding of relatively flat, first-year ice but had little impact on the pond coverage of rougher, multiyear ice. We discuss our results, including model shortcomings and areas of experimental uncertainty.

A three-dimensional open-framework indium selenide: [C7H10N][In9Se14]

An open-framework indium selenide, [C7H10N][In9Se14], has been prepared under solvothermal conditions in the presence of 3,5-dimethylpyridine, and characterized by single crystal diffraction, thermogravimetry, elemental analysis, FTIR spectroscopy and UV-Vis diffuse reflectance. The crystal structure of [C7H10N][In9Se14] contains an unusual building unit, in which corner-linked and edge-linked InSe45- tetrahedra coexist. The presence of one-dimensional circular channels, of ca. 6 Å diameter, results in approximately 25% of solvent accessible void space.

STABILITY ANALYSIS WITH APPLICATIONS OF A TWO-DIMENSIONAL DYNAMICAL SYSTEM ARISING FROM A STOCHASTIC MODEL FOR AN ASSET MARKET

We analyze the stability properties of equilibrium solutions and periodicity of orbits in a two-dimensional dynamical system whose orbits mimic the evolution of the price of an asset and the excess demand for that asset. The construction of the system is grounded upon a heterogeneous interacting agent model for a single risky asset market. An advantage of this construction procedure is that the resulting dynamical system becomes a macroscopic market model which mirrors the market quantities and qualities that would typically be taken into account solely at the microscopic level of modeling. The system`s parameters correspond to: (a) the proportion of speculators in a market; (b) the traders` speculative trend; (c) the degree of heterogeneity of idiosyncratic evaluations of the market agents with respect to the asset`s fundamental value; and (d) the strength of the feedback of the population excess demand on the asset price update increment. This correspondence allows us to employ our results in order to infer plausible causes for the emergence of price and demand fluctuations in a real asset market. The employment of dynamical systems for studying evolution of stochastic models of socio-economic phenomena is quite usual in the area of heterogeneous interacting agent models. However, in the vast majority of the cases present in the literature, these dynamical systems are one-dimensional. Our work is among the few in the area that construct and study analytically a two-dimensional dynamical system and apply it for explanation of socio-economic phenomena.

Thomas-Fermi ground state of dipolar fermions in a circular storage ring

Recent developments in the field of ultracold gases has led to the production of degenerate samples of polar molecules. These have large static electric-dipole moments, which in turn causes the molecules to interact strongly. We investigate the interaction of polar particles in waveguide geometries subject to an applied polarizing field. For circular waveguides, tilting the direction of the polarizing field creates a periodic inhomogeneity of the interparticle interaction. We explore the consequences of geometry and interaction for stability of the ground state within the Thomas-Fermi model. Certain combinations of tilt angles and interaction strengths are found to preclude the existence of a stable Thomas-Fermi ground state. The system is shown to exhibit different behavior for quasi-one-dimensional and three-dimensional trapping geometries.

Spectral element model updating for damage identification using clonal selection algorithm

A spectral element model updating procedure is presented to identify damage in a structure using Guided wave propagation results. Two damage spectral elements (DSE1 and DSE2) are developed to model the local (cracks in reinforcement bar) and global (debonding between reinforcement bar and concrete) damage in one-dimensional homogeneous and composite waveguide, respectively. Transfer matrix method is adopted to assemble the stiffness matrix of multiple spectral elements. In order to solve the inverse problem, clonal selection algorithm is used for the optimization calculations. Two displacement-based functions and two frequency-based functions are used as objective functions in this study. Numerical simulations of wave propagation in a bare steel bar and in a reinforcement bar without and with various assumed damage scenarios are carried out. Numerically simulated data are then used to identify local and global damage of the steel rebar and the concrete-steel interface using the proposed method. Results show that local damage is easy to be identified by using any considered objective function with the proposed method while only using the wavelet energy-based objective function gives reliable identification of global damage. The method is then extended to identify multiple damages in a structure. To further verify the proposed method, experiments of wave propagation in a rectangular steel bar before and after damage are conducted. The proposed method is used to update the structural model for damage identification. The results demonstrate the capability of the proposed method in identifying cracks in steel bars based on measured wave propagation data.

Peak picking and the assessment of separation performance in two-dimensional high performance liquid chromatography

An algorithm was developed for 2DHPLC that automated the process of peak recognition, measuring their retention times, and then subsequently plotting the information in a two-dimensional retention plane. Following the recognition of peaks, the software then performed a series of statistical assessments of the separation performance, measuring for example, correlation between dimensions, peak capacity and the percentage of usage of the separation space. Peak recognition was achieved by interpreting the first and second derivatives of each respective one-dimensional chromatogram to determine the 1D retention times of each solute and then compiling these retention times for each respective fraction ‘cut’. Due to the nature of comprehensive 2DHPLC adjacent cut fractions may contain peaks common to more than one cut fraction. The algorithm determined which components were common in adjacent cuts and subsequently calculated the peak maximum profile by interpolating the space between adjacent peaks. This algorithm was applied to the analysis of a two-dimensional separation of an apple flesh extract separated in a first dimension comprising a cyano stationary phase and an aqueous/THF mobile phase as the first dimension and a second dimension comprising C18-Hydro with an aqueous/MeOH mobile phase. A total of 187 peaks were detected.

Improving peak shapes with counter gradients in two-dimensional high performance liquid chromatography

To achieve the greatest peak capacity in two-dimensional high performance liquid chromatography (2D-HPLC) a gradient should be operated in both separation dimensions. However, it is known that when an injection solvent that is stronger than the initial mobile phase composition is deleterious to peak performance, thus causing problems when cutting a portion from one gradient into another. This was overcome when coupling hydrophilic interaction with reversed phase chromatography by introducing a counter gradient that changed the solvent strength of the second dimension injection. It was found that an injection solvent composition of 20% acetonitrile in water gave acceptable results in one-dimensional simulations with an initial composition of 5% acetonitrile. When this was transferred to a 2D-HPLC separation of standards it was found that a marked improvement in peak shape was gained for the moderately retained analytes (phenol and dimethyl phthalate), some improvement for the weakly retained caffeine and very little change for the strongly retained n-propylbenzene and anthracene which already displayed good chromatographic profiles. This effect was transferred when applied to a 2D-HPLC separation of a coffee extract where the indecipherable retention profile was transformed to a successful application multidimensional chromatography with peaks occupying 71% of the separation space according to the geometric approach to factor analysis.

Bloch-Kohn and Wannier-Kohn functions in one dimension

Bloch and Wannier functions of the Kohn type for a quite general one-dimensional Hamiltonian with inversion symmetry are studied. Important clarifications on null minigaps and the symmetry of those functions are given, with emphasis on the Kronig-Penney model. The lack of a general selection rule on the miniband index for optical transitions between edge states in semiconductor superlattices is discussed. A direct method for the calculation of Wannier-Kohn functions is presented.

Energy-dependent point interactions in one dimension

We consider a new type of point interaction in one-dimensional quantum mechanics. It is characterized by a boundary condition at the origin that involves the second and/or higher order derivatives of the wavefunction. The interaction is effectively energy dependent. It leads to a unitary S-matrix for the transmission-reflection problem. The energy dependence of the interaction can be chosen such that any given unitary S-matrix (or the transmission and reflection coefficients) can be reproduced at all energies. Generalization of the results to coupled-channel cases is discussed.

Gap solitons in a dipolar Bose-Einstein condensate on a three-dimensional optical lattice

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Dipolar Bose-Einstein condensate soliton on a two-dimensional optical lattice

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

SURFACE PERTURBATIONS OF A SHALLOW VISCOUS-FLUID HEATED FROM BELOW AND THE (2+1)-DIMENSIONAL BURGERS-EQUATION

The (2 + 1)-dimensional Burgers equation is obtained as the equation of motion governing the surface perturbations of a shallow viscous fluid heated from below, provided the Rayleigh number of the system satisfies the condition R not-equal 30. A solution to this equation is explicitly exhibited and it is argued that it describes the nonlinear evolution of a nearly one-dimensional kink.