Recent developments in the quantum dynamical characterization of unimolecular resonances

We give a selective review of quantum mechanical methods for calculating and characterizing resonances in small molecular systems, with an emphasis on recent progress in Chebyshev and Lanczos iterative methods. Two archetypal molecular systems are discussed: isolated resonances in HCO, which exhibit regular mode and state specificity, and overlapping resonances in strongly bound HO2, which exhibit irregular and chaotic behavior. Recent progresses for non-zero total angular momentum J calculations of resonances including parallel computing models are also included and future directions in this field are discussed.

Australian data and meta-analysis lend support for alpha-synuclein (NACP-Rep1) as a risk factor for Parkinson's disease

It remains unclear whether genetic variants in SNCA (the alpha-synuclein gene) alter risk for sporadic Parkinson's disease (PD). The polymorphic mixed sequence repeat (NACP-Rep I) in the promoter region of SNCA has been previously examined as a potential susceptibility factor for PD with conflicting results. We report genotype and allele distributions at this locus from 369 PD cases and 370 control subjects of European Australian ancestry, with alleles designated as -1, 0, +1, +2, and +3 as previously described. Allele frequencies designated (0) were less common in Australian cases compared to controls (OR = 0.80, 95% CI 0.62-1.03). Combined analysis including all previously published ancestral European Rep1 data yielded a highly significant association between the 0 allele and a reduced risk for PD (OR = 0.79, 95% CI 0.70-0.89, p = 0.0001). Further study must now proceed to examine in detail this interesting and biologically plausible genetic association. (C) 2004 Elsevier Ireland Ltd. All rights reserved.

A comparison of surface brightness profiles for ultracompact dwarfs and dwarf elliptical nuclei: implications for the “threshing” scenario

Using imaging from the Hubble Space Telescope, we derive surface brightness profiles for ultracompact dwarfs in the Fornax Cluster and for the nuclei of dwarf elliptical galaxies in the Virgo Cluster. Ultracompact dwarfs are more extended and have higher surface brightnesses than typical dwarf nuclei, while the luminosities, colors, and sizes of the nuclei are closer to those of Galactic globular clusters. This calls into question the production of ultracompact dwarfs via threshing, whereby the lower surface brightness envelope of a dwarf elliptical galaxy is removed by tidal processes, leaving behind a bare nucleus. Threshing may still be a viable model if the relatively bright Fornax ultracompact dwarfs considered here are descended from dwarf elliptical galaxies whose nuclei are at the upper end of their luminosity and size distributions.

Afterglows, redshifts, and properties of swift gammaray bursts

We present optical, near-IR, and radio follow-up of 16 Swift bursts, including our discovery of nine afterglows and a redshift determination for three. These observations, supplemented by data from the literature, provide an afterglow recovery rate of 52% in the optical/near-IR, much higher than in previous missions (BeppoSAX, HETE-2, INTEGRAL, and IPN). The optical/near-IR afterglows of Swift events are on average 1.8 mag fainter at t = 12 hr than those of previous missions. The X-ray afterglows are similarly fainter than those of pre-Swift bursts. In the radio the limiting factor is the VLA threshold, and the detection rate for Swift bursts is similar to that for past missions. The redshift distribution of pre-Swift bursts peaked at z similar to 1, whereas the six Swift bursts with measured redshifts are distributed evenly between 0.7 and 3.2. From these results we conclude that ( 1) the pre-Swift distributions were biased in favor of bright events and low-redshift events, ( 2) the higher sensitivity and accurate positions of Swift result in a better representation of the true burst redshift and brightness distributions ( which are higher and dimmer, respectively), and (3) similar to 10% of the bursts are optically dark, as a result of a high redshift and/or dust extinction. We remark that the apparent lack of low-redshift, low-luminosity Swift bursts and the lower event rate than prelaunch estimates ( 90 vs. 150 per year) are the result of a threshold that is similar to that of BATSE. In view of these inferences, afterglow observers may find it advisable to make significant changes in follow-up strategies of Swift events. The faintness of the afterglows means that large telescopes should be employed as soon as the burst is localized. Sensitive observations in RIz and near-IR bands will be needed to discriminate between a typical z similar to 2 burst with modest extinction and a high-redshift event. Radio observations will be profitable for a small fraction (similar to 10%) of events. Finally, we suggest that a search for bright host galaxies in untriggered BAT localizations may increase the chance of finding nearby low-luminosity GRBs.

Moreton Bay, Queensland, Australia: An example of the co-existence of significant marine mammal populations and large-scale coastal development

Recent analyses assert that large marine vertebrates such as marine mammals are now 'functionally or entirely extinct in most coastal ecosystems'. Moreton Bay is a large diverse marine ecosystem bordering the fastest growing area in Australia. The human population is over 1.6 million and increasing yearly by between 10% and 13% with resultant impacts upon the adjoining marine environment. Nonetheless, significant populations of three species of marine mammals are resident within Moreton Bay and a further 14 species are seasonal or occasional visitors. This paper reviews the current and historical distributions and abundance of these species in the context of the current management regime and suggests initiatives to increase the resilience of marine mammal populations to the changes wrought by the burgeoning human population in coastal environments. (C) 2004 Elsevier Ltd. All rights reserved.

The 2dF-SDSS LRG and QSO (2SLAQ) luminous red galaxy survey

We present a spectroscopic survey of almost 15 000 candidate intermediate-redshift luminous red galaxies (LRGs) brighter than i = 19.8, observed with 2dF on the Anglo-Australian Telescope. The targets were selected photometrically from the Sloan Digital Sky Survey (SDSS) and lie along two narrow equatorial strips covering 180 deg(2). Reliable redshifts were obtained for 92 per cent of the targets and the selection is very efficient: over 90 per cent have 0.45 < z < 0.8. More than 80 per cent of the similar to 11 000 red galaxies have pure absorption-line spectra consistent with a passively evolving old stellar population. The redshift, photometric and spatial distributions of the LRGs are described. The 2SLAQ data will be released publicly from mid-2006, providing a powerful resource for observational cosmology and the study of galaxy evolution.

Issues of robustness and high dimensionality in cluster analysis

Finite mixture models are being increasingly used to model the distributions of a wide variety of random phenomena. While normal mixture models are often used to cluster data sets of continuous multivariate data, a more robust clustering can be obtained by considering the t mixture model-based approach. Mixtures of factor analyzers enable model-based density estimation to be undertaken for high-dimensional data where the number of observations n is very large relative to their dimension p. As the approach using the multivariate normal family of distributions is sensitive to outliers, it is more robust to adopt the multivariate t family for the component error and factor distributions. The computational aspects associated with robustness and high dimensionality in these approaches to cluster analysis are discussed and illustrated.

Close packed transitions in slit-shaped pores: Density functional theory study of methane adsorption capacity in carbon

An important feature of improving lattice gas models and classical isotherms is the incorporation of a pore size dependent capacity, which has hitherto been overlooked. In this paper, we develop a model for predicting the temperature dependent variation in capacity with pore size. The model is based on the analysis of a lattice gas model using a density functional theory approach at the close packed limit. Fluid-fluid and solid-fluid interactions are modeled by the Lennard-Jones 12-6 potential and Steele's 10-4-3, potential respectively. The capacity of methane in a slit-shaped carbon pore is calculated from the characteristic parameters of the unit cell, which are extracted by minimizing the grand potential of the unit cell. The capacities predicted by the proposed model are in good agreement with those obtained from grand canonical Monte Carlo simulation, for pores that can accommodate up to three adsorbed layers. Single particle and pair distributions exhibit characteristic features that correspond to the sequence of buckling and rhombic transitions that occur as the slit pore width is increased. The model provides a useful tool to model continuous variation in the microstructure of an adsorbed phase, namely buckling and rhombic transitions, with increasing pore width. (C) 2002 American Institute of Physics.

The estrogen receptor 1 G594A polymorphism is associated with migraine susceptibility in two independent case/control groups

Migraine is a painful and debilitating disorder with a significant genetic component. Steroid hormones, in particular estrogen, have long been considered to play a role in migraine, as variations in hormone levels are associated with migraine onset in many sufferers of the disorder. Steroid hormones mediate their activity via hormone receptors, which have a wide tissue distribution. Estrogen receptors have been localized to the brain in regions considered to be involved in migraine pathogenesis. Hence it is possible that genetic variation in the estrogen receptor gene may play a role in migraine susceptibility. This study thus examined the estrogen receptor 1 (ESRalpha) gene for a potential role in migraine pathogenesis and susceptibility. A population-based cohort of 224 migraine sufferers and 224 matched controls were genotyped for the G594A polymorphism located in exon 8 of the ESR1 gene. Statistical analysis indicated a significant difference between migraineurs and non-migraineurs in both the allele frequencies (P=0.003) and genotype distributions (P=0.008) in this sample. An independent follow-up study was then undertaken using this marker in an additional population-based cohort of 260 migraine sufferers and 260 matched controls. This resulted in a significant association between the two groups with regard to allele frequencies (P=8x10(-6)) and genotype distributions (P=4x10(-5)). Our findings support the hypothesis that genetic variation in hormone receptors, in particular the ESR1 gene, may play a role in migraine.

Rank test of location optimal for hyperbolic secant distribution

There are at least two reasons for a symmetric, unimodal, diffuse tailed hyperbolic secant distribution to be interesting in real-life applications. It displays one of the common types of non normality in natural data and is closely related to the logistic and Cauchy distributions that often arise in practice. To test the difference in location between two hyperbolic secant distributions, we develop a simple linear rank test with trigonometric scores. We investigate the small-sample and asymptotic properties of the test statistic and provide tables of the exact null distribution for small sample sizes. We compare the test to the Wilcoxon two-sample test and show that, although the asymptotic powers of the tests are comparable, the present test has certain practical advantages over the Wilcoxon test.

Projected Gross-Pitaevskii equation for harmonically confined Bose gases at finite temperature

We extend the projected Gross-Pitaevskii equation formalism of Davis [Phys. Rev. Lett. 87, 160402 (2001)] to the experimentally relevant case of thermal Bose gases in harmonic potentials and outline a robust and accurate numerical scheme that can efficiently simulate this system. We apply this method to investigate the equilibrium properties of the harmonically trapped three-dimensional projected Gross-Pitaevskii equation at finite temperature and consider the dependence of condensate fraction, position, and momentum distributions and density fluctuations on temperature. We apply the scheme to simulate an evaporative cooling process in which the preferential removal of high-energy particles leads to the growth of a Bose-Einstein condensate. We show that a condensate fraction can be inferred during the dynamics even in this nonequilibrium situation.

R-estimator of location of the generalized secant hyperbolic distribution

The generalized secant hyperbolic distribution (GSHD) proposed in Vaughan (2002) includes a wide range of unimodal symmetric distributions, with the Cauchy and uniform distributions being the limiting cases, and the logistic and hyperbolic secant distributions being special cases. The current article derives an asymptotically efficient rank estimator of the location parameter of the GSHD and suggests the corresponding one- and two-sample optimal rank tests. The rank estimator derived is compared to the modified MLE of location proposed in Vaughan (2002). By combining these two estimators, a computationally attractive method for constructing an exact confidence interval of the location parameter is developed. The statistical procedures introduced in the current article are illustrated by examples.

Accelerated leap methods for simulating discrete stochastic chemical kinetics

Biologists are increasingly conscious of the critical role that noise plays in cellular functions such as genetic regulation, often in connection with fluctuations in small numbers of key regulatory molecules. This has inspired the development of models that capture this fundamentally discrete and stochastic nature of cellular biology - most notably the Gillespie stochastic simulation algorithm (SSA). The SSA simulates a temporally homogeneous, discrete-state, continuous-time Markov process, and of course the corresponding probabilities and numbers of each molecular species must all remain positive. While accurately serving this purpose, the SSA can be computationally inefficient due to very small time stepping so faster approximations such as the Poisson and Binomial τ-leap methods have been suggested. This work places these leap methods in the context of numerical methods for the solution of stochastic differential equations (SDEs) driven by Poisson noise. This allows analogues of Euler-Maruyuma, Milstein and even higher order methods to be developed through the Itô-Taylor expansions as well as similar derivative-free Runge-Kutta approaches. Numerical results demonstrate that these novel methods compare favourably with existing techniques for simulating biochemical reactions by more accurately capturing crucial properties such as the mean and variance than existing methods.

Determination of pore size distributions by regularization and finite element collocation

The problem of extracting pore size distributions from characterization data is solved here with particular reference to adsorption. The technique developed is based on a finite element collocation discretization of the adsorption integral, with fitting of the isotherm data by least squares using regularization. A rapid and simple technique for ensuring non-negativity of the solutions is also developed which modifies the original solution having some negativity. The technique yields stable and converged solutions, and is implemented in a package RIDFEC. The package is demonstrated to be robust, yielding results which are less sensitive to experimental error than conventional methods, with fitting errors matching the known data error. It is shown that the choice of relative or absolute error norm in the least-squares analysis is best based on the kind of error in the data. (C) 1998 Elsevier Science Ltd. All rights reserved.