Spin coupling around a carbon atom vacancy in graphene

We investigate the details of the electronic structure in the neighborhoods of a carbon atom vacancy in graphene by employing magnetization-constrained density-functional theory on periodic slabs, and spin-exact, multi-reference, second-order perturbation theory on a finite cluster. The picture that emerges is that of two local magnetic moments (one \pi-like and one \sigma-like) decoupled from the \pi- band and coupled to each other. We find that the ground state is a triplet with a planar equilibrium geometry where an apical C atom opposes a pentagonal ring. This state lies ~0.2 eV lower in energy than the open-shell singlet with one spin flipped, which is a bistable system with two equivalent equilibrium lattice configurations (for the apical C atom above or below the lattice plane) and a barrier ~0.1 eV high separating them. Accordingly, a bare carbon-atom vacancy is predicted to be a spin-one paramagnetic species, but spin-half paramagnetism can be accommodated if binding to foreign species, ripples, coupling to a substrate, or doping are taken into account

Local-in-Time Error in Variational Quantum Dynamics

The McLachlan "minimum-distance" principle for optimizing approximate solutions of the time-dependent Schrodinger equation is revisited, with a focus on the local-in-time error accompanying the variational solutions. Simple, exact expressions are provided for this error, which are then evaluated in illustrative cases, notably the widely used mean-field approach and the adiabatic quantum molecular dynamics. Based on these findings, we demonstrate the rigorous formulation of an adaptive scheme that resizes on the fly the underlying variational manifold and, hence, optimizes the overall computational cost of a quantum dynamical simulation. Such adaptive schemes are a crucial requirement for devising and applying direct quantum dynamical methods to molecular and condensed-phase problems

Universal Markovian reduction of Brownian particle dynamics

Non-Markovian processes can often be turned Markovian by enlarging the set of variables. Here we show, by an explicit construction, how this can be done for the dynamics of a Brownian particle obeying the generalized Langevin equation. Given an arbitrary bath spectral density $J_{0}$, we introduce an orthogonal transformation of the bath variables into effective modes, leading stepwise to a semi-infinite chain with nearest-neighbor interactions. The transformation is uniquely determined by $J_{0}$ and defines a sequence $\{J_{n}\}_{n\in\mathbb{N}}$ of residual spectral densities describing the interaction of the terminal chain mode, at each step, with the remaining bath. We derive a simple, one-term recurrence relation for this sequence, and show that its limit is the quasi-Ohmic expression provided by the Rubin model of dissipation. Numerical calculations show that, irrespective of the details of $J_{0}$, convergence is fast enough to be useful in practice for an effective Markovian reduction of quantum dissipative dynamics

Few simple rules governing hydrogenation of graphene dots

We investigated binding of hydrogen atoms to small Polycyclic Aromatic Hydrocarbons (PAHs) - i.e. graphene dots with hydrogen-terminated edges - using density functional theory and correlated wavefunction techniques. We considered a number of PAHs with 3 to 7 hexagonal rings and computed binding energies for most of the symmetry unique sites, along with the minimum energy paths for significant cases. The chosen PAHs are small enough to not present radical character at their edges, yet show a clear preference for adsorption at the edge sites which can be attributed to electronic effects. We show how the results, as obtained at different level of theory, can be rationalized in detail with the help of few simple concepts derivable from a tight-binding model of the $\pi$ electrons

Electron transport in carbon wires in contact with Ag electrodes : a detailed first principles investigation

The structure and electronic properties of carbon atom chains Cn in contact with Ag electrodes are investigated in detail with first principles means. The ideal Ag(100) surface is used as a model for binding, and electron transport through the chains is studied as a function of their length, applied bias voltage, presence of capping atoms (Si, S) and adsorption site. It is found that the metal\u2013molecule bond largely influences electronic coupling to the leads. Without capping atoms the quality of the electric contact improves when increasing the carbon atom coordination number to the metal (1, 2 and 4 for adsorption on a top, bridge and hollow position, respectively) and this finding translates almost unchanged in more realistic tip-like contacts which present one, two or four metal atoms at the contact. Current\u2013voltage characteristics show Ohmic behaviour over a wide range of bias voltages and the resulting conductances change only weakly when increasing the wire length. The effect of a capping species is typically drastic, and either largely reduces (S) or largely increases (Si) the coupling of the wire to the electrodes. Comparison of our findings with recent experimental results highlights the limits of the adopted approach, which can be traced back to the known gap problem of density-functional-theory

Comment on "theoretical study of the dynamics of atomic hydrogen adsorbed on graphene multilayers"

It is shown that the theoretical prediction of a transient magnetization in bilayer and multilayer graphene (M. Moaied, Phys. Rev. B 91, 155419 (2015)PRBMDO1098-012110.1103/PhysRevB.91.155419) relies on an incorrect physical scenario for adsorption, namely, one in which H atoms adsorb barrierless on graphitic substrates and form a random adsorption pattern of monomers. Rather, according to experimental evidence, H atom sticking is an activated process, and adsorption is under kinetic control, largely ruled by a preferential sticking mechanism that leads to stable, nonmagnetic dimers at all but the smallest coverages (<0.004). Theory and experiments are reconciled by reconsidering the hydrogen atom adsorption energetics with the help of van der Waals-inclusive density functional calculations that properly account for the basis set superposition error. It is shown that today van der Waals-density functional theory predicts a shallow physisorption well that nicely agrees with available experimental data and suggests that the hydrogen atom adsorption barrier in graphene is 180 meV high, within 3c5 meV accuracy

Quantum dynamics of hydrogen atoms on graphene. I. System-bath modeling

An accurate system-bath model to investigate the quantum dynamics of hydrogen atoms chemisorbed on graphene is presented. The system comprises a hydrogen atom and the carbon atom from graphene that forms the covalent bond, and it is described by a previously developed 4D potential energy surface based on density functional theory ab initio data. The bath describes the rest of the carbon lattice and is obtained from an empirical force field through inversion of a classical equilibrium correlation function describing the hydrogen motion. By construction, model building easily accommodates improvements coming from the use of higher level electronic structure theory for the system. Further, it is well suited to a determination of the system-environment coupling by means of ab initio molecular dynamics. This paper details the system-bath modeling and shows its application to the quantum dynamics of vibrational relaxation of a chemisorbed hydrogen atom, which is here investigated at T = 0 K with the help of the multi-configuration time-dependent Hartree method. Paper II deals with the sticking dynamics

Classical and quantum dynamics at surfaces : basic concepts from simple models

Elementary processes involving atomic and molecular species at surfaces are reviewed. The emphasis is on simple classical and quantum models that help to single out unifying dynamical themes and to identify the basic physical mechanisms that underlie the rich variety of phenomena of surface chemistry. Starting from an elementary description of the energy transfer between a gas-phase species and a surface\u2014for both classical and quantum lattices\u2014the key processes establishing the formation of an adsorbed phase (sticking, diffusion and vibrational relaxation) are discussed. This is instrumental for introducing the simplest chemical transformations involving adsorbed species and/or scattering of gas-phase molecules: Langmuir\u2013Hinshelwood, Hot-Atom, and Eley\u2013Rideal reactions forming complex molecules from elementary constituents, and dissociative chemisorption of molecules into smaller fragments. Applications are also provided illustrating the ideas developed along the way at work in real-world gas-surface problems

Multi-configurational Ehrenfest simulations of ultrafast nonadiabatic dynamics in a charge-transfer complex

Multi-configurational Ehrenfest (MCE) approaches, which are intended to remedy the lack of correlations in the standard mean-field Ehrenfest method, have been proposed as coherent-state based ansatze for quantum propagation [D. V. Shalashilin, J. Chem. Phys. 130, 244101 (2009)] and also as the classical limit of the variational Gaussian-based multiconfiguration time dependent Hartree (G-MCTDH) method [S. Romer and I. Burghardt, Mol. Phys. 111, 3618 (2013)]. In the present paper, we establish the formal connection between these schemes and assess the performance of MCE for a coherent-state representation of the classical-limit subsystem. As a representative model system, we address the ultrafast, coherent charge transfer dynamics in an oligothiophene-fullerene donor acceptor complex described by a two-state linear vibronic coupling model. MCE calculations are compared with reference calculations performed with the MC IDH method, for 10-40 vibrational modes. Beyond a dimensionality of 10 modes, it is shown that the correct representation of electronic coherence depends crucially on the sampling of initially unoccupied Gaussians