title: Existence of Atoms and Molecules in the Mean-Field Approximation of No-Photon Quantum Electrodynamics

abstract: The Bogoliubov-Dirac-Fock (BDF) model is the mean-field approximation of no-photon Quantum Electrodynamics. The present paper is devoted to the study of the minimization of the BDF energy functional under a charge constraint. An associated minimizer, if it exists, will usually represent the ground state of a system of $N$ electrons interacting with the Dirac sea, in an external electrostatic field generated by one or several fixed nuclei. We prove that such a minimizer exists when a binding (HVZ-type) condition holds. We also derive, study and interpret the equation satisfied by such a minimizer. Finally, we provide two regimes in which the binding condition is fulfilled, obtaining the existence of a minimizer in these cases. The first is the weak coupling regime for which the coupling constant $\alpha$ is small whereas $\alpha Z$ and the particle number $N$ are fixed. The second is the non-relativistic regime in which the speed of light tends to infinity (or equivalently $\alpha$ tends to zero) and $Z$, $N$ are fixed. We also prove that the electronic solution converges in the non-relativistic limit towards a Hartree-Fock ground state.


title: Symmetric symplectic spaces with Ricci-type curvature

abstract: We determine the isomorphism classes of symmetric symplectic manifolds of dimension at least 4 which are connected, simply-connected and have a curvature tensor which has only one non-vanishing irreducible component -- the Ricci tensor.


title: Algebraic Rieffel Induction, Formal Morita Equivalence, and Applications to Deformation Quantization

abstract: In this paper we consider algebras with involution over a ring C which is given by the quadratic extension by i of an ordered ring R. We discuss the *-representation theory of such *-algebras on pre-Hilbert spaces over C and develop the notions of Rieffel induction and formal Morita equivalence for this category analogously to the situation for C^*-algebras. Throughout this paper the notion of positive functionals and positive algebra elements will be crucial for all constructions. As in the case of C^*-algebras, we show that the GNS construction of *-representations can be understood as Rieffel induction and, moreover, that formal Morita equivalence of two *-algebras, which is defined by the existence of a bimodule with certain additional structures, implies the equivalence of the categories of strongly non-degenerate *-representations of the two *-algebras. We discuss various examples like finite rank operators on pre-Hilbert spaces and matrix algebras over *-algebras. Formal Morita equivalence is shown to imply Morita equivalence in the ring-theoretic framework. Finally we apply our considerations to deformation theory and in particular to deformation quantization and discuss the classical limit and the deformation of equivalence bimodules.


title: Realizing coalgebras over the Steenrod algebra

abstract: We describe algebraic obstruction theories for realizing an abstract coalgebra K_* over the mod p Steenrod algebra as the homology of a topological space, and for distinguishing between the p-homotopy types of different realizations. The theories are expressed in terms of the Quillen cohomology of K_*.


title: CW simplicial resolutions of spaces, with an application to loop spaces

abstract: We show how a certain type of CW simplicial resolutions of space by wedges of spheres may be constructed for any topological space, and how such resolutions yield an obstruction theory for a given space X to be a loop space.


title: U_{q}(sl_{2}) at fourth root of unity

abstract: In this paper, we study finite dimensional quotients of the Hopf algebra U_{q}(sl_{2}) at fourth root of unity. We give explicitly all the corresponding automorphisms, real forms and idempotents.


title: The global quantum duality principle

abstract: Let R be an integral domain, let a non-zero h in R be such that k := R/hR is a field, and let HA be the category of torsionless (or flat) Hopf algebras over R. We call H in HA a "quantized function algebra" (=QFA), resp. "quantized restricted universal enveloping algebras" (=QrUEA), at h if - roughly speaking - H/hH is the function algebra of a connected Poisson group, resp. the (restricted, if R/hR has positive characteristic) universal enveloping algebra of a (restricted) Lie bialgebra. Extending a result of Drinfeld, we establish an "inner" Galois' correspondence on HA, via two endofunctors, ( )^\vee and ( )' , of HA such that H^\vee is a QrUEA and H' is a QFA (for all H in HA). In addition: (a) the image of ( )^\vee , resp. of ( )' , is the full subcategory of all QrUEAs, resp. of all QFAs; (b) if p := Char(R/hR) = 0 , the restrictions of ( )^\vee to the QFAs and of ( )' to the QrUEAs yield equivalences inverse to each other; (c) if p = 0 , starting from a QFA over a Poisson group G, resp. from a QrUEA over a Lie bialgebra g, the functor ( )^\vee , resp. ( )' , gives a QrUEA, resp. a QFA, over the dual Lie bialgebra, resp. the dual Poisson group. Several, far-reaching applications are developed in detail in [Ga2-4].


title: Interpolation of subspaces and applications to exponential bases

abstract: This paper has been withdrawn by the author becouse the conjecture presented in this paper is false. The correct study of metrics in interpolation spaces for and applcation to nonharmonic Fourier series will be published in S.Petersburg Math. J. jointly with N. Kalton.


title: Multiple vertex coverings by specified induced subgraphs

abstract: Given graphs H_1,...,H_k, we study the minimum order of a graph G such that for each i, the induced copies of H_i in G cover V(G). We prove a general upper bound of twice the sum of the numbers m_i, where m_i is one less than the order of H_i. When k=2 and one graph is an independent set of size n, we determine the optimum within a constant. When k=2 and the graphs are a star and an independent set, we determine the answer exactly.


title: Linearization Problems for Lie Algebroids and Lie Groupoids

abstract: In this survey, we discuss a series of linearization problems--for Poisson structures, Lie algebroids, and Lie groupoids. The last problem involves a conjecture on the structure of proper groupoids. Attempting to prove this by the method of averaging leads to problems concerning almost actions of compact groups and almost invariant submanifolds for compact group actions. The paper ends with a discussion of possible extensions of the convexity theorems for momentum maps of hamiltonian actions of compact groups.


title: Omni-Lie Algebras

abstract: We show that the space R^n x gl(n,R) with a certain antisymmetric bracket operation contains all n-dimensional Lie algebras. The bracket does not satisfy the Jacobi identity, but it does satisfy it for subalgebras which are isotropic under a certain symmetric bilinear form with values in R^n. We ask what the corresponding "group-like" object should be. The bracket may be obtained by linearizing at a point the bracket on TM + T*M introduced by T. Courant for the definition of Dirac structures, a notion which encompasses Poisson structures, closed 2-forms, and foliations.


title: Large localizations of finite simple groups

abstract: A group homomorphism eta:H-->G is called a localization of H if every homomorphism phi:H-->G can be `extended uniquely' to a homomorphism Phi:G-->G in the sense that Phi eta=phi. Libman showed that a localization of a finite group need not be finite. This is exemplified by a well-known representation A_n-->SO_{n-1}(R) of the alternating group A_n, which turns out to be a localization for n even and n>9. Dror Farjoun asked if there is any upper bound in cardinality for localizations of A_n. In this paper we answer this question and prove, under the generalized continuum hypothesis, that every non abelian finite simple group H, has arbitrarily large localizations. This shows that there is a proper class of distinct homotopy types which are localizations of a given Eilenberg--Mac Lane space K(H,1) for any non abelian finite simple group H.


title: Dual Forms on Supermanifolds and Cartan Calculus

abstract: The complex of "stable forms" on supermanifolds is studied. Stable forms on $M$ are represented by certain Lagrangians of "copaths" (formal systems of equations, which may or may not specify actual surfaces) on $M\times\mathbb R^D$. Changes of $D$ give rise to stability isomorphisms. The Cartan--de Rham complex made of stable forms extends both in positive and negative degree and its positive half is isomorphic to the complex of forms defined as Lagrangians of paths. Considering the negative half is necessary, in particular, for homotopy invariance. We introduce analogs of exterior multiplication by covectors and of contraction with vectors. We find (anti)commutation relations for them. An analog of Cartan's homotopy identity is proved. Before stabilization it contains some "stability operator" $\sigma$.


title: Continuous selections with respect to extension dimension

abstract: The paper is devoted to generalization of well-known Michael's Selection theorem on the case of extension dimension.


title: An arithmetic group associated with a Pisot unit and its symbolic-dynamical representatiom

abstract: To a given Pisot unit $\beta$ we associate a finite abelian group whose size appears to be equal to the discriminant of $\beta$. We call it the Pisot group and find its representation in the two-sided $\beta$-compactum in the case of $\beta$ satisfying the relation $Fin(\beta)=\Bbb Z[\beta]\cap[0,1)$. As a motivation for the definition, we show that the Pisot group is the kernel of some important arithmetic coding of the toral automorphism given by the companion matrix naturally associated with $\beta$.


title: Introduction to the vertex algebra approach to mirror symmetry

abstract: The goal of this paper is to make the vertex operator algebra approach to mirror symmetry accessible to algebraic geometers. Compared to better-known approaches using moduli spaces of stable maps and special Lagrangian fibrations, this approach follows more closely the original line of thinking that lead to the discovery of mirror symmetry by physicists. The ultimate goal of the vertex algebra approach is to give precise mathematical definitions of N=(2,2) superconformal field theories called A and B models associated to any Calabi-Yau variety and then show that thus constructed theories are related by the mirror involution for all known examples of mirror symmetric varieties.


title: Infinite-dimensional algebraic varieties and proof of the Jacobian Conjecture

abstract: This paper has been withdrawn by the author, due an error in the proof of Proposion 2.13.


title: Cayley-Hamilton-Newton identities and quasitriangular Hopf algebras

abstract: In the framework of the Drinfeld theory of twists in Hopf algebras we construct quantum matrix algebras which generalize the Reflection Equation and the RTT algebras. Finite-dimensional representations of these algebras related to the theory of nonultralocal spin chains are presented. The Cayley-Hamilton-Newton identities are demonstrated. These identities allow to define the quantum spectrum for the quantum matrices. We mention possible applications of the new quantum matrix algebras to constructions of noncommutative analogs of Minkowski space and quantum Poincar\'e algebras.


title: Solutions of the KZB equations in genus \geq 1

abstract: We introduce a flat version of the KZB connection. This connection is defined on the complement of the locus of Weierstrass points on the moduli space of genus $g$ complex curves with marked points. We then give integral formulas for flat sections of this connection, using the parametrization of conformal blocks and of the KZB connection of our earlier paper math/9807145.


title: Torus actions, equivariant moment-angle complexes, and coordinate subspace arrangements

abstract: We show that the cohomology algebra of the complement of a coordinate subspace arrangement in m-dimensional complex space is isomorphic to the cohomology algebra of Stanley-Reisner face ring of a certain simplicial complex on m vertices. (The face ring is regarded as a module over the polynomial ring on m generators.) Then we calculate the latter cohomology algebra by means of the standard Koszul resolution of polynomial ring. To prove these facts we construct an equivariant with respect to the torus action homotopy equivalence between the complement of a coordinate subspace arrangement and the moment-angle complex defined by the simplicial complex. The moment-angle complex is a certain subset of a unit poly-disk in m-dimensional complex space invariant with respect to the action of an m-dimensional torus. This complex is a smooth manifold provided that the simplicial complex is a simplicial sphere, but otherwise has more complicated structure. Then we investigate the equivariant topology of the moment-angle complex and apply the Eilenberg-Moore spectral sequence. We also relate our results with well known facts in the theory of toric varieties and symplectic geometry.

categories: math.AT math.AC math.AG math.CO math.RA math.SG


title: The first main theorem on complements: from global to local

abstract: The aim of this paper is to clarify and generalize techniques of works alg-geom/9711024 (see also math.AG/9810097 and math.AG/9901004). Roughly speaking, we prove that for local Fano contractions the existence of complements can be reduced to the existence of complements for lower dimensional projective Fano varieties.


title: On Global existence for nonlinear wave equations outside of convex obstacles

abstract: We prove global existence for semilinear hyperbolic equations that satisfy the null condition of Christodoulou and Klainerman in the exterior of convex domains. We use a combination of the conformal method of Christodoulou and the direct method of Klainerman.


title: Negative results for Nikodym maximal functions and related oscillatory integrals in curved space

abstract: We expand on counterxamples of Bourgain showing how Nikodym maximal estimates and oscillatory integral estimates can break down in the non-Euclidean case. Our examples show the role that the non-existence of totally geodesic submanifolds can play in these problems.


title: Concerning Nikodym-type sets in 3-dimensional curved spaces

abstract: We obtain estimates for maximal functions that arise when one studies Nikodym-type sets. We also formulate a curvature condition that allows favorable estimates for these maximal functions.


title: Global Strichartz estimates for nontrapping perturbations of the Laplacian

abstract: The authors prove global Strichartz estimates for compact perturbations of the wave operator in odd dimensions when a non-trapping assumption is satisfied.


title: Averages over curves with torsion

abstract: We establish $L^p$ Sobolev mapping properties for averages over certain curves in $\R^3$, which improve upon the estimates obtained by $L^2-L^\infty$ interpolation.


title: Weighted Strichartz estimates and global existence for semilinear wave equations

abstract: We prove certain weighted Strichartz estimates and use these to prove a sharp theorem for global existence of small amplitude solutions of $\square u= |u|^p$, thus verifying the so-called "Strauss conjecture".


title: Null form estimates for (1/2,1/2) symbols and local existence for a quasilinear Dirichlet-wave equation

abstract: The authors show that bilinear estimates for null forms hold for Dirichlet-wave equations outside of convex obstacle. This generalizes results for the Euclidean case of Klainerman and Machedon, and of Sogge for the variable coefficient boundaryless case. The estimates are used to prove a local existence theorem for semilinear wave equations satisfying the null condition.


title: Gamma-functions of representations and lifting

abstract: Let G be the group of points of a quasi-split reductive algebraic group over a local field F. It follows from the local Langlands conjectures that to every non-trivial additive character of F and every representation of the Langlands dual group of F one can associate certain meromorphic function on the set of isomorphism classed of irreducible representations of G (which we call gamma-functions). In this paper we describe a general framework for an explicit construction of these gamma-functions. We make this idea precise in certain special cases. As a byproduct we give explicit conjectural formulas for the representation $l_E(\theta)$ of the group GL(n,F) associated to a character $\theta$ of the multiplicative group $E^*$ of a separable extension E of F of degree n. We also describe a conjectural analogue of the Poisson summation formula associated to a representation of the Langlands dual group, which implies the existence of the meromorphic continuation and functional equation for the corresponding automorphic L-functions.


title: The A^{(1)}_M automata related to crystals of symmetric tensors

abstract: A soliton cellular automaton associated with crystals of symmetric tensor representations of the quantum affine algebra U'_q(A^{(1)}_M) is introduced. It is a crystal theoretic formulation of the generalized box-ball system in which capacities of boxes and carriers are arbitrary and inhomogeneous. Scattering matrices of two solitons coincide with the combinatorial R matrices of U'_q(A^{(1)}_{M-1}). A piecewise linear evolution equation of the automaton is identified with an ultradiscrete limit of the nonautonomous discrete KP equation. A class of N soliton solutions is obtained through the ultradiscretization of soliton solutions of the latter.


title: Proof of the volume conjecture for torus knots

abstract: The volume conjecture, formulated recently by H. Murakami and J. Murakami, is proved for the case of torus knots.


title: Cotensor products of modules

abstract: Let C be a coalgebra over a field k and A its dual algebra. The category of C-comodules is equivalent to a category of A-modules. We use this to interpret the cotensor product M \square N of two comodules in terms of the appropriate Hochschild cohomology of the A-bimodule M \otimes N, when A is finite-dimensional, profinite, graded or differential-graded. The main applications are to Galois cohomology, comodules over the Steenrod algebra, and the homology of induced fibrations.


title: Hyperplane Arrangement Cohomology and Monomials in the Exterior Algebra

abstract: We show that if X is the complement of a complex hyperplane arrangement, then the homology of X has linear free resolution as a module over the exterior algebra on the first cohomology of X. We study invariants of X that can be deduced from this resolution. A key ingredient is a result of Aramova, Avramov, and Herzog [2000] on resolutions of monomial ideals in the exterior algebra. We give a new conceptual proof of this result.


title: Isomorphism classes of A-hypergeometric systems

abstract: For a finite set A of integral vectors, Gel'fand, Kapranov and Zelevinskii defined a system of differential equations with a parameter vector as a D-module, which system is called an A-hypergeometric (or a GKZ hypergeometric) system. Classifying the parameters according to the D-isomorphism classes of their corresponding A-hypergeometric systems is one of the most fundamental problems in the theory. In this paper we give a combinatorial answer for the problem under the assumption that the finite set A lies in a hyperplane off the origin, and illustrate it in two particularly simple cases: the normal case and the monomial curve case.


title: On the foundations of nonlinear generalized functions I

abstract: We construct a diffeomorphism invariant (Colombeau-type) differential algebra canonically containing the space of distributions in the sense of L. Schwartz. Employing differential calculus in infinite dimensional (convenient) vector spaces, previous attempts in this direction are unified and completed. Several classification results are achieved and applications to nonlinear differential equations involving singularities are given.


title: On the foundations of nonlinear generalized functions II

abstract: This paper gives a comprehensive analysis of algebras of Colombeau-type generalized functions in the range between the diffeomorphism-invariant quotient algebra $\mathcal{G}^d = \mathcal{E}_M/\mathcal{N}$ introduced in part I and Colombeau's original algebra $\mathcal{G}^e$. Three main results are established: First, a simple criterion describing membership in $\mathcal{N}$ (applicable to all types of Colombeau algebras) is given. Second, two counterexamples demonstrate that $\mathcal{G}^d$ is not injectively included in $\mathcal{G}^e$. Finally, it is shown that in the range ``between'' $\mathcal{G}^d$ and $\mathcal{G}^e$ only one more construction leads to a diffeomorphism invariant algebra. In analyzing the latter, several classification results essential for obtaining an intrinsic description of $\mathcal{G}^d$ on manifolds are derived.


title: A global theory of algebras of generalized functions

abstract: We present a geometric approach to defining an algebra $\hat{\mathcal G}(M)$ (the Colombeau algebra) of generalized functions on a smooth manifold $M$ containing the space ${\mathcal D}'(M)$ of distributions on $M$. Based on differential calculus in convenient vector spaces we achieve an intrinsic construction of $\hat{\mathcal G}(M)$. $\hat{\mathcal G}(M)$ is a{\em differential} algebra, its elements possessing Lie derivatives with respect to arbitrary smooth vector fields. Moreover, we construct a canonical linear embedding of ${\mathcal D}'(M)$ into $\hat{\mathcal G}(M)$ that renders ${\mathcal C}^\infty (M)$ a faithful subalgebra of $\hat{\mathcal G}(M)$. Finally, it is shown that this embedding commutes with Lie derivatives. Thus $\hat{\mathcal G}(M)$ retains all the distinguishing properties of the local theory in a global context.


title: Group analysis of differential equations and generalized functions

abstract: We present an extension of the methods of classical Lie group analysis of differential equations to equations involving generalized functions (in particular: distributions). A suitable framework for such a generalization is provided by Colombeau's theory of algebras of generalized functions. We show that under some mild conditions on the differential equations, symmetries of classical solutions remain symmetries for generalized solutions. Moreover, we introduce a generalization of the infinitesimal methods of group analysis that allows to compute symmetries of linear and nonlinear differential equations containing generalized function terms. Thereby, the group generators and group actions may be given by generalized functions themselves.


title: Commuting difference operators with elliptic coefficients from Baxter's vacuum vestors

abstract: For quantum integrable models with elliptic R-matrix, we construct the Baxter Q-operator in infinite-dimensional representations of the algebra of observables.


title: Regularized derivatives in a 2-dimensional model of self-interacting fields with singular data

abstract: The coupled Maxwell-Lorentz system describes feed-back action of electromagnetic fields in classical electrodynamics. When applied to point-charge sources (viewed as limiting cases of charged fluids) the resulting nonlinear weakly hyperbolic system lies beyond the scope of classical distribution theory. Using regularized derivatives in the framework of Colombeau algebras of generalized functions we analyze a two-dimensional analogue of the Maxwell-Lorentz system. After establishing existence and uniqueness of solutions in this setting we derive some results on distributional limits of solutions with $\delta$-like initial values.


title: Factorization theorem for the transfer function of a 2x2 operator matrix with unbounded couplings

abstract: We consider the analytic continuation of the transfer function associated with a 2x2 operator matrix having unbounded couplings into unphysical sheets of its Riemann surface. We construct a family of non-selfadjoint operators which factorize the transfer function and reproduce certain parts of its spectrum including the nonreal (resonance) spectrum situated in the unphysical sheets neighboring the physical sheet.

categories: math.SP math-ph math.FA math.MP nucl-th quant-ph


title: Classifying subcategories of modules

abstract: In this paper, we classify certain subcategories of modules over a ring R. A wide subcategory of R-modules is an Abelian subcategory of R-Mod that is closed under extensions. We give a complete classification of wide subcategories of finitely presented modules when R is a quotient of a regular commutative coherent ring by a finitely generated ideal. This includes all finitely presented algebras over a principal ideal domain, as well as polynomial rings on infinitely many variable over a PID. The classification is in terms of subsets of Spec R, and depends heavily on Thomason's classification of thick subcategories of small objects in the derived category. We also classify all wide subcategories closed under arbitrary coproducts for any Noetherian commutative ring R. These correspond to arbitrary subsets of Spec R, and this classification depends on Neeman's classification of localizing subcategories of the derived category.


title: A New Greedoid: The Family of Local Maximum Stable Sets of a Forest

abstract: A maximum stable set in a graph G is a stable set of maximum cardinality. S is a local maximum stable set if it is a maximum stable set of the subgraph of G spanned by the union of S and N(S), where N(S) is the neighborhood of S. One theorem of Nemhauser and Trotter Jr., working as a useful sufficient local optimality condition for the weighted maximum stable set problem, ensures that any local maximum stable set of G can be enlarged to a maximum stable set of G. In this paper we demonstrate that an inverse assertion is true for forests. Namely, we show that for any non-empty local maximum stable set S of a forest T there exists a local maximum stable set S1 of T, such that S1 is included in S and |S1| = |S| - 1. Moreover, as a further strengthening of both the theorem of Nemhauser and Trotter Jr. and its inverse, we prove that the family of all local maximum stable sets of a forest forms a greedoid on its vertex set.


title: Adiabatic Limits and Foliations

abstract: We use adiabatic limits to study foliated manifolds. The Bott connection naturally shows up as the adiabatic limit of Levi-Civita connections. As an application, we then construct certain natural elliptic operators associated to the foliation and present a direct geometric proof of a vanshing theorem of Connes[Co], which extends the Lichnerowicz vanishing theorem [L] to foliated manifolds with spin leaves, for what we call almost Riemannian foliations. Several new vanishing theorems are also proved by using our method.


title: Microstates free entropy and cost of equivalence relations

abstract: We define an analog of Voiculescu's free entropy for n-tuples of unitaries (u_{1},...,u_{n}) in a tracial von Neumann algebra M, normalizing a unital diffuse abelian subalgebra B in M. Using this quantity, we define the free dimension \delta_{0}(u_{1},..,u_{n}\btw B). This number depends on (u_{1},... ,u_{n}) only up ``orbit equivalence'' over B. In particular, if R is an measirable equivalence relation on [0,1] generated by n automorphisms \alpha_{1},...,\alpha_{n}, let u_{1},..., u_{n} be the unitaries implementing \alpha_{1},..,\alpha_{n} in the Feldman-Moore crossed product algebra M=W^{*}([0,1],R), and let B be the canonical copy of L^\infty functions on [0,1] inside M. In this way, we obtain an invariant \delta (R)=\delta _{0}(u_{1},...,u_{n},\btw B) of the equivalence relation (R). If R is treeable, \delta (R) coincides with the cost C(R) of R in the sense of Gaboriau. For a general equivalence relation R posessing a finite graphing, \delta(R)\leq C(R). Using the notion of free dimension, we define an dynamical entropy invariant for an automorphism of a measurable equivalence relation (or more generally of an r-discrete measure groupoid), and give examples.


title: Layered Multishift Coupling for use in Perfect Sampling Algorithms (with a primer on CFTP)

abstract: In this article we describe a new coupling technique which is useful in a variety of perfect sampling algorithms. A multishift coupler generates a random function f() so that for each real x, f(x)-x is governed by the same fixed probability distribution, such as a normal distribution. We develop the class of layered multishift couplers, which are simple and have several useful properties. For the standard normal distribution, for instance, the layered multishift coupler generates an f() which (surprisingly) maps an interval of length L to fewer than 2+L/2.35 points --- useful in applications which perform computations on each such image point. The layered multishift coupler improves and simplifies algorithms for generating perfectly random samples from several distributions, including the autogamma distribution, posterior distributions for Bayesian inference, and the steady state distribution for certain storage systems. We also use the layered multishift coupler to develop a Markov-chain based perfect sampling algorithm for the autonormal distribution. At the request of the organizers, we begin by giving a primer on CFTP (coupling from the past); CFTP and Fill's algorithm are the two predominant techniques for generating perfectly random samples using coupled Markov chains.


title: A Duality Theorem for Quantum Groupoids

abstract: We prove a duality theorem for quantum groupoid (weak Hopf algebra) actions that extends the well-known result for usual Hopf algebras.


title: Translated tori in the characteristic varieties of complex hyperplane arrangements

abstract: We give examples of complex hyperplane arrangements for which the top characteristic variety contains positive-dimensional irreducible components that do not pass through the origin of the character torus. These examples answer several questions of Libgober and Yuzvinsky. As an application, we exhibit a pair of arrangements for which the resonance varieties of the Orlik-Solomon algebra are (abstractly) isomorphic, yet whose characteristic varieties are not isomorphic. The difference comes from translated components, which are not detected by the tangent cone at the origin.


title: Seeley's Theory of Pseudodifferential Operators on Orbifolds

abstract: We present the theory of pseudodifferential operators acting on a vector orbibundle over an orbifold, construct the zeta function of an elliptic pseudodifferential operator and show the existence of a meromorphic extension to the complex plane with at most simple poles. We give formulas for generalized densities on the orbifold whose integrals compute the residues of the zeta function.


title: Toric Prevarieties and Subtorus Actions

abstract: Dropping separatedness in the definition of a toric variety, one obtains the more general notion of a toric prevariety. Toric prevarieties occur as ambient spaces in algebraic geometry and moreover they appear naturally as intermediate steps in quotient constructions. We first provide a complete description of the category of toric prevarieties in terms of convex-geometrical data, so-called systems of fans. In a second part, we consider actions of subtori H of the big torus of a toric prevariety X and investigate quotients for such actions. Using our language of systems of fans, we characterize existence of good prequotients for the action of H on X. Moreover, we show by means of an algorithmic construction that there always exists a toric prequotient for the action of H on X, that means an H-invariant toric morphism p from X to a toric prevariety Y such that every H-invariant toric morphism from X to a toric prevariety factors through p. Finally, generalizing a result of D. Cox, we prove that every toric prevariety occurs as the image of a categorical prequotient of an open toric subvariety of some complex affine space.


title: A `replacement sequence' method for finding the largest real root of an integer monic polynomial

abstract: To every integer monic polynomial of degree m can be associated a `replacement rule' that generates a word W* from another word W consisting of symbols belonging to a finite `alphabet' of size 2m. This rule applied iteratively on almost any initial word Wo, yields a sequence of words {Wi}. From acount of different symbols in the word Wi, one can obtain a rational approximate to the largest real root of the polynomial.


title: Q-multilinear Algebra

abstract: The Cayley-Hamilton-Newton theorem - which underlies the Newton identities and the Cayley-Hamilton identity - is reviewed, first, for the classical matrices with commuting entries, second, for two q-matrix algebras, the RTT-algebra and the RLRL-algebra. The Cayley-Hamilton-Newton identities for these q-algebras are related by the factorization map. A class of algebras M(R,F) is presented. The algebras M(R,F) include the RTT-algebra and the RLRL-algebra as particular cases. The algebra M(R,F) is defined by a pair of compatible matrices R and F. The Cayley-Hamilton-Newton theorem for the algebras M(R,F) is stated. A nontrivial example of a compatible pair is given.


title: Bifurcation of relative equilibria in mechanical systems with symmetry

abstract: The relative equilibria of a symmetric Hamiltonian dynamical system are the critical points of the so-called augmented Hamiltonian. The underlying geometric structure of the system is used to decompose the critical point equations and construct a collection of implicitly defined functions and reduced equations describing the set of relative equilibria in a neighborhood of a given relative equilibrium. The structure of the reduced equations is studied in a few relevant situations. In particular, a persistence result of Lerman and Singer [LS98] is generalized to the framework of Abelian proper actions. Also, a Hamiltonian version of the Equivariant Branching Lemma and a study of bifurcations with maximal isotropy are presented. An elementary example is presented to illustrate the use of this approach.


title: Percolation in the Hyperbolic Plane

abstract: This is a study of percolation in the hyperbolic plane and on regular tilings in the hyperbolic plane. The processes discussed include Bernoulli site and bond percolation on planar hyperbolic graphs, invariant dependent percolations on such graphs, and Poisson-Voronoi-Bernoulli percolation. We prove the existence of three distinct nonempty phases for the Bernoulli processes. In the first phase, $p\in(0, p_c]$, there are no unbounded clusters, but there is a unique infinite cluster for the dual process. In the second phase, $p\in(p_c,p_u)$, there are infinitely many unbounded clusters for the process and for the dual process. In the third phase, $p\in [p_u,1)$, there is a unique unbounded cluster, and all the clusters of the dual process are bounded. We also study the dependence of $p_c$ in the Poisson-Voronoi-Bernoulli percolation process on the intensity of the underlying Poisson process.


title: On $\alpha $-Square-Stable Graphs

abstract: The stability number of a graph G, denoted by alpha(G), is the cardinality of a maximum stable set, and mu(G) is the cardinality of a maximum matching in G. If alpha(G) + mu(G) equals its order, then G is a Koenig-Egervary graph. We call G an $\alpha $-square-stable graph, shortly square-stable, if alpha(G) = alpha(G*G), where G*G denotes the second power of G. These graphs were first investigated by Randerath and Wolkmann. In this paper we obtain several new characterizations of square-stable graphs. We also show that G is an square-stable Koenig-Egervary graph if and only if it has a perfect matching consisting of pendant edges. Moreover, we find that well-covered trees are exactly square-stable trees. To verify this result we give a new proof of one Ravindra's theorem describing well-covered trees.


title: Classification of infinite-dimensional simple groups of supersymmetries and quantum field theory

abstract: Talk given at the conference "Visions in Mathematics toward the year 2000", August 1999, Tel-Aviv.


title: Classification of positive definite lattices

abstract: In this paper we describe an algorithm for classifying orbits of vectors in Lorentzian lattices. The main point of this is that isomorphism classes of positive definite lattices in some genus often correspond to orbits of vectors in some Lorentzian lattice, so we can classify some positive definite lattices. As an application we give the classification of the 665 25-dimensional unimodular positive definite lattices and the 121 even 25 dimensional positive definite lattices of determinant 2. We also use this algorithm to show that there is a unique 26 dimensional unimodular positive definite lattice with no roots.


title: Global Stability of McKeithan's Kinetic Proofreading Model for T-Cell Receptor Signal Transduction

abstract: This note presents a self-contained and streamlined exposition of chemical-network results due to Feinberg, Horn, and Jackson. As an application, it shows the global asymptotic stability to equilibria of McKeithan's kinetic proofreading model for T-cell receptor signal transduction.


title: Some fibered and non-fibered links at infinity of hyperbolic complex line arrangements

abstract: Let F be R or C, d the dimension of F over R. Denote by P(F) either the affine plane A(F) or the hyperbolic plane H(F) over F. An arrangement L of k lines in P(F) (pairwise non-parallel in the hyperbolic case) has a link at infinity K(L) comprising k unknotted (d-1)-spheres in the (2d-1)-sphere, whose topology reflects the combinatorics of L `at infinity'. The class of links at infinity of affine F-line arrangements is properly included in the class of links at infinity of hyperbolic F-line arrangements. Many links at infinity of (essentially non-affine) connected hyperbolic C-line arrangements are far from being fibered; the proof is a direct construction, using ``Legendrian inscription''. In contrast, if the (affine or hyperbolic) R-line arrangement L is connected, then its complexification (an affine or hyperbolic C-line arrangement) has a fibered link at infinity; the proof uses A'Campo's divides.


title: Strong Connections and Chern-Connes Pairing in the Hopf-Galois Theory

abstract: We reformulate the concept of connection on a Hopf-Galois extension $B\subseteq P$ in order to apply it in computing the Chern-Connes pairing between the cyclic cohomology $HC^{2n} (B)$ and $K_0 (B)$. This reformulation allows us to show that a Hopf-Galois extension admitting a strong connection is projective and left faithfully flat. It also enables us to conclude that a strong connection is a Cuntz-Quillen-type bimodule connection. To exemplify the theory, we construct a strong connection (super Dirac monopole) to find out the Chern-Connes pairing for the super line bundles associated to super Hopf fibration.


title: A connected sum formula for the SU(3) Casson invariant

abstract: We provide a formula for the SU(3) Casson invariant for 3-manifolds given as the connected sum of two integral homology 3-spheres.


title: On the canonical threefolds with strictly nef anticanonical divisors

abstract: Any ample Cartier divisor D on a projective variety X is strictly nef (i.e. D.C>0 for any effective curve C on X). In general, the converse statement does not hold. But this is conjectured to be true for anticanonical divisors. The present paper establishes this fact for normal complex projective threefolds with canonical singularities. This result extends several previously known special cases. The proof rests mainly on sophisticated techniques of three dimensional birational geometry developed in the last two decades.


title: Computation of some examples of Brown's spectral measure in free probability

abstract: We use free probability techniques for computing spectra and Brown measures of some non hermitian operators in finite von Neumann algebras. Examples include u_n+u_oo where u_n and u_oo are the generators of Z_n and Z respectively, in the free product group Z_n*Z, or elliptic elements, of the form S_a+iS_b where S_a and S_b are free semi-circular elements of variance a and b. We give some pictorial evidence for connections with spectra of random matrices.


title: Quantum homogeneous spaces and quasi-Hopf algebras

abstract: We propose a formulation of the quantization problem of Manin quadruples, and show that a solution to this problem yields a quantization of the corresponding Poisson homogeneous spaces. We then solve both quantization problems in an example related to quantum spheres.


title: Scattering Spaces and a Decomposition of Continuous Spectral Subspace

abstract: We introduce the notion of scattering space $S_b^r$ for $N$-body quantum mechanical systems, where $b$ is a cluster decomposition with $2\le |b|\le N$ and $r$ is a real number $0\le r\le 1$. Utilizing these spaces, we give a decomposition of continuous spectral subspace by $S_b^1$ for $N$-body quantum systems with long-range pair potentials $V_\alpha^L(x_\alpha)=O(|x_\al|^{-\ep})$. This is extended to a decomposition by $S_b^r$ with $0\le r\le 1$ for some long-range case. We also prove a characterization of ranges of wave operators by $S_b^0$.


title: A Semiregularity Map for Modules and Applications to Deformations

abstract: This paper contains the details and complete proofs of our earlier announcement in math.AG/9907004 . We construct a general semiregularity map for algebraic cycles as asked for by S. Bloch in 1972. The existence of such a semiregularity map has well known consequences for the structure of the Hilbert scheme and for the variational Hodge conjecture. Aside from generalizing and extending considerably previously known results in this direction, we give new applications to deformations of modules that encompass, for example, results of Artamkin and Mukai. The formation of the semiregularity map here involves powers of the cotangent complex, Atiyah classes, and trace maps, and is defined not only for subspaces of manifolds but for perfect complexes on arbitrary complex spaces. It generalizes in particular Illusie's treatment of the Chern character to the analytic context and specializes to Bloch's earlier description of the semiregularity map for locally complete intersections as well as to the infinitesimal Abel-Jacobi map for submanifolds.


title: Calibrated embeddings in the special Lagrangian and coassociative cases

abstract: Every closed, oriented, real analytic Riemannian 3-manifold can be isometrically embedded as a special Lagrangian submanifold of a Calabi-Yau 3-fold, even as the real locus of an antiholomorphic, isometric involution. Every closed, oriented, real analytic Riemannian 4-manifold whose bundle of self-dual 2-forms is trivial can be isometrically embedded as a coassociative submanifold in a G_2-manifold, even as the fixed locus of an anti-G_2 involution. These results, when coupled with McLean's analysis of the moduli spaces of such calibrated submanifolds, yield a plentiful supply of examples of compact calibrated submanifolds with nontrivial deformation spaces.


title: Syzygies of Unimodular Lawrence Ideals

abstract: Infinite hyperplane arrangements whose vertices form a lattice are studied from the point of view of commutative algebra. The quotient of such an arrangement modulo the lattice action represents the minimal free resolution of the associated binomial ideal, which defines a toric subvariety in a product of projective lines. Connections to graphic arrangements and to Beilinson's spectral sequence are explored.


title: An elementary approach to the mapping class group of a surface

abstract: We consider an oriented surface S and a cellular complex X of curves on S, defined by Hatcher and Thurston in 1980. We prove by elementary means, without Cerf theory, that the complex X is connected and simply connected. From this we derive an explicit simple presentation of the mapping class group of S, following the ideas of Hatcher-Thurston and Harer.


title: Ranks of elliptic curves in families of quadratic twists

abstract: We show that the unboundedness of the ranks of the quadratic twists of an elliptic curve is equivalent to the divergence of certain infinite series.


title: A two variable Artin conjecture

abstract: Let a and b be non-zero rational numbers that are multiplicatively independent. We study the natural density of the set of primes p for which the subgroup of the multiplicative group of the finite field with p elements generated by (a\mod p) contains (b\mod p). It is shown that, under assumption of the generalized Riemann hypothesis (GRH), this density exists and equals a positive rational multiple of the universal constant S=\prod_{p prime}(1-p/(p^3-1)). An explicit value of the density is given under mild conditions on a and b. This extends and corrects earlier work of P.J. Stephens (1976). Our result, in combination with earlier work of the second author, allows us to deduce that any second order linear recurrence with reducible characteristic polynomial having integer elements, has a positive density of prime divisors (under GRH).


title: A Gelfand-Beurling formula for heights on endomorphism rings

abstract: In this paper we prove a limit formula for heights on the endomorphism ring of a vector space over a number field. This limit formula can be considered as the analogue for both the Gelfand-Beurling formula for the spectral radius on a Banach Algebra and Tate's averaging procedure for constructing canonical heights on Abelian Varieties.


title: On the density of primes in arithmetic progression having a prescribed primitive root

abstract: Let a,f and g be integers, with a and f coprime. Under the generalized Riemann hypothesis it follows from work of Hooley and Lenstra that the set of primes p such that p=a(mod f) and g is primitive root mod p has a natural density. In this note we explicitly evaluate this density and give some applications of this result.


title: Prime divisors of the Lagarias sequence

abstract: For integer a let us consider the sequence X_a={x_0,x_1,x_2,...} defined by x_0=a, x_1=1 and, for n>=1, x_{n+1}=x_n+x_{n-1}. We say that a prime p divides X_a if p divides at least one term of the sequence. It is easy to see that every prime p divides X_1, the sequence of Fibonacci numbers. Lagarias, using a technique involving the computation of degrees of various Kummerian extensions first employed by Hasse, showed in 1985 that X_2, the set of primes dividing some Lucas number has natural density 2/3 and posed as a challenge finding the density of prime divisors of X_3. In this paper we resolve this challenge, assuming GRH, by showing that the density of X_3 equals 1573727S/1569610, with S the so called Stephens constant. This is the first example of a `non-torsion' second order recurrent sequence with irreducible recurrence relation for which we can determine the associated density of prime divisors.


title: Band gap and stability in the ternary intermetallic compounds NiSnM (M = Ti, Zr, Hf): A first principles study

abstract: The structural stability and electronic properties of the ternary intermetallic compounds NiSnM (M = Ti, Zr, Hf) and the closely related Heusler compounds Ni$_2$SnM are discussed using the results of ab initio pseudopotential total energy and band-structure calculations performed with a plane wave basis set using the conjugate gradients algorithm. The results characterize the lowest energy phase of NiSnM compounds, with a SnM rocksalt structure sublattice, as narrow gap semiconductors with indirect gaps near 0.5 eV. Two other atomic arrangements for NiSnM in the MgAgAs structure result in energetically unfavorable compounds that are metallic. The gap formation in the lowest energy structure of NiSnZr and relative stability of the three atomic arrangements are investigated within a tight-binding framework and by considering the decompositions of each ternary compound into a binary substructure plus a third element sublattice. The stabilization of the lowest energy phase of NiSnZr is found to be mainly due to the relative stability of the SnZr rocksalt substructure, while the opening of the gap induced by the addition of the symmetry-breaking Ni sublattice makes a relatively minor contribution. From analysis of structural and chemical trends, CoVSn is predicted to be a new semiconducting intermetallic in the MgAgAs structure with a 0.8 eV indirect gap.


title: Reciprocal Space Analysis of Short Range Order Intensities by the Cluster Variation Method

abstract: A reciprocal space formulation of the Cluster Variation is used in order to extract effective pair interactions in alloys from experimental short-range order diffuse intensities. The method is applied to the analysis of the short range order contribution to the neutron diffuse scattering of a $Fe-19.5 {\%} Al$ single crystal. A detailed comparison with real space methods is carried out for three different levels of approximations. For the highest level of approximation used in this study, effective pair interactions up to fifth neighbors are obtained.


title: Frustrated H-Induced Instability of Mo(110)

abstract: Using helium atom scattering Hulpke and L"udecke recently observed a giant phonon anomaly for the hydrogen covered W(110) and Mo(110) surfaces. An explanation which is able to account for this and other experiments is still lacking. Below we present density-functional theory calculations of the atomic and electronic structure of the clean and hydrogen-covered Mo(110) surfaces. For the full adsorbate monolayer the calculations provide evidence for a strong Fermi surface nesting instability. This explains the observed anomalies and resolves the apparent inconsistencies of different experiments.


title: Three-dimensional electronic instabilities in polymerized solid A1C60

abstract: The low-temperature structure of A1C60 (A=K, Rb) is an ordered array of polymerized C60 chains, with magnetic properties that suggest a non-metallic ground state. We study the paramagnetic state of this phase using first-principles electronic-structure methods, and examine the magnetic fluctuations around this state using a model Hamiltonian. The electronic and magnetic properties of even this polymerized phase remain strongly three dimensional, and the magnetic fluctuations favor an unusual three-dimensional antiferromagnetically ordered structure with a semi-metallic electronic spectrum.


title: Competing structural instabilities in cubic perovskites

abstract: We study the antiferrodistortive instability and its interaction with ferroelectricity in cubic perovskite compounds. Our first-principles calculations show that coexistence of both instabilities is very common. We develop a first-principles scheme to study the thermodynamics of these compounds when both instabilities are present, and apply it to SrTiO$_3$. We find that increased pressure enhances the antiferrodistortive instability while suppressing the ferroelectric one. Moreover, the presence of one instability tends to suppress the other. A very rich $P$--$T$ phase diagram results.


title: Localized Basis for Effective Lattice Hamiltonians: Lattice Wannier Functions

abstract: A systematic method is presented for constructing effective Hamiltonians for general phonon-related structural transitions. The key feature is the application of group theoretical methods to identify the subspace in which the effective Hamiltonian acts and construct for it localized basis vectors, which are the analogue of electronic Wannier functions. The results of the symmetry analysis for the perovskite, rocksalt, fluorite and A15 structures and the forms of effective Hamiltonians for the ferroelectric transition in $PbTiO_3$ and $BaTiO_3$, the oxygen-octahedron rotation transition in $SrTiO_3$, the Jahn-Teller instability in $La_{1-x}(Ca,Sr,Ba)_xMnO_3$ and the antiferroelectric transition in $PbZrO_3$ are discussed. For the oxygen- octahedron rotation transition in $SrTiO_3$, this method provides an alternative to the rotational variable approach which is well behaved throughout the Brillouin zone. The parameters appearing in the Wannier basis vectors and in the effective Hamiltonian, given by the corresponding invariant energy expansion, can be obtained for individual materials using first- principles density-functional-theory total energy and linear response techniques, or any technique that can reliably calculate force constants and distortion energies. A practical approach to the determination of these parameters is presented and the application to ferroelectric $PbTiO_3$ discussed.


title: Density-Functional-Based Determination of the CH3-CH4 Hydrogen Exchange Reaction Barrier

abstract: Due to the overbinding that is inherent in existing {\em local} approximations to the density-functional formalism, certain reaction energies have not been accessible. Since the generalized gradient approximation significantly decreases the overbinding, prospects for density-functional-based reaction dynamics are promising. Results on the generalized-gradient based determination of the CH3-CH4 hydrogen exchange reaction are presented. Including all Born-Oppenheimer effects an energy barrier of 9.5 kcal/Mole is found which is a very significant improvement over the local-density approximation.


title: High-dimensional quantum dynamics of adsorption and desorption of H$_2$ at Cu(111)

abstract: We performed high-dimensional quantum dynamical calculations of the dissociative adsorption and associative desorption of hydrogen on Cu(111). The potential energy surface (PES) is obtained from density functional theory calculations. Two regimes of dynamics are found, at low energies sticking is determined by the minimum energy barrier, at high energies by the distribution of barrier heights. Experimental results are well-reproduced qualitatively, but some quantitative discrepancies are identified as well.


title: Order-N Density-Matrix Electronic-Structure Method for General Potentials

abstract: A new order-N method for calculating the electronic structure of general (non-tight-binding) potentials is presented. The method uses a combination of the ``purification''-based approaches used by Li, Nunes and Vanderbilt, and Daw, and a representation of the density matrix based on ``travelling basis orbitals''. The method is applied to several one-dimensional examples, including the free electron gas, the ``Morse'' bound-state potential, a discontinuous potential that mimics an interface, and an oscillatory potential that mimics a semiconductor. The method is found to contain Friedel oscillations, quantization of charge in bound states, and band gap formation. Quantitatively accurate agreement with exact results is found in most cases. Possible advantages with regard to treating electron-electron interactions and arbitrary boundary conditions are discussed.


title: Born Effective Charges of Barium Titanate: band by band decomposition and sensitivity to structural features

abstract: The Born effective charge tensors of Barium Titanate have been calculated for each of its 4 phases. Large effective charges of Ti and O, also predicted by shell model calculations and made plausible by a simplified model, reflect the partial covalent character of the chemical bond. A band by band decomposition confirms that orbital hybridization is not restricted to Ti and O atoms but also involves Ba which appears more covalent than generally assumed. Our calculations reveal a strong dependence of the effective charges on the atomic positions contrasting with a relative insensitivity on isotropic volume changes.


title: Commensurate and modulated magnetic phases in orthorhombic A1C60

abstract: Competing magnetically ordered structures in polymerized orthorhombic A1C60 are studied. A mean-field theory for the equilibrium phases is developed using an Ising model and a classical Heisenberg model to describe the competition between inter- and intra-chain magnetic order in the solid. In the Ising model, the limiting commensurate one-dimensional and three-dimensional phases are separated by a commensurate three-sublattice state and by two sectors containing higher-order commensurate phases. For the Heisenberg model the quasi-1D phase is never the equilibrium state; instead the 3D commensurate phases exhibits a transition to a continuum of coplanar spiral magnetic phases.


title: Real Space Approach to Electronic-Structure Calculations

abstract: We have applied the Finite Element Method to the self-consistent electronic structure calculations of molecules and solids for the first time. In this approach all the calculations are performed in "real space" and the use of non-uniform mesh is made possible, thus enabling us to deal with localized systems with ease. To illustrate the utility of this method, we perform an all-electron calculation of hydrogen molecule in a supercell with LDA approximation. Our method is also applicable to mesoscopic systems.


title: Theory of adsorption and desorption of H_2/Si(001)

abstract: While the small sticking coefficient for molecular hydrogen on the Si(001) surface apparently requires a large energy barrier of adsorption, no such barrier is observed in desorption experiments. We have calculated the potential-energy surface of an H_2 molecule in front of a Si(001) surface. If we relax the Si substrate, we find an optimum desorption path with a low (\lesssim 0.3 eV) adsorption energy barrier. While molecules impinging on the surface will mostly be reflected at the larger barrier of some frozen-substrate, molecules adsorbed on the surface can desorb along the low-barrier path.


title: Scalable Parallel Numerical Methods and Software Tools for Material Design

abstract: A new method of solution to the local spin density approximation to the electronic Schr\"{o}dinger equation is presented. The method is based on an efficient, parallel, adaptive multigrid eigenvalue solver. It is shown that adaptivity is both necessary and sufficient to accurately solve the eigenvalue problem near the singularities at the atomic centers. While preliminary, these results suggest that direct real space methods may provide a much needed method for efficiently computing the forces in complex materials.


title: Total energy global optimizations using non orthogonal localized orbitals

abstract: An energy functional for orbital based $O(N)$ calculations is proposed, which depends on a number of non orthogonal, localized orbitals larger than the number of occupied states in the system, and on a parameter, the electronic chemical potential, determining the number of electrons. We show that the minimization of the functional with respect to overlapping localized orbitals can be performed so as to attain directly the ground state energy, without being trapped at local minima. The present approach overcomes the multiple minima problem present within the original formulation of orbital based $O(N)$ methods; it therefore makes it possible to perform $O(N)$ calculations for an arbitrary system, without including any information about the system bonding properties in the construction of the input wavefunctions. Furthermore, while retaining the same computational cost as the original approach, our formulation allows one to improve the variational estimate of the ground state energy, and the energy conservation during a molecular dynamics run. Several numerical examples for surfaces, bulk systems and clusters are presented and discussed.


title: Ab-initio Calculations of Energies and Self-Diffusion on Flat and Stepped Surfaces of Al and their Implications on Crystal Growth

abstract: Using density-functional theory we investigate several properties of Al(111), Al(100), Al(110), and stepped Al(111) surfaces. We report results of formation energies of surfaces, steps, adatoms, and vacancies. For the adsorption and diffusion of Al on flat regions of Al(111) surfaces we find the hcp site energetically slightly preferred over the fcc site. The energy barrier for self-diffusion on Al(111) is very low (0.04eV). Coming close to one of the two sorts of close packed, monoatomic steps on Al(111), labeled according to their {111} and {100} micro-facets, Al adatoms experience an attraction of <~ 0.1eV already before direct contact with the edge of the step. This attraction has a range of several atomic spacings and is of electronic origin. Upon arrival at the lower step edge, the adatom attaches with no barrier at a low energy five-fold coordinated site. Coming from the upper terrace, it incorporates into the step by an atomic exchange process, which has a barrier below 0.1eV for both sorts of close packed steps. The barrier for diffusion ...


title: First-Principles Determination of Chain-Structure Instability in KNbO$_3$

abstract: A complete mapping in the Brillouin zone of the structural instability associated with the ferroelectric phase transitions of KNbO$_3$ has been obtained by first-principles calculations using an LAPW linear response approach. The wavevector dependence of the instability reveals pronounced two-dimensional character, which corresponds to chains oriented along $\left<100\right>$\ directions of displaced Nb atoms. The results are discussed in relation to models of the ferroelectric phase transitions.


title: Phase Coexistence of a Stockmayer Fluid in an Applied Field

abstract: We examine two aspects of Stockmayer fluids which consists of point dipoles that additionally interact via an attractive Lennard-Jones potential. We perform Monte Carlo simulations to examine the effect of an applied field on the liquid-gas phase coexistence and show that a magnetic fluid phase does exist in the absence of an applied field. As part of the search for the magnetic fluid phase, we perform Gibbs ensemble simulations to determine phase coexistence curves at large dipole moments, $\mu$. The critical temperature is found to depend linearly on $\mu^2$ for intermediate values of $\mu$ beyond the initial nonlinear behavior near $\mu=0$ and less than the $\mu$ where no liquid-gas phase coexistence has been found. For phase coexistence in an applied field, the critical temperatures as a function of the applied field for two different $\mu$ are mapped onto a single curve. The critical densities hardly change as a function of applied field. We also verify that in an applied field the liquid droplets within the two phase coexistence region become elongated in the direction of the field.


title: Localized Excitons and Breaking of Chemical Bonds at III-V (110) Surfaces

abstract: Electron-hole excitations in the surface bands of GaAs(110) are analyzed using constrained density-functional theory calculations. The results show that Frenkel-type autolocalized excitons are formed. The excitons induce a local surface unrelaxation which results in a strong exciton-exciton attraction and makes complexes of two or three electron-hole pairs more favorable than separate excitons. In such microscopic exciton "droplets" the electron density is mainly concentrated in the dangling orbital of a surface Ga atom whereas the holes are distributed over the bonds of this atom to its As neighbors thus weakening the bonding to the substrate. This finding suggests the microscopic mechanism of a laser-induced emission of neutral Ga atoms from GaAs and GaP (110) surfaces.


title: Theory of alkali metal adsorption on close-packed metal surfaces

abstract: Results of recent density functional theory calculations for alkali metal adsorbates on close-packed metal surfaces are discussed. Single adatoms on the (111) surface of Al and Cu are studied with the self-consistent surface Green-function method by which the pure adsorbate-substrate interaction may be analyzed. Higher coverage ordered adlayers of K on Al(111), Na on Al(111), and Na on Al(001) are treated using the ab-initio pseudopotential plane wave method which affords the prediction of coverage dependent stable and metastable adsorbate geometries and phase transitions of the adsorbate layers. Together, these studies give insight and understanding into current key issues in alkali metal adsorption, namely, the nature of the adsorbate-substrate bond at low coverage and the occurrence of hitherto unanticipated adsorbate geometries, and the associated electronic properties.


title: Electronic structure of intentionally disordered AlAs/GaAs superlattices

abstract: We use realistic pseudopotentials and a plane-wave basis to study the electronic structure of non-periodic, three-dimensional, 2000-atom (AlAs)_n/(GaAs)_m (001) superlattices, where the individual layer thicknesses n,m = {1,2,3} are randomly selected. We find that while the band gap of the equivalent (n = m = 2) ordered superlattice is indirect, random fluctuations in layer thicknesses lead to a direct gap in the planar Brillouin zone, strong wavefunction localization along the growth direction, short radiative lifetimes, and a significant band-gap reduction, in agreement with experiments on such intentionally grown disordered superlattices.


title: A Self-Consistent First-Principles Technique Having Linear Scaling

abstract: An algorithm for first-principles electronic structure calculations having a computational cost which scales linearly with the system size is presented. Our method exploits the real-space localization of the density matrix, and in this respect it is related to the technique of Li, Nunes and Vanderbilt. The density matrix is expressed in terms of localized support functions, and a matrix of variational parameters, L, having a finite spatial range. The total energy is minimized with respect to both the support functions and the elements of the L matrix. The method is variational, and becomes exact as the ranges of the support functions and the L matrix are increased. We have tested the method on crystalline silicon systems containing up to 216 atoms, and we discuss some of these results.


title: Density-Polarization Functional Theory of the response of a periodic insulating solid to an electric field.

abstract: The response of an infinite, periodic, insulating, solid to an infinitesimally small electric field is investigated in the framework of Density Functional Theory. We find that the applied perturbing potential is not a unique functional of the periodic density change~: it depends also on the change in the macroscopic {\em polarization}. Moreover, the dependence of the exchange-correlation energy on polarization induces an exchange-correlation electric field. These effects are exhibited for a model semiconductor. We also show that the scissor-operator technique is an approximate way of bypassing this polarization dependence.


title: A Classification Scheme for Toroidal Molecules

abstract: We construct a class of periodic tilings of the plane, which corresponds to toroidal arrangements of trivalent atoms, with pentagonal, hexagonal and heptagonal rings. Each tiling is characterized by a set of four integers and determines a toroidal molecule. The tiling rules are motivated by geometric considerations and the tiling patterns are rich enough to describe a wide class of toroidal carbon molecules, with a broad range of shapes and numbers of atoms. The molecular dimensions are simply related to the integers that determine the tiling. The configurational energy and the delocalisation energy of several molecules obtained in this way were computed for Tersoff and H\"uckel models. The results indicate that many of these molecules are not strained, and may be expected to be stable. We studied the influence of size on the H\"{u}ckel spectrum: it bears both similarities and differences as compared with the case of tubules.


title: First-principle study of excitonic self-trapping in diamond

abstract: We present a first-principles study of excitonic self-trapping in diamond. Our calculation provides evidence for self-trapping of the 1s core exciton and gives a coherent interpretation of recent experimental X-ray absorption and emission data. Self-trapping does not occur in the case of a single valence exciton. We predict, however, that self-trapping should occur in the case of a valence biexciton. This process is accompanied by a large local relaxation of the lattice which could be observed experimentally.


title: Frustrated H-induced Instability of Mo(110)

abstract: The large interest in the properties of transition metal surfaces has been recently fostered by inelastic helium atom scattering experiments carried out by Hulpke and Luedecke. Sharp and giant anomalies in the surface phonon dispersion curves along Gamma-H and Gamma-S have been detected on W(110) and Mo(110) at a coverage of one monolayer of hydrogen. At the same critical wave vectors a smaller second indentation is present in the experimental phonon branches. Recently we have proposed a possible interpretation which is able to explain this and other experiments. In this paper we discuss results of our recent ab initio calculations of the atomic and electronic properties of the clean and H-covered Mo(110) surface in more details. For the full monolayer coverage the calculated Fermi-surface contours are characterized by strong nesting features which originate the observed anomalies.


title: First-principles theory of ferroelectric phase transitions for perovskites: The case of BaTiO3

abstract: We carry out a completely first-principles study of the ferroelectric phase transitions in BaTiO$_3$. Our approach takes advantage of two features of these transitions: the structural changes are small, and only low-energy distortions are important. Based on these observations, we make systematically improvable approximations which enable the parameterization of the complicated energy surface. The parameters are determined from first-principles total-energy calculations using ultra-soft pseudopotentials and a preconditioned conjugate-gradient scheme. The resulting effective Hamiltonian is then solved by Monte Carlo simulation. The calculated phase sequence, transition temperatures, latent heats, and spontaneous polarizations are all in good agreement with experiment. We find the transitions to be intermediate between order-disorder and displacive character. We find all three phase transitions to be of first order. The roles of different interactions are discussed.


title: Static displacements and chemical correlations in alloys

abstract: Recent experiments in metallic solid solutions have revealed interesting correlations between static pair-displacements and the ordering behavior of these alloys. This paper discusses a simple theoretical model which successfully explains these observations and which provides a natural framework for analyzing experimental measurements of pair-displacements and chemical correlations in solid solutions. The utility and scope of this model is demonstrated by analyzing results of experiments on $Ni-Fe$ and $Cr-Fe$ alloys and results of simulations of $Cu-Au$ and $Cu-Ag$ alloys.


title: The Static, Dynamic and Electronic Properties of Liquid Gallium Studied by First-Principles Simulation

abstract: First-principles molecular dynamics simulations having a duration of 8 ps have been used to study the static, dynamic and electronic properties of l-Ga at the temperatures 702 K and 982 K. The simulations use the density-functional pseudopotential method and the system is maintained on the Born-Oppenheimer surface by conjugate gradients relaxation. The static structure factor and radial distribution function of the simulated system agree very closely with experimental data, but the diffusion coefficient is noticeably lower than measured values. The long simulations allow us to calculate the dynamical structure factor $S(q,\omega)$. A sound-wave peak is clearly visible in $S(q,\omega)$ at small wavevectors, and we present results for the dispersion curve and hence the sound velocity, which is close to the experimental value. The electronic density of states is very close to the free-electron form. Values of the electrical conductivity calculated from Kubo-Greenwood formula are in satisfactory accord with measured data.


title: Ab initio studies of structural instabilities in magnesium silicate perovskite

abstract: Density-functional simulations are used to calculate structural properties and high-symmetry phonons of the hypothetical cubic phase, the stable orthorhombic phase and an intermediate tetragonal phase of magnesium silicate perovskite. We show that the structure of the stable phase is well described by freezing in a small number of unstable phonons into the cubic phase. We use the frequencies of these unstable modes to estimate transition temperatures for cubic--tetragonal and tetragonal--orthorhombic phase transitions. These are investigated further to find that the coupling with the strain suggests that phonons give a better representation than rigid unit modes. The phonons of an intermediate tetragonal phase were found to be stable except for two rotational modes. The eigenvectors of the most unstable mode of each of the cubic and tetragonal phases account for all the positional parameters of the orthorhombic phase. The phase boundary for the orthorhombic--tetragonal transition intersects possible mantle geotherms, suggesting that the tetragonal phase may be present in the lower mantle.


title: Structure and shear response in nanometer thick films

abstract: Simulations of the structure and dynamics of fluid films confined to a thickness of a few molecular diameters are described. Confining walls introduce layering and in-plane order in the adjacent fluid. The latter is essential to transfer of shear stress. As the film thickness is decreased, by increasing pressure or decreasing the number of molecular layers, the entire film may undergo a phase transition. Spherical molecules tend to crystallize, while short chain molecules enter a glassy state with strong local orientational and translational order. These phase transitions lead to dramatic changes in the response of the film to imposed shear velocities $v$. Spherical molecules show an abrupt transition from Newtonian response to a yield stress as they crystallize. Chain molecules exhibit a continuously growing regime of non-Newtonian behavior where the shear viscosity drops as $v^{-2/3}$ at constant normal load. The same power law is found for a wide range of parameters, and extends to lower and lower velocities as a glass transition is approached. Once in the glassy state, chain molecules exhibit a finite yield stress. Shear may occur either within the film or at the film/wall interface. Interfacial shear dominates when films become glassy and when the film viscosity is increased by increasing the chain length.


title: Application of a new Tight-Binding method for transition metals: Manganese

abstract: A new tight-binding total energy method, which has been shown to accurately predict ground state properties of transition and noble metals, is applied to Manganese, the element with the most complex ground state structure among the $d$ metals. We show that the tight-binding method correctly predicts the ground state structure of Mn, and offers some insight into the magnetic properties of this state.


title: Ab initio Molecular Dynamics in Adaptive Coordinates

abstract: We present a new formulation of ab initio molecular dynamics which exploits the efficiency of plane waves in adaptive curvilinear coordinates, and thus provides an accurate treatment of first-row elements. The method is used to perform a molecular dynamics simulation of the CO_2 molecule, and allows to reproduce detailed features of its vibrational spectrum such as the splitting of the Raman sigma+_g mode caused by Fermi resonance. This new approach opens the way to highly accurate ab initio simulations of organic compounds.


title: Atomistic Aspects of Diffusion and Growth on the Si and Ge (111) Surfaces

abstract: The stability of interfaces and the mechanisms of thin film growth on semiconductors are issues of central importance in electronic devices. These issues can only be understood through detailed study of the relevant microscopic processes. Experimental studies are able to provide detailed, atomic scale information for model systems. Theoretical analysis of experimental results is essential in explaining certain surprising observations and in providing guidance for optimizing conditions and methods of growth. We review recent theoretical work on the diffusion of adatoms, the structure of adsorbate monolayers, and their implications for growth on the Si and Ge (111) surfaces. The theoretical analysis consists of first-principles calculations of the total-energy and entropy factors for stable, metastable and saddle-point configurations. These calculations are supplemented by Monte Carlo simulations of simple models that afford direct contact with experimental observations.


title: Electronic-structure-induced deformations of liquid metal clusters

abstract: Ab initio molecular dynamics is used to study deformations of sodium clusters at temperatures $500\cdots 1100$ K. Open-shell Na$_{14}$ cluster has two shape isomers, prolate and oblate, in the liquid state. The deformation is stabilized by opening a gap at the Fermi level. The closed-shell Na$_8$ remains magic also at the liquid state.


title: High Pressure Effects on Thermal Properties of MgO

abstract: Using the non-empirical Variational Induced Breathing (VIB) model, the thermal properties of periclase (MgO) under high pressures and temperatures are investigated using molecular dynamics, which includes all anharmonic effects. Equations of state for temperatures up to 3000K and pressures up to 310 GPa were calculated. Bulk modulus, thermal expansivity, Anderson-Gruneisen parameter, thermal pressure, Gruneisen parameter and their pressure and temperature dependencies are studied in order to better understand high pressure effects on thermal properties. The results agree very well with experiments and show that the thermal expansivity decreases with pressure up to about 100 GPa ($\eta$=0.73), and is almost pressure and temperature independent above this compression. It is also effected by anharmonicity at zero pressure and temperatures above 2500K. The thermal pressure changes very little with increasing pressures and temperatures, and the Gruneisen parameter is temperature independent and decreases slightly with pressure.


title: A direct pathway for sticking/desorption of H$_2$ on Si(100)

abstract: The energetics of H$_2$ interacting with the Si(100) surface is studied by means of {\em ab initio} total energy calculations within the framework of density functional theory. We find a direct desorption pathway from the mono-hydride phase which is compatible with experimental activation energies and demonstrate the importance of substrate relaxation for this process. Both the transition state configuration and barrier height depend crucially on the degree of buckling of the Si dimers on the Si(100) surface. The adsorption barrier height on the clean surface is governed by the buckling via its influence on the surface electronic structure. We discuss the consequences of this coupling for adsorption experiments and the relation between adsorption and desorption.


title: Electrostatic Energy Calculations for Molecular Dynamics

abstract: The evaluation of Coulomb forces is a difficult task. The summations that are involved converge only conditionally and care has to be taken in selecting the appropriate procedure to define the limits. The Ewald method is a standard method for obtaining Coulomb forces, but this method is rather slow, since it depends on the square of the number of atoms in a unit cell. In this paper we have adapted the plane-wise summation method for the evaluation of Coulomb forces. The use of this method allows for larger computational cells in molecular dynamics calculations.


title: Study of Oxygen Diffusion and Clustering in Silicon using an Empirical Interatomic Potential

abstract: The diffusion path and diffusivity of oxygen in crystalline silicon are computed using an empirical interatomic potential which was recently developed for modelling the interactions between oxygen and silicon atoms. The diffusion path is determined by constrained energy minimization, and the diffusivity is computed using jump rate theory. The calculated diffusivity is in excellent agreement with experiemental data. The same interatomic potntial also is used to study the formation of small clusters of oxygen atoms in silicon. The structures of these clusters are found by NPT molecular dynamics simulations, and their free energies are calculated by thermodynamic integration. These free energies are used to predict the temperature dependence of the equilibrium partitioning of oxygen into clusters of different sizes. The calculations show that, for given total oxygen concentration, most oxygen atoms are in clusters at temperature below 1300K, and that the average cluster size increases with decreasing temperature. These results are in qualitative agreement with effects of thermal annealing on oxygen precipitation in silicon crystals.


title: The Energetics and Electronic Structure of Defective and Irregular Surfaces on MgO

abstract: Ab initio calculations based on the density-functional pseudopotential approach have been used to study the fully relaxed structure, the electron distribution and the electronic density of states of (001) terraces, steps, corners and reverse corners, and of F-centers at these surface features on MgO. The calculations confirm earlier predictions of the relaxed structures of surface irregularities based on simple interaction models. A substantial narrowing of the band-gap is found at the surface, which for terraces and steps is due to surface states at the bottom of the conduction band, but for the corner and reverse corner is also due to surface states at the top of the valence band. The F-center formation energy decreases steadily as the coordination of the oxygen site is reduced. The energy of the F-center level shows a tendency to approach the top of the valence band as the coordination of its site decreases.


title: Chemical Hardness, Linear Response, and Pseudopotential Transferability

abstract: We propose a systematic method of analyzing pseudopotential transferability based on linear-response properties of the free atom, including self-consistent chemical hardness and polarizability. Our calculation of hardness extends the approach of Teter\cite{teter} not only by including self-consistency, but also by generalizing to non-diagonal hardness matrices, thereby allowing us to test for transferability to non-spherically symmetric environments. We apply the method to study the transferability of norm-conserving pseudopotentials for a variety of elements in the Periodic Table. We find that the self-consistent corrections are frequently significant, and should not be neglected. We prove that the partial-core correction improves the pseudopotential hardness of alkali metals considerably. We propose a quantity to represent the average hardness error and calculate this quantity for many representative elements as a function of pseudopotential cutoff radii. We find that the atomic polarizabilities are usually well reproduced by the norm-conserving pseudopotentials. Our results provide useful guidelines for making optimal choices in the pseudopotential generation procedure.


title: Monte-Carlo simulations of the recombination dynamics in porous silicon

abstract: A simple lattice model describing the recombination dynamics in visible light emitting porous Silicon is presented. In the model, each occupied lattice site represents a Si crystal of nanometer size. The disordered structure of porous Silicon is modeled by modified random percolation networks in two and three dimensions. Both correlated (excitons) and uncorrelated electron-hole pairs have been studied. Radiative and non-radiative processes as well as hopping between nearest neighbor occupied sites are taken into account. By means of extensive Monte-Carlo simulations, we show that the recombination dynamics in porous Silicon is due to a dispersive diffusion of excitons in a disordered arrangement of interconnected Si quantum dots. The simulated luminescence decay for the excitons shows a stretched exponential lineshape while for uncorrelated electron-hole pairs a power law decay is suggested. Our results successfully account for the recombination dynamics recently observed in the experiments. The present model is a prototype for a larger class of models describing diffusion of particles in a complex disordered system.


title: ELECTRONIC STRUCTURE OF FeSi

abstract: The full set of high-energy spectroscopy measurements including X-ray photoelectron valence band spectra and soft X-ray emission valence band spectra of both components of FeSi (Fe K_beta_5, Fe L_alpha, Si K_beta_1,3 and Si L_2,3) are performed and compared with the results of ab-initio band structure calculations using the linearized muffin-tin orbital method and linearized augmented plane wave method.


title: Fracture driven by a Thermal Gradient

abstract: Motivated by recent experiments by Yuse and Sano (Nature, 362, 329 (1993)), we propose a discrete model of linear springs for studying fracture in thin and elastically isotropic brittle films. The method enables us to draw a map of the stresses in the material. Cracks generated by the model, imposing a moving thermal gradient in the material, can branch or wiggle depending on the driving parameters. The results may be used to compare with other recent theoretical work, or to design future experiments.


title: Periodic Slab LAPW Computations for Ferroelectric BaTiO$_3$

abstract: Linearized augmented plane wave (LAPW) calculations are performed for periodic (001) and (111) slabs of BaTiO$_3$ to understand the effects of surfaces on ferroelectric BaTiO$_3$. The (111) slab is found to be much less stable than the (001) slab. The average surface energies are respectively 3700 erg/cm$^2$ and 1600 erg/cm$^2$. The depolarization field is sufficiently large in the ideal unrelaxated slab to completely inhibit the ferroelectric instability. No mid-gap surface states are evident, but there are surface states in the upper gap in the unrelaxed slab and at the top of the valence band. The dangling surface Ti bonds self-heal making the Ti-O surface very reactive and an excellent epitaxial substrate. The charge density on atoms only one unit cell away from the surface are almost identical to the bulk. Keywords: ferroelectric, surface, slab, electronic structure, depolarization, BaTiO$_3$, thin film


title: A scaling hypothesis for corrections to total energy and stress in plane- wave based ab initio calculations

abstract: We present a new technique aimed at preventing plane-wave based total energy and stress calculations from the effect of abrupt changes in basis set size. This s cheme relies on the interpolation of energy as a function of the number of plane waves, and on a scaling hypothesis that allows to perform the interpolation for a unique reference volume. From a theoretical point of view, the new method is compared to those already pr oposed in the literature, and its more rigorous derivation is emphasized. From a practical point of view, we illustrate the importance of the correction o n different materials (Si, BaTiO3, and He) corresponding to different types of b onding, and to different k-point samplings and cut-off energies. Then, we compare the different approaches for the calculation of the lattice par ameter, the bulk modulus, and its derivative versus pressure in bulk silicon.


title: Polymorphism and metastability in NbN: Structural predictions from first principles

abstract: We use ab initio pseudopotential total energy calculations with a plane wave basis set to investigate the structural energetics of various phases of polymorphic NbN. Particular attention is given to its recently discovered superconducting phase with a $T_c$ of 16.4 K, reported to have the NbO structure type. Results of total energy calculations show that it is in fact energetically unfavorable for NbN to form in this cubic structure, and its predicted theoretical lattice constant is significantly smaller than the experimental value. Various approaches to the identification of an alternative structure are discussed. In preliminary investigations, we have found two new structures that are energetically more favorable than the NbO structure.


title: Ferroelectric Phase Transitions from First Principles

abstract: An effective Hamiltonian for the ferroelectric transition in $PbTiO_3$ is constructed from first-principles density-functional-theory total-energy and linear-response calculations through the use of a localized, symmetrized basis set of ``lattice Wannier functions.'' Preliminary results of Monte Carlo simulations for this system show a first-order cubic-tetragonal transition at 660 K. The involvement of the Pb atom in the lattice instability and the coupling of local distortions to strain are found to be particularly important in producing the behavior characteristic of the $PbTiO_3$ transition. A tentative explanation for the presence of local distortions experimentally observed above $T_c$ is suggested. Further applications of this method to a variety of systems and structures are proposed for first-principles study of finite-temperature structural properties in individual materials.


title: Metastability of NbN in the ordered vacancy NbO phase

abstract: A metastable phase of NbN with superconducting Tc=16.4 K was reported recently by Treece and collaborators. The reported structure of the thin film sample deviates from the rocksalt (B1) NbN structure with 25% ordered vacancies on each sublattice (space group Pm3m) and a lattice constant of 4.442 Angstroms. Using full potential electronic structure methods, we contrast the electronic structure with that of B1 NbN. The calculated energy, 1.00 eV/molecule higher than B1 NbN, and calculated lattice constant of 4.214 Angstroms indicate that the new phase must be something other than the ordered stoichiometric Pm3m phase.


title: The Structure of the Stoichiometric and Reduced SnO2 (110) Surface

abstract: First-principles calculations based on density functional theory (DFT) and the pseudopotential method have been used to study the stoichiometric and reduced SnO2 (110) surface. The ionic relaxations are found to be moderate for both the stoichiometric and reduced surfaces, and are very similar to those found in recent DFT-pseudopotential work on TiO2. Removal of neutral oxygen leaves two electrons per oxygen on the surface, which are distributed in channels passing through bridging oxygen sites. The associated electron density can be attributed to reduction of tin from Sn4+ to Sn2+, but only if the charge distribution on Sn2+ is recognized to be highly asymmetric. Reduction of the surface gives rise to a broad distribution of gap states, in qualitative agreement with spectroscopic measurements.


title: Insulator-metal transition in solid hydrogen: Implication of electronic structure calculations for recent experiments

abstract: Electronic structure calculations for compressed molecular hydrogen are performed to provide more insight into the diversity of phenomena recently observed experimentally. We perform full-potential LAPW calculations and analyze them in terms of a molecular tight-binding model. We show that $\sigma _g$ and $\sigma _u$ bands overlap occurs at rather low pressure, while an insulating state persists at specific orientations up to 200 GPa, in accord with previous work, and is due to opening of hybridization gaps at the Fermi level. We also present calculations of electronic properties such as plasma frequencies and electron-vibron coupling.


title: Complex ion formation in liquid Ag-Se alloys

abstract: Ab initio molecular dynamics simulations are used to investigate the structure and electronic properties of the liquid Ag-Se system at three compositions. The realism of the simulations is demonstrated by comparison with diffraction data for the stoichiometric case Ag$_2$Se. As the Se content is increased beyond the stoichiometric value, short-lived Se$_n$ complexes are formed. The concentration of complexes and the associated changes of electronic structure can be explained using a simple ionic model.


title: Mechanism of Laser Induced Compaction

abstract: Fused Silica is the material of choice for UV optical systems such as the projection optics in microlithography systems. Unfortunately fused silica is not stable under UV irradition and undergoes compaction and color center formation. The UV damage properties of fused silica have been the subject of a number of studies. The detailed microscopic mechanism for color center formation is complex and remains obscure but the fact that compaction must be a generic property of a glassy system seems to have been overlooked by most researchers although Krajnovich, et. al., do briefly discuss compaction with respect to a specific glass model. The purpose of this note is to briefly discuss how the basic mechanism of compaction follows directly from the modern theoretical understanding of glassy systems. This understanding of the basic mechanism of compaction fits the experimental data quite well. It also allows for the prediction of some of qualitative properties of compaction which have not so far been experimentally determined, such as, the absence of a minimum damage threshold for the onset of compaction and the eventual long term cessation of compaction after sufficient shrinkage has occurred.


title: Molecular geometry optimization with a genetic algorithm

abstract: We present a method for reliably determining the lowest energy structure of an atomic cluster in an arbitrary model potential. The method is based on a genetic algorithm, which operates on a population of candidate structures to produce new candidates with lower energies. Our method dramatically outperforms simulated annealing, which we demonstrate by applying the genetic algorithm to a tight-binding model potential for carbon. With this potential, the algorithm efficiently finds fullerene cluster structures up to ${\rm C}_{60}$ starting from random atomic coordinates.


title: Real-Space Adaptive-Coordinate Electronic Structure Calculations

abstract: We present a real-space adaptive-coordinate method, which combines the advantages of the finite-difference approach with the accuracy and flexibility of the adaptive coordinate method. The discretized Kohn-Sham equations are written in generalized curvilinear coordinates and solved self-consistently by means of an iterative approach. The Poisson equation is solved in real space using the Multigrid algorithm. We implemented the method on a massively parallel computer, and applied it to the calculation of the equilibrium geometry and harmonic vibrational frequencies of the CO_2, CO, N_2 and F_2 molecules, yielding excellent agreement with the results of accurate quantum chemistry and Local Density Functional calculations.


title: Large Scale Electronic Structure Calculations with Multigrid Acceleration

abstract: We have developed a set of techniques for performing large scale ab initio calculations using multigrid accelerations and a real-space grid as a basis. The multigrid methods permit efficient calculations on ill-conditioned systems with long length scales or high energy cutoffs. The technique has been applied to systems containing up to 100 atoms, including a highly elongated diamond cell, an isolated C$_{60}$ molecule, and a 32-atom cell of GaN with the Ga d-states in valence. The method is well suited for implementation on both vector and massively parallel architectures.


title: Spin-polarization-induced structural selectivity in Pd$_3X$ and Pt$_3X$ ($X=3d$) compounds

abstract: Spin-polarization is known to lead to important {\it magnetic} and {\it optical} effects in open-shell atoms and elemental solids, but has rarely been implicated in controlling {\it structural} selectivity in compounds and alloys. Here we show that spin-polarized electronic structure calculations are crucial for predicting the correct $T=0$ crystal structures for Pd$_3X$ and Pt$_3X$ compounds. Spin-polarization leads to (i) stabilization of the $L1_2$ structure over the $DO_{22}$ structure in Pt$_3$Cr, Pd$_3$Cr, and Pd$_3$Mn, (ii) to the stabilization of the $DO_{22}$ structure over the $L1_2$ structure in Pd$_3$Co and to (iii) ordering (rather than phase-separation) in Pt$_3$Co and Pd$_3$Cr. The results are analyzed in terms of first-principles local spin density calculations.


title: Ab-initio simulation of high-temperature liquid selenium

abstract: Ab initio molecular dynamics simulation is used to investigate the structure and dynamics of liquid Se at temperatures of 870 and 1370~K. The calculated static structure factor is in excellent agreement with experimental data. The calculated radial distribution function gives a mean coordination number close to 2, but we find a significant fraction of one-fold and three-fold atoms, particularly at 1370~K, so that the chain structure is considerably disrupted. The self-diffusion coefficient has values ($\sim 1 \times 10^{-8}$~m~s$^{-1}$) typical of liquid metals.


title: Structural and Magnetic Properties of Trigonal Iron

abstract: First principles calculations of the electronic structure of trigonal iron were performed using density function theory. The results are used to predict lattice spacings, magnetic moments and elastic properties; these are in good agreement with experiment for both the bcc and fcc structures. We find however, that in extracting these quantities great care must be taken in interpreting numerical fits to the calculated total energies. In addition, the results for bulk iron give insight into the properties of thin iron films. Thin films grown on substrates with mismatched lattice constants often have non-cubic symmetry. If they are thicker than a few monolayers their electronic structure is similar to a bulk material with an appropriately distorted geometry, as in our trigonal calculations. We recast our bulk results in terms of an iron film grown on the (111) surface of an fcc substrate, and find the predicted strain energies and moments accurately reflect the trends for iron growth on a variety of substrates.