Effect of superfluidity on the structure of the inner crust of neutron stars

A self-consistent, completely quantum calculation of the structure of the inner crust of neutron stars is carried out in the Wigner-Seitz approximation with a realistic phenomenological nuclear energy functional, where pair correlations of neutrons and protons are included in the explicit form. It has been shown that the superfluidity of neutrons and protons affects the structure of the ground state of the crust.     

Deformed Quantum Calogero-Moser Problems and Lie Superalgebras

The deformed quantum Calogero-Moser-Sutherland problems related to the root systems of the contragredient Lie superalgebras are introduced. The construction is based on the notion of the generalized root systems suggested by V. Serganova. For the classical series a recurrent formula for the quantum integrals is found, which implies the integrability of these problems. The corresponding algebras of the quantum integrals are investigated, the explicit formulas for their Poincare series for generic values of the deformation parameter are presented.     

Methods for implementing the stream boundary condition in DSMC computations

In the direct simulation Monte‐Carlo (DSMC) method for simulating rarefied gas flows, the velocities of simulator particles that cross a simulation boundary and enter the simulation space are typically generated using the acceptance–rejection procedure that samples the velocities from a truncated theoretical velocity distribution that excludes low and high velocities. This paper analyses an alternative technique, where the velocities of entering particles are obtained by extending the simulation procedures to a region adjacent to the simulation space, and considering the movement of particles generated within that region during the simulation time step. The alternative method may be considered as a form of acceptance–rejection procedure, and permits the generation of all possible velocities, although the population of high velocities is depleted with respect to the theoretical distribution. Nevertheless, this is an improvement over the standard acceptance–rejection method. Previous implementations of the alternative method gave a number flux lower than the theoretical number required. Two methods for obtaining the correct number flux are presented. For upstream boundaries in high‐speed flows, the alternative method is more computationally efficient than the acceptance–rejection method. However, for downstream boundaries, the alternative method is extremely inefficient. The alternative method, with the correct theoretical number flux, should therefore be used in DSMC computations in favour of the acceptance–rejection method for upstream boundaries in high‐speed flows. Copyright © 2003 John Wiley & Sons, Ltd.     

An improved numerical simulator for multiphase flow in porous media

In basin modelling the thermodynamics of a multicomponent multiphase fluid flux are computationally too expensive when derived from an equation of state and the Gibbs equality constraints. In this article we present a novel implicit molar mass formulation technique using binary mixture thermodynamics. The two proposed solution methods, with and without cross derivative terms between components, are based on a preconditioned Newton‐GMRES scheme for each time‐step with analytical computation of the derivatives. These new algorithms reduce significantly the numerical effort for the computation of the molar masses, and we illustrate the behavior of these methods with numerical computations. Copyright © 2004 John Wiley & Sons Ltd.     

Assessment of two‐equation turbulence modelling for high Reynolds number hydrofoil flows

This paper presents an evaluation of the capability of turbulence models available in the commercial CFD code FLUENT 6.0 for their application to hydrofoil turbulent boundary layer separation flow at high Reynolds numbers. Four widely applied two‐equation RANS turbulence models were assessed through comparison with experimental data at Reynolds numbers of 8.284×106 and 1.657×107. They were the standard k–εmodel, the realizable k–εmodel, the standard k–ωmodel and the shear‐stress‐transport (SST) k–ωmodel. It has found that the realizable k–εturbulence model used with enhanced wall functions and near‐wall modelling techniques, consistently provides superior performance in predicting the flow characteristics around the hydrofoil. Copyright © 2004 John Wiley & Sons, Ltd.     

Electron spin lifetimes in Hg0.78Cd0.22 Te and InSb

We have made direct pump–probe measurements of spin lifetimes in long wavelength narrow-gap semiconductors at wavelengths between 4 and 10 μm and from 4 to 300 K. In particular, we measure remarkably long spin lifetimes, τ_s300 ps, even at 300 K for epilayers of degenerate n-type InSb. In this material the mobility is approximately constant between 77 and 300 K, and we find that τ_s is approximately constant in this temperature range. In order to determine the dominant spin relaxation mechanism we have investigated the temperature dependence of τ_s in non-degenerate lightly n-type Hg_0.78Cd_0.22Te of approximately the same band gap as InSb, and find that τ_s varies from 356 ps at 150 K to 24 ps at 300 K. Our results, both in magnitude and temperature dependence of τ_s, imply that the Elliott–Yafet model dominates in these materials.     

Structural and vibrational properties of dilute GaNxAs1−x(P1−x)

We report a comprehensive analyzes of the Fourier transform infrared (FTIR) absorption and Raman scattering data on the structural and vibrational properties of dilute ternary GaAs_1−xN_x,[GaP_1−xN_x] (x<0.03) alloys grown on GaAs [GaP] by metal organic chemical vapor deposition (MOCVD) and solid source molecular beam epitaxy (MBE). By using realistic total energy and lattice dynamical calculations, the origin of experimentally observed N-induced vibrational features are characterized. Useful information is obtained about the structural stability, vibrational frequencies, lattice relaxations and compositional disorder in GaNAs (GaNP) alloys. At lower composition (x<0.015) most of the N atoms occupy the As [P] sublattice {N_As[N_P]}—they prefer moving out of their substitutional sites to more energetically favorable locations at higher x. Our results for the N-isotopic shifts of local mode frequencies compare favorably well with the existing FTIR data.