Electronic Degrees Of Freedom Research Articles

Non-adiabatic molecular dynamics simulations provide new insights into the exciton transfer in the Fenna-Matthews-Olson complex.

A trajectory surface hopping approach, which uses machine learning to speed up the most time-consuming steps, has been adopted to investigate the exciton transfer in light-harvesting systems. The present neural networks achieve high accuracy in predicting both Coulomb couplings and excitation energies. The latter are predicted taking into account the environment of the pigments. Direct simulation of exciton dynamics through light-harvesting complexes on significant time scales is usually challenging due to the coupled motion of nuclear and electronic degrees of freedom in these rather large systems containing several relatively large pigments. In the present approach, however, we are able to evaluate a statistically significant number of non-adiabatic molecular dynamics trajectories with respect to exciton delocalization and exciton paths. The formalism is applied to the Fenna-Matthews-Olson complex of green sulfur bacteria, which transfers energy from the light-harvesting chlorosome to the reaction center with astonishing efficiency. The system has been studied experimentally and theoretically for decades. In total, we were able to simulate non-adiabatically more than 30 ns, sampling also the relevant space of parameters within their uncertainty. Our simulations show that the driving force supplied by the energy landscape resulting from electrostatic tuning is sufficient to funnel the energy towards site 3, from where it can be transferred to the reaction center. However, the high efficiency of transfer within a picosecond timescale can be attributed to the rather unusual properties of the BChl a molecules, resulting in very low inner and outer-sphere reorganization energies, not matched by any other organic molecule, e.g., used in organic electronics. A comparison with electron and exciton transfer in organic materials is particularly illuminating, suggesting a mechanism to explain the comparably high transfer efficiency.

Optimal photoelectron circular dichroism of a model chiral system.

We optimize the internuclear geometry and electronic structure of a model chiral system to achieve a maximal photoelectron circular dichroism (PECD) in its one-photon ionization by circularly polarized light. The electronic structure calculations are performed by the single center method, while the optimization is done using quantum alchemy employing a Taylor series expansion. Thereby, the effect of bond lengths and uncompensated charge distributions on the chiral response of the model is investigated theoretically in some detail. It is demonstrated that manipulating a chiral asymmetry of the ionic potential may enhance the dichroic parameter (i.e., the PECD) of the randomly oriented model system well beyond β1 = 25%. Furthermore, we demonstrate that quantum alchemy is applicable to PECD despite the unusually strong coupling of spatial and electronic degrees of freedom and discuss the relative impact of the individual degrees of freedom in this model system. We define the necessary conditions for the computational design of PECD for real (non-model) chiral molecules using our approach.

Electronic Couplings for Triplet-Triplet Annihilation Upconversion in Crystal Rubrene.

Triplet-triplet annihilation photon upconversion (TTA-UC) is a process able to repackage two low-frequency photons into light of higher energy. This transformation is typically orchestrated by the electronic degrees of freedom within organic compounds possessing suitable singlet and triplet energies and electronic couplings. In this work, we propose a computational protocol for the assessment of electronic couplings crucial to TTA-UC in molecular materials and apply it to the study of crystal rubrene. Our methodology integrates sophisticated yet computationally affordable approaches to quantify couplings in singlet and triplet energy transfer, the binding of triplet pairs, and the fusion to the singlet exciton. Of particular significance is the role played by charge-transfer states along the b-axis of rubrene crystal, acting as both partial quenchers of singlet energy transfer and mediators of triplet fusion. Our calculations identify the π-stacking direction as holding notable triplet energy transfer couplings, consistent with the experimentally observed anisotropic exciton diffusion. Finally, we have characterized the impact of thermally induced structural distortions, revealing their key role in the viability of triplet fusion and singlet fission. We posit that our approaches are transferable to a broad spectrum of organic molecular materials, offering a feasible means to quantify electronic couplings.

Colossal c-Axis Response and Lack of Rotational Symmetry Breaking within the Kagome Planes of the CsV_{3}Sb_{5} Superconductor.

The kagome materials AV_{3}Sb_{5} (A=K, Rb, Cs) host an intriguing interplay between unconventional superconductivity and charge-density waves. Here, we investigate CsV_{3}Sb_{5} by combining high-resolution thermal-expansion, heat-capacity, and electrical resistance under strain measurements. We directly unveil that the superconducting and charge-ordered states strongly compete, and that this competition is dramatically influenced by tuning the crystallographic c axis. In addition, we report the absence of additional bulk phase transitions within the charge-ordered state, notably associated with rotational symmetry breaking within the kagome planes. This suggests that any breaking of the C_{6} invariance occurs via different stacking of C_{6}-symmetric kagome patterns. Finally, we find that the charge-density-wave phase exhibits an enhanced A_{1g}-symmetric elastoresistance coefficient, whose large increase at low temperature is driven by electronic degrees of freedom.

Disentangling the Unusual Magnetic Anisotropy of the Near‐Room‐Temperature Ferromagnet Fe4GeTe2

AbstractIn the quest for 2D conducting materials with high ferromagnetic ordering temperature the new family of the layered FenGeTe2 compounds, especially the near‐room‐temperature ferromagnet Fe4GeTe2, receives a significant attention. Fe4GeTe2 features a peculiar spin reorientation transition at TSR ≈ 110 K suggesting a non‐trivial temperature evolution of the magnetic anisotropy (MA)—one of the main contributors to the stabilization of the magnetic order in the low‐dimensional systems. An electron spin resonance (ESR) spectroscopic study reported here provides quantitative insights into the unusual magnetic anisotropy of Fe4GeTe2. At high temperatures the total MA is mostly given by the demagnetization effect with a small contribution of the counteracting intrinsic magnetic anisotropy of an easy‐axis type, whose growth below a characteristic temperature Tshape ≈ 150 K renders the sample seemingly isotropic at TSR. Below one further temperature Td ≈ 50 K the intrinsic MA becomes even more complex. Importantly, all the characteristic temperatures found in the ESR experiment match those observed in transport measurements, suggesting an inherent coupling between magnetic and electronic degrees of freedom in Fe4GeTe2. This finding together with the observed signatures of the intrinsic two‐dimensionality should facilitate optimization routes for the use of Fe4GeTe2 in the magneto‐electronic devices, potentially even in the monolayer limit.

Diagonalizing the Born-Oppenheimer Hamiltonian via Moyal perturbation theory, nonadiabatic corrections, and translational degrees of freedom.

This article describes a method for calculating higher order or nonadiabatic corrections in Born-Oppenheimer theory and its interaction with the translational degrees of freedom. The method uses the Wigner-Weyl correspondence to map nuclear operators into functions on the classical phase space and the Moyal star product to represent operator multiplication on those functions. These are explained in the body of the paper. The result is a power series in κ2, where κ = (m/M)1/4 is the usual Born-Oppenheimer parameter. The lowest order term is the usual Born-Oppenheimer approximation, while higher order terms are nonadiabatic corrections. These are needed in calculations of electronic currents, momenta, and densities. The separation of nuclear and electronic degrees of freedom takes place in the context of the exact symmetries (for an isolated molecule) of translations and rotations, and these, especially translations, are explicitly incorporated into our discussion. This article presents an independent derivation of the Moyal expansion in molecular Born-Oppenheimer theory. We show how electronic currents and momenta can be calculated within the framework of Moyal perturbation theory; we derive the transformation laws of the electronic Hamiltonian, the electronic eigenstates, and the derivative couplings under translations; we discuss in detail the rectilinear motion of the molecular center of mass in the Born-Oppenheimer representation; and we show how the elimination of the translational components of the derivative couplings leads to a unitary transformation that has the effect of exactly separating the translational degrees of freedom.

Coupled cluster cavity Born-Oppenheimer approximation for electronic strong coupling.

Chemical and photochemical reactivity, as well as supramolecular organization and several other molecular properties, can be modified by strong interactions between light and matter. Theoretical studies of these phenomena require the separation of the Schrödinger equation into different degrees of freedom as in the Born-Oppenheimer approximation. In this paper, we analyze the electron-photon Hamiltonian within the cavity Born-Oppenheimer approximation (CBOA), where the electronic problem is solved for fixed nuclear positions and photonic parameters. In particular, we focus on intermolecular interactions in representative dimer complexes. The CBOA potential energy surfaces are compared with those obtained using a polaritonic approach, where the photonic and electronic degrees of freedom are treated at the same level. This allows us to assess the role of electron-photon correlation and the accuracy of CBOA.

Dynamics and resilience of the unconventional charge density wave in ScV6Sn6 bilayer kagome metal

Long-range electronic ordering descending from a metallic parent state constitutes a rich playground to study the interplay of structural and electronic degrees of freedom. In this framework, kagome metals are in the most interesting regime where both phonon and electronically mediated couplings are significant. Several of these systems undergo a charge density wave transition. However, to date, the origin and the main driving force behind this charge order is elusive. Here, we use the kagome metal ScV6Sn6 as a platform to investigate this problem, since it features both a kagome-derived nested Fermi surface and van-Hove singularities near the Fermi level, and a charge-ordered phase that strongly affects its physical properties. By combining time-resolved reflectivity, first principles calculations and photo-emission experiments, we identify the structural degrees of freedom to play a fundamental role in the stabilization of charge order, indicating that ScV6Sn6 features an instance of charge order predominantly originating from phonons.

Precursor region with full phonon softening above the charge-density-wave phase transition in 2H-TaSe2

Research on charge-density-wave (CDW) ordered transition-metal dichalcogenides continues to unravel new states of quantum matter correlated to the intertwined lattice and electronic degrees of freedom. Here, we report an inelastic x-ray scattering investigation of the lattice dynamics of the canonical CDW compound 2H-TaSe2 complemented by angle-resolved photoemission spectroscopy and density functional perturbation theory. Our results rule out the formation of a central-peak without full phonon softening for the CDW transition in 2H-TaSe2 and provide evidence for a novel precursor region above the CDW transition temperature TCDW, which is characterized by an overdamped phonon mode and not detectable in our photoemission experiments. Thus, 2H-TaSe2 exhibits structural before electronic static order and emphasizes the important lattice contribution to CDW transitions. Our ab-initio calculations explain the interplay of electron-phonon coupling and Fermi surface topology triggering the CDW phase transition and predict that the CDW soft phonon mode promotes emergent superconductivity near the pressure-driven CDW quantum critical point.

Elastocaloric signatures of symmetric and antisymmetric strain-tuning of quadrupolar and magnetic phases in DyB2C2

The adiabatic elastocaloric effect measures the temperature change of a given system with strain and provides a thermodynamic probe of the entropic landscape in the temperature-strain space. Here, we demonstrate that the DC bias strain-dependence of AC elastocaloric effect allows decomposition of the latter into symmetric (rotation-symmetry-preserving) and antisymmetric (rotation-symmetry-breaking) strain channels, using a tetragonal [Formula: see text]-electron intermetallic DyB[Formula: see text]C[Formula: see text]-whose antiferroquadrupolar order breaks local fourfold rotational symmetries while globally remaining tetragonal-as a showcase example. We capture the strain evolution of its quadrupolar and magnetic phase transitions using both singularities in the elastocaloric coefficient and its jumps at the transitions, and the latter we show follows a modified Ehrenfest relation. We find that antisymmetric strain couples to the underlying order parameter in a biquadratic (linear-quadratic) manner in the antiferroquadrupolar (canted antiferromagnetic) phase, which are attributed to a preserved (broken) global tetragonal symmetry, respectively. The broken tetragonal symmetry in the magnetic phase is further evidenced by elastocaloric strain-hysteresis and optical birefringence. Additionally, within the staggered quadrupolar order, the observed elastocaloric response reflects a quadratic increase of entropy with antisymmetric strain, analogous to the role magnetic field plays for Ising antiferromagnetic orders by promoting pseudospin flips. Our results demonstrate AC elastocaloric effect as a compact and incisive thermodynamic probe into the coupling between electronic degrees of freedom and strain in free energy, which holds the potential for investigating and understanding the symmetry of a wide variety of ordered phases in broader classes of quantum materials.

A multi-stage, first-order phase transition in LaFe11.8Si1.2: Interplay between the structural, magnetic, and electronic degrees of freedom

Alloys with a first-order magnetic transition are central to solid-state refrigeration technology, sensors and actuators, or spintronic devices. The discontinuous nature of the transition in these materials is a consequence of the coupling between the magnetic, electronic, and structural subsystems, and such transition can, in principle, cross several metastable states, where at one point, the transition takes place within the magnetic subsystem, while at another, the changes occur in the structural or electronic subsystems. To address this issue, we conducted simultaneous measurements of the macroscopic properties—magnetization, temperature change of the sample, longitudinal, and transversal magnetostrictions—to reveal the rich details of the magneto-structural, first-order transition occurring in the prototypical alloy LaFe11.8Si1.2. We found that the transition does not complete in one but in two distinct stages. The presence of the intermediate state changes the potential-energy landscape, which then impacts strongly on the width of the hysteresis associated with the first-order transition. We complement these findings with experiments on the atomistic scale, i.e., x-ray absorption spectroscopy, x-ray magnetic circular dichroism, and Mössbauer spectroscopy, and then combine them with first-principles calculations to reveal the full complexity and two-stage nature of the transition. This new approach can be successfully extended to a large class of advanced magnetic materials that exhibit analogous transformations.

Theoretical study of nonadiabatic hydrogen atom scattering dynamics on metal surfaces using the hierarchical equations of motion method.

Hydrogen atom scattering on metal surfaces is investigated based on a simplified Newns-Anderson model. Both the nuclear and electronic degrees of freedom are treated quantum mechanically. By partitioning all the surface electronic states as the bath, the hierarchical equations of motion method for the fermionic bath is employed to simulate the scattering dynamics. It is found that, with a reasonable set of parameters, the main features of the recent experimental studies of hydrogen atom scattering on metal surfaces can be reproduced. Vibrational states on the chemisorption state whose energies are close to the incident energy are found to play an important role, and the scattering process is dominated by a single-pass electronic transition forth and back between the diabatic physisorption and chemisorption states. Further study on the effects of the atom-surface coupling strength reveals that, upon increasing the atom-surface coupling strength, the scattering mechanism changes from typical nonadiabatic transitions to dynamics in the electronic friction regime.

Observation of local atomic displacements intrinsic to the double zigzag chain structure of 1T−MTe2 ( M=V, Nb, Ta )

We describe the existence of local distortion discovered in the synchrotron x-ray single-crystal structure analysis of layered ditelluride $1T\text{\ensuremath{-}}M{\mathrm{Te}}_{2} (M=\mathrm{V}$, Nb, Ta). In $1T\text{\ensuremath{-}}{\mathrm{TaTe}}_{2}$, the double zigzag chain structure of Ta is deformed at about 170 K, and heptamer molecules are formed periodically at low temperatures. We found that some of the Ta atoms that compose the double zigzag chain structure appearing at high temperatures are locally displaced, resulting in local dimerization. This tendency weakens when Ta is replaced by V or Nb. Our results indicate that the local distortion persistently survives in these ditellurides, where the electronic degrees of freedom, including orbitals, are weakened. We further discuss the origin of local distortion in these ditellurides, which is different from many usual material systems where molecular formation occurs at low temperatures.

Revisiting the electronic properties of disclinated graphene sheets

The interplay between topological defects, such as dislocations or disclinations, and the electronic degrees of freedom in graphene have been extensively studied. In the literature, for the study of this kind of problems, it is in general used either a gauge theory or a curved spatial Riemannian geometry approach, where, in the geometric case, the information about the defects is contained in the metric and the spin-connection. However, these topological defects can also be associated with a Riemann–Cartan geometry where curvature and torsion plays an important role. In this article, we study the interplay between a wedge disclination in a planar graphene sheet and the properties of its electronic degrees of freedom. Our approach relies on its relation with elasticity theory through the so called elastic-gauge, where their typical coefficients, as for example the Poisson’s ratio, appear directly in the metric, and consequently also in the electronic spectrum.

Isosymmetric Lattice Collapse in Mixed-Valence Rare-Earth Fullerides at High Pressure─Coupling of Lattice and Electronic Degrees of Freedom

The orthorhombic-structured (Sm1/3Ca2/3)2.75C60 is a member of the family of rare-earth metal intercalated C60 fullerides with stoichiometry, RE2.75C60 (RE = rare-earth metal). At ambient conditions, it crystallizes with the formation of a 2a × 2a × 2a (a ∼ 14 Å) supercell of the face-centered cubic (fcc) A3C60 (A = alkali metal) unit cell. The superstructure is stabilized by the long-range ordering of tetrahedral-site Sm/Ca metal defects (O phase). Synchrotron X-ray powder diffraction measurements at high pressure and ambient temperature reveal that (Sm1/3Ca2/3)2.75C60 is a compressible solid with a bulk modulus, K0 = 24(1) GPa, that transforms to a more densely packed isostructural high-pressure O′ phase above ∼4 GPa. The first-order phase transition retains the formation of the superstructure and is accompanied by a discontinuous lattice size decrease (ΔV/V ∼ 2%). The O′ phase is stabilized by the release of the steric crowding that develops upon compression; the Sm/Ca metal ions, which reside in the tetrahedral and octahedral holes of the fulleride sublattice, shift from their off-centered positions at low pressure (O phase) to nest at the centers of the interstices (O′ phase). The nearly exact coincidence of the pressure response of the reversible hysteretic O ↔ O′ phase transformation with that of the samarium valence transition from +2.33(2) to +2.71(3), as established before by synchrotron X-ray absorption spectroscopy, unambiguously establishes an intimate link between the crystal and electronic structures of the (Sm1/3Ca2/3)2.75C60 fulleride.

Valley and spin splittings in a valley layer locked monolayer via atomic adsorption

To date, a number of valley materials have been discovered with spin-valley or valley-layer couplings. It is highly desirable to realize the interplay of more electronic degrees of freedom in a valley material. Based on the first-principles calculations, we demonstrate the valley and spin degeneracy liftings in the band structure of a TiSiCO monolayer, with the help of the selective expression of the layer pseudospin via atomic adsorption. The introduction of the transition-metal adatoms provides an effective electric field and magnetic proximity effect, giving rise to considerable valley and spin splittings. The larger valley splitting and the largest spin splitting for the band-edge states have magnitudes of up to 190 meV and 148 meV, respectively. The valley and spin splittings can be further tuned by applying an external electric field. According to the modified band structure, various interlayer excitons with different combinations of spins and electric dipoles are selectively created under the optical field of appropriate frequencies. The tunable spin and valley splittings in atom-adsorbed TiSiCO monolayer offer opportunities for exploring the interactions between spin, valley, and layer pseudospins and designing advanced optoelectronic devices.

Shadow energy functionals and potentials in Born-Oppenheimer molecular dynamics.

In Born-Oppenheimer molecular dynamics (BOMD) simulations based on the density functional theory (DFT), the potential energy and the interatomic forces are calculated from an electronic ground state density that is determined by an iterative self-consistent field optimization procedure, which, in practice, never is fully converged. The calculated energies and forces are, therefore, only approximate, which may lead to an unphysical energy drift and instabilities. Here, we discuss an alternative shadow BOMD approach that is based on backward error analysis. Instead of calculating approximate solutions for an underlying exact regular Born-Oppenheimer potential, we do the opposite. Instead, we calculate the exact electron density, energies, and forces, but for an underlying approximate shadow Born-Oppenheimer potential energy surface. In this way, the calculated forces are conservative with respect to the approximate shadow potential and generate accurate molecular trajectories with long-term energy stabilities. We show how such shadow Born-Oppenheimer potentials can be constructed at different levels of accuracy as a function of the integration time step, δt, from the constrained minimization of a sequence of systematically improvable, but approximate, shadow energy density functionals. For each energy functional, there is a corresponding ground state Born-Oppenheimer potential. These pairs of shadow energy functionals and potentials are higher-level generalizations of the original "zeroth-level" shadow energy functionals and potentials used in extended Lagrangian BOMD [Niklasson, Eur. Phys. J. B 94, 164 (2021)]. The proposed shadow energy functionals and potentials are useful only within this extended dynamical framework, where also the electronic degrees of freedom are propagated as dynamical field variables together with the atomic positions and velocities. The theory is quite general and can be applied to MD simulations using approximate DFT, Hartree-Fock, or semi-empirical methods, as well as to coarse-grained flexible charge models.

Dynamical Simulations of Carotenoid Photoexcited States Using Density Matrix Renormalization Group Techniques.

We present a dynamical simulation scheme to model the highly correlated excited state dynamics of linear polyenes. We apply it to investigate the internal conversion processes of carotenoids following their photoexcitation. We use the extended Hubbard-Peierls model, , to describe the π-electronic system coupled to nuclear degrees of freedom. This is supplemented by a Hamiltonian, , that explicitly breaks both the particle-hole and two-fold rotation symmetries of idealized carotenoid structures. The electronic degrees of freedom are treated quantum mechanically by solving the time-dependent Schrödinger equation using the adaptive time-dependent DMRG (tDMRG) method, while nuclear dynamics are treated via the Ehrenfest equations of motion. By defining adiabatic excited states as the eigenstates of the full Hamiltonian, , and diabatic excited states as eigenstates of , we present a computational framework to monitor the internal conversion process from the initial photoexcited 11Bu+ state to the singlet triplet-pair states of carotenoids. We further incorporate Lanczos-DMRG to the tDMRG-Ehrenfest method to calculate transient absorption spectra from the evolving photoexcited state. We describe in detail the accuracy and convergence criteria for DMRG, and show that this method accurately describes the dynamical processes of carotenoid excited states. We also discuss the effect of the symmetry-breaking term, , on the internal conversion process, and show that its effect on the extent of internal conversion can be described by a Landau-Zener-type transition. This methodological paper is a companion to our more explanatory discussion of carotenoid excited state dynamics in Manawadu, D.; Georges, T. N.; Barford, W. Photoexcited State Dynamics and Singlet Fission in Carotenoids. J. Phys. Chem. A 2023, 127, 1342.

Non-equilibrium dynamics at the gas-liquid interface: State-resolved studies of NO evaporation from a benzyl alcohol liquid microjet.

First measurements of internal quantum-state distributions for nitric oxide (NO) evaporating from liquid benzyl alcohol are presented over a broad range of temperatures, performed by liquid-microjet techniques in an essentially collision-free regime, with rotational/spin-orbit populations in the 2Π1/2,3/2 manifolds measured by laser-induced fluorescence. The observed rotational distributions exhibit highly linear (i.e., thermal) Boltzmann plots but notably reflect rotational temperatures (Trot) as much as 30K lower than the liquid temperature (Tjet). A comparable lack of equilibrium behavior is also noted in the electronic degrees of freedom but with populations corresponding to spin-orbit temperatures (TSO) consistently higher than Trot by ∼15K. These results unambiguously demonstrate evaporation into a non-equilibrium distribution, which, by detailed-balance considerations, predict quantum-state-dependent sticking coefficients for incident collisions of NO at the gas-liquid interface. Comparison and parallels with previous experimental studies of NO thermal desorption and molecular-beam scattering in other systems are discussed, which suggests the emergence of a self-consistent picture for the non-equilibrium dynamics.

Multipolar quantum electrodynamics of localized charge-current distributions: Spectral theory and renormalization

We formulate a nonrelativistic quantum field theory to model interactions between quantized electromagnetic fields and localized charge-current distributions. The electronic degrees of freedom are encoded in microscopic polarization and magnetization field operators whose moments are identified with the multipole moments of the charge-current distribution. The multipolar Hamiltonian is obtained from the minimal coupling Hamiltonian through a unitary transformation, often referred to as the Power-Zienau-Woolley transformation; we renormalize this Hamiltonian using perturbation theory, the result of which is used to compute the leading-order radiative corrections to the electronic energy levels due to interactions between the electrons and quantum vacuum fluctuations in the electromagnetic field. Our renormalized energy shift constitutes a generalization of the Lamb shift in atomic hydrogen, valid for general localized assemblies of atoms and molecules, possibly with net charge but absent a free current. By expanding the fields in a series of multipole moments, our results can be used to study contributions to this energy shift coming from specific multipole moments of arbitrary order.