Linear fuzzy gene network models obtained from microarray data by exhaustive search

BACKGROUND: Recent technological advances in high-throughput data collection allow for experimental study of increasingly complex systems on the scale of the whole cellular genome and proteome. Gene network models are needed to interpret the resulting large and complex data sets. Rationally designed perturbations (e.g., gene knock-outs) can be used to iteratively refine hypothetical models, suggesting an approach for high-throughput biological system analysis. We introduce an approach to gene network modeling based on a scalable linear variant of fuzzy logic: a framework with greater resolution than Boolean logic models, but which, while still semi-quantitative, does not require the precise parameter measurement needed for chemical kinetics-based modeling. RESULTS: We demonstrated our approach with exhaustive search for fuzzy gene interaction models that best fit transcription measurements by microarray of twelve selected genes regulating the yeast cell cycle. Applying an efficient, universally applicable data normalization and fuzzification scheme, the search converged to a small number of models that individually predict experimental data within an error tolerance. Because only gene transcription levels are used to develop the models, they include both direct and indirect regulation of genes. CONCLUSION: Biological relationships in the best-fitting fuzzy gene network models successfully recover direct and indirect interactions predicted from previous knowledge to result in transcriptional correlation. Fuzzy models fit on one yeast cell cycle data set robustly predict another experimental data set for the same system. Linear fuzzy gene networks and exhaustive rule search are the first steps towards a framework for an integrated modeling and experiment approach to high-throughput "reverse engineering" of complex biological systems

Vertex-corrected tunneling inversion in superconductors: Pb

The McMillan-Rowell tunneling inversion program, which extracts the electron-phonon spectral function $\alpha^2F(\Omega)$ and the Coulomb pseudopotential $\mu^*$ from experimental tunneling data, is generalized to include the lowest-order vertex correction. We neglect the momentum dependence of the electron-phonon matrix elements, which is equivalent to using a local approximation. The perturbation theory is performed on the imaginary axis and then an exact analytic continuation is employed to produce the density of states on the real axis. Comparison is made with the experimental data for Pb.Comment: 14 pages, typeset in ReVTeX, including three encapsulated postscript figure

Electronic interactions in fullerene spheres

The electron-phonon and Coulomb interactions inC$_{60}$, and larger fullerene spheres are analyzed. The coupling between electrons and intramolecular vibrations give corrections $\sim 1 - 10$ meV to the electronic energies for C$_{60}$, and scales as $R^{-4}$ in larger molecules. The energies associated with electrostatic interactions are of order $\sim 1 - 4$ eV, in C$_{60}$ and scale as $R^{-1}$. Charged fullerenes show enhanced electron-phonon coupling, $\sim 10$ meV, which scales as $R^{-2}$. Finally, it is argued that non only C$_{60}^{-}$, but also C$_{60}^{--}$ are highly polarizable molecules. The polarizabilities scale as $R^3$ and $R^4$, respectively. The role of this large polarizability in mediating intermolecular interactions is also discussed.Comment: 12 pages. No figure

Plasmon Lifetime in K: A Case Study of Correlated Electrons in Solids Amenable to Ab Initio Theory

On the basis of a new ab initio, all-electron response scheme, formulated within time-dependent density-functional theory, we solve the puzzle posed by the anomalous dispersion of the plasmon linewidth in K. The key damping mechanism is shown to be decay into particle-hole pairs involving empty states of d-symmetry. While the effect of many-particle correlations is small, the correlations built into the "final-state" -d-bands play an important, and novel, role ---which is related to the phase-space complexity associated with these flat bands. Our case study of plasmon lifetime in K illustrates the importance of ab initio paradigms for the study of excitations in correlated-electron systems.Comment: 12 pages, 4 figures, for html browsing see http://web.utk.edu/~weik

Beyond Eliashberg superconductivity in MgB2: anharmonicity, two-phonon scattering, and multiple gaps

Density-functional calculations of the phonon spectrum and electron-phonon coupling in MgB$_2$ are presented. The $E_{2g}$ phonons, which involve in-plane B displacements, couple strongly to the $p_{x,y}$ electronic bands. The isotropic electron-phonon coupling constant is calculated to be about 0.8. Allowing for different order parameters in different bands, the superconducting $\lambda$ in the clean limit is calculated to be significantly larger. The $E_{2g}$ phonons are strongly anharmonic, and the non-linear contribution to the coupling between the $E_{2g}$ modes and the p$_{x,y}$ bands is significant.Comment: 4 pages, 3 figure

Vibrational signatures for low-energy intermediate-sized Si clusters

We report low-energy locally stable structures for the clusters Si20 and Si21. The structures were obtained by performing geometry optimizations within the local density approximation. Our calculated binding energies for these clusters are larger than any previously reported for this size regime. To aid in the experimental identification of the structures, we have computed the full vibrational spectra of the clusters, along with the Raman and IR activities of the various modes using a recently developed first-principles technique. These represent, to our knowledge, the first calculations of Raman and IR spectra for Si clusters of this size

Reevaluating electron-phonon coupling strengths: Indium as a test case for ab initio and many-body-theory methods

Using indium as a test case, we investigate the accuracy of the electron-phonon coupling calculated with state-of-the-art ab initio and many-body theory methods. The ab initio calculations -- where electrons are treated in the local-density approximation, and phonons and the electron-phonon interaction are treated within linear response -- predict an electron-phonon spectral function alpha^2 F(omega) which translates into a relative tunneling conductance that agrees with experiment to within one part in 1000. The many-body theory calculations -- where alpha^2 F(omega) is extracted from tunneling data by means of the McMillan-Rowell tunneling inversion method -- provide spectral functions that depend strongly on details of the inversion process. For the the most important moment of alpha^2 F(omega), the mass-renormalization parameter lambda, we report 0.9 +/- 0.1, in contrast to the value 0.805 quoted for nearly three decades in the literature. The ab initio calculations also provide the transport electron-phonon spectral function alpha_{tr}^2 F(omega), from which we calculate the resistivity as a function of temperature in good agreement with experiment.Comment: 16 pages, 5 figure

Comparison of structural transformations and superconductivity in compressed Sulfur and Selenium

Density-functional calculations are presented for high-pressure structural phases of S and Se. The structural phase diagrams, phonon spectra, electron-phonon coupling, and superconducting properties of the isovalent elements are compared. We find that with increasing pressure, Se adopts a sequence of ever more closely packed structures (beta-Po, bcc, fcc), while S favors more open structures (beta-Po, simple cubic, bcc). These differences are shown to be attributable to differences in the S and Se core states. All the compressed phases of S and Se considered are calculated to have weak to moderate electron-phonon coupling strengths consistent with superconducting transition temperatures in the range of 1 to 20 K. Our results compare well with experimental data on the beta-Po --> bcc transition pressure in Se and on the superconducting transition temperature in beta-Po S. Further experiments are suggested to search for the other structural phases predicted at higher pressures and to test theoretical results on the electron-phonon interaction and superconducting properties

Exchange-correlation kernels for excited states in solids

The performance of several common approximations for the exchange-correlation kernel within time-dependent density-functional theory is tested for elementary excitations in the homogeneous electron gas. Although the adiabatic local-density approximation gives a reasonably good account of the plasmon dispersion, systematic errors are pointed out and traced to the neglect of the wavevector dependence. Kernels optimized for atoms are found to perform poorly in extended systems due to an incorrect behavior in the long-wavelength limit, leading to quantitative deviations that significantly exceed the experimental error bars for the plasmon dispersion in the alkali metals.Comment: 7 pages including 5 figures, RevTe

Plasmonic excitations in noble metals: The case of Ag

The delicate interplay between plasmonic excitations and interband transitions in noble metals is described by means of {\it ab initio} calculations and a simple model in which the conduction electron plasmon is coupled to the continuum of electron-hole pairs. Band structure effects, specially the energy at which the excitation of the $d$-like bands takes place, determine the existence of a subthreshold plasmonic mode, which manifests itself in Ag as a sharp resonance at 3.8 eV. However, such a resonance is not observed in the other noble metals. Here, this different behavior is also analyzed and an explanation is provided.Comment: 9 pages, 8 figure