Dynamical properties of a strongly correlated model for quarter-filled layered organic molecular crystals

The dynamical properties of an extended Hubbard model, which is relevant to quarter-filled layered organic molecular crystals, are analyzed. We have computed the dynamical charge correlation function, spectral density, and optical conductivity using Lanczos diagonalization and large-N techniques. As the ratio of the nearest-neighbor Coulomb repulsion, V, to the hopping integral, t, increases there is a transition from a metallic phase to a charge-ordered phase. Dynamical properties close to the ordering transition are found to differ from the ones expected in a conventional metal. Large-N calculations display an enhancement of spectral weight at low frequencies as the system is driven closer to the charge-ordering transition in agreement with Lanczos calculations. As V is increased the charge correlation function displays a collective mode which, for wave vectors close to (pi,pi), increases in amplitude and softens as the charge-ordering transition is approached. We propose that inelastic x-ray scattering be used to detect this mode. Large-N calculations predict superconductivity with d(xy) symmetry close to the ordering transition. We find that this is consistent with Lanczos diagonalization calculations, on lattices of 20 sites, which find that the binding energy of two holes becomes negative close to the charge-ordering transition.

Theory of strongly correlated electron systems. I. Intersite coulomb interaction and the approximation of renormalized fermions in total energy calculations

The diagrammatic strong-coupling perturbation theory (SCPT) for correlated electron systems is developed for intersite Coulomb interaction and for a nonorthogonal basis set. The construction is based on iterations of exact closed equations for many - electron Green functions (GFs) for Hubbard operators in terms of functional derivatives with respect to external sources. The graphs, which do not contain the contributions from the fluctuations of the local population numbers of the ion states, play a special role: a one-to-one correspondence is found between the subset of such graphs for the many - electron GFs and the complete set of Feynman graphs of weak-coupling perturbation theory (WCPT) for single-electron GFs. This fact is used for formulation of the approximation of renormalized Fermions (ARF) in which the many-electron quasi-particles behave analogously to normal Fermions. Then, by analyzing: (a) Sham's equation, which connects the self-energy and the exchange- correlation potential in density functional theory (DFT); and (b) the Galitskii and Migdal expressions for the total energy, written within WCPT and within ARF SCPT, a way we suggest a method to improve the description of the systems with correlated electrons within the local density approximation (LDA) to DFT. The formulation, in terms of renormalized Fermions LIDA (RF LDA), is obtained by introducing the spectral weights of the many electron GFs into the definitions of the charge density, the overlap matrices, effective mixing and hopping matrix elements, into existing electronic structure codes, whereas the weights themselves have to be found from an additional set of equations. Compared with LDA+U and self-interaction correction (SIC) methods, RF LDA has the advantage of taking into account the transfer of spectral weights, and, when formulated in terms of GFs, also allows for consideration of excitations and nonzero temperature. Going beyond the ARF SCPT, as well as RF LIDA, and taking into account the fluctuations of ion population numbers would require writing completely new codes for ab initio calculations. The application of RF LDA for ab initio band structure calculations for rare earth metals is presented in part 11 of this study (this issue). (c) 2005 Wiley Periodicals, Inc.

Theory of strongly correlated electron systems. II. Including correlation effects into electronic structure calculations

We have previously shown that a division of the f-shell into two subsystems gives a better understanding of the cohesive properties as well the general behavior of lanthanide systems. In this article, we present numerical computations, using the suggested method. We show that the picture is consistent with most experimental data, e.g., the equilibrium volume and electronic structure in general. Compared with standard energy band calculations and calculations based on the self-interaction correction and LIDA + U, the f-(non-f)-mixing interaction is decreased by spectral weights of the many-body states of the f-ion. (c) 2005 Wiley Periodicals, Inc.

New insights into NO-carbon and N2O-carbon reactions from quantum mechanical calculations

Density functional theory calculations were used to investigate the mechanisms of NO-carbon and N2O-carbon reactions. It was the first time that the importance of surface nitrogen groups was addressed in the kinetic behaviors of the NO-carbon reaction. It was found that the off-plane nitrogen groups that are adjacent to the zigzag edge sites and in-plane nitrogen groups that are located on the armchair sites make the bond energy of oxygen desorption even ca. 20% lower than that of the off-plane epoxy group adjacent to zigzag edge sites and in-plane o-quinone oxygen atoms on armchair sites; this may explain the reason why the experimentally obtained activation energy of the NO-carbon reaction is ca. 20% lower than that of the O-2-carbon reaction over 923 K. A higher ratio of oxygen atoms can be formed in the N2O-carbon reaction, because of the lower dissociation energy of N2O, which results in a higher ratio of off-plane epoxy oxygen atoms. The desorption energy of semiquinone with double adjacent off-plane oxygen groups is ca. 20% less than that of semiquinone with only one adjacent off-plane oxygen group. This may be the reason why the activation energy of N2O is also ca. 20% less than that of the O-2-carbon reaction. The new mechanism can also provide a good qualitative comparison for the relative reaction rates of NO-, N2O-, and O-2-carbon reactions. The anisotropic characters of these gas-carbon reactions can also be well explained.

Characterization of carbon molecular sieves using methane and carbon dioxide as adsorptive probes

Nitrogen adsorption at 77 K is the current standard means for pore size determination of adsorbent materials. However, nitrogen adsorption reaches limitations when dealing with materials such as molecular sieving carbon with a high degree of ultramicroporosity. In this investigation, methane and carbon dioxide adsorption is explored as a possible alternative to the standard nitrogen probe. Methane and carbon dioxide adsorption equilibria and kinetics are measured in a commercially derived carbon molecular sieve over a range of temperatures. The pore size distribution is determined from the adsorption equilibrium, and the kinetics of adsorption is shown to be Fickian for carbon dioxide and non-Fickian for methane. The non-Fickian response is attributed to transport resistance at the pore mouth experienced by the methane molecules but not by the carbon dioxide molecules. Additionally, the change in the rate of adsorption with loading is characterized by the Darken relation in the case of carbon dioxide diffusion but is greater than that predicted by the Darken relation for methane transport. Furthermore, the proposition of inkbottle-shaped micropores in molecular sieving carbon is supported by the determination of the activation energy for the transport of methane and subsequent sizing of the pore-mouth barrier by molecular potential calculations.

Molecular size as the main determinant of solute maximum flux across the skin

One of the most important determinants of dermatological and systemic penetration after topical application is the delivery or flux of solutes into or through the skin. The maximum dose of solute able to be delivered over a given period of time and area of application is defined by its maximum flux (J(max), mol per cm(2) per h) from a given vehicle. In this work, J(max) values from aqueous solution across human skin were acquired or estimated from experimental data and correlated with solute physicochemical properties. Whereas epidermal permeability coefficients (k(p)) are optimally correlated to solute octanol-water partition coefficient (K-ow) and molecular weight (MW) was found to be the dominant determinant of J(max) for this literature data set: log J(max)=-3.90-0.0190MW (n=87, r(2)=0.847, p

Converged quantum calculations of HO2 bound states and resonances for J=6 and 10

Bound and resonance states of HO2 are calculated quantum mechanically using both the Lanczos homogeneous filter diagonalization method and the real Chebyshev filter diagonalization method for nonzero total angular momentum J=6 and 10, using a parallel computing strategy. For bound states, agreement between the two methods is quite satisfactory; for resonances, while the energies are in good agreement, the widths are in general agreement. The quantum nonzero-J specific unimolecular dissociation rates for HO2 are also calculated. (C) 2004 American Institute of Physics.

Transition-state-theory calculations for reactions of O(3P) with halogenated olefins

The potential energy surfaces for the reactions of atomic oxygen in its ground electronic state, O(P-3), with the olefins: CF2=CCl2 and CF2=CF - CF3, have been characterized using ab initio molecular orbital calculations. Geometry optimization and vibrational frequency calculations were performed for reactants, transition states and products at the MP2 and QCISD levels of theory using the 6-31G(d) basis set. This database was then used to calculate the rate constants by means of Transition-State-Theory. To obtain a better reference and to test the reliability of the activation barriers we have also carried out computations using the CCSD(T)(fc)/6-311Gdagger, MP4(SDQ)(fc)/CBSB4 and MP2(fc)/CBSB3 single point energy calculations at both of the above levels of theory, as well as with the composite CBS-RAD procedure ( P. M. Mayer, C. J. Parkinson, D. M. Smith and L. Radom, J. Chem. Phys., 1998, 108, 604) and a modi. cation of this approach, called: CBS-RAD( MP2, MP2). It was found that the kinetic parameters obtained in this work particularly with the CBS-RAD ( MP2, MP2) procedure are in reasonable agreement with the experimental values. For both reactions it is found that the channels leading to the olefin double-bond addition predominates with respect to any other reaction pathway. However, on account of the different substituents in the alkenes we have located, at all levels of theory, two transition states for each reaction. Moreover, we have found that, for the reactions studied, a correlation exists between the activation energies and the electronic structure of the transition states which can explain the influence of the substituent effect on the reactivity of the halo-olefins.

Reduced hepatic extraction of palmitate in steatosis correlated to lower level of liver fatty acid binding protein

Nonalcoholic fatty liver disease is the most common of all liver diseases. The hepatic disposition [H-3]palmitate and its low-molecular-weight metabolites in perfused normal and steatotic rat liver were studied using the multiple indicator dilution technique and a physiologically based slow diffusion/bound pharmacokinetic model. The steatotic rat model was established by administration of 17alpha-ethynylestradiol to female Wistar rats. Serum biochemistry markers and histology of treated and normal animals were assessed and indicated the presence of steatosis in the treatment group. The steatotic group showed a significantly higher alanine aminotransferase-to-aspartate aminotransferase ratio, lower levels of liver fatty acid binding protein and cytochrome P-450, as well as microvesicular steatosis with an enlargement of sinusoidal space. Hepatic extraction for unchanged [H-3]palmitate and production of low-molecular-weight metabolites were found to be significantly decreased in steatotic animals. Pharmacokinetic analysis suggested that the reduced extraction and sequestration for palmitate and its metabolites was mainly attributed to a reduction in liver fatty acid binding protein in steatosis.

Unimolecular Rovibrational Bound and Resonance States for Large Angular Momentum: J=20 Calculations for HO2

We explore the calculation of unimolecular bound states and resonances for deep-well species at large angular momentum using a Chebychev filter diagonalization scheme incorporating doubling of the autocorrelation function as presented recently by Neumaier and Mandelshtam [Phys. Rev. Lett. 86, 5031 (2001)]. The method has been employed to compute the challenging J=20 bound and resonance states for the HO2 system. The methodology has firstly been tested for J=2 in comparison with previous calculations, and then extended to J=20 using a parallel computing strategy. The quantum J-specific unimolecular dissociation rates for HO2-> H+O-2 in the energy range from 2.114 to 2.596 eV have been reported for the first time, and comparisons with the results of Troe and co-workers [J. Chem. Phys. 113, 11019 (2000) Phys. Chem. Chem. Phys. 2, 631 (2000)] from statistical adiabatic channel method/classical trajectory calculations have been made. For most of the energies, the reported statistical adiabatic channel method/classical trajectory rate constants agree well with the average of the fluctuating quantum-mechanical rates. Near the dissociation threshold, quantum rates fluctuate more severely, but their average is still in agreement with the statistical adiabatic channel method/classical trajectory results.

Sociality and the rate of molecular evolution

The molecular clock does not tick at a uniform rate in all taxa but maybe influenced by species characteristics. Eusocial species (those with reproductive division of labor) have been predicted to have faster rates of molecular evolution than their nonsocial relatives because of greatly reduced effective population size; if most individuals in a population are nonreproductive and only one or few queens produce all the offspring, then eusocial animals could have much lower effective population sizes than their solitary relatives, which should increase the rate of substitution of nearly neutral mutations. An earlier study reported faster rates in eusocial honeybees and vespid wasps but failed to correct for phylogenetic nonindependence or to distinguish between potential causes of rate variation. Because sociality has evolved independently in many different lineages, it is possible to conduct a more wide-ranging study to test the generality of the relationship. We have conducted a comparative analysis of 25 phylogenetically independent pairs of social lineages and their nonsocial relatives, including bees, wasps, ants, termites, shrimps, and mole rats, using a range of available DNA sequences (mitochondrial and nuclear DNA coding for proteins and RNAs, and nontranslated sequences). By including a wide range of social taxa, we were able to test whether there is a general influence of sociality on rates of molecular evolution and to test specific predictions of the hypothesis: (1) that social species have faster rates because they have reduced effective population sizes; (2) that mitochondrial genes would show a greater effect of sociality than nuclear genes; and (3) that rates of molecular evolution should be correlated with the degree of sociality. We find no consistent pattern in rates of molecular evolution between social and nonsocial lineages and no evidence that mitochondrial genes show faster rates in social taxa. However, we show that the most highly eusocial Hymenoptera do have faster rates than their nonsocial relatives. We also find that social parasites (that utilize the workers from related species to produce their own offspring) have faster rates than their social relatives, which is consistent with an effect of lower effective population size on rate of molecular evolution. Our results illustrate the importance of allowing for phylogenetic nonindependence when conducting investigations of determinants of variation in rate of molecular evolution.

Gaussian operator bases for correlated fermions

We formulate a general multi-mode Gaussian operator basis for fermions, to enable a positive phase-space representation of correlated Fermi states. The Gaussian basis extends existing bosonic phase-space methods to Fermi systems and thus allows first-principles dynamical or equilibrium calculations in quantum many-body Fermi systems. We prove the completeness of the basis and derive differential forms for products with one- and two-body operators. Because the basis satisfies fermionic superselection rules, the resulting phase space involves only c-numbers, without requiring anticommuting Grassmann variables. Furthermore, because of the overcompleteness of the basis, the phase-space distribution can always be chosen positive. This has important consequences for the sign problem in fermion physics.