We study interaction-driven ultrafast dynamics in finite graphene flakes following an optical pump quench in an interacting tight-binding model. By comparing exact real-time evolution with simulations restricted to particle-hole excitation subspaces, we assess when relaxation can be captured by low-order many-body processes and when this is not sufficient. The single-particle orbital entropy provides a compact diagnostic for dynamic correlation growth. For the systems studied here, periodic graphene flakes are well described by low-order excitations, whereas confined geometries require substantial higher-order contributions even for relatively small interaction strengths. The quench protocol combines simple initial-state preparation with strongly correlated dynamics, identifying a promising benchmark problem for future quantum-computing simulations.

Quantum computers are increasingly accessible, yet demonstrations of physically meaningful simulations for real materials remain scarce. In our work we simulate low-energy magnetic excitations, specifically spin-wave spectra, of chromium tri-halide monolayers. Starting from ab-initio electronic structure calculations for these two-dimensional magnets, we derive an effective spin model and simulate low-energy spin excitations using a real-time propagation of the spin system on the commercial quantum computing cloud platform IQM Resonance. The results for systems with up to 48 qubits are validated against classical benchmarks. While some spectral features remain challenging for today's NISQ devices, our simulation achieves good agreement at quasi-constant wall-time scaling, compared to the exponential scaling of classical methods. Our results demonstrate that, even in the absence of quantum advantage, useful quantum simulations of real materials are becoming possible for domain experts via commercial cloud access to quantum computers.

Magnetism and spin physics are true quantum mechanical effects and their description usually requires multi reference methods and is often hidden in the standard description of molecules in quantum chemistry. In this work we present a twofold approach to the description of spin physics in molecules and solids. First, we present a method that identifies the single-particle basis in which a given subset of the orbitals is equivalent to spin degrees of freedom for models and materials which feature significant spin physics at low energies. We introduce a metric for the spin-like character of a basis orbital, of which the optimization yields the basis containing the optimum spin-like basis orbitals. Second, we demonstrate an extended Schrieffer-Wolff transformation method to derive the effective Hamiltonian acting on the subspace of the Hilbert space in which the charge degree of freedom of electron densities in the spin-like orbitals is integrated out. The method then yields an effective Hamiltonian describing spins coupled to a fermionic environment. This extended Schrieffer-Wolff transformation is applicable to a wide range of Hamiltonians and has been utilized in this work for model Hamiltonians as well as the Hamiltonian describing the active orbital space of molecular chromium bromide. This is achieved by reformulating the highly non-linear Schrieffer-Wolff equations into a linear set of equations corresponding to an operator basis.

Magnetism and spin physics are true quantum mechanical effects and their description usually requires multi-reference methods and is often hidden in the standard description of molecules in quantum chemistry. In this work, we present a twofold approach to the description of spin physics in molecules and solids. First, we present a method that identifies the single-particle basis in which a given subset of the orbitals is equivalent to spin degrees of freedom for models and materials, which feature significant spin physics at low energies. We introduce a metric for the spin-like character of a basis orbital, of which the optimization yields the basis containing the optimum spin-like basis orbitals. Second, we demonstrate an extended Schrieffer–Wolff transformation method to derive the effective Hamiltonian acting on the subspace of the Hilbert space in which the charge degree of freedom of electron densities in the spin-like orbitals is integrated out. The method then yields an effective Hamiltonian describing spins coupled to a fermionic environment. This extended Schrieffer–Wolff transformation is applicable to a wide range of Hamiltonians and has been utilized in this work for model Hamiltonians as well as the Hamiltonian describing the active orbital space of molecular chromium bromide. This is achieved by reformulating the highly non-linear Schrieffer–Wolff equations into a linear set of equations corresponding to an operator basis.

We propose and analyze a method for improving quantum chemical energy calculations on a quantum computer impaired by decoherence and shot noise. The error mitigation approach relies on the fact that the one- and two-particle reduced density matrices (1- and 2-RDM) of a chemical system need to obey so-called N-representability constraints. We post-process the result of an RDM measurement by projecting it into the subspace where certain N-representability conditions are fulfilled. Furthermore, we utilize that such constraints also hold in the hole and particle-hole sector and perform projections in these sectors as well. We expand earlier work by conducting a careful analysis of the method's performance in the context of quantum computing. Specifically, we consider typical decoherence channels (dephasing, damping, and depolarizing noise) as well as shot noise due to a finite number of projective measurements. We provide analytical considerations and examine numerically three example systems, \ch{H2}, \ch{LiH}, and \ch{BeH2}. From these investigations, we derive our own practical yet effective method to best employ the various projection options. Our results show the approach to significantly lower energy errors and measurement variances of (simulated) quantum computations.

In this work, we performed extensive first-principles simulations of high-harmonic generation in the topological Diract semimetal Na 3 Bi using a first-principles time-dependent density functional theory framework, focusing on the effect of spin-orbit coupling (SOC) on the harmonic response. We also derived an analytical model describing the microscopic mechanism of strong-field dynamics in presence of spin-orbit coupling, starting from a locally U (1) × S U (2) gauge-invariant Hamiltonian. Our results reveal that SOC: (i) affects the strong-field excitation of carriers to the conduction bands by modifying the bandstructure of Na 3 Bi, (ii) makes each spin channel reacts differently to the driven laser by modifying the electron velocity (iii) changes the emission timing of the emitted harmonics. Moreover, we show that the SOC affects the harmonic emission by directly coupling the charge current to the spin currents, paving the way to the high-harmonic spectroscopy of spin currents in solids.

Variational algorithms are promising candidates to be implemented on near-term quantum computers. The variational quantum eigensolver (VQE) is a prominent example, where a parametrized trial state of the quantum mechanical wave function is optimized to obtain the ground state energy. In our work, we investigate the variational Hamiltonian Ansatz (VHA), where the trial state is given by a non-interacting reference state modified by unitary rotations using generators that are part of the Hamiltonian describing the system. The lowest energy is obtained by optimizing the angles of those unitary rotations. A standard procedure to optimize the variational parameters is to use gradient-based algorithms. However, shot noise and the intrinsic noise of the quantum device affect the evaluation of the required gradients. We studied how different methods for obtaining the gradient, specifically the finite-difference and the parameter-shift rule, are affected by shot noise and noise of the quantum computer. To this end, we simulated a simple quantum circuit, as well as the 2-site and 6-site Hubbard model.

Cavity modification of material properties and phenomena is a novel research field largely motivated by the advances in strong light-matter interactions. Despite this progress, exact solutions for extended systems strongly coupled to the photon field are not available, and both theory and experiments rely mainly on finite-system models. Therefore, a paradigmatic example of an exactly solvable extended system in a cavity becomes highly desirable. To fill this gap we revisit Sommerfeld's theory of the free electron gas in cavity quantum electrodynamics. We solve this system analytically in the long-wavelength limit for an arbitrary number of noninteracting electrons, and we demonstrate that the electron-photon ground state is a Fermi liquid which contains virtual photons. In contrast to models of finite systems, no ground state exists if the diamagentic A2 term is omitted. Further, by performing linear response we show that the cavity field induces plasmon-polariton excitations and modifies the optical and the DC conductivity of the electron gas. Our exact solution allows us to consider the thermodynamic limit for both electrons and photons by constructing an effective quantum field theory. The continuum of modes leads to a many-body renormalization of the electron mass, which modifies the fermionic quasiparticle excitations of the Fermi liquid and the Wigner-Seitz radius of the interacting electron gas. Last, we show how the matter-modified photon field leads to a repulsive Casimir force and how the continuum of modes introduces dissipation into the light-matter system. Several of the presented findings should be experimentally accessible.

A formalism to describe steady-state electronic and thermal transport in the framework of density functional theory is presented. A one-to-one correspondence is proven between the three basic variables of the theory, i.e., the density on as well as the electrical and heat currents through the junction, and the three basic potentials, i.e., the local potential in as well as the DC bias and thermal gradient across the junction. Consequently, the Kohn-Sham system of the theory requires three exchange-correlations potentials. In linear response, the formalism leads to exact expressions for the many-body transport coefficients (both electrical and thermal conductances and Seebeck coefficient) in terms of both the corresponding Kohn-Sham coefficients and derivatives of the exchange-correlations potentials. The theory is applied to the single impurity Anderson model and an accurate analytic parametrization for these derivatives in the Coulomb blockade regime is constructed through reverse engineering.

We develop a trajectory-based approach for excited-state molecular dynamics simulations of systems subject to an external periodic drive. We combine the exact-factorization formalism, allowing us to treat electron-nuclear systems in nonadiabatic regimes, with the Floquet formalism for time-periodic processes. The theory is developed starting with the molecular time-dependent Schrödinger equation with the inclusion of an external periodic drive that couples to the system dipole moment. With the support of the Floquet formalism, quantum dynamics is approximated by combining classical-like, trajectory-based, nuclear evolution with electronic dynamics represented in the Floquet basis. The resulting algorithm, which is an extension of the coupled-trajectory mixed quantum-classical scheme for periodically driven systems, is applied to a model study, exactly solvable, with different field intensities.

We investigate non-adiabatic effects on the magnon frequency in an interacting system of localized spins and itinerant electrons. Including the lowest order corrections to the adiabatic dynamics in an analytically solvable model, applicable to simple ferromagnets like Fe, Co and Ni, we find that the magnon frequency is renormalized by a geometrical torque arising from the electronic spin Berry curvature. Comparison to exact numerical simulations reveals that our analytical solution captures essential low-energy features, and provides a mechanism for the magnon frequency hardening observed in recent first principles calculations for Fe, provided the geometrical torque is taken into account.

We propose an efficient and non-perturbative scheme to compute magnetic excitations for extended systems employing the framework of time-dependent density functional theory. Within our approach, we drive the system out of equilibrium using an ultrashort magnetic kick perpendicular to the ground-state magnetization of the material. The dynamical properties of the system are obtained by propagating the time-dependent Kohn-Sham equations in real time, and the analysis of the time-dependent magnetization reveals the transverse magnetic excitation spectrum of the magnet. We illustrate the performance of the method by computing the magnetization dynamics, obtained from a real-time propagation, for iron, cobalt, and nickel and compare them to known results obtained using the linear-response formulation of time-dependent density functional theory. Moreover, we point out that our time-dependent approach is not limited to the linear-response regime, and we present the first results for nonlinear magnetic excitations from first principles in iron.

Over the last few years, extraordinary advances in experimental and theoretical tools have allowed us to monitor and control matter at short time and atomic scales with a high degree of precision. An appealing and challenging route toward engineering materials with tailored properties is to find ways to design or selectively manipulate materials, especially at the quantum level. To this end, having a state-of-the-art ab initio computer simulation tool that enables a reliable and accurate simulation of light-induced changes in the physical and chemical properties of complex systems is of utmost importance. The first principles real-space-based Octopus project was born with that idea in mind, i.e., to provide a unique framework that allows us to describe non-equilibrium phenomena in molecular complexes, low dimensional materials, and extended systems by accounting for electronic, ionic, and photon quantum mechanical effects within a generalized time-dependent density functional theory. This article aims to present the new features that have been implemented over the last few years, including technical developments related to performance and massive parallelism. We also describe the major theoretical developments to address ultrafast light-driven processes, such as the new theoretical framework of quantum electrodynamics density-functional formalism for the description of novel light-matter hybrid states. Those advances, and others being released soon as part of the Octopus package, will allow the scientific community to simulate and characterize spatial and time-resolved spectroscopies, ultrafast phenomena in molecules and materials, and new emergent states of matter (quantum electrodynamical-materials).

The search for exchange-correlation functionals going beyond the adiabatic approximation has always been a challenging task for time-dependent density-functional theory. Starting from known results and using symmetry properties, we put forward a nonadiabatic exchange-correlation functional for lattice models describing a generic transport setup. We show that this functional reduces to known results for a single quantum dot connected to one or two reservoirs and furthermore yields the adiabatic local-density approximation in the static limit. Finally, we analyze the features of the exchange-correlation potential and the physics it describes in a linear chain connected to two reservoirs where the transport is induced by a bias voltage applied to the reservoirs. We find that the Coulomb blockade is correctly described for a half-filled chain, while additional effects arise as the doping of the chain changes.

We present a kinetic-energy density-functional theory and the corresponding kinetic-energy Kohn-Sham (keKS) scheme on a lattice and show that, by including more observables explicitly in a density-functional approach, already simple approximation strategies lead to very accurate results. Here, we promote the kinetic-energy density to a fundamental variable alongside the density and show for specific cases (analytically and numerically) that there is a one-to-one correspondence between the external pair of on-site potential and site-dependent hopping and the internal pair of density and kinetic-energy density. On the basis of this mapping, we establish two unknown effective fields, the mean-field exchange-correlation potential and the mean-field exchange-correlation hopping, which force the keKS system to generate the same kinetic-energy density and density as the fully interacting one. We show, by a decomposition based on the equations of motions for the density and the kinetic-energy density, that we can construct simple orbital-dependent functionals that outperform the corresponding exact-exchange Kohn-Sham (KS) approximation of standard density-functional theory. We do so by considering the exact KS and keKS systems and comparing the unknown correlation contributions as well as by comparing self-consistent calculations based on the mean-field exchange (for the effective potential) and a uniform (for the effective hopping) approximation for the keKS and the exact-exchange approximation for the KS system, respectively.

The wide-band limit is a commonly used approximation to analyze transport through nanoscale devices. In this work we investigate its applicability to the study of charge and heat transport through molecular break junctions exposed to voltage biases and temperature gradients. We find by comparative simulations that while the wide-band-limit approximation faithfully describes the long-time charge and heat transport, it fails to characterize the short-time behavior of the junction. In particular, we show that the charge current flowing through the device shows a discontinuity when a temperature gradient is applied, while the energy flow is discontinuous when a voltage bias is switched on and even diverges when the junction is exposed to both a temperature gradient and a voltage bias. We provide an explanation for this pathological behavior and propose two possible solutions to this problem.

Ab initio calculations of the magnon dispersion in ferromagnetic materials typically rely on the adiabatic local density approximation (ALDA) in which the effective exchange-correlation field is everywhere parallel to the magnetization. These calculations, however, tend to overestimate the “magnon stiffness”, defined as the curvature of the magnon frequency vs. wave vector relation evaluated at zero wave vector. Here we suggest a simple procedure to improve the magnon dispersion by taking into account gradient corrections to the ALDA at the exchange-only level. We find that this gradient correction always reduces the magnon stiffness. The surprisingly large size of these corrections (~30%) greatly improves the agreement between the calculated and the observed magnon stiffness for cobalt and nickel, which are known to be overestimated within the ALDA.

A semirelativistic density-functional theory that includes spin-orbit couplings and Zeeman fields on equal footing with the electromagnetic potentials, is an appealing framework to develop a unified first-principles computational approach for noncollinear magnetism, spintronics, orbitronics, and topological states. The basic variables of this theory include the paramagnetic current and the spin-current density, besides the particle and the spin density, and the corresponding exchange-correlation (xc) energy functional is invariant under local U(1) x SU(2) gauge transformations. The xc-energy functional must be approximated to enable practical applications, but, contrary to the case of the standard density functional theory, finding simple approximations suited to deal with realistic atomistic inhomogeneities has been a long-standing challenge. Here we propose a way out of this impasse by showing that approximate gauge-invariant functionals can be easily generated from existing approximate functionals of ordinary density-functional theory by applying a simple minimal substitution on the kinetic energy density, which controls the short-range behavior of the exchange hole. Our proposal opens the way to the construction of approximate, yet nonempirical functionals, which do not assume weak inhomogeneity and therefore may have a wide range of applicability in atomic, molecular, and condensed matter physics.

We discuss the potential advantages of calculating the effectivemass of quasiparticles in the interacting electron liquid from the low-temperature free energy vis-a-vis the conventional approach, in which the effective mass is obtained from approximate calculations of the self-energy, or from a quantum Monte Carlo evaluation of the energy of a variational "quasiparticle wave function". While raw quantum Monte Carlo data are presently too sparse to allow for an accurate determination of the effective mass, the values estimated by this method are numerically close to the ones obtained in previous calculations using diagrammatic many-body theory. In contrast to this, a recently published parametrization of quantum Monte Carlo data for the free energy of the homogeneous electron liquid yields effective masses that considerably deviate from previous calculations and even change sign for low densities, reflecting an unphysical negative entropy. We suggest that this anomaly is related to the treatment of the exchange energy at finite temperature.