Statistical theory of finite Fermi systems with chaotic excited eigenstates

A theory of strongly interacting Fermi systems of a few particles is developed. At high excit at ion energies (a few times the single-parti cle level spacing) these systems are characterized by an extreme degree of complexity due to strong mixing of the shell-model-based many-part icle basis st at es by the residual two- body interaction. This regime can be described as many-body quantum chaos. Practically, it occurs when the excitation energy of the system is greater than a few single-particle level spacings near the Fermi energy. Physical examples of such systems are compound nuclei, heavy open shell atoms (e.g. rare earths) and multicharged ions, molecules, clusters and quantum dots in solids. The main quantity of the theory is the strength function which describes spreading of the eigenstates over many-part icle basis states (determinants) constructed using the shell-model orbital basis. A nonlinear equation for the strength function is derived, which enables one to describe the eigenstates without diagonalization of the Hamiltonian matrix. We show how to use this approach to calculate mean orbital occupation numbers and matrix elements between chaotic eigenstates and introduce typically statistical variable s such as t emperature in an isolated microscopic Fermi system of a few particles.

A low-energy core-collapse supernova without a hydrogen envelope

The final fate of massive stars depends on many factors. Theory suggests that some with initial masses greater than 25 to 30 solar masses end up as Wolf-Rayet stars, which are deficient in hydrogen in their outer layers because of mass loss through strong stellar winds. The most massive of these stars have cores which may form a black hole and theory predicts that the resulting explosion of some of them produces ejecta of low kinetic energy, a faint optical luminosity and a small mass fraction of radioactive nickel. An alternative origin for low-energy supernovae is the collapse of the oxygen-neon core of a star of 7-9 solar masses. No weak, hydrogen-deficient, core-collapse supernovae have hitherto been seen. Here we report that SN 2008ha is a faint hydrogen-poor supernova. We propose that other similar events have been observed but have been misclassified as peculiar thermonuclear supernovae (sometimes labelled SN 2002cx-like events). This discovery could link these faint supernovae to some long-duration gamma-ray bursts, because extremely faint, hydrogen-stripped core-collapse supernovae have been proposed to produce such long gamma-ray bursts, the afterglows of which do not show evidence of associated supernovae.

Collective excitations of trapped Bose-Einstein condensates in the energy and time domains

A time-dependent method for calculating the collective excitation frequencies and densities of a trapped, inhomogeneous Bose-Einstein condensate with circulation is presented. The results are compared with time-independent solutions of the Bogoliubov-de Gennes equations. The method is based on time-dependent linear-response theory combined with spectral analysis of moments of the excitation modes of interest. The technique is straightforward to apply, extremely efficient in our implementation with parallel fast Fourier transform methods, and produces highly accurate results. For high dimensionality or low symmetry the time-dependent approach is a more practical computational scheme and produces accurate and reliable data. The method is suitable for general trap geometries, condensate flows and condensates permeated with defects and vortex structures.

Approximating the Set of Local Minima in Partial RNA Folding Landscapes

Motivation: We study a stochastic method for approximating the set of local minima in partial RNA folding landscapes associated with a bounded-distance neighbourhood of folding conformations. The conformations are limited to RNA secondary structures without pseudoknots. The method aims at exploring partial energy landscapes pL induced by folding simulations and their underlying neighbourhood relations. It combines an approximation of the number of local optima devised by Garnier and Kallel (2002) with a run-time estimation for identifying sets of local optima established by Reeves and Eremeev (2004).

Results: The method is tested on nine sequences of length between 50 nt and 400 nt, which allows us to compare the results with data generated by RNAsubopt and subsequent barrier tree calculations. On the nine sequences, the method captures on average 92% of local minima with settings designed for a target of 95%. The run-time of the heuristic can be estimated by O(n2D?ln?), where n is the sequence length, ? is the number of local minima in the partial landscape pL under consideration and D is the maximum number of steepest descent steps in attraction basins associated with pL.

Acetaldehyde Production in the Direct Ethanol Fuel Cell: Mechanistic Elucidation by Density Functional Theory

This study employs density functional theory (DFT) calculations to examine the mechanism by which acetaldehyde is formed on platinum in a typical direct ethanol fuel cell (DEFC). A pathway is found involving the formation of a strongly hydrogen-bonded complex between adsorbed ethanol and the surface hydroxyl (OH) species, followed by the facile alpha-dehydrogenation of ethanol, with spontaneous weakening of the hydrogen bond in favor of adsorbed acetaldehyde and water. This mechanism is found to be comparably viable on both the close-packed surface and the monatomic steps. Comparison of further reactions on these two sites strongly indicates that the steps act as net removers of acetaldehyde from the product stream, while the flat surface acts as a net producer.

B-spline methods in R-matrix theory for scattering in two-electron systems

The use of B-spline basis sets in R-matrix theory for scattering processes has been investigated. In the present approach a B-spline basis is used for the description of the inner region, which is matched to the physical outgoing wavefunctions by the R-matrix. Using B-splines, continuum basis functions can be determined easily, while pseudostates can be included naturally. The accuracy for low-energy scattering processes is demonstrated by calculating inelastic scattering cross sections for e colliding on H. Very good agreement with other calculations has been obtained. Further extensions of the codes to quasi two-electron systems and general atoms are discussed as well as the application to (multi) photoionization.

Construction of Drug–Polymer Thermodynamic Phase Diagrams Using Flory–Huggins Interaction Theory: Identifying the Relevance of Temperature and Drug Weight Fraction to Phase Separation within Solid Dispersions

Amorphous drug-polymer solid dispersions have the potential to enhance the dissolution performance and thus bioavailability of BCS class II drug compounds. The principle drawback of this approach is the limited physical stability of amorphous drug within the dispersion. Accurate determination of the solubility and miscibility of drug in the polymer matrix is the key to the successful design and development of such systems. In this paper, we propose a novel method, based on Flory-Huggins theory, to predict and compare the solubility and miscibility of drug in polymeric systems. The systems chosen for this study are (1) hydroxypropyl methylcellulose acetate succinate HF grade (HPMCAS-HF)-felodipine (FD) and (2) Soluplus (a graft copolymer of polyvinyl caprolactam-polyvinyl acetate-polyethylene glycol)-FD. Samples containing different drug compositions were mixed, ball milled, and then analyzed by differential scanning calorimetry (DSC). The value of the drug-polymer interaction parameter ? was calculated from the crystalline drug melting depression data and extrapolated to lower temperatures. The interaction parameter ? was also calculated at 25 °C for both systems using the van Krevelen solubility parameter method. The rank order of interaction parameters of the two systems obtained at this temperature was comparable. Diagrams of drug-polymer temperature-composition and free energy of mixing (?G mix) were constructed for both systems. The maximum crystalline drug solubility and amorphous drug miscibility may be predicted based on the phase diagrams. Hyper-DSC was used to assess the validity of constructed phase diagrams by annealing solid dispersions at specific drug loadings. Three different samples for each polymer were selected to represent different regions within the phase diagram

Density Functional Theory Study on the Cleavage Mechanism of the Carbonyl Bond in Amides on Flat and Stepped Ru Surfaces: Hydrogen-Induced or Direct C-O Bond Breaking?

We have performed density functional theory (DFT) calculations to investigate the reaction mechanism of the cleavage of the carbonyl bond in amides on both flat and stepped Ru surfaces. The simplest amide molecule, N,N-dimethylacetamide (DMA), was used as the exemplar model molecule. Through the calculations, the most stable transition states (TSs) in all the pathways on both flat and stepped Ru surfaces are identified. Comparing the energy profiles of different reaction pathways, we find that a direct cleavage mechanism is always energetically favored as compared with an alternative hydrogen-induced mechanism on either the flat or stepped Ru surface. It is easier for the dissociation process to occur on the stepped surface than on the flat surface. However, as compared with the terrace, the superiority of step sites boosting the C-O bond dissociation is not as evident as that on CO dissociation.

Time-dependent density functional theory study of charge transfer in collisions

We study the charge transfer between colliding ions, atoms, or molecules, within time-dependent density functional theory. Two particular cases are presented, the collision between a proton and a Helium atom, and between a gold atom and a butane molecule. In the first case, proton kinetic energies between 16 keV and 1.2 MeV are considered, with impact parameters between 0.31 and 1.9 angstrom. The partial transfer of charge is monitored with time. The total cross-section is obtained as a function of the proton kinetic energy. In the second case, we analyze one trajectory and discuss spin-dependent charge transfer between the different fragments.

Electronic properties of interfaces and defects from many-body perturbation theory: Recent developments and applications

We review some recent developments in many body perturbation theory (MBPT) calculations that have enabled the study of interfaces and defects. Starting from the theoretical basis of MBPT, Hedin's equations are presented, leading to the CW and CWI' approximations. We introduce the perturbative approach, that is the one most commonly used for obtaining quasiparticle (QP) energies. The practical strategy presented for dealing with the frequency dependence of the self energy operator is based on either plasmon-pole models (PPM) or the contour deformation technique, with the latter being more accurate. We also discuss the extrapolar method for reducing the number of unoccupied states which need to be included explicity in the calculations. The use of the PAW method in the framework of MBPT is also described. Finally, results which have been obtained using, MBPT for band offsets a interfaces and for defects presented, with companies on the main difficulties and cancels.

Schematic representation of the QP corrections (marked with ) to the band edges (E and E-v) and a defect level (F) for a Si/SiO2 interface (Si and O atoms are represented in blue and red, respectively, in the ball and stick model) with an oxygen vacancy leading to a Si-Si bond (the Si atoms involved in this bond are colored light blue).

Electronic properties of zircon and hafnon from many-body perturbation theory

The electronic properties of zircon and hafnon, two wide-gap high-kappa materials, are investigated using many-body perturbation theory (MBPT) combined with the Wannier interpolation technique. For both materials, the calculated band structures differ from those obtained within density-functional theory and MBPT by (i) a slight displacement of the highest valence-band maximum from the Gamma point and (ii) an opening of the indirect band gap to 7.6 and 8.0 eV for zircon and hafnon, respectively. The introduction of vertex corrections in the many-body self-energy does not modify the results except for a global rigid shift of the many-body corrections.

Density functionals from many-body perturbation theory: The band gap for semiconductors and insulators

Theoretically the Kohn-Sham band gap differs from the exact quasiparticle energy gap by the derivative discontinuity of the exchange-correlation functional. In practice for semiconductors and insulators the band gap calculated within any local or semilocal density approximations underestimates severely the experimental energy gap. On the other hand, calculations with an "exact" exchange potential derived from many-body perturbation theory via the optimized effective potential suggest that improving the exchange-correlation potential approximation can yield a reasonable agreement between the Kohn-Sham band gap and the experimental gap. The results in this work show that this is not the case. In fact, we add to the exact exchange the correlation that corresponds to the dynamical (random phase approximation) screening in the GW approximation. This accurate exchange-correlation potential provides band structures similar to the local density approximation with the corresponding derivative discontinuity that contributes 30%-50% to the energy gap. Our self-consistent results confirm substantially the results for Si and other semiconductors obtained perturbatively [R. W. Godby , Phys. Rev. B 36, 6497 (1987)] and extend the conclusion to LiF and Ar, a wide-gap insulator and a noble-gas solid. (c) 2006 American Institute of Physics.

Surface energy and surface proton order of the ice Ih basal and prism surfaces

Density-functional theory (DFT) is used to examine the basal and prism surfaces of ice Ih. Similar surface energies are obtained for the two surfaces; however, in each case a strong dependence of the surface energy on surface proton order is identified. This dependence, which can be as much as 50% of the absolute surface energy, is significantly larger than the bulk dependence (< 1%) on proton order, suggesting that the thermodynamic ground state of the ice surface will remain proton ordered well above the bulk order-disorder temperature of about 72 K. On the basal surface this suggestion is supported by Monte Carlo simulations with an empirical potential and solution of a 2D Ising model with nearest neighbor interactions taken from DFT. Order parameters that define the surface energy of each surface in terms of nearest neighbor interactions between dangling OH bonds (those which point out of the surface into vacuum) have been identified and are discussed. Overall, these results suggest that proton order-disorder effects have a profound impact on the stability of ice surfaces and will most likely have an effect on ice surface reactivity as well as ice crystal growth and morphology. S Supplementary data are available from stacks.iop.org/JPhysCM/22/074209/mmedia

Testing metabolic scaling theory using intraspecific allometries in Antarctic microarthropods

Quantitative scaling relationships among body mass, temperature and metabolic rate of organisms are still controversial, while resolution may be further complicated through the use of different and possibly inappropriate approaches to statistical analysis. We propose the application of a modelling strategy based on the theoretical approach of Akaike's information criteria and non-linear model fitting (nlm). Accordingly, we collated and modelled available data at intraspecific level on the individual standard metabolic rate of Antarctic microarthropods as a function of body mass (M), temperature (T), species identity (S) and high rank taxa to which species belong (G) and tested predictions from metabolic scaling theory (mass-metabolism allometric exponent b = 0.75, activation energy range 0.2-1.2 eV). We also performed allometric analysis based on logarithmic transformations (lm). Conclusions from lm and nlm approaches were different. Best-supported models from lm incorporated T, M and S. The estimates of the allometric scaling exponent linking body mass and metabolic rate resulted in a value of 0.696 +/- 0.105 (mean +/- 95% CI). In contrast, the four best-supported nlm models suggested that both the scaling exponent and activation energy significantly vary across the high rank taxa (Collembola, Cryptostigmata, Mesostigmata and Prostigmata) to which species belong, with mean values of b ranging from about 0.6 to 0.8. We therefore reached two conclusions: 1, published analyses of arthropod metabolism based on logarithmic data may be biased by data transformation; 2, non-linear models applied to Antarctic microarthropod metabolic rate suggest that intraspecific scaling of standard metabolic rate in Antarctic microarthropods is highly variable and can be characterised by scaling exponents that greatly vary within taxa, which may have biased previous interspecific comparisons that neglected intraspecific variability.