Local order and magnetic field effects on the electronic properties of disordered binary alloys in the quantum site percolation limit

Electronic properties of disordered binary alloys are studied via the calculation of the average Density of States (DOS) in two and three dimensions. We propose a new approximate scheme that allows for the inclusion of local order effects in finite geometries and extrapolates the behavior of infinite systems following finite-size scaling ideas. We particularly investigate the limit of the Quantum Site Percolation regime described by a tight-binding Hamiltonian. This limit was chosen to probe the role of short range order (SRO) properties under extreme conditions. The method is numerically highly efficient and asymptotically exact in important limits, predicting the correct DOS structure as a function of the SRO parameters. Magnetic field effects can also be included in our model to study the interplay of local order and the shifted quantum interference driven by the field. The average DOS is highly sensitive to changes in the SRO properties and striking effects are observed when a magnetic field is applied near the segregated regime. The new effects observed are twofold: there is a reduction of the band width and the formation of a gap in the middle of the band, both as a consequence of destructive interference of electronic paths and the loss of coherence for particular values of the magnetic field. The above phenomena are periodic in the magnetic flux. For other limits that imply strong localization, the magnetic field produces minor changes in the structure of the average DOS. © World Scientific Publishing Company.

Ordered phases in the extended Hubbard model

A novel method to probe the diverse phases for the extended Hubbard model (EHM), including the correlated hopping term, is presented. We extend an effective medium approach [1] to a bipartite lattice, allowing for charge- and/or spin-ordered phases. We calculate the necessary correlation functions to build the EHM phase diagram.

Efficient intrinsic spin relaxation in graphene due to flexural distortions

We propose an intrinsic spin scattering mechanism in graphene originated by the interplay of atomic spin-orbit interaction and the local curvature induced by flexural distortions of the atomic lattice. Starting from a multiorbital tight-binding Hamiltonian with spin-orbit coupling considered non-perturbatively, we derive an effective Hamiltonian for the spin scattering of the Dirac electrons due to flexural distortions. We compute the spin lifetime due to both flexural phonons and ripples and we find values in the 1-10 ns range at room temperature. The proposed mechanism dominates the spin relaxation in high mobility graphene samples and should also apply to other planar aromatic compounds.

Anisotropic intrinsic spin relaxation in graphene due to flexural distortions

We propose an intrinsic spin scattering mechanism in graphene originated by the interplay of atomic spin-orbit interaction and the local curvature induced by flexural distortions of the atomic lattice. Starting from a multiorbital tight-binding Hamiltonian with spin-orbit coupling considered nonperturbatively, we derive an effective Hamiltonian for the spin scattering of the Dirac electrons due to flexural distortions. We compute the spin lifetime due to both flexural phonons and ripples and we find values in the microsecond range at room temperature. Interestingly, this mechanism is anisotropic on two counts. First, the relaxation rate is different for off-plane and in-plane spin quantization axis. Second, the spin relaxation rate depends on the angle formed by the crystal momentum with the carbon-carbon bond. In addition, the spin lifetime is also valley dependent. The proposed mechanism sets an upper limit for spin lifetimes in graphene and will be relevant when samples of high quality can be fabricated free of extrinsic sources of spin relaxation.

Molecular modelling of organic functional materials

The investigation of the mechanisms lying behind the (photo-)chemical processes is fundamental to address and improve the design of new organic functional materials. In many cases, dynamics simulations represent the only tool to capture the system properties emerging from complex interactions between many molecules. Despite the outstanding progresses in calculation power, the only way to carry out such computational studies is to introduce several approximations with respect to a fully quantum mechanical (QM) description. This thesis presents an approach that combines QM calculations with a classical Molecular Dynamics (MD) approach by means of accurate QM-derived force fields. It is based on a careful selection of the most relevant molecular degrees of freedom, whose potential energy surface is calculated at QM level and reproduced by the analytic functions of the force field, as well as by an accurate tuning of the approximations introduced in the model of the process to be simulated. This is made possible by some tools developed purposely, that allow to obtain and test the FF parameters through comparison with the QM frequencies and normal modes. These tools were applied in the modelling of three processes: the npi* photoisomerisation of azobenzene, where the FF description was extended to the excited state too and the non-adiabatic events were treated stochastically with Tully fewest switching algorithm; the charge separation in donors-acceptors bulk heterojunction organic solar cells, where a tight-binding Hamiltonian was carefully parametrised and solved by means of a code, also written specifically; the effect of the protonation state on the photoisomerisation quantum yield of the aryl-azoimidazolium unit of the axle molecule of a rotaxane molecular shuttle. In each case, the QM-based MD models that were specifically developed gave noteworthy information about the investigated phenomena, proving to be a fundamental key for a deeper comprehension of several experimental evidences.

Charge transfer in DNA: hole charge is confined to a single base pair due to solvation effects

We include solvation effects in tight-binding Hamiltonians for hole states in DNA. The corresponding linear-response parameters are derived from accurate estimates of solvation energy calculated for several hole charge distributions in DNA stacks. Two models are considered: (A) the correction to a diagonal Hamiltonian matrix element depends only on the charge localized on the corresponding site and (B) in addition to this term, the reaction field due to adjacent base pairs is accounted for. We show that both schemes give very similar results. The effects of the polar medium on the hole distribution in DNA are studied. We conclude that the effects of polar surroundings essentially suppress charge delocalization in DNA, and hole states in (GC)n sequences are localized on individual guanines

Charge transfer in DNA: hole charge is confined to a single base pair due to solvation effects

We include solvation effects in tight-binding Hamiltonians for hole states in DNA. The corresponding linear-response parameters are derived from accurate estimates of solvation energy calculated for several hole charge distributions in DNA stacks. Two models are considered: (A) the correction to a diagonal Hamiltonian matrix element depends only on the charge localized on the corresponding site and (B) in addition to this term, the reaction field due to adjacent base pairs is accounted for. We show that both schemes give very similar results. The effects of the polar medium on the hole distribution in DNA are studied. We conclude that the effects of polar surroundings essentially suppress charge delocalization in DNA, and hole states in (GC)n sequences are localized on individual guanines

Three-dimensional structure of ribonuclease T-1 complexed with an isosteric phosphonate substrate analogue of GpU: Alternate substrate binding modes and catalysis

The X-ray crystal structure of a complex between ribonuclease T-1 and guanylyl(3'-6')-6'-deoxyhomouridine (GpcU) has been determined at 2.0 Angstrom resolution. This Ligand is an isosteric analogue of the minimal RNA substrate, guanylyl(3'-5')uridine (GpU), where a methylene is substituted for the uridine 5'-oxygen atom. Two protein molecules are part of the asymmetric unit and both have a GpcU bound at the active site in the same manner. The protein-protein interface reveals an extended aromatic stack involving both guanines and three enzyme phenolic groups. A third GpcU has its guanine moiety stacked on His92 at the active site on enzyme molecule A and interacts with GpcU on molecule B in a neighboring unit via hydrogen bonding between uridine ribose 2'- and 3'-OH groups. None of the uridine moieties of the three GpcU molecules in the asymmetric unit interacts directly with the protein. GpcU-active-site interactions involve extensive hydrogen bonding of the guanine moiety at the primary recognition site and of the guanosine 2'-hydroxyl group with His40 and Glu58. on the other hand, the phosphonate group is weakly bound only by a single hydrogen bond with Tyr38, unlike ligand phosphate groups of other substrate analogues and 3'-GMP, which hydrogen-bonded with three additional active-site residues. Hydrogen bonding of the guanylyl 2'-OH group and the phosphonate moiety is essentially the same as that recently observed for a novel structure of a RNase T-1-3'-GMP complex obtained immediately after in situ hydrolysis of exo-(S-p)-guanosine 2',3'-cyclophosphorothioate [Zegers et al. (1998) Nature Struct. Biol. 5, 280-283]. It is likely that GpcU at the active site represents a nonproductive binding mode for GpU [:Steyaert, J., and Engleborghs (1995) fur. J. Biochem. 233, 140-144]. The results suggest that the active site of ribonuclease T-1 is adapted for optimal tight binding of both the guanylyl 2'-OH and phosphate groups (of GpU) only in the transition state for catalytic transesterification, which is stabilized by adjacent binding of the leaving nucleoside (U) group.

Temperature dependence of structural and electronic properties of the spin-liquid candidate κ-(BEDT-TTF) 2Cu 2(CN) 3

We investigate the effect that the temperature dependence of the crystal structure of a two-dimensional organic charge-transfer salt has on the low-energy Hamiltonian representation of the electronic structure. For that, we determine the crystal structure of κ-(BEDT-TTF) 2Cu 2(CN) 3 for a series of temperatures between T=5 and 300 K by single crystal X-ray diffraction and analyze the evolution of the electronic structure with temperature by using density functional theory and tight binding methods. We find a considerable temperature dependence of the corresponding triangular lattice Hubbard Hamiltonian parameters. We conclude that even in the absence of a change of symmetry, the temperature dependence of quantities like frustration and interaction strength can be significant and should be taken into account. © 2012 American Physical Society.

Numerical Modelling of Graphene Nanoribbon-fets for Analog and Digital Applications

Graphene, that is a monolayer of carbon atoms arranged in a honeycomb lattice, has been isolated only recently from graphite. This material shows very attractive physical properties, like superior carrier mobility, current carrying capability and thermal conductivity. In consideration of that, graphene has been the object of large investigation as a promising candidate to be used in nanometer-scale devices for electronic applications. In this work, graphene nanoribbons (GNRs), that are narrow strips of graphene, for which a band-gap is induced by the quantum confinement of carriers in the transverse direction, have been studied. As experimental GNR-FETs are still far from being ideal, mainly due to the large width and edge roughness, an accurate description of the physical phenomena occurring in these devices is required to have valuable predictions about the performance of these novel structures. A code has been developed to this purpose and used to investigate the performance of 1 to 15-nm wide GNR-FETs. Due to the importance of an accurate description of the quantum effects in the operation of graphene devices, a full-quantum transport model has been adopted: the electron dynamics has been described by a tight-binding (TB) Hamiltonian model and transport has been solved within the formalism of the non-equilibrium Green's functions (NEGF). Both ballistic and dissipative transport are considered. The inclusion of the electron-phonon interaction has been taken into account in the self-consistent Born approximation. In consideration of their different energy band-gap, narrow GNRs are expected to be suitable for logic applications, while wider ones could be promising candidates as channel material for radio-frequency applications.

Sperm-binding fibronectin type II-module proteins are genetically linked and functionally related

Fibronectin type II (Fn2) module-containing proteins in the male genital tract are characterized by different numbers of Fn2 modules. Predominantly two classes exist which are distinct by having either two or four Fn2 modules. Minor variants with three Fn2 modules were also found in the human and the porcine epididymis. To reveal their relationship, mRNAs and proteins of representatives of these classes were studied in human, in Sus scrofa, and in rodents. Adult boars expressed members of both classes, i.e. ELSPBP1 and pB1, in subsequent regions of the epididymis, and both were under androgenic control. Human and rodent epididymides, on the other hand, alternatively contained only representatives of one of these two classes, i.e. ELSPBP1 in the human and two different pB1-related counterparts in rodents. ELSPBP1 and pB1-related genomic sequences were closely linked in chromosomal regions HSA 19q and SSC 6 q11-q21; conserved synteny between these regions is well established. On the other hand, in a syntenic region on mouse chromosome 7, ELSPBP1-related sequences were lacking. Tight binding to the sperm membrane via a choline-mediated mechanism was a common feature of the two classes of Fn2-module proteins, suggesting related function(s). However, differences in their regionalized expression patterns along the male genital tract as well as in association sites on the sperm surface suggested a species-specific sequential order in sperm binding.

Cellulose-binding domains promote hydrolysis of different sites on crystalline cellulose

The cohesin-dockerin interaction in Clostridium thermocellum cellulosome mediates the tight binding of cellulolytic enzymes to the cellulosome-integrating protein CipA. Here, this interaction was used to study the effect of different cellulose-binding domains (CBDs) on the enzymatic activity of C. thermocellum endoglucanase CelD (1,4-β-d endoglucanase, EC3.2.1.4) toward various cellulosic substrates. The seventh cohesin domain of CipA was fused to CBDs originating from the Trichoderma reesei cellobiohydrolases I and II (CBDCBH1 and CBDCBH2) (1,4-β-d glucan-cellobiohydrolase, EC3.2.1.91), from the Cellulomonas fimi xylanase/exoglucanase Cex (CBDCex) (β-1,4-d glucanase, EC3.2.1.8), and from C. thermocellum CipA (CBDCipA). The CBD-cohesin hybrids interacted with the dockerin domain of CelD, leading to the formation of CelD-CBD complexes. Each of the CBDs increased the fraction of cellulose accessible to hydrolysis by CelD in the order CBDCBH1 < CBDCBH2 ≈ CBDCex < CBDCipA. In all cases, the extent of hydrolysis was limited by the disappearance of sites accessible to CelD. Addition of a batch of fresh cellulose after completion of the reaction resulted in a new burst of activity, proving the reversible binding of the intact complexes despite the apparent binding irreversibility of some CBDs. Furthermore, burst of activity also was observed upon adding new batches of CelD–CBD complexes that contained a CBD differing from the first one. This complementation between different CBDs suggests that the sites made available for hydrolysis by each of the CBDs are at least partially nonoverlapping. The only exception was CBDCipA, whose sites appeared to overlap all of the other sites.

Identification of a second aryl phosphate-binding site in protein-tyrosine phosphatase 1B: A paradigm for inhibitor design

The structure of the catalytically inactive mutant (C215S) of the human protein-tyrosine phosphatase 1B (PTP1B) has been solved to high resolution in two complexes. In the first, crystals were grown in the presence of bis-(para-phosphophenyl) methane (BPPM), a synthetic high-affinity low-molecular weight nonpeptidic substrate (Km = 16 μM), and the structure was refined to an R-factor of 18.2% at 1.9 Å resolution. In the second, crystals were grown in a saturating concentration of phosphotyrosine (pTyr), and the structure was refined to an R-factor of 18.1% at 1.85 Å. Difference Fourier maps showed that BPPM binds PTP1B in two mutually exclusive modes, one in which it occupies the canonical pTyr-binding site (the active site), and another in which a phosphophenyl moiety interacts with a set of residues not previously observed to bind aryl phosphates. The identification of a second pTyr molecule at the same site in the PTP1B/C215S–pTyr complex confirms that these residues constitute a low-affinity noncatalytic aryl phosphate-binding site. Identification of a second aryl phosphate binding site adjacent to the active site provides a paradigm for the design of tight-binding, highly specific PTP1B inhibitors that can span both the active site and the adjacent noncatalytic site. This design can be achieved by tethering together two small ligands that are individually targeted to the active site and the proximal noncatalytic site.

Large spin splitting in the conduction band of transition metal dichalcogenide monolayers

We study the conduction band spin splitting that arises in transition metal dichalcogenide (TMD) semiconductor monolayers such as MoS2, MoSe2, WS2, and WSe2 due to the combination of spin-orbit coupling and lack of inversion symmetry. Two types of calculation are done. First, density functional theory (DFT) calculations based on plane waves that yield large splittings, between 3 and 30 meV. Second, we derive a tight-binding model that permits to address the atomic origin of the splitting. The basis set of the model is provided by the maximally localized Wannier orbitals, obtained from the DFT calculation, and formed by 11 atomiclike orbitals corresponding to d and p orbitals of the transition metal (W, Mo) and chalcogenide (S, Se) atoms respectively. In the resulting Hamiltonian, we can independently change the atomic spin-orbit coupling constant of the two atomic species at the unit cell, which permits to analyze their contribution to the spin splitting at the high symmetry points. We find that—in contrast to the valence band—both atoms give comparable contributions to the conduction band splittings. Given that these materials are most often n-doped, our findings are important for developments in TMD spintronics.

Derivation of the spin Hamiltonians for Fe in MgO

A method to calculate the effective spin Hamiltonian for a transition metal impurity in a non-magnetic insulating host is presented and applied to the paradigmatic case of Fe in MgO. In the first step we calculate the electronic structure employing standard density functional theory (DFT), based on generalized gradient approximation (GGA), using plane waves as a basis set. The corresponding basis of atomic-like maximally localized Wannier functions is derived and used to represent the DFT Hamiltonian, resulting in a tight-binding model for the atomic orbitals of the magnetic impurity. The third step is to solve, by exact numerical diagonalization, the N electron problem in the open shell of the magnetic atom, including both effects of spin–orbit and Coulomb repulsion. Finally, the low energy sector of this multi-electron Hamiltonian is mapped into effective spin models that, in addition to the spin matrices S, can also include the orbital angular momentum L when appropriate. We successfully apply the method to Fe in MgO, considering both the undistorted and Jahn–Teller (JT) distorted cases. Implications for the influence of Fe impurities on the performance of magnetic tunnel junctions based on MgO are discussed.