Electronic properties and hyperfine fields of nickel-related complexes in diamond

We carried out a first-principles investigation on the microscopic properties of nickel-related defect centers in diamond. Several configurations, involving substitutional and interstitial nickel impurities, have been considered either in isolated configurations or forming complexes with other defects, such as vacancies and boron and nitrogen dopants. The results, in terms of spin, symmetry, and hyperfine fields, were compared with the available experimental data on electrically active centers in synthetic diamond. Several microscopic models, previously proposed to explain those data, have been confirmed by this investigation, while some models could be discarded. We also provided insights into the microscopic structure of several of those centers.

Quasiharmonic elastic constants corrected for deviatoric thermal stresses

The quasiharmonic approximation (QHA), in its simplest form also called the statically constrained (SC) QHA, has been shown to be a straightforward method to compute thermoelastic properties of crystals. Recently we showed that for noncubic solids SC-QHA calculations develop deviatoric thermal stresses at high temperatures. Relaxation of these stresses leads to a series of corrections to the free energy that may be taken to any desired order, up to self-consistency. Here we show how to correct the elastic constants obtained using the SC-QHA. We exemplify the procedure by correcting to first order the elastic constants of MgSiO(3) perovskite and MgSiO(3) postperovskite, the major phases of the Earth's lower mantle. We show that this first-order correction is quite satisfactory for obtaining the aggregated elastic averages of these minerals and their velocities in the lower mantle. This type of correction is also shown to be applicable to experimental measurements of elastic constants in situations where deviatoric stresses can develop, such as in diamond-anvil cells.

Structural basis of GC-1 selectivity for thyroid hormone receptor isoforms

Background: Thyroid receptors, TRa and TR beta, are involved in important physiological functions such as metabolism, cholesterol level and heart activities. Whereas metabolism increase and cholesterol level lowering could be achieved by TR beta isoform activation, TRa activation affects heart rates. Therefore, beta-selective thyromimetics have been developed as promising drug-candidates for treatment of obesity and elevated cholesterol level. GC-1 [ 3,5-dimethyl-4-(4'-hydroxy- 3'-isopropylbenzyl)-phenoxy acetic acid] has ability to lower LDL cholesterol with 600-to 1400-fold more potency and approximately two-to threefold more efficacy than atorvastatin (Lipitor(C)) in studies in rats, mice and monkeys. Results: To investigate GC-1 specificity, we solved crystal structures and performed molecular dynamics simulations of both isoforms complexed with GC-1. Crystal structures reveal that, in TRa Arg228 is observed in multiple conformations, an effect triggered by the differences in the interactions between GC-1 and Ser277 or the corresponding asparagine (Asn331) of TR beta. The corresponding Arg282 of TR beta is observed in only one single stable conformation, interacting effectively with the ligand. Molecular dynamics support this model: our simulations show that the multiple conformations can be observed for the Arg228 in TR alpha, in which the ligand interacts either strongly with the ligand or with the Ser277 residue. In contrast, a single stable Arg282 conformation is observed for TR beta, in which it strongly interacts with both GC-1 and the Asn331. Conclusion: Our analysis suggests that the key factors for GC-1 selectivity are the presence of an oxyacetic acid ester oxygen and the absence of the amino group relative to T(3). These results shed light into the beta-selectivity of GC-1 and may assist the development of new compounds with potential as drug candidates to the treatment of hypercholesterolemia and obesity.

Phase transitions and spatially ordered counterion association in ionic-lipid membranes: A statistical model

We propose a statistical model to account for the gel-fluid anomalous phase transitions in charged bilayer- or lamellae-forming ionic lipids. The model Hamiltonian comprises effective attractive interactions to describe neutral-lipid membranes as well as the effect of electrostatic repulsions of the discrete ionic charges on the lipid headgroups. The latter can be counterion dissociated (charged) or counterion associated (neutral), while the lipid acyl chains may be in gel (low-temperature or high-lateral-pressure) or fluid (high-temperature or low-lateral-pressure) states. The system is modeled as a lattice gas with two distinct particle types-each one associated, respectively, with the polar-headgroup and the acyl-chain states-which can be mapped onto an Ashkin-Teller model with the inclusion of cubic terms. The model displays a rich thermodynamic behavior in terms of the chemical potential of counterions (related to added salt concentration) and lateral pressure. In particular, we show the existence of semidissociated thermodynamic phases related to the onset of charge order in the system. This type of order stems from spatially ordered counterion association to the lipid headgroups, in which charged and neutral lipids alternate in a checkerboard-like order. Within the mean-field approximation, we predict that the acyl-chain order-disorder transition is discontinuous, with the first-order line ending at a critical point, as in the neutral case. Moreover, the charge order gives rise to continuous transitions, with the associated second-order lines joining the aforementioned first-order line at critical end points. We explore the thermodynamic behavior of some physical quantities, like the specific heat at constant lateral pressure and the degree of ionization, associated with the fraction of charged lipid headgroups.

Hydrogen bond interactions between acetone and supercritical water

Hydrogen bond interactions between acetone and supercritical water are investigated using a combined and sequential Monte Carlo/quantum mechanics (S-MC/QM) approach. Simulation results show a dominant presence of con. gurations with one hydrogen bond for different supercritical states, indicating that this specific interaction plays an important role on the solvation properties of acetone in supercritical water. Using QM MP2/aug-cc-pVDZ the calculated average interaction energy reveals that the hydrogen-bonded acetone-water complex is energetically more stable under supercritical conditions than ambient conditions and its stability is little affected by variations of temperature and/or pressure. All average results reported here are statistically converged.

The isotropic nuclear magnetic shielding constants of acetone in supercritical water: A sequential Monte Carlo/quantum mechanics study including solute polarization

The nuclear isotropic shielding constants sigma((17)O) and sigma((13)C) of the carbonyl bond of acetone in water at supercritical (P=340.2 atm and T=673 K) and normal water conditions have been studied theoretically using Monte Carlo simulation and quantum mechanics calculations based on the B3LYP/6-311++G(2d,2p) method. Statistically uncorrelated configurations have been obtained from Monte Carlo simulations with unpolarized and in-solution polarized solute. The results show that solvent effects on the shielding constants have a significant contribution of the electrostatic interactions and that quantitative estimates for solvent shifts of shielding constants can be obtained modeling the water molecules by point charges (electrostatic embedding). In supercritical water, there is a decrease in the magnitude of sigma((13)C) but a sizable increase in the magnitude of sigma((17)O) when compared with the results obtained in normal water. It is found that the influence of the solute polarization is mild in the supercritical regime but it is particularly important for sigma((17)O) in normal water and its shielding effect reflects the increase in the average number of hydrogen bonds between acetone and water. Changing the solvent environment from normal to supercritical water condition, the B3LYP/6-311++G(2d,2p) calculations on the statistically uncorrelated configurations sampled from the Monte Carlo simulation give a (13)C chemical shift of 11.7 +/- 0.6 ppm for polarized acetone in good agreement with the experimentally inferred result of 9-11 ppm. (C) 2008 American Institute of Physics.

Piezomagnetic behavior of Co-doped ZnO nanoribbons

Using ab initio total energy calculations, we show that bilayer systems of ZnO nanoribbons, (ZnO)(2)NR, doped with Co atoms exhibit a piezomagnetic behavior. We find the formation of energetically stable zigzag chains of Co atoms along the edge sites of (ZnO)(2)NR's, Co(Zn(chain))-(ZnO)(2)NR. At the ground state, the antiferromagnetic and the ferromagnetic states are very close in energy, whereas upon longitudinal stretch, parallel to the nanoribbon growth direction, it becomes ferromagnetic. Further electronic structure calculations indicate that not only the magnetic state but also the electronic structure of CoZn(chain)-(ZnO)(2)NR can be tuned by the mechanical stretch. In this case, we find that stretched NR's exhibit dispersive unpaired electronic states within the (ZnO)(2)NR band gap.

Splitting of the zero-energy edge states in bilayer graphene

Bilayer graphene nanoribbons with zigzag termination are studied within the tight-binding model. We also include single-site electron-electron interactions via the Hubbard model within the unrestricted Hartree-Fock approach. We show that either the interactions between the outermost edge atoms or the presence of a magnetic order can cause a splitting of the zero-energy edge states. Two kinds of edge alignments are considered. For one kind of edge alignment (?) the system is nonmagnetic unless the Hubbard parameter U becomes greater than a critical value Uc. For the other kind of edge alignment (?) the system is magnetic for any U>0. Our results agree very well with ab initio density functional theory calculations.

Anomalous lattice parameter of magnetic semiconductor alloys

The addition of transition metals to III-V semiconductors radically changes their electronic, magnetic, and structural properties. We show by ab initio calculations that in contrast to the conventional semiconductor alloys, the lattice parameter in magnetic semiconductor alloys, including those with diluted concentration, strongly deviates from Vegard's law. We find a direct correlation between the magnetic moment and the anion-transition metal bond lengths and derive a simple and general formula that determines the lattice parameter of a particular magnetic semiconductor by considering both the composition and magnetic moment. This dependence can explain some experimentally observed anomalies and stimulate other kind of investigations.

Electronic, structural, and transport properties of Ni-doped graphene nanoribbons

We have investigated the electronic and transport properties of zigzag Ni-adsorbed graphene nanoribbons (Ni/GNRs) using ab initio calculations. We find that the Ni adatoms lying along the edge of zigzag GNRs represent the energetically most stable configuration, with an energy difference of approximately 0.3 eV when compared to the adsorption in the middle of the ribbon. The carbon atoms at the ribbon edges still present nonzero magnetic moments as in the pristine GNR even though there is a quenching by a factor of almost five in the value of the local magnetic moments at the C atoms bonded to the Ni. This quenching decays relatively fast and at approximately 9 A from the Ni adsorption site the magnetic moments have already values close to the pristine ribbon. At the opposite edge and at the central carbon atoms the changes in the magnetic moments are negligible. The energetic preference for the antiparallel alignment between the magnetization at the opposite edges of the ribbon is still maintained upon Ni adsorption. We find many Ni d-related states within an energy window of 1 eV above and below the Fermi energy, which gives rise to a spin-dependent charge transport. These results suggest the possibility of manufacturing spin devices based on GNRs doped with Ni atoms.

Defect-induced effects on carrier migration through one-dimensional poly(para-phenylenevinylene) chains

Defects in one-dimensional (1D) systems can be intrinsically distinct from its three-dimensional counterparts, and polymer films are good candidates for showing both extremes that are difficult to individuate in the experimental data. We study theoretically the impact of simple hydrogen and oxygen defects on the electron transport properties of one-dimensional poly(para-phenylenevinylene) chains through a multiscale technique, starting from classical structural simulations for crystalline films to extensive ab initio calculations within density functional theory for the defects in single crystalline-constrained chains. The most disruptive effect on carrier transport comes from conjugation breaking imposed by the overcoordination of a carbon atom in the vinyl group independently from the chemical nature of the defect. The particular case of the [C=O] (keto-defect) shows in addition unexpected electron-hole separation, suggesting that the experimentally detected photoluminescence bleaching and photoconductivity enhancement could be due to exciton dissociation caused by the 1D characteristics of the defect.

Optical absorption and electronic band structure first-principles calculations of alpha-glycine crystals

Light absorption of alpha-glycine crystals grown by slow evaporation at room temperature was measured, indicating a 5.11 +/- 0.02 eV energy band gap. Structural, electronic, and optical absorption properties of alpha-glycine crystals were obtained by first-principles quantum mechanical calculations using density functional theory within the generalized gradient approximation in order to understand this result. To take into account the contribution of core electrons, ultrasoft and norm-conserving pseudopotentials, as well as an all electron approach were considered to compute the electronic density of states and band structure of alpha-glycine crystals. They exhibit three indirect energy band gaps and one direct Gamma-Gamma energy gap around 4.95 eV. The optical absorption related to transitions between the top of the valence band and the bottom of the conduction band involves O 2p valence states and C, O 2p conduction states, with the carboxyl group contributing significantly to the origin of the energy band gap. The calculated optical absorption is highly dependent on the polarization of the incident radiation due to the spatial arrangement of the dipolar glycine molecules; in the case of a polycrystalline sample, the first-principles calculated optical absorption is in good agreement with the measurement when a rigid energy shift is applied.

Elementary Excitations at Magnetic Surfaces and Their Spin Dependence

The elementary surface excitations are studied by spin-polarized electron energy loss spectroscopy on a prototype oxide surface [an oxygen passivated Fe(001)-p(1 x 1) surface], where the various excitations coexist. For the first time, the surface phonons and magnons are measured simultaneously and are distinguished based on their different spin nature. The dispersion relation of all excitations is probed over the entire Brillouin zone. The different phonon modes observed in our experiment are described by means of ab initio calculations.

Contrasting LH-HH subband splitting of strained quantum wells grown along [001] and [113] directions

Contrasting responses for the temperature tuning of the electronic structure in semiconductor quantum wells are discussed for heterolayered structures grown along (001) and (113) directions. The temperature affects the strain modulation of the deformation potentials and the effective optical gap is tuned along with the intersub-band splitting in the valence band. A multiband theoretical model accounts for the characterization of the electronic structure, highlighting the main qualitative and quantitative differences between the two systems under study. The microscopic source of strain fields and the detailed mapping of their distribution are provided by a simulation using classical molecular-dynamics technics.

Intermediate motions and dipolar couplings as studied by Lee-Goldburg cross-polarization NMR: Hartmann-Hahn matching profiles

In this article, we evaluate the use of simple Lee-Goldburg cross-polarization (LG-CP) NMR experiments for obtaining quantitative information of molecular motion in the intermediate regime. In particular, we introduce the measurement of Hartmann-Hahn matching profiles for the assessment of heteronuclear dipolar couplings as well as dynamics as a reliable and robust alternative to the more common analysis of build-up curves. We have carried out dynamic spin dynamics simulations in order to test the method's sensitivity to intermediate motion and address its limitations concerning possible experimental imperfections. We further demonstrate the successful use of simple theoretical concepts, most prominently Anderson-Weiss (AW) theory, to analyze the data. We further propose an alternative way to estimate activation energies of molecular motions, based upon the acquisition of only two LG-CP spectra per temperature at different temperatures. As experimental tests, molecular jumps in imidazole methyl sulfonate, trimethylsulfoxonium iodide, and bisphenol A polycarbonate were investigated with the new method.