Functional evaluation of noncovalent interactions in neuroreceptors and progress toward the expansion of unnatural amino acid methodology

This dissertation primarily describes chemical-scale studies of G protein-coupled receptors and Cys-loop ligand-gated ion channels to better understand ligand binding interactions and the mechanism of channel activation using recently published crystal structures as a guide. These studies employ the use of unnatural amino acid mutagenesis and electrophysiology to measure subtle changes in receptor function.

In chapter 2, the role of a conserved aromatic microdomain predicted in the D3 dopamine receptor is probed in the closely related D2 and D4 dopamine receptors. This domain was found to act as a structural unit near the ligand binding site that is important for receptor function. The domain consists of several functionally important noncovalent interactions including hydrogen bond, aromatic-aromatic, and sulfur-π interactions that show strong couplings by mutant cycle analysis. We also assign an alternate interpretation for the linear fluorination plot observed at W6.48, a residue previously thought to participate in a cation-π interaction with dopamine.

Chapter 3 outlines attempts to incorporate chemically synthesized and in vitro acylated unnatural amino acids into mammalian cells. While our attempts were not successful, method optimizations and data for nonsense suppression with an in vivo acylated tRNA are included. This chapter is aimed to aid future researchers attempting unnatural amino acid mutagenesis in mammalian cells.

Chapter 4 identifies a cation-π interaction between glutamate and a tyrosine residue on loop C in the GluClβ receptor. Using the recently published crystal structure of the homologous GluClα receptor, other ligand-binding and protein-protein interactions are probed to determine the similarity between this invertebrate receptor and other more distantly related vertebrate Cys-loop receptors. We find that many of the interactions previously observed are conserved in the GluCl receptors, however care must be taken when extrapolating structural data.

Chapter 5 examines inherent properties of the GluClα receptor that are responsible for the observed glutamate insensitivity of the receptor. Chimera synthesis and mutagenesis reveal the C-terminal portion of the M4 helix and the C-terminus as contributing to formation of the decoupled state, where ligand binding is incapable of triggering channel gating. Receptor mutagenesis was unable to identify single residue mismatches or impaired protein-protein interactions within this domain. We conclude that M4 helix structure and/or membrane dynamics are likely the cause of ligand insensitivity in this receptor and that the M4 helix has an role important in the activation process.

Design, synthesis, and biological activity of rhodium metalloinsertors

Deficiencies in the mismatch repair (MMR) pathway are associated with several types of cancers, as well as resistance to commonly used chemotherapeutics. Rhodium metalloinsertors have been found to bind DNA mismatches with high affinity and specificity in vitro, and also exhibit cell-selective cytotoxicity, targeting MMR-deficient cells over MMR-proficient cells.

Here we examine the biological fate of rhodium metalloinsertors bearing dipyridylamine ancillary ligands. These complexes are shown to exhibit accelerated cellular uptake which permits the observation of various cellular responses, including disruption of the cell cycle and induction of necrosis, which occur preferentially in the MMR-deficient cell line. These cellular responses provide insight into the mechanisms underlying the selective activity of this novel class of targeted anti-cancer agents.

In addition, ten distinct metalloinsertors with varying lipophilicities are synthesized and their mismatch binding affinities and biological activities studied. While they are found to have similar binding affinities, their cell-selective antiproliferative and cytotoxic activities vary significantly. Inductively coupled plasma mass spectrometry (ICP-MS) experiments show that all of these metalloinsertors localize in the nucleus at sufficient concentrations for binding to DNA mismatches. Furthermore, metalloinsertors with high rhodium localization in the mitochondria show toxicity that is not selective for MMR-deficient cells. This work supports the notion that specific targeting of the metalloinsertors to nuclear DNA gives rise to their cytotoxic and antiproliferative activities that are selective for cells deficient in MMR.

To explore further the basis of the unique selectivity of the metlloinsertors in targeting MMR-deficient cells, experiments were conducted using engineered NCI-H23 lung adenocarcinoma cells that contain a doxycycline-inducible shRNA which suppresses the expression of the MMR gene MLH1. Here we use this new cell line to further validate rhodium metalloinsertors as compounds capable of differentially inhibiting the proliferation of MMR-deficient cancer cells over isogenic MMR-proficient cells. General DNA damaging agents, such as cisplatin and etoposide, in contrast, are less effective in the induced cell line defective in MMR.

Finally, we describe a new subclass of metalloinsertors with enhanced potency and selectivity, in which the complexes show Rh-O coordination. In particular, it has been found that both Δ and Λ enantiomers of [Rh(chrysi)(phen)(DPE)]²⁺ bind to DNA with similar affinities, suggesting a possible different binding conformation than previous metalloinsertors. Remarkably, all members of this new family of compounds have significantly increased potency in a range of cellular assays; indeed, all are more potent than the FDA-approved anticancer drugs cisplatin and MNNG. Moreover, these activities are coupled with high levels of selectivity for MMR-deficient cells.

A simple approach for the estimation of food consumption from growth rates at different environmental conditions and its application to juvenile cod (Gadus morhua L.) of a fjordic sea loch on the west coast of Scotland

Abstract Growth and condition of fish are functions of available food and environmental conditions. This led to the idea of using fish as a “consumption sensor” for the measurement of food intake over a defined period of time. A bio-physical model for the estimation of food consumption was developed based on the von Bertalanffy model. Whereas some of the input variables of the model, the initial and final lengths and masses of a fish and the temperature within the time period considered can easily be measured, internal characteristics of the species have to be determined indirectly. Three internal parameters are used in the model: the maintenance consumption at 0°C, the temperature dependence of this consumption and the food efficiency, the percentage of the ingested food utilized. Estimates of the parameters for a given species can be determined by feeding experiments. Here, data from published feeding experiments on juvenile cod, Gadus morhua L., were used to validate the model. The average of the relative error for the food intake predicted by the model for individual fish was about 24 %, indicating that fish used the food with different efficiencies. However, grouping the fish according to size classes and temperature lowered the relative error of the predicted food intake for the group to typically 5 %. For a group containing all fish of the feeding experiment the relative prediction error was about 2 %. Zusammenfassung Wachstum und Kondition der Fische sind von der verfügbaren Nahrung und von Umweltbedingungen abhängig. Dies führte zur Idee, Fisch als „Konsum-Sensor“ für die Messung der Nahrungsaufnahme über einen definierten Zeitraum zu verwenden. Auf Grundlage des von Bertalanffy-Modells wurde ein bio-physikalisches Modell zur Schätzung der Futteraufnahme entwickelt. Während einige der Eingangsgrößen des Modells leicht gemessen werden können (Anfangs- und Endlänge und -körpermasse der Fische und die Temperatur innerhalb des betrachteten Zeitraum), können interne Parameter der betrachteten Art nur indirekt bestimmt werden. Drei interne Parameter werden in dem Modell verwendet: Die Erhaltungskonsumtion bei 0° C, die Temperaturabhängigkeit dieser Rate und der Wirkungsgrad der Nahrung (der Anteil der Nahrung ,der aufgenommen und verwendet und nicht ungenutzt wieder ausgeschieden wird). Die Modellparameter für eine bestimmte Art können durch Fütterungsversuche bestimmt werden. Um das Modell zu validieren wurden Daten aus veröffentlichten Fütterungsversuchen mit juvenilen Kabeljau (Gadus morhua L.) verwendet. Modell und Wirklichkeit weichen in der Regel voneinander ab. Der durchschnittliche relative Fehler der durch das Modell vorhergesagten Nahrungsaufnahme betrug für Einzelfische etwa 24%, was darauf hinweist, dass einzelne Fisch die Nahrung mit unterschiedlichen Wirkungsgraden verwerten. Allerdings senkte die Gruppierung der Fische nach Größenklassen und Temperatur den relativen Vorhersagefehler für die Nahrungsaufnahme der Gruppe auf etwa 5%. Für alle Fische im Fütterungsversuch ist der relative Vorhersagefehler etwa 2%.

Strength of tantalum at high pressures through Richtmyer-Meshkov laser compression experiments and simulations

Strength at extreme pressures (>1 Mbar or 100 GPa) and high strain rates (106-108 s-1) of materials is not well characterized. The goal of the research outlined in this thesis is to study the strength of tantalum (Ta) at these conditions. The Omega Laser in the Laboratory for Laser Energetics in Rochester, New York is used to create such extreme conditions. Targets are designed with ripples or waves on the surface, and these samples are subjected to high pressures using Omega’s high energy laser beams. In these experiments, the observational parameter is the Richtmyer-Meshkov (RM) instability in the form of ripple growth on single-mode ripples. The experimental platform used for these experiments is the “ride-along” laser compression recovery experiments, which provide a way to recover the specimens having been subjected to high pressures. Six different experiments are performed on the Omega laser using single-mode tantalum targets at different laser energies. The energy indicates the amount of laser energy that impinges the target. For each target, values for growth factor are obtained by comparing the profile of ripples before and after the experiment. With increasing energy, the growth factor increased.

Engineering simulations are used to interpret and correlate the measurements of growth factor to a measure of strength. In order to validate the engineering constitutive model for tantalum, a series of simulations are performed using the code Eureka, based on the Optimal Transportation Meshfree (OTM) method. Two different configurations are studied in the simulations: RM instabilities in single and multimode ripples. Six different simulations are performed for the single ripple configuration of the RM instability experiment, with drives corresponding to laser energies used in the experiments. Each successive simulation is performed at higher drive energy, and it is observed that with increasing energy, the growth factor increases. Overall, there is favorable agreement between the data from the simulations and the experiments. The peak growth factors from the simulations and the experiments are within 10% agreement. For the multimode simulations, the goal is to assist in the design of the laser driven experiments using the Omega laser. A series of three-mode and four-mode patterns are simulated at various energies and the resulting growth of the RM instability is computed. Based on the results of the simulations, a configuration is selected for the multimode experiments. These simulations also serve as validation for the constitutive model and the material parameters for tantalum that are used in the simulations.

By designing samples with initial perturbations in the form of single-mode and multimode ripples and subjecting these samples to high pressures, the Richtmyer-Meshkov instability is investigated in both laser compression experiments and simulations. By correlating the growth of these ripples to measures of strength, a better understanding of the strength of tantalum at high pressures is achieved.

A quantitative investigation of the solar modulation of cosmic-ray protons and helium nuclei

The differential energy spectra of cosmic-ray protons and He nuclei have been measured at energies up to 315 MeV/nucleon using balloon- and satellite-borne instruments. These spectra are presented for solar quiet times for the years 1966 through 1970. The data analysis is verified by extensive accelerator calibrations of the detector systems and by calculations and measurements of the production of secondary protons in the atmosphere.

The spectra of protons and He nuclei in this energy range are dominated by the solar modulation of the local interstellar spectra. The transport equation governing this process includes as parameters the solar-wind velocity, V, and a diffusion coefficient, K(r,R), which is assumed to be a scalar function of heliocentric radius, r, and magnetic rigidity, R. The interstellar spectra, j_D, enter as boundary conditions on the solutions to the transport equation. Solutions to the transport equation have been calculated for a broad range of assumed values for K(r,R) and j_D and have been compared with the measured spectra.

It is found that the solutions may be characterized in terms of a dimensionless parameter, ψ(r,R) = ^∞∫_r V dr'/(K(r',R). The amount of modulation is roughly proportional to ψ. At high energies or far from the Sun, where the modulation is weak, the solution is determined primarily by the value of ψ (and the interstellar spectrum) and is not sensitive to the radial dependence of the diffusion coefficient. At low energies and for small r, where the effects of adiabatic deceleration are found to be large, the spectra are largely determined by the radial dependence of the diffusion coefficient and are not very sensitive to the magnitude of ψ or to the interstellar spectra. This lack of sensitivity to j_D implies that the shape of the spectra at Earth cannot be used to determine the interstellar intensities at low energies.

Values of ψ determined from electron data were used to calculate the spectra of protons and He nuclei near Earth. Interstellar spectra of the form j_D α (W - 0.25m)^-2.65 for both protons and He nuclei were found to yield the best fits to the measured spectra for these values of ψ, where W is the total energy and m is the rest energy. A simple model for the diffusion coefficient was used in which the radial and rigidity dependence are separable and K is independent of radius inside a modulation region which has a boundary at a distance D. Good agreement was found between the measured and calculated spectra for the years 1965 through 1968, using typical boundary distances of 2.7 and 6.1 A.U. The proton spectra observed in 1969 and 1970 were flatter than in previous years. This flattening could be explained in part by an increase in D, but also seemed to require that a noticeable fraction of the observed protons at energies as high at 50 to 100 MeV be attributed to quiet-time solar emission. The turnup in the spectra at low energies observed in all years was also attributed to solar emission. The diffusion coefficient used to fit the 1965 spectra is in reasonable agreement with that determined from the power spectra of the interplanetary magnetic field (Jokipii and Coleman, 1968). We find a factor of roughly 3 increase in ψ from 1965 to 1970, corresponding to the roughly order of magnitude decrease in the proton intensity at 250 MeV. The change in ψ might be attributed to a decrease in the diffusion coefficient, or, if the diffusion coefficient is essentially unchanged over that period (Mathews et al., 1971), might be attributed to an increase in the boundary distance, D.

First principles based multiparadigm modeling of electronic structures and dynamics

Electronic structures and dynamics are the key to linking the material composition and structure to functionality and performance.

An essential issue in developing semiconductor devices for photovoltaics is to design materials with optimal band gaps and relative positioning of band levels. Approximate DFT methods have been justified to predict band gaps from KS/GKS eigenvalues, but the accuracy is decisively dependent on the choice of XC functionals. We show here for CuInSe₂ and CuGaSe₂, the parent compounds of the promising CIGS solar cells, conventional LDA and GGA obtain gaps of 0.0-0.01 and 0.02-0.24 eV (versus experimental values of 1.04 and 1.67 eV), while the historically first global hybrid functional, B3PW91, is surprisingly the best, with band gaps of 1.07 and 1.58 eV. Furthermore, we show that for 27 related binary and ternary semiconductors, B3PW91 predicts gaps with a MAD of only 0.09 eV, which is substantially better than all modern hybrid functionals, including B3LYP (MAD of 0.19 eV) and screened hybrid functional HSE06 (MAD of 0.18 eV).

The laboratory performance of CIGS solar cells (> 20% efficiency) makes them promising candidate photovoltaic devices. However, there remains little understanding of how defects at the CIGS/CdS interface affect the band offsets and interfacial energies, and hence the performance of manufactured devices. To determine these relationships, we use the B3PW91 hybrid functional of DFT with the AEP method that we validate to provide very accurate descriptions of both band gaps and band offsets. This confirms the weak dependence of band offsets on surface orientation observed experimentally. We predict that the CBO of perfect CuInSe₂/CdS interface is large, 0.79 eV, which would dramatically degrade performance. Moreover we show that band gap widening induced by Ga adjusts only the VBO, and we find that Cd impurities do not significantly affect the CBO. Thus we show that Cu vacancies at the interface play the key role in enabling the tunability of CBO. We predict that Na further improves the CBO through electrostatically elevating the valence levels to decrease the CBO, explaining the observed essential role of Na for high performance. Moreover we find that K leads to a dramatic decrease in the CBO to 0.05 eV, much better than Na. We suggest that the efficiency of CIGS devices might be improved substantially by tuning the ratio of Na to K, with the improved phase stability of Na balancing phase instability from K. All these defects reduce interfacial stability slightly, but not significantly.

A number of exotic structures have been formed through high pressure chemistry, but applications have been hindered by difficulties in recovering the high pressure phase to ambient conditions (i.e., one atmosphere and room temperature). Here we use dispersion-corrected DFT (PBE-ulg flavor) to predict that above 60 GPa the most stable form of N₂O (the laughing gas in its molecular form) is a 1D polymer with an all-nitrogen backbone analogous to cis-polyacetylene in which alternate N are bonded (ionic covalent) to O. The analogous trans-polymer is only 0.03-0.10 eV/molecular unit less stable. Upon relaxation to ambient conditions both polymers relax below 14 GPa to the same stable non-planar trans-polymer, accompanied by possible electronic structure transitions. The predicted phonon spectrum and dissociation kinetics validate the stability of this trans-poly-NNO at ambient conditions, which has potential applications as a new type of conducting polymer with all-nitrogen chains and as a high-energy oxidizer for rocket propulsion. This work illustrates in silico materials discovery particularly in the realm of extreme conditions.

Modeling non-adiabatic electron dynamics has been a long-standing challenge for computational chemistry and materials science, and the eFF method presents a cost-efficient alternative. However, due to the deficiency of FSG representation, eFF is limited to low-Z elements with electrons of predominant s-character. To overcome this, we introduce a formal set of ECP extensions that enable accurate description of p-block elements. The extensions consist of a model representing the core electrons with the nucleus as a single pseudo particle represented by FSG, interacting with valence electrons through ECPs. We demonstrate and validate the ECP extensions for complex bonding structures, geometries, and energetics of systems with p-block character (C, O, Al, Si) and apply them to study materials under extreme mechanical loading conditions.

Despite its success, the eFF framework has some limitations, originated from both the design of Pauli potentials and the FSG representation. To overcome these, we develop a new framework of two-level hierarchy that is a more rigorous and accurate successor to the eFF method. The fundamental level, GHA-QM, is based on a new set of Pauli potentials that renders exact QM level of accuracy for any FSG represented electron systems. To achieve this, we start with using exactly derived energy expressions for the same spin electron pair, and fitting a simple functional form, inspired by DFT, against open singlet electron pair curves (H₂ systems). Symmetric and asymmetric scaling factors are then introduced at this level to recover the QM total energies of multiple electron pair systems from the sum of local interactions. To complement the imperfect FSG representation, the AMPERE extension is implemented, and aims at embedding the interactions associated with both the cusp condition and explicit nodal structures. The whole GHA-QM+AMPERE framework is tested on H element, and the preliminary results are promising.

Biological Activity of Rhodium Metalloinsertors and the Design of Bifunctional Conjugates

The Barton laboratory has established that octahedral rhodium complexes bearing the sterically expansive 5,6-chrysene diimine ligand can target thermodynamically destabilized sites, such as base pair mismatches, in DNA with high affinity and selectivity. These complexes approach DNA from the minor groove, ejecting the mismatched base pairs from the duplex in a binding mode termed metalloinsertion. In recent years, we have shown that these metalloinsertor complexes also exhibit cytotoxicity preferentially in cancer cells that are deficient in the mismatch repair (MMR) machinery.

Here, we establish that a sensitive structure-activity relationship exists for rhodium metalloinsertors. We studied the relationship between the chemical structures of metalloinsertors and their effect on biological activity for ten complexes with similar DNA binding affinities, but wide variation in their lipophilicity. Drastic differences were observed in the selectivities of the complexes for MMR-deficient cells. Compounds with hydrophilic ligands were highly selective, exhibiting preferential cytotoxicity in MMR-deficient cells at low concentrations and short incubation periods, whereas complexes with lipophilic ligands displayed poor cell-selectivity. It was discovered that all of the complexes localized to the nucleus in concentrations sufficient for mismatch binding; however, highly lipophilic complexes also exhibited high mitochondrial uptake. Significantly, these results support the notion that mitochondrial DNA is not the desired target for our metalloinsertor complexes; instead, selectivity stems from targeting mismatches in genomic DNA.

We have also explored the potential for metalloinsertors to be developed into more complex structures with multiple functionalities that could either enhance their overall potency or impart mismatch selectivity onto other therapeutic cargo. We have constructed a family of bifunctional metalloinsertor conjugates incorporating cis-platinum, each unique in its chemical structure, DNA binding interactions, and biological activity. The study of these complexes in MMR-deficient cells has established that the cell-selective biological activity of rhodium metalloinsertors proceeds through a critical cellular pathway leading to necrosis.

We further explored the underlying mechanisms surrounding the biological response to mismatch recognition by metalloinsertors in the genome. Immunofluorescence assays of MMR-deficient and MMR-proficient cells revealed that a critical biomarker for DNA damage, phosphorylation of histone H2AX (γH2AX) rapidly accumulates in response to metalloinsertor treatment, signifying the induction of double strand breaks in the genome. Significantly, we have discovered that our metalloinsertor complexes selectively inhibit transcription in MMR-deficient cells, which may be a crucial checkpoint in the eventual breakdown of the cell via necrosis. Additionally, preliminary in vivo studies have revealed the capability of these compounds to traverse the complex environments of multicellular organisms and accumulate in MMR-deficient tumors. Our ever-increasing understanding of metalloinsertors, as well as the development of new generations of complexes both monofunctional and bifunctional, enables their continued progress into the clinic as promising new chemotherapeutic agents.

A study of T=1, T=3/2, and T=2 states in some light nuclei using (He^3,n) reactions

The (He³, n) reactions on B¹¹, N¹⁵, O¹⁶, and O¹⁸ targets have been studied using a pulsed-beam time-of-flight spectrometer. Special emphasis was placed upon the determination of the excitation energies and properties of states with T = 1 (in Ne¹⁸), T = 3/2 (in N¹³ and F¹⁷) and T = 2 (in Ne²⁰). The identification of the T = 3/2 and T = 2 levels is based on the structure of these states as revealed by intensities and shapes of angular distributions. The reactions are interpreted in terms of double stripping theory. Angular distributions have been compared with plane and distorted wave stripping theories. Results for the four reactions are summarized below:

1) O¹⁶ (He³, n). The reaction has been studied at incident energies up to 13.5 MeV and two previously unreported levels in Ne¹⁸ were observed at E_x = 4.55 ± .015 MeV (Γ = 70 ± 30 keV) and E_x = 5.14 ± .018 MeV (Γ = 100 ± 40 keV).

2) B¹¹ (He³, n). The reaction has been studied at incident energies up to 13.5 MeV. Three T = 3/2 levels in N¹³ have been identified at E_x = 15.068 ± .008 MeV (Γ ˂ 15 keV), E_x = 18.44 ± .04, and E_x 18.98 ± .02 MeV (Γ = 40 ± 20 keV).

3) N¹⁵ (He³, n). The reaction has been studied at incident energies up to 11.88 MeV. T = 3/2 levels in F¹⁷ have been identified at E_x = 11.195 ± .007 MeV (Γ ˂ 20 keV), E_x = 12.540 ± .010 MeV (Γ ˂ 25 keV), and E_x = 13.095 ± .009 MeV (Γ ˂ 25 keV).

4) O¹⁸ (He³, n). The reaction has been studied at incident energies up to 9.0 MeV. The excitation energy of the lowest T = 2 level in Ne²⁰ has been found to be 16.730 ± .006 MeV (Γ ˂ 20 keV).

Angular distributions of the transitions leading to the above higher isospin states are well described by double stripping theory. Analog correspondences are established by comparing the present results with recent studies (t, p) and (He³, p) reactions on the same targets.

Beta and gamma vibrational bands in deformed nuclei

Energies and relative intensities of gamma transitions in ¹⁵²Sm, ¹⁵²Gd, ¹⁵⁴Gd, ¹⁶⁶Er, and ²³²U following radioactive decay have been measured with a Ge(Li) spectrometer. A peak fitting program has been developed to determine gamma ray energies and relative intensities with precision sufficient to give a meaningful test of nuclear models. Several previously unobserved gamma rays were placed in the nuclear level schemes. Particular attention has been paid to transitions from the beta and gamma vibrational bands, since the gamma ray branching ratios are sensitive tests of configuration mixing in the nuclear levels. As the reduced branching ratios depend on the multipolarity of the gamma transitions, experiments were performed to measure multipole mixing ratios for transitions from the gamma vibrational band. In ¹⁵⁴Gd, angular correlation experiments showed that transitions from the gamma band to the ground state band were predominantly electric quadrupole, in agreement with the rotational model. In ²³²U, the internal conversion spectrum has been studied with a Si(Li) spectrometer constructed for electron spectroscopy. The strength of electric monopole transitions and the multipolarity of some gamma transitions have been determined from the measured relative electron intensities.

The results of the experiments have been compared with the rotational model and several microscopic models. Relative B(E2) strengths for transitions from the gamma band in ²³²U and ¹⁶⁶Er are in good agreement with a single parameter band mixing model, with values of z₂= 0.025(10) and 0.046(2), respectively. Neither the beta nor the gamma band transition strengths in ¹⁵²Sm and ¹⁵⁴Gd can be accounted for by a single parameter theory, nor can agreement be found by considering the large mixing found between the beta and gamma bands. The relative B(E2) strength for transitions from the gamma band to the beta band in ²³²U is found to be five times greater than the strength to the ground state band, indicating collective transitions with strength approximately 15 single particle units.

1. Ultrasonic studies of binary liquid structure in the critical region. Theory and experiment for the 2,6-lutidine/water system. 2. Hartree-Fock calculations of electric polarizabilities of some simple atoms and molecules, and their practicality. 3. Calculation of vibrational transition probabilities in collinear atom-diatom and diatom-diatom collisions with Lennard-Jones interaction

Part 1. Many interesting visual and mechanical phenomena occur in the critical region of fluids, both for the gas-liquid and liquid-liquid transitions. The precise thermodynamic and transport behavior here has some broad consequences for the molecular theory of liquids. Previous studies in this laboratory on a liquid-liquid critical mixture via ultrasonics supported a basically classical analysis of fluid behavior by M. Fixman (e. g., the free energy is assumed analytic in intensive variables in the thermodynamics)--at least when the fluid is not too close to critical. A breakdown in classical concepts is evidenced close to critical, in some well-defined ways. We have studied herein a liquid-liquid critical system of complementary nature (possessing a lower critical mixing or consolute temperature) to all previous mixtures, to look for new qualitative critical behavior. We did not find such new behavior in the ultrasonic absorption ascribable to the critical fluctuations, but we did find extra absorption due to chemical processes (yet these are related to the mixing behavior generating the lower consolute point). We rederived, corrected, and extended Fixman's analysis to interpret our experimental results in these more complex circumstances. The entire account of theory and experiment is prefaced by an extensive introduction recounting the general status of liquid state theory. The introduction provides a context for our present work, and also points out problems deserving attention. Interest in these problems was stimulated by this work but also by work in Part 3.

Part 2. Among variational theories of electronic structure, the Hartree-Fock theory has proved particularly valuable for a practical understanding of such properties as chemical binding, electric multipole moments, and X-ray scattering intensity. It also provides the most tractable method of calculating first-order properties under external or internal one-electron perturbations, either developed explicitly in orders of perturbation theory or in the fully self-consistent method. The accuracy and consistency of first-order properties are poorer than those of zero-order properties, but this is most often due to the use of explicit approximations in solving the perturbed equations, or to inadequacy of the variational basis in size or composition. We have calculated the electric polarizabilities of H₂, He, Li, Be, LiH, and N₂ by Hartree-Fock theory, using exact perturbation theory or the fully self-consistent method, as dictated by convenience. By careful studies on total basis set composition, we obtained good approximations to limiting Hartree-Fock values of polarizabilities with bases of reasonable size. The values for all species, and for each direction in the molecular cases, are within 8% of experiment, or of best theoretical values in the absence of the former. Our results support the use of unadorned Hartree-Pock theory for static polarizabilities needed in interpreting electron-molecule scattering data, collision-induced light scattering experiments, and other phenomena involving experimentally inaccessible polarizabilities.

Part 3. Numerical integration of the close-coupled scattering equations has been carried out to obtain vibrational transition probabilities for some models of the electronically adiabatic H₂-H₂ collision. All the models use a Lennard-Jones interaction potential between nearest atoms in the collision partners. We have analyzed the results for some insight into the vibrational excitation process in its dependence on the energy of collision, the nature of the vibrational binding potential, and other factors. We conclude also that replacement of earlier, simpler models of the interaction potential by the Lennard-Jones form adds very little realism for all the complication it introduces. A brief introduction precedes the presentation of our work and places it in the context of attempts to understand the collisional activation process in chemical reactions as well as some other chemical dynamics.

One-dimensional potential for image-potential states on graphene

In the framework of dielectric theory, the static non-local self-energy of an electron near an ultra-thin polarizable layer has been calculated and applied to study binding energies of image-potential states near free-standing graphene. The corresponding series of eigenvalues and eigenfunctions have been obtained by numerically solving the one-dimensional Schrodinger equation. The imagepotential state wave functions accumulate most of their probability outside the slab. We find that the random phase approximation (RPA) for the nonlocal dielectric function yields a superior description for the potential inside the slab, but a simple Fermi-Thomas theory can be used to get a reasonable quasi-analytical approximation to the full RPA result that can be computed very economically. Binding energies of the image-potential states follow a pattern close to the Rydberg series for a perfect metal with the addition of intermediate states due to the added symmetry of the potential. The formalism only requires a minimal set of free parameters: the slab width and the electronic density. The theoretical calculations are compared with experimental results for the work function and image-potential states obtained by two-photon photoemission.

Correlations between the performance characteristics of a liquid crystal laser and the macroscopic material properties.

In order to understand how the performance of a liquid-crystal laser depends on the physical properties of the low molar mass nematic host, we have studied the energy threshold and slope efficiency of ten optically pumped liquid-crystal lasers based on different hosts. Specifically, this leads to a variation in the birefringence, the orientational order parameter, and the order parameter of the transition dipole moment of the dye. It is found that low threshold energies and high slope efficiencies correlate with high order parameters and large birefringences. To a first approximation this can be understood by considering analytical expressions for the threshold and slope efficiency, which are derived from the space-independent rate equations for a two-level system, in terms of the macroscopic liquid crystal properties.

Probabilistic fault identification using vibration data and neural networks

Bayesian formulated neural networks are implemented using hybrid Monte Carlo method for probabilistic fault identification in cylindrical shells. Each of the 20 nominally identical cylindrical shells is divided into three substructures. Holes of (12±2) mm in diameter are introduced in each of the substructures and vibration data are measured. Modal properties and the Coordinate Modal Assurance Criterion (COMAC) are utilized to train the two modal-property-neural-networks. These COMAC are calculated by taking the natural-frequency-vector to be an additional mode. Modal energies are calculated by determining the integrals of the real and imaginary components of the frequency response functions over bandwidths of 12% of the natural frequencies. The modal energies and the Coordinate Modal Energy Assurance Criterion (COMEAC) are used to train the two frequency-response-function-neural-networks. The averages of the two sets of trained-networks (COMAC and COMEAC as well as modal properties and modal energies) form two committees of networks. The COMEAC and the COMAC are found to be better identification data than using modal properties and modal energies directly. The committee approach is observed to give lower standard deviations than the individual methods. The main advantage of the Bayesian formulation is that it gives identities of damage and their respective confidence intervals.

The identification and quantification of highly stable 'common hairpin' in the dynamic process of co-transcriptional mRNA folding

In our studies, 88 human mRNA samples were collected from the Integrated Sequence-Structure database and then the dynamic process in co-transcriptional mRNA folding was simulated using the RNAstructure version 4.1 program. Through statistical analyses of the frequencies of occurrence of hairpins, a group of special folding structures-the 'common hairpins'-were identified. These 'common hairpins' have lower energies and occur in all the subsequent folding units that formed in the dynamic folding process. By applying the formulas (1)-(4) of the 'common hairpins' statistical model, 163 'common hairpins' were found, to make up about 7% of the total of 2286 hairpins. Classified studies further show that the 'common hairpins' that were studied may oscillate in the dynamic folding process. However, the hairpin loops of the 'common hairpins' and stems proximal to those 'common hairpins' loops maintain topologically stable structures, while other loops and stems distal to the 'common hairpins' loops are shown to be alterable structures. Strikingly, further studies indicate that the stable structures of these 'common hairpins' may have unbeknown effects on controlling the formation of protein structures in the translation process (unpublished results). (c) 2005 Elsevier B.V. All rights reserved.

Observing graphene grow: catalyst-graphene interactions during scalable graphene growth on polycrystalline copper.

Complementary in situ X-ray photoelectron spectroscopy (XPS), X-ray diffractometry, and environmental scanning electron microscopy are used to fingerprint the entire graphene chemical vapor deposition process on technologically important polycrystalline Cu catalysts to address the current lack of understanding of the underlying fundamental growth mechanisms and catalyst interactions. Graphene forms directly on metallic Cu during the high-temperature hydrocarbon exposure, whereby an upshift in the binding energies of the corresponding C1s XPS core level signatures is indicative of coupling between the Cu catalyst and the growing graphene. Minor carbon uptake into Cu can under certain conditions manifest itself as carbon precipitation upon cooling. Postgrowth, ambient air exposure even at room temperature decouples the graphene from Cu by (reversible) oxygen intercalation. The importance of these dynamic interactions is discussed for graphene growth, processing, and device integration.