Molecular Computational Elements Encode Large Populations of Small Objects

Since the introduction of molecular computation1, 2, experimental molecular computational elements have grown3, 4, 5 to encompass small-scale integration6, arithmetic7 and games8, among others. However, the need for a practical application has been pressing. Here we present molecular computational identification (MCID), a demonstration that molecular logic and computation can be applied to a widely relevant issue. Examples of populations that need encoding in the microscopic world are cells in diagnostics or beads in combinatorial chemistry (tags). Taking advantage of the small size9 (about 1 nm) and large 'on/off' output ratios of molecular logic gates and using the great variety of logic types, input chemical combinations, switching thresholds and even gate arrays in addition to colours, we produce unique identifiers for members of populations of small polymer beads (about 100 m) used for synthesis of combinatorial libraries10, 11. Many millions of distinguishable tags become available. This method should be extensible to far smaller objects, with the only requirement being a 'wash and watch' protocol12. Our focus on converting molecular science into technology concerning analog sensors13, 14, turns to digital logic devices in the present work.

Accurate computational methods for two-electron atom-laser interactions

We discuss the application of quantitatively accurate computational methods to the study of laser-driven two-electron atoms in short intense laser pulses. The fundamental importance of such calculations to the subject area is emphasized. Calculations of single- and double-electron ionization rates at 390 nm are presented. (C) 2001 Optical Society of America.

Computational challenges in atomic, molecular and optical physics

Six challenges are discussed. These are the laser-driven helium atom; the laser-driven hydrogen molecule and hydrogen molecular ion: electron scattering (with ionization) from one-electron atoms; the vibrational and rotational structure of molecules such as H-3(+) and water at their dissociation limits; laser- heated clusters; and quantum degeneracy and Bose-Einstein condensation. The first four concern fundamental few-body systems where use of high-performance computing (HPC) is currently making possible accurate modelling from first principles. This leads to reliable predictions and support for laboratory experiment as well as true understanding of the dynamics. Important aspects of these challenges addressable only via a terascale facility are set out. Such a facility makes the last two challenges in the above list meaningfully accessible for the first time, and the scientific interest together with the prospective role for HPC in these is emphasized.

High-pressure deformation mechanism in scolecite: A combined computational-experimental study

Pressure-induced structural modifications in scolecite were studied by means of in situ synchrotron X-ray powder diffraction and density functional computations. The experimental cell parameters were refined up to 8.5 GPa. Discontinuities in the slope of the unit-cell parameters vs. pressure dependence were observed; as a consequence, an increase in the slope of the linear pressure-volume dependence is observed at about 6 GPa, suggesting an enhanced compressibility at higher pressures. Weakening and broadening of the diffraction peaks reveals increasing structural disorder with pressure, preventing refinement of the lattice parameters above 8.5 GPa. Diffraction patterns collected during decompression show that the disorder is irreversible. Atomic coordinates within unit cells of different dimensions were determined by means of Car-Parrinello simulations. The discontinuous rise in compressibility at about 6 GPa is reproduced by the computation, allowing us to attribute it to re-organization of the hydrogen bonding network, with the formation of water dimers. Moreover we found that, with increasing pressure, the tetrahedral chains parallel to c rotate along their elongation axis and display an increasing twisting along a direction perpendicular to c. At the same time, we observed the compression of the channels. We discuss the modification of the Ca polyhedra under pressure, and the increase in coordination number (from 4 to 5) of one of the two Al atoms, resulting from the approach of a water molecule. We speculate that this last transformation triggers the irreversible disordering of the system.

Zipf and Type-Token rules for the English, Spanish, Irish and Latin languages

The Zipf curves of log of frequency against log of rank for a large English corpus of 500 million word tokens, 689,000 word types and for a large Spanish corpus of 16 million word tokens, 139,000 word types are shown to have the usual slope close to –1 for rank less than 5,000, but then for a higher rank they turn to give a slope close to –2. This is apparently mainly due to foreign words and place names. Other Zipf curves for highlyinflected Indo-European languages, Irish and ancient Latin, are also given. Because of the larger number of word types per lemma, they remain flatter than the English curve maintaining a slope of –1 until turning points of about ranks 30,000 for Irish and 10,000 for Latin. A formula which calculates the number of tokens given the number of types is derived in terms of the rank at the turning point, 5,000 for both English and Spanish, 30,000 for Irish and 10,000 for Latin.

Electron collisions with iron peak elements: a computational grand challenge

This paper discusses the calculation of electron impact collision strengths and effective collision strengths for iron peak elements of importance in the analysis of many astronomical and laboratory spectra. It commences with a brief overview of R-matrix theory which is the basis of computer programs which have been widely used to calculate the relevant atomic data used in this analysis. A summary is then given of calculations carried out over the last 20 y for electron collisions with Fe II. The grand challenge, represented by the calculation of accurate collision strengths and effective collision strengths for this ion, is then discussed. A new parallel R-matrix program PRMAT, which is being developed to meet this challenge, is then described and results of recent calculations, using this program to determine optically forbidden transitions in e- – Ni IV on a Cray T3E-1200 parallel supercomputer, are presented. The implications of this e- – Ni IV calculation for the determination of accurate data from an isoelectronic e- – Fe II calculation are discussed and finally some future directions of research are reviewed.

Ion Association in [bmim][PF6]/Naphthalene Mixtures: An Experimental and Computational Study

Mixtures of room temperature ionic liquids (IL) with neutral organic molecules provide a valuable testing ground to investigate the interplay of the ionic and molecular-dipolar state in dense Coulomb systems at near ambient conditions. In the present study, the viscosity eta and the ionic conductivity a of 1-n-butyl-3-methylimidazolium hexafluorophosphate ([bmim][PF6])/naphthalene mixtures at T = 80 degrees C have been measured at 10 stoichiometries spanning the composition range from pure naphthalene to pure [bmim][PF6]. The viscosity grows nearly monotonically with increasing IL mole fraction (x), whereas the conductivity per ion displays a clear peak at x approximate to 15%. The origin of this maximum has been investigated using molecular dynamics simulations based on a classical force field. Snapshots of the simulated samples show that the conductivity maximum is due to the gradual transition in the IL component from an ionic state at high x to a dipolar fluid made of neutral ion pairs at low x. At concentrations x <0.20 the ion pairs condense into molecular-thin filaments bound by dipolar forces and extending in between nanometric droplets of IL. These results are confirmed and complemented by the computation of dynamic and transport properties in [bmim][PF6]/naphthalene mixtures at low IL concentration.

A computational model for the cooling phase of injection moulding

This paper discusses the approaches and techniques used to build a realistic numerical model to analyse the cooling phase of the injection moulding process. The procedures employed to select an appropriate mesh and the boundary and initial conditions for the problem are discussed and justified. The final model is validated using direct comparisons with experimental results generated in an earlier study. The model is shown to be a useful tool for further studies aimed at optimising the cooling phase of the injection moulding process. Using the numerical model provides additional information relating to changes in conditions throughout the process, which otherwise could not be deduced or assessed experimentally. These results, and other benefits related to the use of the model, are also discussed in the paper. © 2007 Elsevier B.V. All rights reserved.

The energetics of tetrahydrocarbazole aromatization over Pd(111): A computational analysis

The carbazole moiety is a component of many important pharmaceuticals including anticancer and anti-HIV agents and is commonly utilized in the production of modern polymeric materials with novel photophysical and electronic properties. Simple carbazoles are generally produced via the aromatization of the respective tetrahydrocarbazole (THCZ). In this work, density functional theory calculations are used to model the reaction pathway of tetrahydrocarbazole aromatization over Pd(111). The geometry of each of the intermediate surface species has been determined and how each structure interacts with the metal surface addressed. The reaction energies and barriers of each of the elementary surface reactions have also been calculated, and a detailed analysis of the energetic trends performed. Our calculations have shown that the surface intermediates remain fixed to the surface via the aromatic ring in a manner similar to that of THCZ. Moreover, the aliphatic ring becomes progressively more planer with the dissociation of each subsequent hydrogen atom. Analysis of the reaction energy profile has revealed that the trend in reaction barriers is determined by the two factors: (i) the strength of the dissociating ring-H bond and (ii) the subsequent gain in energy due to the geometric relaxation of the aliphatic ring. (c) 2008 American Institute of Physics.