Optical absorption of non-interacting tight-binding electrons in a Peierls-distorted chain at half band-filling

In this first of three articles on the optical absorption of electrons in half-filled Peierls-distorted chains we present analytical results for non-interacting tight-binding electrons. We carefully derive explicit expressions for the current operator, the dipole transition matrix elements, and the optical absorption for electrons with a cosine dispersion relation of band width $W$ and dimerization parameter $\delta$. New correction (``$\eta$''-)terms to the current operator are identified. A broad band-to-band transition is found in the frequency range $W\delta < \omega < W$ whose shape is determined by the joint density of states for the upper and lower Peierls subbands and the strong momentum dependence of the transition matrix elements.Comment: 17 pages REVTEX 3.0, 2 postscript figures; hardcopy versions before May 96 are obsolete; accepted for publication in The Philosophical Magazine

Dominant particle-hole contributions to the phonon dynamics in the spinless one-dimensional Holstein model

In the spinless Holstein model at half-filling the coupling of electrons to phonons is responsible for a phase transition from a metallic state at small coupling to a Peierls distorted insulated state when the electron-phonon coupling exceeds a critical value. For the adiabatic case of small phonon frequencies, the transition is accompanied by a phonon softening at the Brillouin zone boundary whereas a hardening of the phonon mode occurs in the anti-adiabatic case. The phonon dynamics studied in this letter do not only reveal the expected renormalization of the phonon modes but also show remarkable additional contributions due to electronic particle-hole excitations.Comment: 7 pages, 4 figures and 1 table included; v2: discussion of Luttinger liquid parameters adde

Perturbation theory for optical excitations in the one-dimensional extended Peierls--Hubbard model

For the one-dimensional, extended Peierls--Hubbard model we calculate analytically the ground-state energy and the single-particle gap to second order in the Coulomb interaction for a given lattice dimerization. The comparison with numerically exact data from the Density-Matrix Renormalization Group shows that the ground-state energy is quantitatively reliable for Coulomb parameters as large as the band width. The single-particle gap can almost triple from its bare Peierls value before substantial deviations appear. For the calculation of the dominant optical excitations, we follow two approaches. In Wannier theory, we perturb the Wannier exciton states to second order. In two-step perturbation theory, similar in spirit to the GW-BSE approach, we form excitons from dressed electron-hole excitations. We find the Wannier approach to be superior to the two-step perturbation theory. For singlet excitons, Wannier theory is applicable up to Coulomb parameters as large as half band width. For triplet excitons, second-order perturbation theory quickly fails completely.Comment: 32 pages, 12 figures, submtted to JSTA

Localization Properties of the Periodic Random Anderson Model

We consider diagonal disordered one-dimensional Anderson models with an underlying periodicity. We assume the simplest periodicity, i.e., we have essentially two lattices, one that is composed of the random potentials and the other of non-random potentials. Due to the periodicity special resonance energies appear, which are related to the lattice constant of the non-random lattice. Further on two different types of behaviors are observed at the resonance energies. When a random site is surrounded by non-random sites, this model exhibits extended states at the resonance energies, whereas otherwise all states are localized with, however, an increase of the localization length at these resonance energies. We study these resonance energies and evaluate the localization length and the density of states around these energies.Comment: 4 page

Exact results for the optical absorption of strongly correlated electrons in a half-filled Peierls-distorted chain

In this second of three articles on the optical absorption of electrons in a half-filled Peierls-distorted chain we present exact results for strongly correlated tight-binding electrons. In the limit of a strong on-site interaction $U$ we map the Hubbard model onto the Harris-Lange model which can be solved exactly in one dimension in terms of spinless fermions for the charge excitations. The exact solution allows for an interpretation of the charge dynamics in terms of parallel Hubbard bands with a free-electron dispersion of band-width $W$, separated by the Hubbard interaction $U$. The spin degrees of freedom enter the expressions for the optical absorption only via a momentum dependent but static ground state expectation value. The remaining spin problem can be traced out exactly since the eigenstates of the Harris-Lange model are spin-degenerate. This corresponds to the Hubbard model at temperatures large compared to the spin exchange energy. Explicit results are given for the optical absorption in the presence of a lattice distortion $\delta$ and a nearest-neighbor interaction $V$. We find that the optical absorption for $V=0$ is dominated by a peak at $\omega=U$ and broad but weak absorption bands for $| \omega -U | \leq W$. For an appreciable nearest-neighbor interaction, $V>W/2$, almost all spectral weight is transferred to Simpson's exciton band which is eventually Peierls-split.Comment: 50 pages REVTEX 3.0, 6 postscript figures; hardcopy versions before May 96 are obsolete; accepted for publication in The Philosophical Magazine

Molecular dynamics simulations of lead clusters

Molecular dynamics simulations of nanometer-sized lead clusters have been performed using the Lim, Ong and Ercolessi glue potential (Surf. Sci. {\bf 269/270}, 1109 (1992)). The binding energies of clusters forming crystalline (fcc), decahedron and icosahedron structures are compared, showing that fcc cuboctahedra are the most energetically favoured of these polyhedral model structures. However, simulations of the freezing of liquid droplets produced a characteristic form of ``shaved'' icosahedron, in which atoms are absent at the edges and apexes of the polyhedron. This arrangement is energetically favoured for 600-4000 atom clusters. Larger clusters favour crystalline structures. Indeed, simulated freezing of a 6525-atom liquid droplet produced an imperfect fcc Wulff particle, containing a number of parallel stacking faults. The effects of temperature on the preferred structure of crystalline clusters below the melting point have been considered. The implications of these results for the interpretation of experimental data is discussed.Comment: 11 pages, 18 figues, new section added and one figure added, other minor changes for publicatio

New Tetrahedral Global Minimum for the 98-atom Lennard-Jones Cluster

A new atomic cluster structure corresponding to the global minimum of the 98-atom Lennard-Jones cluster has been found using a variant of the basin-hopping global optimization algorithm. The new structure has an unusual tetrahedral symmetry with an energy of -543.665361, which is 0.022404 lower than the previous putative global minimum. The new LJ_98 structure is of particular interest because its tetrahedral symmetry establishes it as one of only three types of exceptions to the general pattern of icosahedral structural motifs for optimal LJ microclusters. Similar to the other exceptions the global minimum is difficult to find because it is at the bottom of a narrow funnel which only becomes thermodynamically most stable at low temperature.Comment: 3 pages, 2 figures, revte

Entropic effects on the Size Evolution of Cluster Structure

We show that the vibrational entropy can play a crucial role in determining the equilibrium structure of clusters by constructing structural phase diagrams showing how the structure depends upon both size and temperature. These phase diagrams are obtained for example rare gas and metal clusters.Comment: 5 pages, 3 figure

Global Optimization by Basin-Hopping and the Lowest Energy Structures of Lennard-Jones Clusters Containing up to 110 Atoms

We describe a global optimization technique using `basin-hopping' in which the potential energy surface is transformed into a collection of interpenetrating staircases. This method has been designed to exploit the features which recent work suggests must be present in an energy landscape for efficient relaxation to the global minimum. The transformation associates any point in configuration space with the local minimum obtained by a geometry optimization started from that point, effectively removing transition state regions from the problem. However, unlike other methods based upon hypersurface deformation, this transformation does not change the global minimum. The lowest known structures are located for all Lennard-Jones clusters up to 110 atoms, including a number that have never been found before in unbiased searches.Comment: 8 pages, 3 figures, revte

Unbiased Global Optimization of Lennard-Jones Clusters for N <= 201 by Conformational Space Annealing Method

We apply the conformational space annealing (CSA) method to the Lennard-Jones clusters and find all known lowest energy configurations up to 201 atoms, without using extra information of the problem such as the structures of the known global energy minima. In addition, the robustness of the algorithm with respect to the randomness of initial conditions of the problem is demonstrated by ten successful independent runs up to 183 atoms. Our results indicate that the CSA method is a general and yet efficient global optimization algorithm applicable to many systems.Comment: revtex, 4 pages, 2 figures. Physical Review Letters, in pres