Vibronic Research Articles

Vibrational Electronic-Thermofield Coupled Cluster (VE-TFCC) Theory for Quantum Simulations of Vibronic Coupling Systems at Thermal Equilibrium.

A vibrational electronic-thermofield coupled cluster (VE-TFCC) approach is developed to calculate thermal properties of non-adiabatic vibronic coupling systems. The thermofield (TF) theory and a mixed linear exponential ansatz based on second-quantized Bosonic construction operators are introduced to propagate the thermal density operator as a "pure state" in the Bogoliubov representation. Through this compact representation of the thermal density operator, the approach is basis-set-free and scales classically (polynomial) as the number of degrees of freedoms (DoF) in the system increases. The VE-TFCC approach is benchmarked with small test models and a real molecular compound (CoF4- anion) against the conventional sum over states (SOS) method and applied to calculate thermochemistry properties of a gas-phase reaction: CoF3 + F- → CoF4-. Results shows that the VE-TFCC approach, in conjunction with vibronic models, provides an effective protocol for calculating thermodynamic properties of vibronic coupling systems.

Healthy Horndeski cosmologies with torsion

We show that the full Horndeski theory with both curvature and torsion can support nonsingular, stable and subluminal cosmological solutions at all times. Thus, with torsion, the usual No-Go theorem that holds in a curved spacetime is avoided. In particular, it is essential to include the nonminimal derivative couplings of the ℒ5 part of the Horndeski action (Gμν ∇ μ ∇ νϕ, and (∇2 ϕ)3). Without the latter a No-Go already impedes the eternal subluminality of nonsingular, stable cosmologies.

Generalized second law for non-minimally coupled matter theories

We establish the generalized second law (GSL) within the framework of higher curvature gravity theories, considering non-minimal couplings in the matter sector. Our proof pertains to the regime of linearized fluctuations around equilibrium black holes, aligning with previous works by Wall and Sarkar. Notably, while prior proofs addressed various gravity theories such as Lovelock theory and higher curvature gravity, they uniformly assumed minimally coupled matter sectors. In this work, we extend the proof of the linearized semi-classical GSL to encompass scenarios involving non-minimal couplings in the matter sector. Our approach involves a proposal for evaluation of the matter path integral in the expectation value of the stress tensor, adopting an effective field theory treatment for the higher derivative couplings. We leverage the recently established outcome regarding the linearized second law in such theories to substantiate our argument.

On the Topological Phase around Conical Intersections with Tamm-Dancoff Linear-Response Time-Dependent Density Functional Theory.

Regions of nuclear-configuration space away from the Franck-Condon geometry can prove problematic for some electronic structure methods, given the propensity of such regions to possess conical intersections, i.e., (highly connected) points of degeneracy between potential energy surfaces. With the likelihood (perhaps even inevitability) for nonadiabatic dynamics simulations to explore molecular geometries in close proximity to conical intersections, it is vital that the performance of electronic structure methods is routinely examined in this context. In a recent paper [Taylor, J. T. J. Chem. Phys. 2023, 159, 214115.], the ability of linear-response time-dependent density functional theory within the adiabatic approximation (AA LR-TDDFT) to provide a proper description of conical intersections, in terms of their topology and topography, was investigated, with particular attention paid to conical intersections between two excited electronic states. For the same prototypical molecules, protonated formaldimine and pyrazine, we herein consider whether AA LR-TDDFT can correctly reproduce the topological phase accumulated by the adiabatic electronic wave function upon traversing a closed path around an excited-to-excited state conical intersection despite not using the appropriate quadratic-response nonadiabatic coupling vectors. Equally, we probe the ability of the ground-to-excited state intersection ring exhibited by AA LR-TDDFT in protonated formaldimine to give rise to a similar topological phase in spite of its incorrect dimensionality.

Detailed Complementary Consistency: Wave Function Tells Particle How to Hop, Particle Tells Wave Function How to Collapse.

In mixed quantum-classical dynamics, the quantum subsystem can have both wave function and particle-like descriptions. However, they may yield inconsistent results for the expectation value of the same physical quantity. We here propose a novel detailed complementary consistency (DCC) method based on the principle of detailed internal consistency. Namely, the wave function along each trajectory tells the particle how to hop, while the particle tells the wave function how to collapse based on active states in the trajectory ensemble. As benchmarked in a diverse array of representative models with localized nonadiabatic couplings, DCC not only achieves fully consistent results (i.e., identical populations calculated based on wave functions and active states) but also closely reproduces the exact quantum results. Due to the high performance, our new DCC method has great potential to give a consistent and accurate mixed quantum-classical description of general nonadiabatic dynamics after further development.

Light harvesting in purple bacteria does not rely on resonance fine-tuning in peripheral antenna complexes

The ring-like peripheral light-harvesting complex 2 (LH2) expressed by many phototrophic purple bacteria is a popular model system in biological light-harvesting research due to its robustness, small size, and known crystal structure. Furthermore, the availability of structural variants with distinct electronic structures and optical properties has made this group of light harvesters an attractive testing ground for studies of structure–function relationships in biological systems. LH2 is one of several pigment-protein complexes for which a link between functionality and effects such as excitonic coherence and vibronic coupling has been proposed. While a direct connection has not yet been demonstrated, many such interactions are highly sensitive to resonance conditions, and a dependence of intra-complex dynamics on detailed electronic structure might be expected. To gauge the sensitivity of energy-level structure and relaxation dynamics to naturally occurring structural changes, we compare the photo-induced dynamics in two structurally distinct LH2 variants. Using polarization-controlled 2D electronic spectroscopy at cryogenic temperatures, we directly access information on dynamic and static disorder in the complexes. The simultaneous optimal spectral and temporal resolution of these experiments further allows us to characterize the ultrafast energy relaxation, including exciton transport within the complexes. Despite the variations in PPC molecular structure manifesting as clear differences in electronic structure and disorder, the energy-transport and—relaxation dynamics remain remarkably similar. This indicates that the light-harvesting functionality of purple bacteria within a single LH2 complex is highly robust to structural perturbations and likely does not rely on finely tuned electronic- or electron-vibrational resonance conditions.

Nonadiabatic molecular dynamics with subsystem density functional theory: application to crystalline pentacene

In this work, we report the development and assessment of the nonadiabatic molecular dynamics approach with the electronic structure calculations based on the linearly scaling subsystem density functional method. The approach is implemented in an open-source embedded Quantum Espresso/Libra software specially designed for nonadiabatic dynamics simulations in extended systems. As proof of the applicability of this method to large condensed-matter systems, we examine the dynamics of nonradiative relaxation of excess excitation energy in pentacene crystals with the simulation supercells containing more than 600 atoms. We find that increased structural disorder observed in larger supercell models induces larger nonadiabatic couplings of electronic states and accelerates the relaxation dynamics of excited states. We conduct a comparative analysis of several quantum-classical trajectory surface hopping schemes, including two new methods proposed in this work (revised decoherence-induced surface hopping and instantaneous decoherence at frustrated hops). Most of the tested schemes suggest fast energy relaxation occurring with the timescales in the 0.7–2.0 ps range, but they significantly overestimate the ground state recovery rates. Only the modified simplified decay of mixing approach yields a notably slower relaxation timescales of 8–14 ps, with a significantly inhibited ground state recovery.

Energy-Conserving Surface Hopping for Auger Processes.

Auger-type processes are ubiquitous in nanoscale materials because quantum confinement enhances Coulomb interactions, and there exist large densities of states. Modeling Auger processes requires the modification of nonadiabatic (NA) molecular dynamics algorithms to include transitions caused by both NA and Coulomb couplings. The system is split into quantum and classical subsystems, e.g., electrons and vibrations, and as a result, energy conservation becomes nontrivial. In surface hopping, an electronic transition induced by NA coupling is accompanied by a classical velocity readjustment to ensure conservation of the total quantum-classical energy. A different treatment is needed for Auger transitions driven by Coulomb interactions. We develop a nonadiabatic molecular dynamics methodology that meticulously differentiates the energy redistribution accompanying hops induced by the NA coupling and the Coulomb interaction and correctly conserves the total energy at each transition. If the transition is driven by a Coulomb interaction, the hop energy is redistributed within the quantum electronic subsystem only. If the transition is NA, the energy is redistributed between the quantum and classical subsystems. Properly maintaining energy conservation for both types of transitions is crucial to generate a correct order of events, obtain accurate transition times, maintain a proper statistical distribution of state populations, and reach thermodynamic equilibrium. We test the method with biexciton annihilation and Auger-assisted hot electron relaxation in a CdSe quantum dot. The sequence of Auger and phonon-driven processes and the calculated time scales are in excellent agreement with the experimental results. The developed approach can be coupled with any surface-hopping method and provides a crucial practical advance to study charge-carrier dynamics in the nanoscale and condensed matter systems.

Optimizing Conical Intersections without Explicit Use of Non-Adiabatic Couplings.

We present two alternative methods for optimizing minimum energy conical intersection (MECI) molecular geometries without knowledge of the derivative coupling (DC). These methods are based on the utilization of Lagrange multipliers: (i) one method uses an approximate calculation of the DC, while the other (ii) do not require the DC. Both methods use the fact that information on the DC is contained in the Hessian of the squared energy difference. Tests done on a set of small molecular systems, in comparison with other methods, show the ability of the proposed methods to optimize MECIs. Finally, we apply the methods to the furimamide molecule, to optimize and characterize its S1/S2 MECI, and to optimizing the S0/S1 MECI of the silver trimer.

An Effective Model for Capturing the Role of Excitonic Interactions in the Wave-Packet Dynamics of DNA Nucleobases

Investigating exciton dynamics within DNA nucleobases is essential for comprehensively understanding how inherent photostability mechanisms function at the molecular level, particularly in the context of life’s resilience to solar radiation. In this paper, we introduce a mathematical model that effectively simulates the photoexcitation and deactivation dynamics of nucleobases within an ultrafast timeframe, particularly focusing on wave-packet dynamics under conditions of strong nonadiabatic coupling. Employing the hierarchy equation of motion, we simulate two-dimensional electronic spectra (2DES) and calibrate our model by comparing it with experimentally obtained spectra. This study also explores the effects of base stacking on the photo-deactivation dynamics in DNA. Our results demonstrate that, while strong excitonic interactions between nucleobases are present, they have a minimal impact on the deactivation dynamics of the wave packet in the electronic excited states. We further observe that the longevity of electronic excited states increases with additional base stacking and pairing, a phenomenon accurately depicted by our excitonic model. This model enables a detailed examination of the wave packet’s motion on electronic excited states and its rapid transition to the ground state. Additionally, using this model, we studied base stacks in DNA hairpins to effectively capture the primary exciton dynamics at a reasonable computational scale. Overall, this work provides a valuable framework for studying exciton dynamics from single nucleobases to complex structures such as DNA hairpins.

Nuclear Quantum Effects Accelerate Charge Separation and Recombination in g-C3N4/TiO2 Heterojunctions.

We combined ring-polymer molecular dynamics (MD) and ab initio MD with nonadiabatic MD to study the effects of nuclear quantum effects (NQEs) on interlayer electron transfer and electron-hole recombination at the g-C3N4/TiO2 interface. Our simulations indicate that NQEs significantly affect electron transfer and electron-hole recombination dynamics, accelerating both processes. NQEs deform the g-C3N4 layer and expedite the movement of carbon and nitrogen atoms, thus, enhancing charge delocalization and interlayer coupling. This improved overlap between electronic state wave functions enhances nonadiabatic couplings, facilitating electron transfer and recombination. In addition to the enhanced nonadiabatic couplings accelerating electron transfer, the presence of NQEs narrows the energy gap and delays decoherence by mitigating overall fluctuations, because of restricted TiO2 movements overwhelming enhanced g-C3N4 fluctuations, thereby making the recombination faster. This work provides valuable insights into NQEs in light-element systems and contributes to guiding the development of highly efficient photocatalysts.

Nonadiabatic Coupling-Induced Quantum Coherence in Two-Dimensional Materials.

Two-dimensional materials provide a rich platform demonstrating quantum effects, and the process of electron-hole recombination occurring in them has significant applications in the fields of the photocatalytic and optoelectronic community. Here, we present nonadiabatic coupling-induced quantum coherence and quantum beats in Al-doped blue phosphorene. The work improves our understanding and utilization of nonadiabatic coupling in low-dimensional materials from a new perspective. In addition, our investigations provide meaningful guidance for manipulating quantum coherence in low-dimensional materials and promoting their novel optoelectronic properties.

Vibrationally Mode-Specific Molecular Energy Transfer to Surface Electrons in Metastable Formaldehyde Scattering from Cesium-Covered Au(111).

Nonadiabatic interaction of adsorbate nuclear motion with the continuum of electronic states is known to affect the dynamics of chemical reactions at metal surfaces. A large body of work has probed the fundamental mechanisms of such interactions for atomic and diatomic molecules at surfaces. In polyatomic molecules, the possibility of mode-specific damping of vibrational motion due to the effects of electronic friction raises the question of whether such interactions could profoundly affect the outcome of chemistry at surfaces by selectively removing energy from a particular intramolecular adsorbate mode. However, to date, there have not been any fundamental experiments demonstrating nonadiabatic electron-vibration coupling in a polyatomic molecule at a surface. In this work, we scatter excited metastable formaldehyde and formaldehyde-d2 from a low work function surface and detect ejected exoelectrons that accompany molecular relaxation. The exoelectron ejection efficiency exhibits a strong dependence on the vibrational mode that is excited: out-of-plane bending excitation (ν4) leads to significantly more exoelectrons than does CO stretching excitation (ν2). The results provide clear evidence for mode-specific energy transfer from vibration to surface electrons.

The anion photoelectron spectrum and diabatization of tetrazolyl.

The potential energy surface of tetrazolyl [cyclic (N4CH)] has a conical intersection seam between the two lowest-energy electronic states near the ground state minimum geometry. This work treats that molecule. The potential energy surfaces used in this study are based on a least-squares fitting procedure that includes abinitio energies, energy gradients, and derivative couplings described using polynomials up to fourth-order and abinitio data obtained from multireference configuration interaction wave functions. A five-electronic-state description was generated with a root mean square absolute energy error of 9.6cm-1, compared to 326.8cm-1 when only second-order terms were used. The time-independent multimode vibronic coupling in the KDC approximation was used to simulate and analyze the anion ultraviolet photoelectron spectrum of tetrazolide.

Axial H-Bonding Solvent Controls Inhomogeneous Spectral Broadening, While Peripheral H-Bonding Solvent Controls Vibronic Broadening: Cresyl Violet in Methanol.

The dynamics of the nuclei of both a chromophore and its condensed-phase environment control many spectral features, including the vibronic and inhomogeneous broadening present in spectral line shapes. For the cresyl violet chromophore in methanol, we here analyze and isolate the effect of specific chromophore-solvent interactions on simulated spectral densities, reorganization energies, and linear absorption spectra. Employing both chromophore and its condensed-phase environment control many spectral features, including the vibronic and inhomogeneous broadening present in spectral line shapes. For the cresyl violet chromophore in methanol, we here analyze and isolate the effect of specific chromophore-solvent interactions on simulated spectral densities, reorganization energies, and linear absorption spectra. Employing both force field and ab initio molecular dynamics trajectories along with the inclusion of only certain solvent molecules in the excited-state calculations, we determine that the methanol molecules axial to the chromophore are responsible for the majority of inhomogeneous broadening, with a single methanol molecule that forms an axial hydrogen bond dominating the response. The strong peripheral hydrogen bonds do not contribute to spectral broadening, as they are very stable throughout the dynamics and do not lead to increased energy-gap fluctuations. We also find that treating the strong peripheral hydrogen bonds as molecular mechanical point charges during the molecular dynamics simulation underestimates the vibronic coupling. Including these peripheral hydrogen bonding methanol molecules in the quantum-mechanical region in a geometry optimization increases the vibronic coupling, suggesting that a more advanced treatment of these strongly interacting solvent molecules during the molecular dynamics trajectory may be necessary to capture the full vibronic spectral broadening.

Troubles mounting for multipolar dark matter

In this paper, we revisit the experimental constraints on the multipolar dark matter that has derivative coupling to the visible sector mediated by the Standard Model photon. The momentum dependent interaction enables them to be captured efficiently within massive celestial bodies boosted by their steep gravitational potential. This phenomena makes compact celestial bodies as an efficient target to probe such type of dark matter candidates. We demonstrate that a synergy of the updated direct detection results from DarkSide-50 and LUX-ZEPLIN together with IceCube bounds on high energy solar neutrinos from dark matter capture disfavour the viable parameter space of the dipolar dark matter scenario. Whereas, for the anapole dark matter scenario, a narrow window survives that lies within the reach of prospective heating signals due to the capture of dark matter at cold neutron stars.

R‐NG●‐Regulated Electron‐Phonon Coupling to Enhance Hot‐Hole Injection for Inverted Perovskite Solar Cells

AbstractExtracting hot‐carriers before they relax back to the band edge will reduce thermal dissipation above the bandgap, which is one of the primary sources of efficiency loss in perovskite solar cells (PSCs). It requires slow cooling of hot‐carriers together with the efficient extraction of hot‐carriers. Here, a novel interface engineering is demonstrated by embedding radical‐containing nanographene (r‐NG●) between perovskite and poly[bis(4‐phenyl)(2,4‐dimethoxyphenyl)‐amine] (PTOAA) to harvest this excess energy. This strategy accelerates the extraction rate of the hot‐hole (3 times that of control at 400 K) by augmenting the relaxation channel, which is related to the “quasi‐SOMO” energy state of r‐NG●. Nonadiabatic molecular dynamics calculation further revealed that radical character promotes the rotation of the cation in perovskite, which contributes to generating strong nonadiabatic couplings and multiple phonon modes. In addition, the band bending effect of buried perovskite and improved conductivity of PTOAA also facilitate exciton dissociation and hole collection efficiency (97.5%). Consequently, the resultant r‐NG● decorated PSCs delivered an impressive efficiency of 24.20%, along with better thermal and humidity stability. Moreover, it can maintain 93% of its initial efficiency after 300 hours of illumination.

Inflation and reheating in quadratic metric-affine gravity with derivative couplings

Within the framework of metric-affine theories of gravity, where both the metric and connection are treated as independent variables, we consider actions quadratic in the Ricci scalar curvature coupled non-minimally to a scalar field through derivative couplings. Our analysis delves into the inflationary predictions, revealing their consistency with the latest observational constraints across a wide range of parameters. This compatibility permits adjustments such as an increase in the spectral index and a reduction in the tensor-to-scalar ratio. While we do not propose a specific reheating mechanism, our analysis demonstrates that within the quadratic model of inflation, the maximum reheating temperature can reach ∼ 3 × 1015 GeV.

Gravitating kinks with asymptotically flat metrics

In this work, we consider a two-dimensional (2D) dilaton gravity model where the dilaton kinetic term is modified by an additional derivative coupling term . In the case with a canonical scalar matter field, the field equations of this model have a simple first-order formalism, from which exact static kink solutions can be constructed. The novelty of these solutions is that the corresponding metric can be asymptotically flat rather than asymptotically anti-de Sitter. The linear stability and the localization of scalar matter fields are also studied. It was found that the solutions are stable against small linear perturbations, and the localization of scalar matter fields can be realized by introducing scalar-kink interactions.

Generalized Semiclassical Ehrenfest Method: A Route to Wave Function-Free Photochemistry and Nonadiabatic Dynamics with Only Potential Energies and Gradients.

We reconsider recent methods by which direct dynamics calculations of electronically nonadiabatic processes can be carried out while requiring only adiabatic potential energies and their gradients. We show that these methods can be understood in terms of a new generalization of the well-known semiclassical Ehrenfest method. This is convenient because it eliminates the need to evaluate electronic wave functions and their matrix elements along the mixed quantum-classical trajectories. The new approximations and procedures enabling this advance are the curvature-driven approximation to the time-derivative coupling, the generalized semiclassical Ehrenfest method, and a new gradient correction scheme called the time-derivative matrix (TDM) scheme. When spin-orbit coupling is present, one can carry out dynamics calculations in the fully adiabatic basis using potential energies and gradients calculated without spin-orbit coupling plus the spin-orbit coupling matrix elements. Even when spin-orbit coupling is neglected, the method is useful because it allows calculations by electronic structure methods for which nonadiabatic coupling vectors are unavailable. In order to place the new considerations in context, the article starts out with a review of background material on trajectory surface hopping, the semiclassical Ehrenfest scheme, and methods for incorporating decoherence. We consider both internal conversion and intersystem crossing. We also review several examples from our group of successful applications of the curvature-driven approximation.