Generalization of the density-matrix method to a non-orthogonal basis

We present a generalization of the Li, Nunes and Vanderbilt density-matrix method to the case of a non-orthogonal set of basis functions. A representation of the real-space density matrix is chosen in such a way that only the overlap matrix, and not its inverse, appears in the energy functional. The generalized energy functional is shown to be variational with respect to the elements of the density matrix, which typically remains well localized.Comment: 11 pages + 2 postcript figures at the end (search for -cut here

Effect of dimensionality on the charge-density-wave in few-layers 2H-NbSe2_2

We investigate the charge density wave (CDW) instability in single and double layers, as well as in the bulk 2H-NbSe$_{2}$. We demonstrate that the density functional theory correctly describes the metallic CDW state in the bulk 2H-NbSe$_{2}$. We predict that both mono- and bilayer NbSe$_{2}$ undergo a CDW instability. However, while in the bulk the instability occurs at a momentum $\mathbf{q}_{CDW}\approx{2/3}\mathbf{\Gamma M}$, in free-standing layers it occurs at $\mathbf{q}_{CDW}\approx{1/2}\mathbf{\Gamma M}$. Furthermore, while in the bulk the CDW leads to a metallic state, in a monolayer the ground state becomes semimetallic, in agreement with recent experimental data. We elucidate the key role that an enhancement of the electron-phonon matrix element at $\mathbf{q}\approx\mathbf{q}_{CDW}$ plays in forming the CDW ground state.Comment: 4 pages 5 figure

All-electron magnetic response with pseudopotentials: NMR chemical shifts

A theory for the ab initio calculation of all-electron NMR chemical shifts in insulators using pseudopotentials is presented. It is formulated for both finite and infinitely periodic systems and is based on an extension to the Projector Augmented Wave approach of Bloechl [P. E. Bloechl, Phys. Rev. B 50, 17953 (1994)] and the method of Mauri et al [F. Mauri, B.G. Pfrommer, and S.G. Louie, Phys. Rev. Lett. 77, 5300 (1996)]. The theory is successfully validated for molecules by comparison with a selection of quantum chemical results, and in periodic systems by comparison with plane-wave all-electron results for diamond.Comment: 25 pages, 4 tables, submitted to Physical Review

Theory of double-resonant Raman spectra in graphene: intensity and line shape of defect-induced and two-phonon bands

We calculate the double resonant (DR) Raman spectrum of graphene, and determine the lines associated to both phonon-defect processes, and two-phonons ones. Phonon and electronic dispersions reproduce calculations based on density functional theory corrected with GW. Electron-light, -phonon, and -defect scattering matrix elements and the electronic linewidth are explicitly calculated. Defect-induced processes are simulated by considering different kind of idealized defects. For an excitation energy of $\epsilon_L=2.4$ eV, the agreement with measurements is very good and calculations reproduce: the relative intensities among phonon-defect or among two-phonon lines; the measured small widths of the D, $D'$, 2D and $2D'$ lines; the line shapes; the presence of small intensity lines in the 1800, 2000 cm$^{-1}$ range. We determine how the spectra depend on the excitation energy, on the light polarization, on the electronic linewidth, on the kind of defects and on their concentration. According to the present findings, the intensity ratio between the $2D'$ and 2D lines can be used to determine experimentally the electronic linewidth. The intensity ratio between the $D$ and $D'$ lines depends on the kind of model defect, suggesting that this ratio could possibly be used to identify the kind of defects present in actual samples. Charged impurities outside the graphene plane provide an almost undetectable contribution to the Raman signal

Hygrothermal damage mechanisms in graphite-epoxy composites

T300/5209 and T300/5208 graphite epoxy laminates were studied experimentally and analytically in order to: (1) determine the coupling between applied stress, internal residual stress, and moisture sorption kinetics; (2) examine the microscopic damage mechanisms due to hygrothermal cycling; (3) evaluate the effect of absorbed moisture and hygrothermal cycling on inplane shear response; (4) determine the permanent loss of interfacial bond strength after moisture absorption and drying; and (5) evaluate the three dimensional stress state in laminates under a combination of hygroscopic, thermal, and mechanical loads. Specimens were conditioned to equilibrium moisture content under steady exposure to 55% or 95% RH at 70 C or 93 C. Some specimens were tested subsequent to moisture conditioning and 100 cycles between -54 C and either 70 C or 93 C

Acceleration Schemes for Ab-Initio Molecular Dynamics and Electronic Structure Calculations

We study the convergence and the stability of fictitious dynamical methods for electrons. First, we show that a particular damped second-order dynamics has a much faster rate of convergence to the ground-state than first-order steepest descent algorithms while retaining their numerical cost per time step. Our damped dynamics has efficiency comparable to that of conjugate gradient methods in typical electronic minimization problems. Then, we analyse the factors that limit the size of the integration time step in approaches based on plane-wave expansions. The maximum allowed time step is dictated by the highest frequency components of the fictitious electronic dynamics. These can result either from the large wavevector components of the kinetic energy or from the small wavevector components of the Coulomb potential giving rise to the so called {\it charge sloshing} problem. We show how to eliminate large wavevector instabilities by adopting a preconditioning scheme that is implemented here for the first-time in the context of Car-Parrinello ab-initio molecular dynamics simulations of the ionic motion. We also show how to solve the charge-sloshing problem when this is present. We substantiate our theoretical analysis with numerical tests on a number of different silicon and carbon systems having both insulating and metallic character.Comment: RevTex, 9 figures available upon request, to appear in Phys. Rev.

Branching and leaf initiation in the erect aerial system of Stromatopteris moniliformis (Gleicheniaceae)

La partie aérienne dressée du #Stromatopteris moniliformis$ Met., une fougère Gleicheniacée endémique de la Nouvelle-Calédonie, a fait l'objet d'une étude histomorphogénique. Partant d'un système souterrain ramifié (non étudié), des axes aériens dressés traversent la litière. Ils portent des bourgeons latéraux dormants dans les parties profondes, puis plusieurs feuilles au sommet, avant que le méristème terminal ne dégénère. Les bourgeons latéraux se développent à leur tout successivement, dans le sens basipète, pour atteindre la surface de la litière; ils portent aussi des bourgeons latéraux et quelques feuilles de la même façon que l'axe principal, etc. On aboutit ainsi à un buissonnement très dense d'axes ramifiés prolongés par des feuilles, port original dans ce groupe. (Résumé d'auteur