The Application of Computational Chemistry and Chemometrics to Developing a Method for Online Monitoring of Polymer Degradation in the Manufacture of Bioresorbable Medical Implants

Bioresorbable polymers such as PLA have an important role to play in the development of temporary implantable medical devices with significant benefits over traditional therapies. However, development of new devices is hindered by high manufacturing costs associated with difficulties in processing the material. A major problem is the lack of insight on material degradation during processing. In this work, a method of quantifying degradation of PLA using IR spectroscopy coupled with computational chemistry and chemometric modeling is examined. It is shown that the method can predict the quantity of degradation products in solid-state samples with reasonably good accuracy, indicating the potential to adapt the method to developing an on-line sensor for monitoring PLA degradation in real-time during processing.

Electron shell contributions to gamma-ray spectra of positron annihilation in noble gases

Gamma-ray positron annihilation spectra of the noble gases are simulated using computational chemistry tools for the bound electron wavefunctions and plane-wave approximation for the low-energy positron. The present annihilation line shapes, i.e. the full width at half maximum, Delta epsilon, of the gamma-ray annihilation spectra for He and Ar (valence) agree well with available independent atomic calculations using a different algorithm. For other noble gases they achieve moderate agreement with the experimental measurements. It is found that the contributions of various atomic electron shells to the spectra depend significantly on their principal quantum number n and orbital angular momentum quantum number l. The present study further reveals that the outermost ns electrons of the noble gases exhibit spectral line shapes in close agreement with those measured, indicating (as expected) that the measurements are not due to a simple sum over the momentum densities for all atomic electrons. The robust nature of the present approach makes it possible for us to proceed to more complex molecular systems using the tools of modern computational chemistry.

The structural information content of chemical networks

We present an information-theoretic method to measure the structural information content of networks and apply it to chemical graphs. As a result, we find that our entropy measure is more general than classical information indices known in mathematical and computational chemistry. Further, we demonstrate that our measure reflects the essence of molecular branching meaningfully by determining the structural information content of some chemical graphs numerically.

Electronic states and spin-forbidden cooling transitions of AlH and AlF

The feasibility of laser cooling AlH and AlF is investigated using ab initio quantum chemistry. All the electronic states corresponding to the ground and lowest two excited states of the Al atom are calculated using multi-reference configuration interaction (MRCI) and the large AV6Z basis set for AlH. The smaller AVQZ basis set is used to calculate the valence electronic states of AlF. Theoretical Franck-Condon factors are determined for the A(1)Pi -> X(1)Sigma(+) transitions in both radicals and found to agree with the highly diagonal factors found experimentally, suggesting computational chemistry is an effective method for screening suitable laser cooling candidates. AlH does not appear to have a transition quite as diagonal as that in SrF (which has been laser cooled) but the A(1)Pi -> X(1)Sigma(+) transition transition of AlF is a strong candidate for cooling with just a single laser, though the cooling frequency is deep in the UV. Furthermore, the a (3)Pi -> X(1)Sigma(+) transitions are also strongly diagonal and in AlF is a practical method for obtaining very low final temperatures around 3 mu K.

STO and GTO field-induced polarization functions for H to Kr

Field-induced polarization (FIP) functions were proposed over two decades ago to improve the accuracy of calculated response properties, and the FIP functions in GTO form for H and C to F were tested on small molecules, with encouraging results. The concept of FIP,is now extended to all atoms up to Kr. New simplifying approximations for the description of asymptotic highest occupied atomic orbitals. (HOAOs) are introduced in this study. They provide the basis for STO and GTO exponents of a complete set of FIP functions from H to Kr, which are both listed for the convenience of the users. Tests on the polarizabilities of a series of atoms and molecules demonstrate that addition of the FIP basis functions to a series' of standard basis sets drastically improves the performance of all these basis sets compared to converged results. Moreover, the byproduct of this study (approximate asymptotic HOAOs) provides information for the construction of accurate basis sets for long-range ground state properties. (C) 2003 Wiley Periodicals, Inc.

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.

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.

Solid-state analysis of low-melting 1,3-dialkylimidazolium hexafluorophosphate salts (ionic liquids) by combined x-ray crystallographic and computational analyses

The interactions of ions in the solid state for a series of representative 1,3-dialkylimidazolium hexafluorophosphate salts (either ionic liquids or closely related) have been examined by crystallographic analysis, combined with the theoretical estimation of crystal-packing densities and lattice-interaction energies. Efficient close-packing of the ions in the crystalline states is observed, but there was no compelling evidence for specific directional hydrogen-bonding to the hexafluorophosphate anions or the formation of interstitial voids. The close-packing efficiency is supported by the theoretical calculation of ion volumes, crystal lattice energies, and packing densities, which correlated well with experimental data. The crystal density of the salts can be predicted accurately from the summation of free ion volumes and lattice energies calculated. Of even more importance for future work, on these and related salts, the solid-state density of 1,3-dialkylimidazolium hexafluorophosphate salts can be predicted with reasonable accuracy purely on the basis of on ab initio free ion volumes, and this allows prediction of lattice energies without necessarily requiring the crystal structures.

Computational Verification of Two Universal Relations for Simple Ionic Liquids. Kinetic Properties of a Model 2:1 Molten Salt

Two semianalytical relations [Nature, 1996, 381, 137 and Phys. Rev. Lett. 2001, 87, 245901] predicting dynamical coefficients of simple liquids on the basis of structural properties have been tested by extensive molecular dynamics simulations for an idealized 2:1 model molten salt. In agreement with previous simulation studies, our results support the validity of the relation expressing the self-diffusion coefficient as a Function of the radial distribution functions for all thermodynamic conditions such that the system is in the ionic (ie., fully dissociated) liquid state. Deviations are apparent for high-density samples in the amorphous state and in the low-density, low-temperature range, when ions condense into AB(2) molecules. A similar relation predicting the ionic conductivity is only partially validated by our data. The simulation results, covering 210 distinct thermodynamic states, represent an extended database to tune and validate semianalytical theories of dynamical properties and provide a baseline for the interpretation of properties of more complex systems such as the room-temperature ionic liquids.

A computational study of the mechanism of CO oxidation by a ceria supported surface rhodium oxide layer

A mechanism of CO oxidation by a thin surface oxide of Rh supported on ceria is proposed: CO is oxidized by the Rh-oxide film, which is subsequently reoxidized by a ceria surface O atom. The proposed mechanism is supported by in situ Raman spectroscopic investigations.

Comparison of dynamics of wildtype and V94M human UDP-galactose 4-epimerase-A computational perspective on severe epimerase-deficiency galactosemia

UDP-galactose 4'-epimerase (GALE) catalyzes the interconversion of UDP-galactose and UDP-glucose, an important step in galactose catabolism. Type III galactosemia, an inherited metabolic disease, is associated with mutations in human GALE. The V94M mutation has been associated with a very severe form of type III galactosemia. While a variety of structural and biochemical studies have been reported that elucidate differences between the wildtype and this mutant form of human GALE, little is known about the dynamics of the protein and how mutations influence structure and function. We performed molecular dynamics simulations on the wildtype and V94M enzyme in different states of substrate and cofactor binding. In the mutant, the average distance between the substrate and both a key catalytic residue (Tyr157) and the enzyme-bound NAD(+) cofactor and the active site dynamics are altered making substrate binding slightly less stable. However, overall stability or dynamics of the protein is not altered. This is consistent with experimental findings that the impact is largely on the turnover number (kcat), with less substantial effects on Km. Active site fluctuations were found to be correlated in enzyme with substrate bound to just one of the subunits in the homodimer suggesting inter-subunit communication. Greater active site loop mobility in human GALE compared to the equivalent loop in Escherichia coli GALE explains why the former can catalyze the interconversion of UDP-N-acetylgalactosamine and UDP-N-acetylglucosamine while the bacterial enzyme cannot. This work illuminates molecular mechanisms of disease and may inform the design of small molecule therapies for type III galactosemia.