Robust algorithms for solving stochastic partial differential equations

A robust semi-implicit central partial difference algorithm for the numerical solution of coupled stochastic parabolic partial differential equations (PDEs) is described. This can be used for calculating correlation functions of systems of interacting stochastic fields. Such field equations can arise in the description of Hamiltonian and open systems in the physics of nonlinear processes, and may include multiplicative noise sources. The algorithm can be used for studying the properties of nonlinear quantum or classical field theories. The general approach is outlined and applied to a specific example, namely the quantum statistical fluctuations of ultra-short optical pulses in chi((2)) parametric waveguides. This example uses a non-diagonal coherent state representation, and correctly predicts the sub-shot noise level spectral fluctuations observed in homodyne detection measurements. It is expected that the methods used wilt be applicable for higher-order correlation functions and other physical problems as well. A stochastic differencing technique for reducing sampling errors is also introduced. This involves solving nonlinear stochastic parabolic PDEs in combination with a reference process, which uses the Wigner representation in the example presented here. A computer implementation on MIMD parallel architectures is discussed. (C) 1997 Academic Press.

Free field and parafermionic realizations of twisted su(3)(k)((2)) current algebra

Free field and twisted parafermionic representations of twisted su(3)(k)((2)) current algebra are obtained. The corresponding twisted Sugawara energy-momentum tensor is given in terms of three (beta, gamma) pairs and two scalar fields and also in terms of twisted parafermionic currents and one scalar field. Two screening currents of the first kind are presented in terms of the free fields.

Three-dimensional solitons in coupled atomic-molecular Bose-Einstein condensates

We present a theoretical analysis of three-dimensional (3D) matter-wave solitons and their stability properties in coupled atomic and molecular Bose-Einstein condensates (BECs). The soliton solutions to the mean-field equations are obtained in an approximate analytical form by means of a variational approach. We investigate soliton stability within the parameter space described by the atom-molecule conversion coupling, the atom-atom s-wave scattering, and the bare formation energy of the molecular species. In terms of ordinary optics, this is analogous to the process of sub- or second-harmonic generation in a quadratic nonlinear medium modified by a cubic nonlinearity, together with a phase mismatch term between the fields. While the possibility of formation of multidimensional spatiotemporal solitons in pure quadratic media has been theoretically demonstrated previously, here we extend this prediction to matter-wave interactions in BEC systems where higher-order nonlinear processes due to interparticle collisions are unavoidable and may not be neglected. The stability of the solitons predicted for repulsive atom-atom interactions is investigated by direct numerical simulations of the equations of motion in a full 3D lattice. Our analysis also leads to a possible technique for demonstrating the ground state of the Schrodinger-Newton and related equations that describe Bose-Einstein condensates with nonlocal interparticle forces.

Quantum computing and polynomial equations over the finite field Z(2)

What is the computational power of a quantum computer? We show that determining the output of a quantum computation is equivalent to counting the number of solutions to an easily computed set of polynomials defined over the finite field Z(2). This connection allows simple proofs to be given for two known relationships between quantum and classical complexity classes, namely BQP subset of P-#P and BQP subset of PP.

Heat flow for the Yang-Mills-Higgs field and the Hermitian Yang-Mills-Higgs metric

For a parameter lambda > 0, we study a type of vortex equations, which generalize the well-known Hermitian-Einstein equation, for a connection A and a section phi of a holomorphic vector bundle E over a Kahler manifold X. We establish a global existence of smooth solutions to heat flow for a self-dual Yang-Mills-Higgs field on E. Assuming the lambda -stability of (E, phi), we prove the existence of the Hermitian Yang-Mills-Higgs metric on the holomorphic bundle E by studying the limiting behaviour of the gauge flow.

Twisted sl(3, C)(k)((2)) current algebra: free field representation and screening currents

Motivated by application of twisted current algebra in description of the entropy of Ads(3) black hole, we investigate the simplest twisted current algebra sl(3, c)(k)((2)). Free field representation of the twisted algebra, and the corresponding twisted Sugawara energy-momentum tensor are obtained by using three (beta, gamma) pairs and two scalar fields. Primary fields and two screening currents of the first kind are presented. (C) 2001 Published by Elsevier Science B.V.

Investigation of the flow field of a highly heated jet of air

Measurements of mean and fluctuating velocity and temperature and their self- and cross-products to the third-order are presented for a heated axisymmetric air jet. Froude numbers in the range of 3500 13,190, Reynolds numbers in the range of 3470-8500 and non-dimensional streamwise distances. X*, from 0.27 to 1.98 are covered by the data. It was found that turbulence intensity decreases for the heated jet in the region between the inertia dominated and the buoyancy dominated regions which is contrary to findings with helium jets mixing with ambient air to produce density fluctuations. The effects of heating on the turbulent kinetic energy budget and the temperature variance budget show small differences for the inertia dominated region and the intermediate region which help to explain the transition process to the far field plume region. Constants are evaluated for the isotropic eddy diffusivity and generalised gradient hypothesis models as well as the scalar variance model. No significant effect of heating on the dissipation time-scale ratio was found. A novel wire array with an inclined cold wire was used. Measurements obtained with this probe are found to lead to asymmetries in some of the higher-order products. Further investigation suggested that the asymmetries are attributable to an as yet unreported interference effect produced by the leading prong of the inclined temperature wire, The effect may also have implications for inclined velocity wires which contain a temperature component when used in heated flows. (C) 2002 Elsevier Science Inc. All rights reserved.

Microcanonical temperature for a classical field: Application to Bose-Einstein condensation

We show that the projected Gross-Pitaevskii equation (PGPE) can be mapped exactly onto Hamilton's equations of motion for classical position and momentum variables. Making use of this mapping, we adapt techniques developed in statistical mechanics to calculate the temperature and chemical potential of a classical Bose field in the microcanonical ensemble. We apply the method to simulations of the PGPE, which can be used to represent the highly occupied modes of Bose condensed gases at finite temperature. The method is rigorous, valid beyond the realms of perturbation theory, and agrees with an earlier method of temperature measurement for the same system. Using this method we show that the critical temperature for condensation in a homogeneous Bose gas on a lattice with a uv cutoff increases with the interaction strength. We discuss how to determine the temperature shift for the Bose gas in the continuum limit using this type of calculation, and obtain a result in agreement with more sophisticated Monte Carlo simulations. We also consider the behavior of the specific heat.

The gap equations for spin singlet and triplet ferromagnetic superconductors

We derive gap equations for superconductivity in coexistence with ferromagnetism. We treat singlet and triplet states With either equal spin pairing (ESP) or opposite spin pairing (OSP) states, and study the behaviour of these states as a function of exchange splitting. For the s-wave singlet state we find that our gap equations correctly reproduce the Clogston-Chandrasekhar limiting behaviour and the phase diagram of the Baltensperger-Sarma equation (excluding the FFLO region). The singlet superconducting order parameter is shown to be independent of exchange splitting at zero temperature, as is assumed in the derivation of the Clogston-Chandrasekhar limit. P-wave triplet states of the OSP type behave similarly to the singlet state as a function of exchange splitting. On the other hand, ESP triplet states show a very different behaviour. In particular, there is no Clogston-Chandrasekhar limiting and the superconducting critical temperature, T-C, is actually increased by exchange splitting.

Relationship between thirty post-thaw spermatozoal characteristics and the field fertility of 11 high-use Australian dairy AI sires

This study determined the relationship between two measures of field fertility of I I high-use Australian artificial insemination (AI) dairy bulls and thirty standard laboratory assessments of spermatozoal post-thaw viability. The two measures of field fertility used, conception rates (cCR) and non-return rates (cNRR), were both corrected for all major non-bull variables. Sperm viability assessments were conducted on semen collected within the same season as that used to derive the field fertility estimates. These assessments measured sperm concentration, motility, morphology and membrane integrity at thawing, after 2 h incubation and after the swim-up sperm selection procedure. Derivations of these measures and in vitro embryo fertilizing and developmental capacity were also determined. The Genstat Statistical Package [Genstat 5 Release 4.2 Reference Manual, VSN International, Oxford, 20001 was used to conduct an analysis of variance on the viability parameters across semen straws and bulls, and to calculate the strength of correlation between each semen parameter, cNRR and cCR in a correlation matrix. Step forward multiple regression identified the combination of semen parameters that were most highly correlated with cCR and with cNRR. The sperm parameters identified as being most predictive of cCR were the percentage of morphologically normal sperm immediately post-thaw (zeroNorm), the number of morphologically normal sperm after the swim-up procedure (nSuNorm), and the rate of zygote cleavage in vitro (Clv); the predictive equation formed by these parameters accounted for 70% of variance. The predictive equation produced for cNRR contained the variables zeroNorm, the proportion of membrane intact sperm after 2 h incubation at 37 degreesC (twoMem) and Clv and accounted for 76.5% of the variation. ZeroNorm was found to be consistent across straws and semen batches within-bull and the sperm parameter with the strongest individual predictive capacity for both cCR (P = 0.1) and cNRR (P = 0.001). Post-thaw sperm parameters can be used to predict field fertility of Australian dairy sires; the calculated predictive equations are particularly useful for identifying and monitoring bulls of very high and very low potential fertility within a group. (C) 2003 Elsevier B.V. All rights reserved.

Mine waste beach profile and flow resistance equations

Four mine waste beach longitudinal profile equations are compared theoretically and in statistical analyses of profile data from 64 field and laboratory beaches formed by mine tailings, co-disposed coal mine wastes, and sand. All four equations fit the profile data well. The best performing equation both accounts for particle sorting and satisfies hydraulic constraints, and the combination of assumptions underlying it is considered to best represent the processes occurring on mine waste beaches. Combining these assumptions with the Lacey normal equation leads to a variant of the Manning resistance equation. Features that it is desirable to incorporate in theoretical and numerical models of mine waste beaches are listed.

On the Schrodinger equation involving a critical Sobolev exponent and magnetic field

We consider the semilinear Schrodinger equation -Delta(A)u + V(x)u = Q(x)vertical bar u vertical bar(2* -2) u. Assuming that V changes sign, we establish the existence of a solution u not equal 0 in the Sobolev space H-A,V(1) + (R-N). The solution is obtained by a min-max type argument based on a topological linking. We also establish certain regularity properties of solutions for a rather general class of equations involving the operator -Delta(A).

Field observations of instantaneous water slopes and horizontal pressure gradients in the swash-zone

Field observations of instantaneous water surface slopes in the swash zone are presented. For free-surface flows with a hydrostatic pressure distribution the surface slope is equivalent to the horizontal pressure gradient. Observations were made using a novel technique which in its simplest form consists of a horizontal stringline extending seaward from the beach face. Visual observation, still photography or video photography is then sufficient to determine the surface slope where the free-surface cuts the line or between reference points in the image. The method resolves the mean surface gradient over a cross-shore distance of 5 m or more to within +/- 0.001, or 1/20th -1/100th of typical beach gradients. In addition, at selected points and at any instant in time during the swash cycle, the water surface slope can be determined exactly to be dipping either seaward or landward. Close to the location of bore collapse landward dipping water surface slopes of order 0.05-0.1 occur over a very small region (order 0.5 m) at the blunt or convex leading edge of the swash. In the middle and upper swash the water surface slope at this leading edge is usually very close to horizontal or slightly seaward. Behind the leading edge, the water surface slope was observed to be very close to horizontal or dipping seaward at all times throughout the swash uprush. During the backwash the water surface slope was observed to be always dipping seaward, approaching the beach slope, and remained seaward until a new uprush edge or incident bore passed any particular cross-shore location of interest. The observations strongly Suggest that the swash boundary layer is subject to an adverse pressure gradient during uprush and a favourable pressure gradient during the backwash. Furthermore, assuming Euler's equations are a good approximation in the swash, the observations also show that the total fluid acceleration is negative (offshore) for almost the whole of the uprush and for the entire backwash. The observations are contrary to recent work suggesting significant shoreward directed accelerations and pressure gradients occur in the swash (i.e., delta u/delta t > 0 similar to delta p/delta x < 0), but consistent with analytical and numerical solutions for swash uprush and backwash. The results have important implications for sediment transport modelling in the swash zone.