Search

Photo-dynamics can proceed differently at the water/air interface compared to in the respective bulk phases. Second-order non-linear spectroscopy is capable of selectively probing the dynamics of species in such an environment. However, certain conclusions drawn from vibrational and electronic sum-frequency generation spectroscopies do not agree as is the case for the formation and structure of hydrated electrons at the interface. This Perspective aims to highlight these apparent discrepancies, how they can be reconciled, suggests how the two techniques complement one another, and outline the value of performing both techniques on the same system. [ABSTRACT FROM AUTHOR]

We investigate the dynamics and site-selectivity in the dissociation of transient anions formed upon attachment of low energy electrons to acetic acid by anion fragment momentum imaging experiments. The resonances at 6.7 and 7.7 eV are confirmed to dissociate exclusively by the O–H bond, while a third resonance at 9.1 eV dissociates primarily by both C–H break and O–H break. A fourth resonance near 10 eV is found to dissociate by O–H break. For each resonance, the measured kinetic energy release indicates two-body dissociation produces a neutral radical in the ground electronic state, for all four resonances. The measured angular distributions are consistent with all four resonances having A′ symmetry. [ABSTRACT FROM AUTHOR]

The electron-impact ionization of water molecules at low impact energies is investigated using a theoretical approach named M3CWZ. In this model, which considers exchange effects and post-collision interaction, the continuum electrons (incident, scattered, and ejected) are all described by a Coulomb wave that corresponds to distance-dependent charges generated from the molecular target properties. Triple differential cross-sections for low impact energy ionization of either the 1b1 or 3a1 orbitals are calculated for several geometrical and kinematical configurations, all in the dipole regime. The M3CWZ model is thoroughly tested with an extensive comparison with available theoretical results and COLTRIMS measurements performed at projectile energies of Ei = 81 eV [Ren et al., Phys. Rev. A 95, 022701 (2017)] and Ei = 65 eV [Zhou et al., Phys. Rev. A 104, 012817 (2021)]. Similar to other theoretical models, an overall good agreement with both sets of measured data is observed for the angular distributions. Our calculated cross-sections' magnitudes are also satisfactory when compared to the other theoretical results, as well as to the cross-normalized relative scale data at 81 eV impact energy. The 65 eV set of data, measured on an absolute scale, offers a further challenging task for theoretical descriptions, and globally the M3CWZ performs fairly well and comparably to other theories. The proposed approach with variable charges somehow allows to capture the main multicenter distortion effects while avoiding high computational costs. [ABSTRACT FROM AUTHOR]

The accuracy of density-functional theory (DFT) calculations is ultimately determined by the quality of the underlying approximate functionals, namely the exchange-correlation functional in electronic DFT and the excess functional in the classical DFT formalism of fluids. For both electrons and fluids, the exact functional is highly nonlocal, yet most calculations employ approximate functionals that are semi-local or nonlocal in a limited weighted-density form. Machine-learned (ML) nonlocal density-functional approximations show promise in advancing applications of both electronic and classical DFTs, but so far these two distinct research areas have implemented disparate approaches with limited generality. Here, we formulate a universal ML framework and training protocol to learn nonlocal functionals that combine features of equivariant convolutional neural networks and the weighted-density approximation. We prototype this new approach for several 1D and quasi-1D problems and demonstrate that functionals with exactly the same hyperparameters achieve excellent accuracy for a diverse set of systems, including the hard-rod fluid, the inhomogeneous Ising model, the exact exchange energy of electrons, the electron kinetic energy for orbital-free DFT, as well as for liquid water with 1D inhomogeneities. These results lay the foundation for a universal ML approach to approximate exact 3D functionals spanning electronic and classical DFTs. [ABSTRACT FROM AUTHOR]

Dissociative electron attachment rate constants have been measured for Cr(CO)6 under thermal conditions, 296–400 K, yielding Cr(CO)5− product. At 296 K, 2.92 ± 0.70 cm3 s−1 was measured and a small decrease with temperature was observed (2.72 ± 0.70 cm3 s−1 at 400 K). We additionally determined the cation products of Ar+ reacting with Cr(CO)6. [ABSTRACT FROM AUTHOR]

We apply the analytically solvable model of two electrons in two orbitals to diradical molecules, characterized by two unpaired electrons. The effect of doubly occupied and empty orbitals is taken into account by means of random phase approximation (RPA). We show that in the static limit, the direct RPA leads to the renormalization of the parameters of the two-orbital model. We test our model by comparing its predictions for singlet–triplet splitting with the results of several multi-reference methods for a set of thirteen molecules. We find that for this set, the static RPA results are close to those of the NEVPT2 method with two orbitals and two electrons in the active space. [ABSTRACT FROM AUTHOR]

This second paper in a series of two describes the chirped-pulse ice apparatus that permits the detection of buffer gas cooled molecules desorbed from an energetically processed ice using broadband mm-wave rotational spectroscopy. Here, we detail the lower ice stage developed to generate ices at 4 K, which can then undergo energetic processing via UV/VUV photons or high-energy electrons and which ultimately enter the gas phase via temperature-programmed desorption (TPD). Over the course of TPD, the lower ice stage is interfaced with a buffer gas cooling cell that allows for sensitive detection via chirped-pulse rotational spectroscopy in the 60–90 GHz regime. In addition to a detailed description of the ice component of this apparatus, we show proof-of-principle experiments demonstrating the detection of H2CO products formed through irradiation of neat methanol ices or 1:1 CO + CH4 mixed ices. [ABSTRACT FROM AUTHOR]

Dipole-bound states in anions exist when a polar neutral core binds an electron in a diffuse orbital through charge–dipole interaction. Electronically excited polar neutral cores can also bind an electron in a diffuse orbital to form Core-Excited Dipole-Bound States (CE-DBSs), which are difficult to observe because they usually lie above the electron detachment threshold, leading to very short lifetimes and, thus, unstructured transitions. We report here the photodetachment spectroscopy of cryogenically cooled acetylacetonate anion (C5H7O2−) recorded by detecting the neutral radical produced upon photodetachment and the infrared spectroscopy in He-nanodroplets. Two DBSs were identified in this anion. One of them lies close to the electron detachment threshold (∼2.74 eV) and is associated with the ground state of the radical (D0-DBS). Surprisingly, the other DBS appears as resonant transitions at 3.69 eV and is assigned to the CE-DBS associated with the first excited state of the radical (D1-DBS). It is proposed that the resonant transitions of the D1-DBS are observed ∼1 eV above the detachment threshold because its lifetime is determined by the internal conversion to the D0-DBS, after which the fast electron detachment takes place. [ABSTRACT FROM AUTHOR]

We employ the exact-factorization formalism to study the coupled dynamics of photons, electrons, and nuclei at the quantum mechanical level, proposing illustrative examples of model situations of nonadiabatic dynamics and spontaneous emission of electron–nuclear systems in the regime of strong light–matter coupling. We make a particular choice of factorization for such a multi-component system, where the full wavefunction is factored as a conditional electronic amplitude and a marginal photon–nuclear amplitude. Then, we apply the coupled-trajectory mixed quantum–classical (CTMQC) algorithm to perform trajectory-based simulations, by treating photonic and nuclear degrees of freedom on equal footing in terms of classical-like trajectories. The analysis of the time-dependent potentials of the theory along with the assessment of the performance of CTMQC allows us to point out some limitations of the current approximations used in CTMQC. Meanwhile, comparing CTMQC with other trajectory-based algorithms, namely multi-trajectory Ehrenfest and Tully surface hopping, demonstrates the better quality of CTMQC predictions. [ABSTRACT FROM AUTHOR]

Fully relativistic calculations, primarily at the 4-component coupled-cluster singles and doubles with perturbative triples [CCSD(T)] level of theory with the Dirac–Coulomb (DC) Hamiltonian, have been carried out for the superheavy (E125)F molecule using large Gaussian basis sets. The electronic ground state is determined to have an [Og]8s25g16f3 configuration on E125 with an Ω = 6 ground state and an 8p electron largely donated to F. A Mulliken population analysis indicates that the ground state is mainly ionic with a partial charge of +0.79 on E125 and a single sigma bond involving the F 2p and E125 8p spinors. The occupied g spinor is not involved in the bonding. With the largest basis set used in this work, the (0 K) dissociation energy was calculated at the DC-CCSD(T) level of theory to be 7.02 eV. Analogous calculations were also carried out for the E125 atom, both the neutral and its cation. The lowest energy electron configuration of E125+, [Og] 8 s 1 / 2 2 5 g 7 / 2 1 6 f 5 / 2 3 with a J = 6 ground state, was found to be similar to that in (E125)F, while the neutral E125 atom has an [Og] 8 s 1 / 2 2 5 g 7 / 2 1 6 f 5 / 2 2 7 d 3 / 2 1 8 p 1 / 2 1 ground state electron configuration with a J = 17/2 ground state. The ionization energy (IE) of E125 is reported for the first time and is calculated to be 4.70 eV at the DC-CCSD(T) level of theory. Non-relativistic calculations were also carried out on the E125 atom and the (E125)F molecule. The non-relativistic ground state of the E125 atom was calculated to have a 5g5 ground state with an IE of just 3.4 eV. The net effect of relativity on (E125)F is to stabilize its bonding. [ABSTRACT FROM AUTHOR]

The quantum nuclear kinematic interactions with electrons (or nonadiabatic interactions) are the inherent driving force that possibly causes a mixture of the adiabatic electronic states in molecules. Particularly in systems whose electron wavepackets lie in a densely quasi-degenerate electronic-state manifold where many-dimensional and many-state nonadiabatic interactions last continually, we have found before that those extensive mixings can lead to a quantum electronic-state chaos [K. Takatsuka and Y. Arasaki, J. Chem. Phys. 159, 074110 (2023)]. This chaos of electron dynamics is a new kind yet generic. This Communication identifies the mathematical/physical mechanism of this class of chaos by means of the collective coordinate analysis of the nonadiabatic interactions, along with the numerical applications to excited states of boron clusters. Some physical consequences of the present chaos are also discussed. [ABSTRACT FROM AUTHOR]

Electron transfer plays a crucial role in living systems, including the generation of reactive oxygen species (ROS). Oxygen acts as the terminal electron acceptor in the respiratory chains of aerobic organisms as well as in some photoinduced processes followed by the formation of ROS. This is why the participation of exogenous antioxidants in electron transfer processes in living systems is of particular interest. In the present study, using chemically induced dynamic nuclear polarization (CIDNP) and dissociative electron attachment (DEA) techniques, we have elucidated the affinity of solvated and free electrons to glycyrrhetinic acid (GA)—the aglicon of glycyrrhizin (the main active component of Licorice root). CIDNP is a powerful instrument to study the mechanisms of electron transfer reactions in solution, but the DEA technique shows its effectiveness in gas phase processes. For CIDNP experiments, the photoionization of the dianion of 5-sulfosalicylic acid (HSSA2−) was used as a model reaction of solvated electron generation. DEA experiments testify that GA molecules are even better electron acceptors than molecular oxygen, at least under gas-phase conditions. In addition, the effect of the solvent on the energetics of the reactants is discussed. [ABSTRACT FROM AUTHOR]

The efficient and stable electroluminescence of quantum dots (QDs) is of great importance in their applications in new display technologies. The short service life of blue QDs, however, hinders their development and commercialization. Different mechanisms have been proposed for the destabilization of QDs in electroluminescent processes. Based on real-time time-dependent density functional theory studies on the QD models covered by Z-type ligands (XAc2, X = Cd, Zn, Mg), the structural evolution is simulated to reveal the mechanism of the reduction reactions induced by electron injection. Our simulations reproduce the experimental observations that the reduction reactions occur at the QD–ligand interface, and the reduced Cd atom is almost in a zero valence state. However, different sites are predicted for the reactions in which the surface metal atom of the QD instead of the metal atom in the ligands is reduced. As a result, one of the arms of the chelate ligand leaves the QD, which tends to cause damage to its electroluminescent performance. Our findings contribute to a mechanistic understanding of the reduction reactions that occurred at the QD–ligand interface. [ABSTRACT FROM AUTHOR]

We propose an efficient version of ensemble Hartree–Fock/density functional theory to calculate a set of two charge-transfer states for systems with odd-numbers of electrons. The approach does require definitions of donor/acceptor fragments; however, the approach is not very sensitive to such definitions—even in the limit of very strong electronic coupling. The key ansatz is that, by mandating that the vector space spanned by the active orbitals projects equally onto the donor and acceptor fragments, such a constraint eliminates all intra-molecular local excitations and makes it far easier to generate potential energy surfaces that are smooth over a wide region of configuration space. The method is fast, working with only two electron configurations, and should be useful for ab initio non-adiabatic dynamics in the near future. [ABSTRACT FROM AUTHOR]

The dissociative double ionization of CH3I and CH2I2 irradiated with extreme ultraviolet light at hv = 100 eV is investigated by multi-electron–ion coincidence spectroscopy using a magnetic bottle type electron spectrometer. The spin–orbit state-resolved Auger electron spectra for the I 4d core-hole states, (I 4d3/2)−1 and (I 4d5/2)−1, provide clear identifications of electronic states of CH3I2+ and CH2I22+. The dominant ion species produced after the double ionization correlate with the Auger electron energy, showing that different fragmentation pathways are open depending on the electronic states populated by the Auger decay. Theoretical calculations are performed to understand the fragmentation from the doubly charged states and the observed spin–orbit specificity in the Auger decay. [ABSTRACT FROM AUTHOR]

We propose a new formulation of time-dependent coupled cluster with adaptive basis functions and division of the one-particle space into active and secondary subspaces. The formalism is fully bivariational in the sense of a real-valued time-dependent bivariational principle and converges to the complete-active-space solution, a property that is obtained by the use of biorthogonal basis functions. A key and distinguishing feature of the theory is that the active bra and ket functions span the same space by construction. This ensures numerical stability and is achieved by employing a split unitary/non-unitary basis set transformation: the unitary part changes the active space itself, while the non-unitary part transforms the active basis. The formulation covers vibrational as well as electron dynamics. Detailed equations of motion are derived and implemented in the context of vibrational dynamics, and the numerical behavior is studied and compared to related methods. [ABSTRACT FROM AUTHOR]

The photoexcitation dynamics of molecular materials on the 10–100 nm length scale depend on complex interactions between electronic and vibrational degrees of freedom, rendering exact calculations difficult or intractable. The adaptive Hierarchy of Pure States (adHOPS) is a formally exact method that leverages the locality imposed by interactions between thermal environments and electronic excitations to achieve size-invariant scaling calculations for single-excitation processes in systems described by a Frenkel–Holstein Hamiltonian. Here, we extend adHOPS to account for arbitrary couplings between thermal environments and vertical excitation energies, enabling formally exact, size-invariant calculations that involve multiple excitations or states with shared thermal environments. In addition, we introduce a low-temperature correction and an effective integration of the noise to reduce the computational expense of including ultrafast vibrational relaxation in Hierarchy of Pure States (HOPS) simulations. We present these advances in the latest version of the open-source MesoHOPS library and use MesoHOPS to characterize charge separation at a one-dimensional organic heterojunction when both the electron and hole are mobile. [ABSTRACT FROM AUTHOR]

Despite the importance of the one-particle picture provided by the orbital paradigm, a rigorous understanding of the spatial distribution of electrons in molecules is still of paramount importance to chemistry. Considerable progress has been made following the introduction of topological approaches, capable of partitioning space into chemically meaningful regions. They usually provide atomic partitions, for example, through the attraction basins of the electron density in the quantum theory of atoms in molecules (QTAIM) or electron-pair decompositions, as in the case of the electron localization function (ELF). In both cases, the so-called electron distribution functions (EDFs) provide a rich statistical description of the electron distribution in these spatial domains. Here, we take the EDF concept to a new fine-grained limit by calculating EDFs in the QTAIM ∩ ELF intersection domains. As shown in AHn systems based on main group elements, as well as in the CO, NO, and BeO molecules, this approach provides an exquisitely detailed picture of the electron distribution in molecules, allowing for an insightful combination of the distribution of electrons between Lewis entities (such as bonds and lone pairs) and atoms at the same time. Besides mean-field calculations, we also explore the impact of electron correlation through Hartree–Fock (HF), density functional theory (DFT) (B3LYP), and CASSCF calculations. [ABSTRACT FROM AUTHOR]

2D lateral heterostructures possess atomically sharp lateral interfaces, while understanding of their ultrafast photocarrier dynamics from a spatiotemporal viewpoint is rather elusive. In this study, we have investigated the spatiotemporal evolution of photocarrier dynamics across the 1D lateral interface of a WS2–ReS2 2D lateral heterostructure utilizing femtosecond laser pump–probe. The nontrivial band offset across the 1D lateral interface markedly mediates the spatiotemporal photocarrier transfer and transport processes. Subsequently, a hole accumulation region on the WS2 side and an electron accumulation region (1DEG) on the ReS2 side have been spatially identified by correlating ultrafast photocarrier signals. The measured width of the unilateral depletion region is 1360 ± 160 nm. Our work has provided substantial insights into mediated photocarrier dynamics in the 2D lateral heterostructure, which will benefit explorations in 2D interfacial physics and 2D lateral optoelectronic devices. [ABSTRACT FROM AUTHOR]

Within the context of fewest-switch surface hopping (FSSH) dynamics, one often wishes to remove the angular component of the derivative coupling between states J and K . In a previous set of papers, Shu et al. [J. Phys. Chem. Lett. 11, 1135–1140 (2020)] posited one approach for such a removal based on direct projection, while we isolated a second approach by constructing and differentiating a rotationally invariant basis. Unfortunately, neither approach was able to demonstrate a one-electron operator O ̂ whose matrix element J O ̂ K was the angular component of the derivative coupling. Here, we show that a one-electron operator can, in fact, be constructed efficiently in a semi-local fashion. The present results yield physical insight into designing new surface hopping algorithms and are of immediate use for FSSH calculations. [ABSTRACT FROM AUTHOR]

Recently developed locally scaled self-interaction correction (LSIC) is a one-electron SIC method that, when used with a ratio of kinetic energy densities (zσ) as iso-orbital indicator, performs remarkably well for both thermochemical properties as well as for barrier heights overcoming the paradoxical behavior of the well-known Perdew–Zunger self-interaction correction (PZSIC) method. In this work, we examine how well the LSIC method performs for the delocalization error. Our results show that both LSIC and PZSIC methods correctly describe the dissociation of H 2 + and He 2 + but LSIC is overall more accurate than the PZSIC method. Likewise, in the case of the vertical ionization energy of an ensemble of isolated He atoms, the LSIC and PZSIC methods do not exhibit delocalization errors. For the fractional charges, both LSIC and PZSIC significantly reduce the deviation from linearity in the energy vs number of electrons curve, with PZSIC performing superior for C, Ne, and Ar atoms while for Kr they perform similarly. The LSIC performs well at the endpoints (integer occupations) while substantially reducing the deviation. The dissociation of LiF shows both LSIC and PZSIC dissociate into neutral Li and F but only LSIC exhibits charge transfer from Li+ to F− at the expected distance from the experimental data and accurate ab initio data. Overall, both the PZSIC and LSIC methods reduce the delocalization errors substantially. [ABSTRACT FROM AUTHOR]

Computational study of electronic resonances is still a very challenging topic, with the phenomenon of dissociative electron attachment (DEA) being one of the multiple features worth investigating. Recently, we extended the charge stabilization method from energies to properties of conceptual density functional theory and applied this to metastable anionic states of ethene and chlorinated ethene derivatives to study the DEA mechanism present in these compounds. We now present an extension to spatial functions, namely, the electronic Fukui function and the electron localization function. The results of our analysis show that extrapolated spatial functions are relevant and useful for more precise localization of the unbound electron. Furthermore, we report for the first time the combination of the electron localization function with Berlin's binding function for these challenging electronic states. This promising methodology allows for accurate predictions of when and where DEA will happen in the molecules studied and provides more insight into the process. [ABSTRACT FROM AUTHOR]

The photophysics of thiobases—nucleobases in which one or more oxygen atoms are replaced with sulfur atoms— vary greatly depending on the location of sulfonation. Not only are direct dynamics of a neutral thiobase impacted, but also the dynamics of excess electron accommodation. In this work, time-resolved photoelectron spectroscopy is used to measure binary anionic clusters of iodide and 4-thiouracil, I− · 4TU. We investigate charge transfer dynamics driven by excitation at 3.88 eV, corresponding to the lowest ππ* transition of the thiouracil, and at 4.16 eV, near the cluster vertical detachment energy. The photoexcited state dynamics are probed by photodetachment with 1.55 and 3.14 eV pulses. Excitation at 3.88 eV leads to a signal from a valence-bound ion only, indicating a charge accommodation mechanism that does not involve a dipole-bound anion as an intermediate. Excitation at 4.16 eV rapidly gives rise to dipole-bound and valence-bound ion signals, with a second rise in the valence-bound signal corresponding to the decay of the dipole-bound signal. The dynamics associated with the low energy ππ* excitation of 4-thiouracil provide a clear experimental proof for the importance of localized excitation and electron backfilling in halide–nucleobase clusters. [ABSTRACT FROM AUTHOR]

This paper reports the energies and charge and spin distributions of the low-lying excited states in singlet and triplet N2V defects in diamond from direct Δ-SCF calculations based on Gaussian orbitals within the B3LYP, PBE0, and HSE06 functionals. They assign the observed absorption at 2.463 eV, first reported by Davies et al. [Proc. R. Soc. London 351, 245 (1976)], to the excitation of a N(sp3) lone-pair electron in the singlet and triplet states, respectively, with estimates of ∼1.1 eV for that of the unpaired electrons, C(sp3). In both cases, the excited states are predicted to be highly local and strongly excitonic with 81% of the C(sp3) and 87% of the N(sp3) excited charges localized at the three C atoms nearest neighbor (nn) to the excitation sites. Also reported are the higher excited gap states of both the N lone pair and C unpaired electron. Calculated excitation energies of the bonding sp3 hybrids of the C atoms nn to the four inner atoms are close to that of the bulk, which indicates that the N2V defect is largely a local defect. The present results are in broad agreement with those reported by Udvarhelyi et al. [Phys. Rev. B 96, 155211 (2017)] from plane wave HSE06 calculations, notably for the N lone pair excitation energy, for which both predict an energy of ∼2.7 eV but with a difference of ∼0.5 eV for the excitation of the unpaired electron. [ABSTRACT FROM AUTHOR]

Although it plays a critical role in the photophysics and catalysis of lanthanides, spin–orbit coupling of electrons on individual lanthanide atoms in small clusters is not well understood. The major objective of this work is to probe such coupling of the praseodymium (Pr) 4f and 6s electrons in Pr2O2 and Pr2O2+. The approach combines mass-analyzed threshold ionization spectroscopy and spin–orbit multiconfiguration second-order quasi-degenerate perturbation theory. The energies of six ionization transitions are precisely measured; the adiabatic ionization energy of the neutral cluster is 38 045 (5) cm−1. Most of the electronic states involved in these transitions are identified as spin–orbit coupled states consisting of two or more electron spins. The electron configurations of these states are 4f46s2 for the neutral cluster and 4f46s for the singly charged cation, both in planar rhombus-type structures. The spin–orbit splitting due to the coupling of the electrons on the separate Pr atoms is on the order of hundreds of wavenumbers. [ABSTRACT FROM AUTHOR]

The adiabatic connection interaction strength interpolation (ISI)-like method provides a high-level expression for the correlation energy, being, in principle, exact not only in the weak-interaction limit, where it recovers the second-order Görling–Levy perturbation term, but also in the strong-interaction limit that is described by the strictly correlated electron approach. In this work, we construct a genISI functional made accurate for the uniform electron gas, a solid-state physics paradigm that is a very difficult test for ISI-like correlation functionals. We assess the genISI functional for various jellium spheres with the number of electrons Z ≤ 912 and for the non-relativistic noble atoms with Z ≤ 290. For the jellium clusters, the genISI is remarkably accurate, while for the noble atoms, it shows a good performance, similar to other ISI-like methods. Then, the genISI functional can open the path using the ISI-like method in solid-state calculations. [ABSTRACT FROM AUTHOR]

We investigate two limits in open-shell diradical systems: O3, in which the interesting orbitals are in close proximity to one another, and (C21H13)2, where there is a significant spatial separation between the two orbitals. In accord with earlier calculations, we find that standard density-functional approximations do not predict the open-shell character for the former case but uniformly predict the open-shell character for the latter case. We trace the qualitatively incorrect behavior in O3 predicted by these standard density functional approximations to self-interaction error and use the Fermi–Löwdin-orbital-self-interaction-corrected formalism to determine accurate triplet, closed-shell singlet, and open-shell broken-spin-symmetry electronic configurations. Analysis of the resulting many-electron overlap matrices allows us to unambiguously show that the broken-spin-symmetry configurations do not participate in the representation of the Ms = 0 triplet states and allows us to reliably extract the singlet–triplet splitting in O3 by analyzing the energy as a function of Fermi-orbital-descriptor permutations. The results of these analyses predict the percentage of open-shell character in O3, which agrees well with conventional wavefunction-based methods. While these techniques are expected to be required in cases near the Coulson–Fischer point, we find that they will be less necessary in diradical systems with well-separated electrons, such as (C21H13)2. Results based on energies from self-interaction-corrected generalized gradient, local density, and Hartree–Fock approximations and experimental results are in generally good agreement for O3. These results help form the basis for deriving extended Heisenberg-like Hamiltonians that are needed for descriptions of molecular magnets when there are competing low-energy electronic configurations. [ABSTRACT FROM AUTHOR]

The electronic properties of SO2 at the surface of a water cluster were investigated by employing a combination of Born–Oppenheimer molecular dynamics and electron propagator theory (EPT). In our work, we utilized a revised version of the Perdew–Burke–Ernzerhof (PBE) exchange-correlation functional, which incorporates empirical corrections for dispersion interactions in line with a recent study of the air–water interface conducted by Ohto et al. [J. Phys. Chem. Lett. 10(17), 4914–4919 (2019)]. Polarization effects induce a significant broadening of the electron binding energy distribution, as predicted by EPT. This broadening can result in a substantial increase in electron affinity, impacting the chemical reactivity of SO2 at the air–water interface, a topic of significant and recent research interest. We discuss the relationship between electron binding energies (EBEs) and the specific connections of SO2 to water. The results indicate that configurations involving an OS⋯H bond tend to yield higher electron affinities compared to complex formation through S⋯OW bonds. Surprisingly, SO2 molecules not bound to water molecules according to a specific criterion may also exhibit higher electron affinities. This feature can be explained by the role played by the polarization field from water molecules. Our best estimate for the HOMO–LUMO (H–L) gap of SO2 at the surface of a water cluster is 11.6 eV. Very similar H–L gaps are predicted for isolated and micro-solvated SO2. Fukui functions for the gas phase, and the micro-solvated SO2–H2O complex supports the view that the LUMO is predominantly localized on the SO2 moiety. [ABSTRACT FROM AUTHOR]

We numerically study the strong-interaction limit of the exchange–correlation functional for neutral atoms and Bohr atoms as the number of electrons increases. Using a compact representation, we analyze the second-order gradient expansion, comparing it with the one for exchange (weak interaction limit). The two gradient expansions, at strong and weak interaction, turn out to be very similar in magnitude but with opposite signs. We find that the point-charge plus continuum model is surprisingly accurate for the gradient expansion coefficient at strong coupling, while generalized gradient approximations, such as Perdew–Burke–Ernzerhof (PBE) and PBEsol, severely underestimate it. We then use our results to analyze the Lieb–Oxford bound from the point of view of slowly varying densities, clarifying some aspects on the bound at a fixed number of electrons. [ABSTRACT FROM AUTHOR]

A novel implementation of the coupled-cluster singles and doubles (CCSD) approach is presented that is specifically tailored for the treatment of large symmetric systems. It fully exploits Abelian point-group symmetry and the use of the Cholesky decomposition of the two-electron repulsion integrals. In accordance with modern CCSD algorithms, we propose two alternative strategies for the computation of the so-called particle–particle ladder term. The code is driven toward the optimal choice depending on the available hardware resources. As a large-scale application, we computed the frozen-core correlation energy of buckminsterfullerene (C60) with a polarized valence triple-zeta basis set (240 correlated electrons in 1740 orbitals). [ABSTRACT FROM AUTHOR]

In the present work, two quasi-molecular compounds each involving one antiproton and one electron ( p ̄) , H e + − p ̄ and H − p ̄ , are investigated. Using completely relativistic calculations within the finite-basis method adapted to systems with axial symmetry, the adiabatic potential curves are constructed by numerically solving the two-center Dirac equation. The binding energies of electron are obtained as a function of the inter-nuclear distance and compared with the corresponding nonrelativistic values and relativistic leading-order corrections calculated in the framework of other approaches. A semantic analysis of antiproton quasi-molecular ions with compounds containing a proton (p) instead of an antiproton is given. The advantages of the A-DKB method are demonstrated. [ABSTRACT FROM AUTHOR]

A quasi-atomic orbital analysis of the halogen bonded NH3⋯XF complexes (X = F, Cl, Br, and I) is performed to gain insight into the electronic properties associated with these σ-hole interactions. It is shown that significant sharing of electrons between the nitrogen lone pair of the ammonia molecule and the XF molecule occurs, resulting in a weakening of the X–F bond. In addition, the N–X bond shows increasing covalent character as the size of the halogen atom X increases. While the Mulliken outer complex NH3⋯XF appears to be overall the main species, the strength of the covalent interaction of the N–X bond becomes increasingly similar to that of the N–X bond in the [NH3X]+ cation as the size of X increases. [ABSTRACT FROM AUTHOR]

The damage mechanism and exposure tolerance of epoxy resins to fast electrons remain unclear. We quantitatively investigated the effects of electron irradiation on a common epoxy resin by dose-dependent electron energy loss spectroscopy. The results show that sp3 states of nitrogen, oxygen, and their adjacent carbon atoms were converted to sp2 states, forming imine (C=N) and carbonyl (C=O) as the total electron dose increased. The sp3 to sp2 conversion mechanism was proposed. The epoxy resin was very sensitive to fast electrons and the original electronic states were maintained up to a total dose of ∼103e− nm−2 at a low temperature of 103 K. Dose-dependent electron diffraction revealed that the intra- and intermolecular geometries changed below and around the total dose of ∼103e− nm−2. [ABSTRACT FROM AUTHOR]

The adiabatic connection formalism, usually based on the first-order perturbation theory, has been generalized to an arbitrary order. The generalization stems from the observation that the formalism can be derived from a properly arranged Taylor expansion. The second-order theory is developed in detail and applied to the description of two electrons in a parabolic confinement (harmonium). A significant improvement relative to the first-order theory has been obtained. [ABSTRACT FROM AUTHOR]

In a recent study, we developed a kinetic-energy density functional that can be utilized in orbital-free quantum mechanical/molecular mechanical (OF-QM/MM) simulations. The functional includes the nonlocal term constructed from the response function of the reference system of the QM solute. The present work provides a method to combine the OF-QM/MM with a theory of solutions based on the energy representation to compute the solvation free energy of the QM solute in solution. The method is applied to the calculation of the solvation free energy Δμ of a QM water solute in an MM water solvent. It is demonstrated that Δμ is computed as −7.7 kcal/mol, in good agreement with an experimental value of −6.3 kcal/mol. We also develop a theory to map the free energy δμ due to electron density polarization onto the coordinate space of electrons. The free energy density obtained by the free-energy mapping for the QM water clarifies that each hydrogen atom makes a positive contribution (+34.7 kcal/mol) to δμ, and the oxygen atom gives the negative free energy (−71.7 kcal/mol). It is shown that the small polarization free energy −2.4 kcal/mol is generated as a result of the cancellation of these counteracting energies. These analyses are made possible by the OF-QM/MM approach combined with a statistical theory of solutions. [ABSTRACT FROM AUTHOR]

A new generation of diagonal self-energies for the calculation of electron removal energies of molecules and molecular ions that has superseded its predecessors with respect to accuracy, efficiency, and interpretability is extended to include non-diagonal self-energies that permit Dyson orbitals to be expressed as linear combinations of canonical Hartree–Fock orbitals. In addition, an improved algorithm for renormalized methods eliminates the convergence difficulties encountered in the first studies of the new, diagonal self-energies. A dataset of outer-valence, vertical ionization energies with almost full-configuration-interaction quality serves as a standard of comparison in numerical tests. The new non-diagonal, renormalized methods are slightly more accurate than their diagonal counterparts, with mean absolute errors between 0.10 and 0.06 eV for outer-valence final states. This advantage is procured at the cost of an increase in the scaling of arithmetic bottlenecks that accompany the inclusion of non-diagonal self-energy terms. The new, non-diagonal, renormalized self-energies are also more accurate and efficient than their non-diagonal predecessors. [ABSTRACT FROM AUTHOR]

The fragmentation of fulminic acid, HCNO, after excitation and ionization of core electrons was investigated using Auger-electron–photoion coincidence spectroscopy. A considerable degree of site-selectivity is observed. Ionization of the carbon and oxygen 1s electron leads to around 70% CH+ + NO+, while ionization at the central N-atom produces only 37% CH+ + NO+, but preferentially forms O+ + HCN+ and O+ + CN+. The mass-selected Auger-electron spectra show that these fragments are associated with higher binding energy final states. Furthermore, ionization of the C 1s electron leads to a higher propensity for C–H bond fission compared to O 1s ionization. Following resonant Auger–Meitner decay after 1s → 3π excitation, 12 different ionic products are formed. At the C 1s edge, the parent ion HCNO+ is significantly more stable compared to the other two edges, which we also attribute to the higher contribution of final states with low binding energies in the C 1s resonant Auger electron spectra. [ABSTRACT FROM AUTHOR]

We report the first ever measurements of the thermalization length of low-energy electrons injected into solid para-hydrogen at a temperature T ≈ 2.8 K. The use of the pulsed Townsend photoinjection technique has allowed us to investigate the behavior of quasi-free electrons rather than of massive, slow negative charges, as reported in all previous literature. We have found an average thermalization length ⟨z0⟩ = 26.1 nm, which is three to five times longer than that in liquid helium at the same temperature. [ABSTRACT FROM AUTHOR]

When a vacancy is created in an inner-valence orbital of a dimer of atoms or molecules, the resulting species can undergo interatomic/intermolecular Coulombic decay (ICD): the hole is filled through a relaxation process that leads to a doubly ionized cluster with two positively charged atoms or molecules. Since they are subject to electronic decay, inner-valence ionized states are not bound states but electronic resonances whose transient nature can only be described with special quantum-chemical methods. In this work, we explore the capacity of equation-of-motion coupled-cluster theory with two techniques from non-Hermitian quantum mechanics, complex basis functions and Feshbach–Fano projection with a plane wave description of the outgoing electron, to describe ICD. To this end, we compute the decay rates of several dimers: Ne2, NeAr, NeMg, and (HF)2, among which the energy of the outgoing electron varies between 0.3 and 16 eV. We observe that both methods deliver better results when the outgoing electron is fast, but the characteristic R−6 distance dependence of the ICD width is captured much better with complex basis functions. [ABSTRACT FROM AUTHOR]

In the previous study Wang and Dou [J. Chem. Phys. 158, 224109 (2023)], we have derived a Floquet classical master equation (FCME) to treat nonadiabatic dynamics near metal surfaces under Floquet engineering. We have also proposed a trajectory surface hopping algorithm to solve the FCME. In this study, we map the FCME into a Floquet Fokker–Planck equation in the limit of fast Floquet driving and fast electron motion as compared to nuclear motion. The Fokker–Planck equation is then being solved using Langevin dynamics with explicit friction and random force from the nonadiabatic effects of hybridized electrons and Floquet states. We benchmark the Floquet electronic friction dynamics against Floquet quantum master equation and Floquet surface hopping. We find that Floquet driving results in a violation of the second fluctuation–dissipation theorem, which further gives rise to heating effects. [ABSTRACT FROM AUTHOR]

Controlling matter at the level of electrons using ultrafast laser sources represents an important challenge for science and technology. Recently, we introduced a general laser control scheme (the Stark control of electrons at interfaces or SCELI) based on the Stark effect that uses the subcycle structure of light to manipulate electron dynamics at semiconductor interfaces [A. Garzón-Ramírez and I. Franco, Phys. Rev. B 98, 121305 (2018)]. Here, we demonstrate that SCELI is also of general applicability in molecule–semiconductor interfaces. We do so by following the quantum dynamics induced by non-resonant few-cycle laser pulses of intermediate intensity (non-perturbative but non-ionizing) across model molecule–semiconductor interfaces of varying level alignments. We show that SCELI induces interfacial charge transfer regardless of the energy level alignment of the interface and even in situations where charge exchange is forbidden via resonant photoexcitation. We further show that the SCELI rate of charge transfer is faster than those offered by resonant photoexcitation routes as it is controlled by the subcycle structure of light. The results underscore the general applicability of SCELI to manipulate electron dynamics at interfaces on ultrafast timescales. [ABSTRACT FROM AUTHOR]

The rate constant for electron attachment to Mo(CO)6 was determined to be ka = 2.4 ± 0.6 × 10−7 cm3 s−1 at 297 K in a flowing-afterglow Langmuir-probe experiment. The sole anion product is Mo(CO)5−. A small decline in ka was observed up to 450 K, and decomposition was apparent at higher temperatures. The charge transfer reaction of Ar+ with Mo(CO)6 is exothermic by 7.59 ± 0.03 eV, which appears to be sufficient to remove the first 5 ligands from Mo(CO)6+. [ABSTRACT FROM AUTHOR]

We here report theoretical triply differential cross sections (TDCS) for 250 eV electron and positron impact ionization of the methane molecule calculated within the second-order distorted-wave Born approximation (DWBA2) for various momentum transfer conditions. The experimental data taken from Işık et al. [J. Phys. B: At., Mol. Opt. Phys. 49, 065203 (2016)] will be compared with the current theoretical predictions as well as molecular three body distorted wave (M3DW) approximation and generalized Sturmian function (GSF) theoretical models in a non-coplanar geometry. In the low analyzer scattering plane, the results obtained within the DWBA2 theory show better agreement with the experimental results compared to the GSF results. The M3DW results also exhibit agreement with the experimental results, in particular in the perpendicular plane geometry. Furthermore, significant differences between electron and positron TDCS were observed. [ABSTRACT FROM AUTHOR]

The rate constant for electron attachment to Mo(CO)6 was determined to be ka = 2.4 ± 0.6 × 10−7 cm3 s−1 at 297 K in a flowing-afterglow Langmuir-probe experiment. The sole anion product is Mo(CO)5−. A small decline in ka was observed up to 450 K, and decomposition was apparent at higher temperatures. The charge transfer reaction of Ar+ with Mo(CO)6 is exothermic by 7.59 ± 0.03 eV, which appears to be sufficient to remove the first 5 ligands from Mo(CO)6+. [ABSTRACT FROM AUTHOR]

Photoionization spectroscopy and mass spectrometry of doped helium (He) nanodroplets rely on the ability to efficiently detect ions and/or electrons. Using a commercial quadrupole mass spectrometer and a photoelectron–photoion coincidence spectrometer, we systematically measure yields of ions and electrons created in pure and doped He nanodroplets in a wide size range and in two ionization regimes—direct ionization and secondary ionization after resonant photoexcitation of the droplets. For two different types of dopants (oxygen molecules, O2, and lithium atoms, Li), we infer the optimal droplet size to maximize the yield of ejected ions. When dopants are ionized by charge-transfer to photoionized He nanodroplets, the highest yield of O2 and Li ions is detected for a mean size of ∼ 5 × 1 0 4 He atoms per nanodroplet. When dopants are Penning ionized via photoexcitation of the He droplets, the highest yield of O2 and Li ions is detected for ∼ 1 0 3 and ∼ 1 0 5 He atoms per droplet, respectively. At optimum droplet sizes, the detection efficiency of dopant ions in proportion to the number of primary photoabsorption events is up to 20% for charge-transfer ionization of O2 and 2% for Li, whereas for Penning ionization it is 1% for O2 and 4% for Li. Our results are instrumental in determining optimal conditions for mass spectrometric studies and photoionization spectroscopy of molecules and complexes isolated in He nanodroplets. [ABSTRACT FROM AUTHOR]

We report a Kohn–Sham density functional theory calculation of a system with more than 200 000 atoms and 800 000 electrons using a real-space high-order finite-difference method to investigate the electronic structure of large spherical silicon nanoclusters. Our system of choice was a 20 nm large spherical nanocluster with 202 617 silicon atoms and 13 836 hydrogen atoms used to passivate the dangling surface bonds. To speed up the convergence of the eigenspace, we utilized Chebyshev-filtered subspace iteration, and for sparse matrix–vector multiplications, we used blockwise Hilbert space-filling curves, implemented in the PARSEC code. For this calculation, we also replaced our orthonormalization + Rayleigh–Ritz step with a generalized eigenvalue problem step. We utilized all of the 8192 nodes (458 752 processors) on the Frontera machine at the Texas Advanced Computing Center. We achieved two Chebyshev-filtered subspace iterations, yielding a good approximation of the electronic density of states. Our work pushes the limits on the capabilities of the current electronic structure solvers to nearly 106 electrons and demonstrates the potential of the real-space approach to efficiently parallelize large calculations on modern high-performance computing platforms. [ABSTRACT FROM AUTHOR]

The metastability and the dissociation processes of the OCS3+ states produced by the S 2p double Auger decay of OCS are investigated by multi-electron-ion coincidence spectroscopy using a magnetic bottle electron spectrometer. The spectra of the OCS3+ states filtered for the production of individual ions are derived by four-fold (or five-fold) coincidence among three electrons and a product ion (or two ions). The ground OCS3+ state is confirmed to be metastable in the 10 µs regime. The OCS3+ states relevant to the individual channels of two- and three-body dissociations are clarified. [ABSTRACT FROM AUTHOR]

Recombination of vibrationally cold N 2 + ions with electrons was studied in the temperature range of 140–250 K. A cryogenic stationary afterglow apparatus equipped with cavity ring-down spectrometer and microwave diagnostics was utilized to probe in situ the time evolutions of number densities of particular rotational and vibrational states of N 2 + ions and of electrons. The obtained value of the recombination rate coefficient for the recombination of the vibrational ground state of N 2 + with electrons is αv=0 = (2.95 ± 0.50) × 10−7(300/T)(0.28±0.07) cm3 s−1, while that for the first vibrationally excited state was inferred as αv=1 = (4 ± 4) × 10−8 cm3 s−1 at 250 K. [ABSTRACT FROM AUTHOR]

While the dissociative recombination (DR) of ground-state molecular ions with low-energy free electrons is generally known to be exothermic, it has been predicted to be endothermic for a class of transition-metal oxide ions. To understand this unusual case, the electron recombination of titanium oxide ions (TiO+) with electrons has been experimentally investigated using the Cryogenic Storage Ring. In its low radiation field, the TiO+ ions relax internally to low rotational excitation (≲100 K). Under controlled collision energies down to ∼ 2 meV within the merged electron and ion beam configuration, fragment imaging has been applied to determine the kinetic energy released to Ti and O neutral reaction products. Detailed analysis of the fragment imaging data considering the reactant and product excitation channels reveals an endothermicity for the TiO+ dissociative electron recombination of (+4 ± 10) meV. This result improves the accuracy of the energy balance by a factor of 7 compared to that found indirectly from hitherto known molecular properties. Conversely, the present endothermicity yields improved dissociation energy values for D0(TiO) = (6.824 ± 0.010) eV and D0(TiO+) = (6.832 ± 0.010) eV. All thermochemistry values were compared to new coupled-cluster calculations and found to be in good agreement. Moreover, absolute rate coefficients for the electron recombination of rotationally relaxed ions have been measured, yielding an upper limit of 1 × 10−7 cm3 s−1 for typical conditions of cold astrophysical media. Strong variation of the DR rate with the TiO+ internal excitation is predicted. Furthermore, potential energy curves for TiO+ and TiO have been calculated using a multi-reference configuration interaction method to constrain quantum-dynamical paths driving the observed TiO+ electron recombination. [ABSTRACT FROM AUTHOR]

Characterization of the inner-shell decay processes in molecules containing heavy elements is key to understanding x-ray damage of molecules and materials and for medical applications with Auger-electron-emitting radionuclides. The 1s hole states of heavy atoms can be produced by absorption of tunable x rays and the resulting vacancy decays characterized by recording emitted photons, electrons, and ions. The 1s hole states in heavy elements have large x-ray fluorescence yields that transfer the hole to intermediate electron shells that then decay by sequential Auger-electron transitions that increase the ion's charge state until the final state is reached. In molecules, the charge is spread across the atomic sites, resulting in dissociation to energetic atomic ions. We have used x-ray/ion coincidence spectroscopy to measure charge states and energies of Iq+ and Brq′+ atomic ions following 1s ionization at the I and Br K-edges of IBr. We present the charge states and kinetic energies of the two correlated fragment ions associated with core-excited states produced during the various steps of the cascades. To understand the dynamics leading to the ion data, we develop a computational model that combines Monte-Carlo/Molecular-Dynamics (MC/MD) simulations with a classical over-the-barrier model to track inner-shell cascades and redistribution of electrons in valence orbitals and nuclear motion of fragments. [ABSTRACT FROM AUTHOR]