Dynamical aspects of inextensible chains

In the present work the dynamics of a continuous inextensible chain is studied. The chain is regarded as a system of small particles subjected to constraints on their reciprocal distances. It is proposed a treatment of systems of this kind based on a set Langevin equations in which the noise is characterized by a non-gaussian probability distribution. The method is explained in the case of a freely hinged chain. In particular, the generating functional of the correlation functions of the relevant degrees of freedom which describe the conformations of this chain is derived. It is shown that in the continuous limit this generating functional coincides with a model of an inextensible chain previously discussed by one of the authors of this work. Next, the approach developed here is applied to a inextensible chain, called the freely jointed bar chain, in which the basic units are small extended objects. The generating functional of the freely jointed bar chain is constructed. It is shown that it differs profoundly from that of the freely hinged chain. Despite the differences, it is verified that in the continuous limit both generating functionals coincide as it is expected.Comment: 15 pages, LaTeX 2e + various packages, 3 figures. The title has been changed and three references have been added. A large part of the manuscript has been rewritten to improve readability. Chapter 4 has been added. It contains the construction of the generating functional without the shish-kebab approximation and a new derivation of the continuous limit of the freely jointed bar chai

Van der Waals forces in density functional theory: perturbational long-range electron interaction corrections

Long-range exchange and correlation effects, responsible for the failure of currently used approximate density functionals in describing van der Waals forces, are taken into account explicitly after a separation of the electron-electron interaction in the Hamiltonian into short- and long-range components. We propose a "range-separated hybrid" functional based on a local density approximation for the short-range exchange-correlation energy, combined with a long-range exact exchange energy. Long-range correlation effects are added by a second-order perturbational treatment. The resulting scheme is general and is particularly well-adapted to describe van der Waals complexes, like rare gas dimers.Comment: 8 pages, 1 figure, submitted to Phys. Rev.

Theoretical study of O adlayers on Ru(0001)

Recent experiments performed at high pressures indicate that ruthenium can support unusually high concentrations of oxygen at the surface. To investigate the structure and stability of high coverage oxygen structures, we performed density functional theory calculations, within the generalized gradient approximation, for O adlayers on Ru(0001) from low coverage up to a full monolayer. We achieve quantitative agreement with previous low energy electron diffraction intensity analyses for the (2x2) and (2x1) phases and predict that an O adlayer with a (1x1) periodicity and coverage of 1 monolayer can form on Ru(0001), where the O adatoms occupy hcp-hollow sites.Comment: RevTeX, 6 pages, 4 figure

Thermochemistry of Alane Complexes for Hydrogen Storage: A Theoretical and Experimental Comparison

Knowledge of the relative stabilities of alane (AlH3) complexes with electron donors is essential for identifying hydrogen storage materials for vehicular applications that can be regenerated by off-board methods; however, almost no thermodynamic data are available to make this assessment. To fill this gap, we employed the G4(MP2) method to determine heats of formation, entropies, and Gibbs free energies of formation for thirty-eight alane complexes with NH3-nRn (R = Me, Et; n = 0-3), pyridine, pyrazine, triethylenediamine (TEDA), quinuclidine, OH2-nRn (R = Me, Et; n = 0-2), dioxane, and tetrahydrofuran (THF). Monomer, bis, and selected dimer complex geometries were considered. Using these data, we computed the thermodynamics of the key formation and dehydrogenation reactions that would occur during hydrogen delivery and alane regeneration, from which trends in complex stability were identified. These predictions were tested by synthesizing six amine-alane complexes involving trimethylamine, triethylamine, dimethylethylamine, TEDA, quinuclidine, and hexamine, and obtaining upper limits of delta G for their formation from metallic aluminum. Combining these computational and experimental results, we establish a criterion for complex stability relevant to hydrogen storage that can be used to assess potential ligands prior to attempting synthesis of the alane complex. Based on this, we conclude that only a subset of the tertiary amine complexes considered and none of the ether complexes can be successfully formed by direct reaction with aluminum and regenerated in an alane-based hydrogen storage system.Comment: Accepted by the Journal of Physical Chemistry

Warming Up Density Functional Theory

Density functional theory (DFT) has become the most popular approach to electronic structure across disciplines, especially in material and chemical sciences. Last year, at least 30,000 papers used DFT to make useful predictions or give insight into an enormous diversity of scientific problems, ranging from battery development to solar cell efficiency and far beyond. The success of this field has been driven by usefully accurate approximations based on known exact conditions and careful testing and validation. In the last decade, applications of DFT in a new area, warm dense matter, have exploded. DFT is revolutionizing simulations of warm dense matter including applications in controlled fusion, planetary interiors, and other areas of high energy density physics. Over the past decade or so, molecular dynamics calculations driven by modern density functional theory have played a crucial role in bringing chemical realism to these applications, often (but not always) with excellent agreement with experiment. This chapter summarizes recent work from our group on density functional theory at non-zero temperatures, which we call thermal DFT. We explain the relevance of this work in the context of warm dense matter, and the importance of quantum chemistry to this regime. We illustrate many basic concepts on a simple model system, the asymmetric Hubbard dimer

Surface reconstruction induced geometries of Si clusters

We discuss a generalization of the surface reconstruction arguments for the structure of intermediate size Si clusters, which leads to model geometries for the sizes 33, 39 (two isomers), 45 (two isomers), 49 (two isomers), 57 and 61 (two isomers). The common feature in all these models is a structure that closely resembles the most stable reconstruction of Si surfaces, surrounding a core of bulk-like tetrahedrally bonded atoms. We investigate the energetics and the electronic structure of these models through first-principles density functional theory calculations. These models may be useful in understanding experimental results on the reactivity of Si clusters and their shape as inferred from mobility measurements.Comment: 9 figures (available from the author upon request) Submitted to Phys. Rev.

Metal worker’s lung: spatial association with Mycobacterium avium

Background Outbreaks of hypersensitivity pneumonitis (HP) are not uncommon in workplaces where metal working fluid (MWF) is used to facilitate metal turning. Inhalation of microbe-contaminated MWF has been assumed to be the cause, but previous investigations have failed to establish a spatial relationship between a contaminated source and an outbreak. Objectives After an outbreak of five cases of HP in a UK factory, we carried out blinded, molecular-based microbiological investigation of MWF samples in order to identify potential links between specific microbial taxa and machines in the outbreak zone. Methods Custom-quantitative PCR assays, microscopy and phylogenetic analyses were performed on blinded MWF samples to quantify microbial burden and identify potential aetiological agents of HP in metal workers. Measurements and main results MWF from machines fed by a central sump, but not those with an isolated supply, was contaminated by mycobacteria. The factory sump and a single linked machine at the centre of the outbreak zone, known to be the workstation of the index cases, had very high levels of detectable organisms. Phylogenetic placement of mycobacterial taxonomic marker genes generated from these samples indicated that the contaminating organisms were closely related to Mycobacterium avium. Conclusions We describe, for the first time, a close spatial relationship between the abundance of a mycobacterium-like organism, most probably M. avium, and a localised outbreak of MWF-associated HP. The further development of sequence-based analytic techniques should assist in the prevention of this important occupational disease

Tomographic Separation of Composite Spectra. VIII. The Physical Properties of the Massive Compact Binary in the Triple Star System HD 36486 (delta Orionis A)

Double-lined spectroscopic orbital elements have recently been found for the central binary in the massive triple, delta Orionis A based on radial velocities from cross-correlation techniques applied to IUE high dispersion spectra and He I 6678 spectra obtained at Kitt Peak. The primary and secondary velocity amplitudes were found to be 94.9 +/- 0.6 km/s and 186 +/- 9 km/s respectively. Tomographic reconstructions of the primary and secondary stars' spectra confirm the O9.5 II classification of the primary and indicate a B0.5 III type for the secondary. The widths of the UV cross-correlation functions are used to estimate the projected rotational velocities, Vsin i = 157 +/- 6 km/s and 138 +/- 16 km/s for the primary and secondary, respectively implying that both stars rotate faster than their orbital motion. We used the spectroscopic results to make a constrained fit of the Hipparcos light curve of this eclipsing binary, and the model fits limit the inclination to the range between 67 and 77 degrees. The i = 67 degrees solution, which corresponds to a near Roche-filling configuration, results in a primary mass of 11.2 solar masses and a secondary mass of 5.6 solar masses, both of which are substantially below the expected masses for stars of their luminosity. This binary may have experienced a mass ratio reversal caused by Case A Roche lobe overflow, or the system may have suffered extensive mass loss through a binary interaction, perhaps during a common envelope phase, in which most of the primary's mass was lost from the system rather than transferred to the secondary.Comment: 27 pages, 15 figures in press, the Astrophysical Journal, February 1, 200

Electron affinities of the first- and second- row atoms: benchmark ab initio and density functional calculations

A benchmark ab initio and density functional (DFT) study has been carried out on the electron affinities of the first- and second-row atoms. The ab initio study involves basis sets of $spdfgh$ and $spdfghi$ quality, extrapolations to the 1-particle basis set limit, and a combination of the CCSD(T), CCSDT, and full CI electron correlation methods. Scalar relativistic and spin-orbit coupling effects were taken into account. On average, the best ab initio results agree to better than 0.001 eV with the most recent experimental results. Correcting for imperfections in the CCSD(T) method improves the mean absolute error by an order of magnitude, while for accurate results on the second-row atoms inclusion of relativistic corrections is essential. The latter are significantly overestimated at the SCF level; for accurate spin-orbit splitting constants of second-row atoms inclusion of (2s,2p) correlation is essential. In the DFT calculations it is found that results for the 1st-row atoms are very sensitive to the exchange functional, while those for second-row atoms are rather more sensitive to the correlation functional. While the LYP correlation functional works best for first-row atoms, its PW91 counterpart appears to be preferable for second-row atoms. Among ``pure DFT'' (nonhybrid) functionals, G96PW91 (Gill 1996 exchange combined with Perdew-Wang 1991 correlation) puts in the best overall performance. The best results overall are obtained with the 1-parameter hybrid modified Perdew-Wang (mPW1) exchange functionals of Adamo and Barone [J. Chem. Phys. {\bf 108}, 664 (1998)], with mPW1LYP yielding the best results for first-row, and mPW1PW91 for second-row atoms. Indications exist that a hybrid of the type $a$ mPW1LYP + $(1-a)$ mPW1PW91 yields better results than either of the constituent functionals.Comment: Phys. Rev. A, in press (revised version, review of issues concerning DFT and electron affinities added