Partition function for hindered, one-and multi-dimensional rotors

Writing the hindered rotor (hr) partition function as the trace of (rho) over cap = e(-beta(H) over cap hr), we approximate it by the sum of contributions from a set of points in position space. The contribution of the density matrix from each point is approximated by performing a local harmonic expansion around it. The highlight of this method is that it can be easily extended to multidimensional systems. Local harmonic expansion leads to a breakdown of the method a low temperatures. In order to calculate the partition function at low temperatures, we suggest a matrix multiplication procedure. The results obtained using these methods closely agree with the exact partition function at all temperature ranges. Our method bypasses the evaluation of eigenvalues and eigenfunctions and evaluates the density matrix for internal rotation directly. We also suggest a procedure to account for the antisymmetry of the total wavefunction in the same. (C) 2012 Elsevier B.V. All rights reserved.

Universal behaviour of shock precursors in the presence of efficient cosmic ray acceleration

The self-consistent interaction between energetic particles and self-generated hydromagnetic waves in a cosmic ray pressure dominated plasma is considered. Using a three-dimensional hybrid magnetohydrodynamics (MHD)-kinetic code, which utilizes a spherical harmonic expansion of the Vlasov-Fokker-Planck equation, high-resolution simulations of the magnetic field growth including feedback on the cosmic rays are carried out. It is found that for shocks with high cosmic ray acceleration efficiency, the magnetic fields become highly disorganized, resulting in near isotropic diffusion, independent of the initial orientation of the ambient magnetic field. The possibility of sub-Bohm diffusion is demonstrated for parallel shocks, while the diffusion coefficient approaches the Bohm limit from below for oblique shocks. This universal behaviour suggests that Bohm diffusion in the root-mean-squared field inferred from observation may provide a realistic estimate for the maximum energy acceleration time-scale in young supernova remnants. Although disordered, the magnetic field is not self-similar suggesting a non-uniform energy-dependent behaviour of the energetic particle transport in the precursor. Possible indirect radiative signatures of cosmic ray driven magnetic field amplification are discussed.

Vibrational intensities in methane

The absorption intensities of the two infra-red active vibrations in methane have been obtained from a perturbation calculation on the equilibrium wave functions derived in the preceding paper. The perturbation field is the change in the potential field due to the nuclei which results from moving the nuclei in the vibrational coordinate concerned, and a simplified form of second order perturbation theory, developed by Pople and Schofield, is used for the calculation. The main approximation involved is the neglect of f and higher harmonics in the spherical harmonic expansion of the nuclear field. The resulting dipole moment derivatives are approximately three times larger than the experimental values, but they show qualitative features and sign relationships which are significant.

Low-frequency absorption cross section of the electromagnetic waves for extreme Reissner-Nordstrom black holes in higher dimensions

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

Structural properties and energetics of diffuse Rb-87 clusters in three-dimension

A correlated two-body basis function is used to describe the three-dimensional bosonic clusters interacting via two-body van der Waals potential. We calculate the ground state and the zero orbital angular momentum excited states for Rb-N clusters with up to N = 40. We solve the many-particle Schrodinger equation by potential harmonics expansion method, which keeps all possible two-body correlations in the calculation and determines the lowest effective many-body potential. We study energetics and structural properties for such diffuse clusters both at dimer and tuned scattering length. The motivation of the present study is to investigate the possibility of formation of N-body clusters interacting through the van der Waals interaction. We also compare the system with the well studied He, Ne, and Ar clusters. We also calculate correlation properties and observe the generalised Tjon line for large cluster. We test the validity of the shape-independent potential in the calculation of the ground state energy of such diffuse cluster. These are the first such calculations reported for Rb clusters. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4730972]

Code and readme of a MATLAB-based graphical user interface program for computing functionals of the geopotential

We present a novel graphical user interface program GrafLab (GRAvity Field LABoratory) for spherical harmonic synthesis (SHS) created in MATLAB®. This program allows to comfortably compute 38 various functionals of the geopotential up to ultra-high degrees and orders of spherical harmonic expansion. For the most difficult part of the SHS, namely the evaluation of the fully normalized associated Legendre functions (fnALFs), we used three different approaches according to required maximum degree: (i) the standard forward column method (up to maximum degree 1800, in some cases up to degree 2190); (ii) the modified forward column method combined with Horner's scheme (up to maximum degree 2700); (iii) the extended-range arithmetic (up to an arbitrary maximum degree). For the maximum degree 2190, the SHS with fnALFs evaluated using the extended-range arithmetic approach takes only approximately 2-3 times longer than its standard arithmetic counterpart, i.e. the standard forward column method. In the GrafLab, the functionals of the geopotential can be evaluated on a regular grid or point-wise, while the input coordinates can either be read from a data file or entered manually. For the computation on a regular grid we decided to apply the lumped coefficients approach due to significant time-efficiency of this method. Furthermore, if a full variance-covariance matrix of spherical harmonic coefficients is available, it is possible to compute the commission errors of the functionals. When computing on a regular grid, the output functionals or their commission errors may be depicted on a map using automatically selected cartographic projection.

POSITIVE FILTERED P_N METHOD FOR LINEAR TRANSPORT EQUATIONS AND THE ASSOCIATED OPTIMIZATION ALGORITHM

We propose a positive, accurate moment closure for linear kinetic transport equations based on a filtered spherical harmonic (FP_N) expansion in the angular variable. The FP_N moment equations are accurate approximations to linear kinetic equations, but they are known to suffer from the occurrence of unphysical, negative particle concentrations. The new positive filtered P_N (FP_N+) closure is developed to address this issue. The FP_N+ closure approximates the kinetic distribution by a spherical harmonic expansion that is non-negative on a finite, predetermined set of quadrature points. With an appropriate numerical PDE solver, the FP_N+ closure generates particle concentrations that are guaranteed to be non-negative. Under an additional, mild regularity assumption, we prove that as the moment order tends to infinity, the FP_N+ approximation converges, in the L2 sense, at the same rate as the FP_N approximation; numerical tests suggest that this assumption may not be necessary. By numerical experiments on the challenging line source benchmark problem, we confirm that the FP_N+ method indeed produces accurate and non-negative solutions. To apply the FP_N+ closure on problems at large temporal-spatial scales, we develop a positive asymptotic preserving (AP) numerical PDE solver. We prove that the propose AP scheme maintains stability and accuracy with standard mesh sizes at large temporal-spatial scales, while, for generic numerical schemes, excessive refinements on temporal-spatial meshes are required. We also show that the proposed scheme preserves positivity of the particle concentration, under some time step restriction. Numerical results confirm that the proposed AP scheme is capable for solving linear transport equations at large temporal-spatial scales, for which a generic scheme could fail. Constrained optimization problems are involved in the formulation of the FP_N+ closure to enforce non-negativity of the FP_N+ approximation on the set of quadrature points. These optimization problems can be written as strictly convex quadratic programs (CQPs) with a large number of inequality constraints. To efficiently solve the CQPs, we propose a constraint-reduced variant of a Mehrotra-predictor-corrector algorithm, with a novel constraint selection rule. We prove that, under appropriate assumptions, the proposed optimization algorithm converges globally to the solution at a locally q-quadratic rate. We test the algorithm on randomly generated problems, and the numerical results indicate that the combination of the proposed algorithm and the constraint selection rule outperforms other compared constraint-reduced algorithms, especially for problems with many more inequality constraints than variables.

Modeling non-harmonic behavior of materials from experimental inelastic neutron scattering and thermal expansion measurements

Based on thermodynamic principles, we derive expressions quantifying the non-harmonic vibrational behavior of materials, which are rigorous yet easily evaluated from experimentally available data for the thermal expansion coefficient and the phonon density of states. These experimentally- derived quantities are valuable to benchmark first-principles theoretical predictions of harmonic and non-harmonic thermal behaviors using perturbation theory, ab initio molecular-dynamics, or Monte-Carlo simulations. We illustrate this analysis by computing the harmonic, dilational, and anharmonic contributions to the entropy, internal energy, and free energy of elemental aluminum and the ordered compound FeSi over a wide range of temperature. Results agree well with previous data in the literature and provide an efficient approach to estimate anharmonic effects in materials.

Expansion of a Bose-Einstein condensate formed on a joint harmonic and one-dimensional optical-lattice potential

We study the expansion of a Bose-Einstein condensate trapped in a combined optical-lattice and axially-symmetric harmonic potential using the numerical solution of the mean-field Gross-Pitaevskii equation. First, we consider the expansion of such a condensate under the action of the optical-lattice potential alone. In this case the result of numerical simulation for the axial and radial sizes during expansion is in agreement with two experiments by Morsch et al (2002 Phys. Rev. A 66 021601(R) and 2003 Laser Phys. 13 594). Finally, we consider the expansion under the action of the harmonic potential alone. In this case the oscillation, and the disappearance and revival of the resultant interference pattern is in agreement with the experiment by Muller et al (2003 J. Opt. B: Quantum Semiclass. Opt. 5 S38).

The virial expansion of a classical interacting system

We consider N particles interacting pairwise by an inverse square potential in one dimension (Calogero-Sutherland-Moser model). For a system placed in a harmonic trap, its classical partition function for the repulsive regime is recognised in the literature. We start by presenting a concise re-derivation of this result. The equation of state is then calculated both for the trapped and the homogeneous gas. Finally, the classical limit of Wu's distribution function for fractional exclusion statistics is obtained and we re-derive the classical virial expansion of the homogeneous gas using this distribution function.

Strong-coupling expansion for ultracold bosons in an optical lattice at finite temperatures in the presence of superfluidity

We develop a strong-coupling (t << U) expansion technique for calculating the density profile for bosonic atoms trapped in an optical lattice with an overall harmonic trap at finite temperature and finite on-site interaction in the presence of superfluid regions. Our results match well with quantum Monte Carlo simulations at finite temperature. We also show that the superfluid order parameter never vanishes in the trap due to the proximity effect. Our calculations for the scaled density in the vacuum-to-superfluid transition agree well with the experimental data for appropriate temperatures. We present calculations for the entropy per particle as a function of temperature which can be used to calibrate the temperature in experiments. We also discuss issues connected with the demonstration of universal quantum critical scaling in the experiments.

Harmonic generation by noble-gas atoms in the near-IR regime using ab initio time-dependent R-matrix theory

We demonstrate the capability of ab initio time-dependent R-matrix theory to obtain accurate harmonic generation spectra of noble-gas atoms at near-IR wavelengths between 1200 and 1800 nm and peak intensities up to 1.8 × 10^(14) W/cm^(2). To accommodate the excursion length of the ejected electron, we use an angular-momentum expansion up to Lmax=279. The harmonic spectra show evidence of atomic structure through the presence of a Cooper minimum in harmonic generation for Kr, and of multielectron interaction through the giant resonance for Xe. The theoretical spectra agree well with those obtained experimentally.

Extension of explicit formulas in Poissonian white noise analysis using harmonic analysis on configuration spaces

Harmonic analysis on configuration spaces is used in order to extend explicit expressions for the images of creation, annihilation, and second quantization operators in L2-spaces with respect to Poisson point processes to a set of functions larger than the space obtained by directly using chaos expansion. This permits, in particular, to derive an explicit expression for the generator of the second quantization of a sub-Markovian contraction semigroup on a set of functions which forms a core of the generator.

A self-similar expansion model for use in solar wind transient propogation studies

Since the advent of wide-angle imaging of the inner heliosphere, a plethora of techniques have been developed to investigate the three-dimensional structure and kinematics of solar wind transients, such as coronal mass ejections, from their signatures in single- and multi-spacecraft imaging observations. These techniques, which range from the highly complex and computationally intensive to methods based on simple curve fitting, all have their inherent advantages and limitations. In the analysis of single-spacecraft imaging observations, much use has been made of the fixed φ fitting (FPF) and harmonic mean fitting (HMF) techniques, in which the solar wind transient is considered to be a radially propagating point source (fixed φ, FP, model) and a radially expanding circle anchored at Sun centre (harmonic mean, HM, model), respectively. Initially, we compare the radial speeds and propagation directions derived from application of the FPF and HMF techniques to a large set of STEREO/Heliospheric Imager (HI) observations. As the geometries on which these two techniques are founded constitute extreme descriptions of solar wind transients in terms of their extent along the line of sight, we describe a single-spacecraft fitting technique based on a more generalized model for which the FP and HM geometries form the limiting cases. In addition to providing estimates of a transient’s speed and propagation direction, the self-similar expansion fitting (SSEF) technique provides, in theory, the capability to estimate the transient’s angular extent in the plane orthogonal to the field of view. Using the HI observations, and also by performing a Monte Carlo simulation, we assess the potential of the SSEF technique.

Free expansion of fermionic dark solitons in a boson-fermion mixture

We use a time-dependent dynamical mean-field-hydrodynamic model to study the formation of fermionic dark solitons in a trapped degenerate Fermi gas mixed with a Bose-Einstein condensate in a harmonic as well as a periodic optical-lattice potential. The dark soliton with a 'notch' in the probability density with a zero at the minimum is simulated numerically as a nonlinear continuation of the first vibrational excitation of the linear mean-field-hydrodynamic equations, as suggested recently for pure bosons. We study the free expansion of these dark solitons as well as the consequent increase in the size of their central notch and discuss the possibility of experimental observation of the notch after free expansion.