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.

On the use of a time sequence of surface pressures in four-dimensional data assimilation

The possibility of using a time sequence of surface pressure observations in four-dimensional data assimilation is being investigated. It is shown that a linear multilevel quasi-geostrophic model can be updated successfully with surface data alone, provided the number of time levels are at least as many as the number of vertical levels. It is further demonstrated that current statistical analysis procedures are very inefficient to assimilate surface observations, and it is shown by numerical experiments that the vertical interpolation must be carried out using the structure of the most dominating baroclinic mode in order to obtain a satisfactory updating. Different possible ways towards finding a practical solution are being discussed.

Non-ergodicity of inviscid two-dimensional flow on a beta-plane and on the surface of a rotating sphere

It is shown that, for a sufficiently large value of β, two-dimensional flow on a doubly-periodic beta-plane cannot be ergodic (phase-space filling) on the phase-space surface of constant energy and enstrophy. A corresponding result holds for flow on the surface of a rotating sphere, for a sufficiently rapid rotation rate Ω. This implies that the higher-order, non-quadratic invariants are exerting a significant influence on the statistical evolution of the flow. The proof relies on the existence of a finite-amplitude Liapunov stability theorem for zonally symmetric basic states with a non-vanishing absolute-vorticity gradient. When the domain size is much larger than the size of a typical eddy, then a sufficient condition for non-ergodicity is that the wave steepness ε < 1, where ε = 2[surd radical]2Z/βU in the planar case and $\epsilon = 2^{\frac{1}{4}} a^{\frac{5}{2}}Z^{\frac{7}{4}}/\Omega U^{\frac{5}{2}}$ in the spherical case, and where Z is the enstrophy, U the r.m.s. velocity, and a the radius of the sphere. This result may help to explain why numerical simulations of unforced beta-plane turbulence (in which ε decreases in time) seem to evolve into a non-ergodic regime at large scales.

Humans use predictive kinematic models to calibrate visual cues to three-dimensional surface slant

When the sensory consequences of an action are systematically altered our brain can recalibrate the mappings between sensory cues and properties of our environment. This recalibration can be driven by both cue conflicts and altered sensory statistics, but neither mechanism offers a way for cues to be calibrated so they provide accurate information about the world, as sensory cues carry no information as to their own accuracy. Here, we explored whether sensory predictions based on internal physical models could be used to accurately calibrate visual cues to 3D surface slant. Human observers played a 3D kinematic game in which they adjusted the slant of a surface so that a moving ball would bounce off the surface and through a target hoop. In one group, the ball’s bounce was manipulated so that the surface behaved as if it had a different slant to that signaled by visual cues. With experience of this altered bounce, observers recalibrated their perception of slant so that it was more consistent with the assumed laws of kinematics and physical behavior of the surface. In another group, making the ball spin in a way that could physically explain its altered bounce eliminated this pattern of recalibration. Importantly, both groups adjusted their behavior in the kinematic game in the same way, experienced the same set of slants and were not presented with low-level cue conflicts that could drive the recalibration. We conclude that observers use predictive kinematic models to accurately calibrate visual cues to 3D properties of world.

SURFACE ALFVEN WAVE DAMPING IN A THREE-DIMENSIONAL SIMULATION OF THE SOLAR WIND

Here we investigate the contribution of surface Alfven wave damping to the heating of the solar wind in minima conditions. These waves are present in the regions of strong inhomogeneities in density or magnetic field (e.g., the border between open and closed magnetic field lines). Using a three-dimensional (3D) magnetohydrodynamics (MHD) model, we calculate the surface Alfven wave damping contribution between 1 and 4 R(circle dot) (solar radii), the region of interest for both acceleration and coronal heating. We consider waves with frequencies lower than those that are damped in the chromosphere and on the order of those dominating the heliosphere: 3 x 10(-6) to 10(-1) Hz. In the region between open and closed field lines, within a few R(circle dot) of the surface, no other major source of damping has been suggested for the low frequency waves we consider here. This work is the first to study surface Alfven waves in a 3D environment without assuming a priori a geometry of field lines or magnetic and density profiles. We demonstrate that projection effects from the plane of the sky to 3D are significant in the calculation of field line expansion. We determine that waves with frequencies >2.8 x 10(-4) Hz are damped between 1 and 4 R(circle dot). In quiet-Sun regions, surface Alfven waves are damped at further distances compared to active regions, thus carrying additional wave energy into the corona. We compare the surface Alfven wave contribution to the heating by a variable polytropic index and find it as an order of magnitude larger than needed for quiet-Sun regions. For active regions, the contribution to the heating is 20%. As it has been argued that a variable gamma acts as turbulence, our results indicate that surface Alfven wave damping is comparable to turbulence in the lower corona. This damping mechanism should be included self-consistently as an energy driver for the wind in global MHD models.

Suppressing Fermi acceleration in a two-dimensional non-integrable time-dependent oval-shaped billiard with inelastic collisions

Some dynamical properties of a classical particle confined inside a closed region with an oval-shaped boundary are studied. We have considered both the static and time-dependent boundaries. For the static case, the condition that destroys the invariant spanning curves in the phase space was obtained. For the time-dependent perturbation, two situations were considered: (i) non-dissipative and (ii) dissipative. For the non-dissipative case, our results show that Fermi acceleration is observed. When dissipation, via inelastic collisions, is introduced Fermi acceleration is suppressed. The behaviour of the average velocity for both the dissipative as well as the non-dissipative dynamics is described using the scaling approach. (C) 2009 Elsevier B.V. All rights reserved.

Boundary crisis and suppression of Fermi acceleration in a dissipative two-dimensional non-integrable time-dependent billiard

Some dynamical properties for a dissipative time-dependent oval-shaped billiard are studied. The system is described in terms of a four-dimensional nonlinear mapping. Dissipation is introduced via inelastic collisions of the particle with the boundary, thus implying that the particle has a fractional loss of energy upon collision. The dissipation causes profound modifications in the dynamics of the particle as well as in the phase space of the non-dissipative system. In particular, inelastic collisions can be assumed as an efficient mechanism to suppress Fermi acceleration of the particle. The dissipation also creates attractors in the system, including chaotic. We show that a slightly modification of the intensity of the damping coefficient yields a drastic and sudden destruction of the chaotic attractor, thus leading the system to experience a boundary crisis. We have characterized such a boundary crisis via a collision of the chaotic attractor with its own basin of attraction and confirmed that inelastic collisions do indeed suppress Fermi acceleration in two-dimensional time-dependent billiards. (C) 2010 Elsevier B.V. All rights reserved.

Suppressing Fermi acceleration in two-dimensional driven billiards

We consider a dissipative oval-like shaped billiard with a periodically moving boundary. The dissipation considered is proportional to a power of the velocity V of the particle. The three specific types of power laws used are: (i) F proportional to-V; (ii) F proportional to-V-2 and (iii) F proportional to-V-delta with 1 < delta < 2. In the course of the dynamics of the particle, if a large initial velocity is considered, case (i) shows that the decay of the particle's velocity is a linear function of the number of collisions with the boundary. For case (ii), an exponential decay is observed, and for 1 < delta < 2, an powerlike decay is observed. Scaling laws were used to characterize a phase transition from limited to unlimited energy gain for cases (ii) and (iii). The critical exponents obtained for the phase transition in the case (ii) are the same as those obtained for the dissipative bouncer model. Therefore near this phase transition, these two rather different models belong to the same class of universality. For all types of dissipation, the results obtained allow us to conclude that suppression of the unlimited energy growth is indeed observed.

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.

One-dimensional superfluid Bose-Fermi mixture: Mixing, demixing, and bright solitons

We study an ultracold and dilute superfluid Bose-Fermi mixture confined in a strictly one-dimensional (1D) atomic waveguide by using a set of coupled nonlinear mean-field equations obtained from the Lieb-Liniger energy density for bosons and the Gaudin-Yang energy density for fermions. We consider a finite Bose-Fermi interatomic strength gbf and both periodic and open boundary conditions. We find that with periodic boundary conditions-i.e., in a quasi-1D ring-a uniform Bose-Fermi mixture is stable only with a large fermionic density. We predict that at small fermionic densities the ground state of the system displays demixing if gbf >0 and may become a localized Bose-Fermi bright soliton for gbf <0. Finally, we show, using variational and numerical solutions of the mean-field equations, that with open boundary conditions-i.e., in a quasi-1D cylinder-the Bose-Fermi bright soliton is the unique ground state of the system with a finite number of particles, which could exhibit a partial mixing-demixing transition. In this case the bright solitons are demonstrated to be dynamically stable. The experimental realization of these Bose-Fermi bright solitons seems possible with present setups. © 2007 The American Physical Society.

Faraday waves in quasi-one-dimensional superfluid Fermi-Bose mixtures

The generation of Faraday waves in superfluid Fermi-Bose mixtures in elongated traps is investigated. The generation of waves is achieved by periodically changing a parameter of the system in time. Two types of modulations of parameters are considered: a variation of the fermion-boson scattering length and the boson-boson scattering length. We predict the properties of the generated Faraday patterns and study the parameter regions where they can be excited. © 2013 American Physical Society.

Three-dimensional transient complex free surface flows: Numerical simulation of XPP fluid

In this paper we present a finite difference MAC-type approach for solving three-dimensional viscoelastic incompressible free surface flows governed by the eXtended Pom-Pom (XPP) model, considering a wide range of parameters. The numerical formulation presented in this work is an extension to three-dimensions of our implicit technique [Journal of Non-Newtonian Fluid Mechanics 166 (2011) 165-179] for solving two-dimensional viscoelastic free surface flows. To enhance the stability of the numerical method, we employ a combination of the projection method with an implicit technique for treating the pressure on the free surfaces. The differential constitutive equation of the fluid is solved using a second-order Runge-Kutta scheme. The numerical technique is validated by performing a mesh refinement study on a pipe flow, and the numerical results presented include the simulation of two complex viscoelastic free surface flows: extrudate-swell problem and jet buckling phenomenon. © 2013 Elsevier B.V.

Modelling the surface circulation and thermal structure of a tropical reservoir using three-dimensional hydrodynamic lake model and remote-sensing data

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Surface-acoustic-wave counterflow micropumps for on-chip liquid motion control in two-dimensional microchannel arrays

Fully controlled liquid injection and flow in hydrophobic polydimethylsiloxane (PDMS) two-dimensional microchannel arrays based on on-chip integrated, low-voltage-driven micropumps are demonstrated. Our architecture exploits the surface-acoustic-wave (SAW) induced counterflow mechanism and the effect of nebulization anisotropies at crossing areas owing to lateral propagating SAWs. We show that by selectively exciting single or multiple SAWs, fluids can be drawn from their reservoirs and moved towards selected positions of a microchannel grid. Splitting of the main liquid flow is also demonstrated by exploiting multiple SAW beams. As a demonstrator, we show simultaneous filling of two orthogonal microchannels. The present results show that SAW micropumps are good candidates for truly integrated on-chip fluidic networks allowing liquid control in arbitrarily shaped two-dimensional microchannel arrays.