Canonical Bose gas simulations with stochastic gauges

A technique to simulate the grand canonical ensembles of interacting Bose gases is presented. Results are generated for many temperatures by averaging over energy-weighted stochastic paths, each corresponding to a solution of coupled Gross-Pitaevskii equations with phase noise. The stochastic gauge method used relies on an off-diagonal coherent-state expansion, thus taking into account all quantum correlations. As an example, the second-order spatial correlation function and momentum distribution for an interacting 1D Bose gas are calculated.

Theory manual to OctVCE – a cartesian cell CFD code with special application to blast wave problems, Department of Mechanical Engineering Report 2007/12

OctVCE is a cartesian cell CFD code produced especially for numerical simulations of shock and blast wave interactions with complex geometries, in particular, from explosions. Virtual Cell Embedding (VCE) was chosen as its cartesian cell kernel for its simplicity and sufficiency for practical engineering design problems. The code uses a finite-volume formulation of the unsteady Euler equations with a second order explicit Runge-Kutta Godonov (MUSCL) scheme. Gradients are calculated using a least-squares method with a minmod limiter. Flux solvers used are AUSM, AUSMDV and EFM. No fluid-structure coupling or chemical reactions are allowed, but gas models can be perfect gas and JWL or JWLB for the explosive products. This report also describes the code’s ‘octree’ mesh adaptive capability and point-inclusion query procedures for the VCE geometry engine. Finally, some space will also be devoted to describing code parallelization using the shared-memory OpenMP paradigm. The user manual to the code is to be found in the companion report 2007/13.

Truncation errors in finite difference models for solute transport equation with first-order reaction

The truncation errors associated with finite difference solutions of the advection-dispersion equation with first-order reaction are formulated from a Taylor analysis. The error expressions are based on a general form of the corresponding difference equation and a temporally and spatially weighted parametric approach is used for differentiating among the various finite difference schemes. The numerical truncation errors are defined using Peclet and Courant numbers and a new Sink/Source dimensionless number. It is shown that all of the finite difference schemes suffer from truncation errors. Tn particular it is shown that the Crank-Nicolson approximation scheme does not have second order accuracy for this case. The effects of these truncation errors on the solution of an advection-dispersion equation with a first order reaction term are demonstrated by comparison with an analytical solution. The results show that these errors are not negligible and that correcting the finite difference scheme for them results in a more accurate solution. (C) 1999 Elsevier Science B.V. All rights reserved.

A spectral-based cusum test of evolutionary change

We develop a test of evolutionary change that incorporates a null hypothesis of homogeneity, which encompasses time invariance in the variance and autocovariance structure of residuals from estimated econometric relationships. The test framework is based on examining whether shifts in spectral decomposition between two frames of data are significant. Rejection of the null hypothesis will point not only to weak nonstationarity but to shifts in the structure of the second-order moments of the limiting distribution of the random process. This would indicate that the second-order properties of any underlying attractor set has changed in a statistically significant way, pointing to the presence of evolutionary change. A demonstration of the test's applicability to a real-world macroeconomic problem is accomplished by applying the test to the Australian Building Society Deposits (ABSD) model.

A director theory for visco-elastic folding instabilities in multilayered rock

A model for finely layered visco-elastic rock proposed by us in previous papers is revisited and generalized to include couple stresses. We begin with an outline of the governing equations for the standard continuum case and apply a computational simulation scheme suitable for problems involving very large deformations. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered beam under compression. We analyse folding up to 40% shortening. The standard continuum solution becomes unstable for extreme values of the shear/normal viscosity ratio. The instability is a consequence of the neglect of the bending stiffness/viscosity in the standard continuum model. We suggest considering these effects within the framework of a couple stress theory. Couple stress theories involve second order spatial derivatives of the velocities/displacements in the virtual work principle. To avoid C-1 continuity in the finite element formulation we introduce the spin of the cross sections of the individual layers as an independent variable and enforce equality to the spin of the unit normal vector to the layers (-the director of the layer system-) by means of a penalty method. We illustrate the convergence of the penalty method by means of numerical solutions of simple shears of an infinite layer for increasing values of the penalty parameter. For the shear problem we present solutions assuming that the internal layering is oriented orthogonal to the surfaces of the shear layer initially. For high values of the ratio of the normal-to the shear viscosity the deformation concentrates in thin bands around to the layer surfaces. The effect of couple stresses on the evolution of folds in layered structures is also investigated. (C) 2002 Elsevier Science Ltd. All rights reserved.

Quantum interference in optical fields and atomic radiation

We discuss the connection between quantum interference effects in optical beams and radiation fields emitted from atomic systems. We illustrate this connection by a study of the first- and second-order correlation functions of optical fields and atomic dipole moments. We explore the role of correlations between the emitting systems and present examples of practical methods to implement two systems with non-orthogonal dipole moments. We also derive general conditions for quantum interference in a two-atom system and for a control of spontaneous emission. The relation between population trapping and dark states is also discussed. Moreover, we present quantum dressed-atom models of cancellation of spontaneous emission, amplification on dark transitions, fluorescence quenching and coherent population trapping.

On learning context-free and context-sensitive languages

The long short-term memory (LSTM) is not the only neural network which learns a context sensitive language. Second-order sequential cascaded networks (SCNs) are able to induce means from a finite fragment of a context-sensitive language for processing strings outside the training set. The dynamical behavior of the SCN is qualitatively distinct from that observed in LSTM networks. Differences in performance and dynamics are discussed.

The nonexistence of spurious solutions to discrete, two-point boundary value problems

We investigate difference equations which arise as discrete approximations to two-point boundary value problems for systems of second-order, ordinary differential equations. We formulate conditions under which all solutions to the discrete problem satisfy certain a priori bounds which axe independent of the step-size. As a result, the nonexistence of spurious solutions are guaranteed. Some existence and convergence theorems for solutions to the discrete problem are also presented. (C) 2002 Elsevier Science Ltd. All rights reserved.

The uniqueness of solutions to discrete, victor, two-point boundary value problems

Difference equations which may arise as discrete approximations to two-point boundary value problems for systems of second-order, ordinary differential equations are investigated and conditions are formulated under which solutions to the discrete problem are unique. Some existence, uniqueness implies existence, and convergence theorems for solutions to the discrete problem are also presented.

A geometric approach to quantum circuit lower bounds

What is the minimal size quantum circuit required to exactly implement a specified n-qubit unitary operation, U, without the use of ancilla qubits? We show that a lower bound on the minimal size is provided by the length of the minimal geodesic between U and the identity, I, where length is defined by a suitable Finsler metric on the manifold SU(2(n)). The geodesic curves on these manifolds have the striking property that once an initial position and velocity are set, the remainder of the geodesic is completely determined by a second order differential equation known as the geodesic equation. This is in contrast with the usual case in circuit design, either classical or quantum, where being given part of an optimal circuit does not obviously assist in the design of the rest of the circuit. Geodesic analysis thus offers a potentially powerful approach to the problem of proving quantum circuit lower bounds. In this paper we construct several Finsler metrics whose minimal length geodesics provide lower bounds on quantum circuit size. For each Finsler metric we give a procedure to compute the corresponding geodesic equation. We also construct a large class of solutions to the geodesic equation, which we call Pauli geodesics, since they arise from isometries generated by the Pauli group. For any unitary U diagonal in the computational basis, we show that: (a) provided the minimal length geodesic is unique, it must be a Pauli geodesic; (b) finding the length of the minimal Pauli geodesic passing from I to U is equivalent to solving an exponential size instance of the closest vector in a lattice problem (CVP); and (c) all but a doubly exponentially small fraction of such unitaries have minimal Pauli geodesics of exponential length.

Coal ash conversion into effective adsorbents for removal of heavy metals and dyes from wastewater

Fly ash was modified by hydrothermal treatment using NaOH solutions under various conditions for zeolite synthesis. The XRD patterns are presented. The results indicated that the samples obtained after treatment are much different. The XRD profiles revealed a number of new reflexes, suggesting a phase transformation probably occurred. Both heat treatment and chemical treatment increased the surface area and pore volume. It was found that zeolite P would be formed at the conditions of higher NaOH concentration and temperature. The treated fly ash was tested for adsorption of heavy metal ions and dyes in aqueous solution. It was shown that fly ash and the modified forms could effectively absorb heavy metals and methylene blue but not effectively adsorb rhodamine B. Modifying fly ash with NaOH solution would significantly enhance the adsorption capacity depending on the treatment temperature, time, and base concentration. The adsorption capacity of methylene blue would increases with pH of the dye solution and the sorption capacity of FA-NaOH could reach 5 x 10(-5) mol/g. The adsorption isotherm could be described by the Langmuir and Freundlich isotherm equations. Removal of copper and nickel ions could also be achieved on those treated fly ash. The removal efficiency for copper and nickel ions could be from 30% to 90% depending on the initial concentrations. The increase in adsorption temperature will enhance the adsorption efficiency for both heavy metals. The pseudo second-order kinetics would be better for fitting the dynamic adsorption of Cu and Ni ions. (c) 2005 Elsevier B.V. All rights reserved.

Simulations of Bose fields at finite temperature

We introduce a time-dependent projected Gross-Pitaevskii equation to describe a partially condensed homogeneous Bose gas, and find that this equation will evolve randomized initial wave functions to equilibrium. We compare our numerical data to the predictions of a gapless, second order theory of Bose-Einstein condensation [S. A. Morgan, J. Phys. B 33, 3847 (2000)], and find that we can determine a temperature when the theory is valid. As the Gross-Pitaevskii equation is nonperturbative, we expect that it can describe the correct thermal behavior of a Bose gas as long as all relevant modes are highly occupied. Our method could be applied to other boson fields.

Molecular cocrystals of carboxylic acids. XXVII - An investigation into the use of triphenylphosphine oxide cocrystals for non-linear optics: the crystal structure of the adduct of triphenylphosphine oxide with N-methylpyrrole-2-carboxylic acid

Four adducts of triphenylphosphine oxide with aromatic carboxylic acids have been synthesized and tested for second-order non-linear optical properties. These were with N-methylpyrrole-2-carboxylic acid (I), indole-2-carboxylic acid (2), 3-dimethylaminobenzoic acid (3), and thiophen-2-carboxylic acid (4). Compound (1) produced clear, colourless crystals (space group P2(1)2(1)2(1) With a 9.892(1), b 14.033(1), c 15.305(1) Angstrom, Z 4) which allowed the structure to be determined by X-ray diffraction.

Adolescent reputation enhancement: Differentiating delinquent, nondelinquent, and at-risk youths

This research reports the findings of two studies conducted to measure and then investigate differences between delinquent, nondelinquent, and at-risk youths' orientations towards reputation enhancement. In the first study, concerning item selection and scale development, the factor structure and content validity of a potential Reputation Enhancement Scale were tested by examining the item responses of the scale completed by 230 high-school students. In the second study, the scale was validated by comparing the item responses of 80 delinquent, 90 at-risk, and 90 nondelinquent adolescents with the responses of the original students. The instrument was found to be reliable (alphas from .64 to .92), indicating that the factors are dependable across different samples, and the coefficients of congruence were sufficiently high to investigate meaningful group differences. Three second-order factors (Conforming Reputation, Nonconforming Reputation, Self-presentation) were derived from the 15 first-order factors. Although multivariate analyses revealed significant differences between the reputational orientations of delinquent, at-risk, and nondelinquent participants, the self-presentation second-order factor did not differentiate the three groups.