Nicola N. Marzari

Maximally localized Wannier functions for entangled energy bands

We present a method for obtaining well-localized Wannier-like functions (WFs) for energy bands that are attached to or mixed with other bands. The present scheme removes the limitation of the usual maximally localized WFs method [N. Marzari and D. Vanderbilt, Phys. Rev. B 56, 12 847 (1997)] that the bands of interest should form an isolated group, separated by gaps from higher and lower bands everywhere in the Brillouin zone. An energy window encompassing N bands of interest is specified by the user, and the algorithm then proceeds to disentangle these from the remaining bands inside the window by filtering out an optimally connected N-dimensional subspace. This is achieved by minimizing a functional that measures the subspace dispersion across the Brillouin zone. The maximally localized WFs for the optimal subspace are then obtained via the algorithm of Marzari and Vanderbilt. The method, which functions as a postprocessing step using the output of conventional electronic-structure codes, is applied to the s and d bands of copper, and to the valence and low-lying conduction bands of silicon. For the low-lying nearly-free-electron bands of copper we find WFs which are centered at the tetrahedral-interstitial sites, suggesting an alternative tight-binding parametrization.

A classical and ab initio study of the interaction of the myosin triphosphate binding domain with ATP.

We used classical molecular mechanics (MM) simulations and quantum mechanical (QM) structural relaxations to examine the active site of myosin when bound to ATP. Two conformations of myosin have been determined by x-ray crystallography. In one, there is no direct interaction between switch 2 and the nucleotide (open state). In the other (closed state), the universally conserved switch 2 glycine forms a hydrogen bond with a gamma-phosphate oxygen. MM simulations indicate that the two states are thermodynamically stable and allow us to investigate the extent to which the P-loop, switch 1, and switch 2 are involved in hydrolysis. We find that the open structure has a higher affinity for ATP than the closed structure, and that ATP is distorted toward a transition state by interactions with the protein. We also examine how the structure of the binding site changes with either MgATP or CaATP as the nucleotide in myosin in the open conformer. Our analyses suggest that higher CaATPase rates occur because the leaving phosphate (P(i)) group is more weakly bound and dissociation occurs faster. Finally, we validate the use of a particular formulation of a QM methodology (Car-Parrinello) to further refine the structures of the active site.