Bifurcations and Complete Chaos for the Diamagnetic Kepler Problem

We describe the structure of bifurcations in the unbounded classical Diamagnetic Kepler problem. We conjecture that this system does not have any stable orbits and that the non-wandering set is described by a complete trinary symbolic dynamics for scaled energies larger then $\epsilon_c=0.328782\ldots$.Comment: 15 pages PostScript uuencoded with figure

Alas, the dark matter structures were not that trivial

The radial density profile of dark matter structures has been observed to have an almost universal behaviour in numerical simulations, however, the physical reason for this behaviour remains unclear. It has previously been shown that if the pseudo phase-space density, rho/sigma_d^epsilon, is a beautifully simple power-law in radius, with the "golden values" epsilon=3 and d=r (i.e., the phase-space density is only dependent on the radial component of the velocity dispersion), then one can analytically derive the radial variation of the mass profile, dispersion profile etc. That would imply, if correct, that we just have to explain why rho/sigma^3_r ~r^{-alpha}, and then we would understand everything about equilibrated DM structures. Here we use a set of simulated galaxies and clusters of galaxies to demonstrate that there are no such golden values, but that each structure instead has its own set of values. Considering the same structure at different redshifts shows no evolution of the phase-space parameters towards fixed points. There is also no clear connection between the halo virialized mass and these parameters. This implies that we still do not understand the origin of the profiles of dark matter structures.Comment: 4 pages, 3 figures, accepted for publication in ApJ

On the nonlocal viscosity kernel of mixtures

In this report we investigate the multiscale hydrodynamical response of a liquid as a function of mixture composition. This is done via a series of molecular dynamics simulations where the wave vector dependent viscosity kernel is computed for three mixtures each with 7-15 different compositions. We observe that the nonlocal viscosity kernel is dependent on composition for simple atomic mixtures for all the wave vectors studied here, however, for a model polymer melt mixture the kernel is independent of composition for large wave vectors. The deviation from ideal mixing is also studied. Here it is shown that a Lennard-Jones mixture follows the ideal mixing rule surprisingly well for a large range of wave vectors, whereas for both the Kob-Andersen mixture and the polymer melt large deviations are found. Furthermore, for the polymer melt the deviation is wave vector dependent such that there exists a critical length scale at which the ideal mixing goes from under-estimating to over-estimating the viscosity

Infall near clusters of galaxies: comparing gas and dark matter velocity profiles

We consider the dynamics in and near galaxy clusters. Gas, dark matter and galaxies are presently falling into the clusters between approximately 1 and 5 virial radii. At very large distances, beyond 10 virial radii, all matter is following the Hubble flow, and inside the virial radius the matter particles have on average zero radial velocity. The cosmological parameters are imprinted on the infall profile of the gas, however, no method exists, which allows a measurement of it. We consider the results of two cosmological simulations (using the numerical codes RAMSES and Gadget) and find that the gas and dark matter radial velocities are very similar. We derive the relevant dynamical equations, in particular the generalized hydrostatic equilibrium equation, including both the expansion of the Universe and the cosmological background. This generalized gas equation is the main new contribution of this paper. We combine these generalized equations with the results of the numerical simulations to estimate the contribution to the measured cluster masses from the radial velocity: inside the virial radius it is negligible, and inside two virial radii the effect is below 40%, in agreement the earlier analyses for DM. We point out how the infall velocity in principle may be observable, by measuring the gas properties to distance of about two virial radii, however, this is practically not possible today.Comment: 7 pages, 3 figures, to appear in MNRA

An Iterative Receiver for OFDM With Sparsity-Based Parametric Channel Estimation

In this work we design a receiver that iteratively passes soft information between the channel estimation and data decoding stages. The receiver incorporates sparsity-based parametric channel estimation. State-of-the-art sparsity-based iterative receivers simplify the channel estimation problem by restricting the multipath delays to a grid. Our receiver does not impose such a restriction. As a result it does not suffer from the leakage effect, which destroys sparsity. Communication at near capacity rates in high SNR requires a large modulation order. Due to the close proximity of modulation symbols in such systems, the grid-based approximation is of insufficient accuracy. We show numerically that a state-of-the-art iterative receiver with grid-based sparse channel estimation exhibits a bit-error-rate floor in the high SNR regime. On the contrary, our receiver performs very close to the perfect channel state information bound for all SNR values. We also demonstrate both theoretically and numerically that parametric channel estimation works well in dense channels, i.e., when the number of multipath components is large and each individual component cannot be resolved.Comment: Major revision, accepted for IEEE Transactions on Signal Processin

Method for locating low-energy solutions within DFT+U

The widely employed DFT+U formalism is known to give rise to many self-consistent yet energetically distinct solutions in correlated systems, which can be highly problematic for reliably predicting the thermodynamic and physical properties of such materials. Here we study this phenomenon in the bulk materials UO_2, CoO, and NiO, and in a CeO_2 surface. We show that the following factors affect which self-consistent solution a DFT+U calculation reaches: (i) the magnitude of U; (ii) initial correlated orbital occupations; (iii) lattice geometry; (iv) whether lattice symmetry is enforced on the charge density; and (v) even electronic mixing parameters. These various solutions may differ in total energy by hundreds of meV per atom, so identifying or approximating the ground state is critical in the DFT+U scheme. We propose an efficient U-ramping method for locating low-energy solutions, which we validate in a range of test cases. We also suggest that this method may be applicable to hybrid functional calculations

Rates of induced abortion in Denmark according to age, previous births and previous abortions

Background: Whereas the effects of various socio-demographic determinants on a woman’s risk of having an abortion are relatively well-documented, less attention has been given to the effect of previous abortions and births. Objective: To study the effect of previous abortions and births on Danish women’s risk of an abortion, in addition to a number of demographic and personal characteristics. Data and methods: From the Fertility of Women and Couples Dataset we obtained data on the number of live births and induced abortions by year (1981-2001), age (16-39), county of residence and marital status. Logistic regression analysis was used to estimate the influence of the explanatory variables on the probability of having an abortion in a relevant year. Main findings and conclusion: A woman’s risk of having an abortion increases with the number of previous births and previous abortions. Some interactions were was found in the way a woman’s risk of abortion varies with calendar year, age and parity. The risk of an abortion for women with no children decreases while the risk of an abortion for women with children increases over time. Furthermore, the risk of an abortion decreases with age, but relatively more so for women with children compared to childless women. Trends for teenagers are discussed in a separate section.abortion parity, abortion rates, family situation, fertility parity, urbanization

Asymmetries in the CMB anisotropy field

We report on the results from two independent but complementary statistical analyses of the WMAP first-year data, based on the power spectrum and N-point correlation functions. We focus on large and intermediate scales (larger than about 3 degrees) and compare the observed data against Monte Carlo ensembles with WMAP-like properties. In both analyses, we measure the amplitudes of the large-scale fluctuations on opposing hemispheres and study the ratio of the two amplitudes. The power-spectrum analysis shows that this ratio for WMAP, as measured along the axis of maximum asymmetry, is high at the 95%-99% level (depending on the particular multipole range included). The axis of maximum asymmetry of the WMAP data is weakly dependent on the multipole range under consideration but tends to lie close to the ecliptic axis. In the N-point correlation function analysis we focus on the northern and southern hemispheres defined in ecliptic coordinates, and we find that the ratio of the large-scale fluctuation amplitudes is high at the 98%-99% level. Furthermore, the results are stable with respect to choice of Galactic cut and also with respect to frequency band. A similar asymmetry is found in the COBE-DMR map, and the axis of maximum asymmetry is close to the one found in the WMAP data.Comment: 6 pages, 5 figures; version to appear in ApJ, textual improvements, added reference