Search

We investigate whether making the friction spatially dependent on the reaction coordinate introduces quantum effects into the thermal reaction rates for dissipative reactions. Quantum rates are calculated using the numerically exact multi-configuration time-dependent Hartree (MCTDH) method, as well as the approximate ring-polymer molecular dynamics (RPMD), ring-polymer instanton (RPI) methods, and classical mechanics. By conducting simulations across a wide range of temperatures and friction strengths, we can identify the various regimes that govern the reactive dynamics. At high temperatures, in addition to the spatial-diffusion and energy-diffusion regimes predicted by Kramer's rate theory, a (coherent) tunnelling-dominated regime is identified at low friction. At low temperatures, incoherent tunnelling dominates most of Kramer's curve, except at very low friction when coherent tunnelling becomes dominant. Unlike in classical mechanics, the bath's influence changes the equilibrium time-independent properties of the system, leading to a complex interplay between spatially dependent friction and nuclear quantum effects even at high temperatures. More specifically, a realistic friction profile can lead to an increase (decrease) of the quantum (classical) rates with friction within the spatial-diffusion regime, showing that classical and quantum rates display qualitatively different behaviours. Except at very low frictions, we find that RPMD captures most of the quantum effects in the thermal reaction rates.

Geometric phases play a crucial role in diverse fields. In chemistry they appear when a reaction path encircles an intersection between adiabatic potential energy surfaces and the molecular wavefunction experiences quantum-mechanical interference effects. This intriguing effect, closely resembling the magnetic Aharonov-Bohm effect, crucially relies on the adiabatic description of the dynamics, and it is uncertain whether and how it persists in an exact quantum dynamical framework. Recent works have shown that the geometric phase is an artifact of the adiabatic approximation, thereby challenging the perceived utility of the geometric phase concept in molecules. Here, we investigate this issue in an exact dynamical framework. We introduce instantaneous, gauge invariant phases separately for the electrons and for the nuclei, and use them to monitor the phase difference between the trailing edges of a wavepacket encircling a conical intersection. In this way we unambiguosly assess the role of the geometric phase in the interference process and shed light on its persistence in molecular systems., Comment: arXiv admin note: text overlap with arXiv:2312.02823

The fate of the molecular geometric phase in an exact dynamical framework is investigated with the help of the exact factorization of the wavefunction and a recently proposed quantum hydrodynamical description of its dynamics. An instantaneous, gauge invariant phase is introduced for arbitrary paths in nuclear configuration space in terms of hydrodynamical variables, and shown to reduce to the adiabatic geometric phase when the state is adiabatic and the path is closed. The evolution of the closed-path phase over time is shown to adhere to a Maxwell-Faraday induction law, with non-conservative forces arising from the electron dynamics that play the role of electromotive forces. We identify the pivotal forces that are able to change the value of the phase, thereby challenging any topological argument. Nonetheless, negligible changes in the phase occur when the local dynamics along the probe loop is approximately adiabatic. In other words, the adiabatic idealization of geometric phase effects may remain suitable for effectively describing certain dynamic observables.

The quantum dynamics of electron-nuclear systems is analyzed from the perspective of the exact factorization of the wavefunction, with the aim of defining gauge invariant equations of motion for both the nuclei and the electrons. For pure states this is accomplished with a quantum hydrodynamical description of the nuclear dynamics and electronic density operators tied to the fluid elements. For statistical mixtures of states the exact factorization approach is extended to two limiting situations that we call "type-n" and "type-e" mixtures, depending on whether the nuclei or the electrons are, respectively, in an intrinsically mixed state. In both cases a fully gauge invariant formulation of the dynamics is obtained again in hydrodynamic form with the help of mechanical momentum moments (MMMs). Nuclear MMMs extend in a gauge invariant way the ordinary momentum moments of the Wigner distribution associated with a density matrix of positional variables, electron MMMs are operator-valued and represent a generalization of the (conditional) density operators used for pure states. The theory presented here bridges exact quantum dynamics with several mixed quantum-classical approaches currently in use to tackle non-adiabatic molecular problems, offering a foundation for systematic improvements. It further connects to non-adiabatic theories in condensed-phase systems. As an example, we re-derive the finite-temperature theory of electronic friction of Dou, Miao \& Subotnik (Phys. Rev. Lett. 119, 046001 (2017)) from the dynamics of "type-e" mixtures and discuss possible improvements.

Carbon atom vacancies in graphene give rise to a local magnetic moment of $\sigma+\pi$ origin, whose magnitude is yet uncertain and debated. Partial quenching of $\pi$ magnetism has been ubiquitously reported in periodic $first-principles$ calculations, with magnetic moments scattered in the range 1.0 - 2.0 $\mu_{B}$, slowly converging to the lower or the upper end, depending on how the diluted limit is approached. By contrast, (ensemble) density functional theory calculations on cluster models neatly converge to the value of $2$ $\mu_{B}$ when increasing the system size. This stunning discrepancy has sparked an ongoing debate about the role of defect-defect interactions and self-doping, and about the importance of the self-interaction-error in the density-functional-theory description of the vacancy-induced states. Here, we settle this puzzle by showing that the problem has a fundamental, mono-electronic origin which is related to the special (periodic) arrangement of defects that results when using the slab-supercell approach. Specifically, we report the existence of resonant states that are $anomalously$ delocalized over the lattice and that make the $\pi$ midgap band $unphysically$ dispersive, hence prone to self-doping and quenching of the $\pi$ magnetism. Hybrid functionals fix the problem by widening the gap between the spin-resolved $\pi$ midgap bands, without reducing their artificial widths. As a consequence, while reconciling the magnetic moment with expectations, they predict a spin-splitting which is one order of magnitude larger than found in experiments.

The simulation of quantum dynamics calls for quantum algorithms working in first quantized grid encodings. Here, we propose a variational quantum algorithm for performing quantum dynamics in first quantization. In addition to the usual reduction in circuit depth conferred by variational approaches, this algorithm also enjoys several advantages compared to previously proposed ones. For instance, variational approaches suffer from the need for a large number of measurements. However, the grid encoding of first quantized Hamiltonians only requires measuring in position and momentum bases, irrespective of the system size. Their combination with variational approaches is therefore particularly attractive. Moreover, heuristic variational forms can be employed to overcome the limitation of the hard decomposition of Trotterized first quantized Hamiltonians into quantum gates. We apply this quantum algorithm to the dynamics of several systems in one and two dimensions. Our simulations exhibit the previously observed numerical instabilities of variational time propagation approaches. We show how they can be significantly attenuated through subspace diagonalization at a cost of an additional $\mathcal{O}(MN^2)$ 2-qubit gates where $M$ is the number of dimensions and $N^M$ is the total number of grid points.

A theory of electronic friction is developed using the exact factorization of the electron-nuclear wavefunction. No assumption is made regarding the electronic bath, which can be made of independent or interacting electrons, and the nuclei are treated quantally. The ensuing equation of motion for the nuclear wavefunction is a non-linear Schr\"{o}dinger equation including a friction term. The resulting friction kernel agrees with a previously derived mixed quantum-classical result by Dou, Miao \& Subotnik (\emph{Phys. Rev. Lett.} \textbf{119}, 046001 (2017)), except for a \emph{pseudo}-magnetic contribution in the latter that is here removed. More specifically, it is shown that the electron dynamics generally washes out the\emph{ gauge} fields appearing in the adiabatic dynamics. However, at T=0 K, the \emph{pseudo}-magnetic force is fully re-established in the typical situation where the electrons respond rapidy on the slow time-scale of the nuclear dynamics (Markov limit). Hence, we predict Berry's phase effects to be observable also in the presence of electronic friction, and non-trivial geometric phases should be attainable for molecules on metallic magnetic surfaces., Comment: 20 pages (with Supplemental Material)

The purpose of the present Comment is to point out the connection between an approach to spawning and regularization that was recently introduced by D. Mendive-Tapia and H.-D. Meyer [J. Chem. Phys. 153, 234114 (2020)] in the context of the Multiconfiguration Time-Dependent Hartree (MCTDH) method, and earlier work where adaptive variational quantum propagation based on the Local-in-Time Error (LITE) was introduced [R. Martinazzo and I. Burghardt, Phys. Rev. Lett. 124, 150601 (2020); arXiv:1907.00841 [quant-ph] (2019)]. Furthermore, we show that the LITE represents a gauge-invariant distance which provides a natural, physically sound tool for adaptive 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. These findings pave the way for the rigorous development of adaptive schemes that re-size on-the-fly the underlying variational manifold and thus optimize the overall computational cost of a quantum dynamical simulation.

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.

Graphene, being one-atom thick, is extremely sensitive to the presence of adsorbed atoms and molecules and, more generally, to defects such as vacancies, holes and/or substitutional dopants. This property, apart from being directly usable in molecular sensor devices, can also be employed to tune graphene electronic properties. Here we briefly review the basic features of atomic-scale defects that can be useful for material design. After a brief introduction on isolated $p_z$ defects, we analyse the electronic structure of multiple defective graphene substrates, and show how to predict the presence of microscopically ordered magnetic structures. Subsequently, we analyse the more complicated situation where the electronic structure, as modified by the presence of some defects, affects chemical reactivity of the substrate towards adsorption (chemisorption) of atomic/molecular species, leading to preferential sticking on specific lattice positions. Then, we consider the reverse problem, that is how to use defects to engineer graphene electronic properties. In this context, we show that arranging defects to form honeycomb-shaped superlattices (what we may call "supergraphenes") a sizeable gap opens in the band structure and new Dirac cones are created right close to the gapped region. Similarly, we show that substitutional dopants such as group IIIA/VA elements may have gapped quasi-conical structures corresponding to massive Dirac carriers. All these possible structures might find important technological applications in the development of graphene-based logic transistors., Comment: 16 pages, 14 figures, "Physics and Applications of Graphene - Theory" - Chapter 3, http://www.intechweb.org/books/show/title/physics-and-applications-of-graphene-theory

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.

We investigate graphene superlattices of nitrogen and boron substitutional defects and by using symmetry arguments and electronic structure calculations we show how such superlattices can be used to modify graphene band structure. Specifically, depending on the superlattice symmetry, the structures considered here can either preserve the Dirac cones (D_{6h} superlattices) or open a band gap (D_{3h}). Relevant band parameters (carriers effective masses, group velocities and gaps, when present) are found to depend on the superlattice constant n as 1/n^{p} where p is in the range 1-2, depending on the case considered. Overall, the results presented here show how one can tune the graphene band structure to a great extent by modifying few superlattice parameters., Comment: accepted, J. Phys. Chem. C

We study nxn honeycomb superlattices of defects in graphene. The considered defects are missing p_z orbitals and can be realized by either introducing C atom vacancies or chemically binding simple atomic species at the given sites. Using symmetry arguments we show how it is possible to open a gap when n=3m+1,3m+2 (m integer), and estimate its value to have an approximate square-root dependence on the defect concentration x=1/n^2. Tight-binding calculations confirm these findings and show that the induced-gaps can be quite large, e.g. ~100 meV for x~10^{-3}. Gradient-corrected density functional theory calculations on a number of superlattices made by H atoms adsorbed on graphene are in good agreement with tight-binding results, thereby suggesting that the proposed structures may be used in practice to open a gap in graphene., Comment: 5 pages, 4 figures

Adsorption of hydrogen atoms on a single graphite sheet (graphene) has been investigated by first-principles electronic structure means, employing plane-wave based, periodic density functional theory. A reasonably large 5x5 surface unit cell has been employed to study single and multiple adsorption of H atoms. Binding and barrier energies for sequential sticking have been computed for a number of configurations involving adsorption on top of carbon atoms. We find that binding energies per atom range from ~0.8 eV to ~1.9 eV, with barriers to sticking in the range 0.0-0.2 eV. In addition, depending on the number and location of adsorbed hydrogen atoms, we find that magnetic structures may form in which spin density localizes on a $\sqrt{3}{x}\sqrt{3}{R}30^{\circ}$ sublattice, and that binding (barrier) energies for sequential adsorption increase (decrease) linearly with the site-integrated magnetization. These results can be rationalized with the help of the valence-bond resonance theory of planar $\pi$ conjugated systems, and suggest that preferential sticking due to barrierless adsorption is limited to formation of hydrogen pairs., Comment: 12 pages, 8 figures and 4 tables