Search

Accurate theoretical prediction of the band offsets at interfaces of semiconductor heterostructures can often be quite challenging. Although density functional theory has been reasonably successful to carry out such calculations, efficient, accurate semilocal functionals are desirable to reduce the computational cost. In general, the semilocal functionals based on the generalized gradient approximation (GGA) significantly underestimate the bulk bandgaps. This, in turn, results in inaccurate estimates of the band offsets at the heterointerfaces. In this paper, we investigate the performance of several advanced meta-GGA functionals in the computational prediction of band offsets at semiconductor heterojunctions. In particular, we investigate the performance of r2SCAN (two times revised strongly constrained and appropriately normed functional), rMGGAC (revised semilocal functional based on cuspless hydrogen model and Pauli kinetic energy density functional), mTASK (modified Aschebrock and Kümmel meta-GGA functional), and local modified Becke–Johnson exchange-correlation functionals. Our results strongly suggest that these meta-GGA functionals for supercell calculations perform quite well, especially, when compared to computationally more demanding GW calculations. We also present band offsets calculated using ionization potentials and electron affinities, as well as band alignment via the branch point energies. Overall, our study shows that the aforementioned meta-GGA functionals can be used within the density functional theory framework to estimate the band offsets in semiconductor heterostructures with predictive accuracy. [ABSTRACT FROM AUTHOR]

The revised, regularized Tao–Mo (rregTM) exchange-correlation density functional approximation (DFA) [A. Patra, S. Jana, and P. Samal, J. Chem. Phys. 153, 184112 (2020) and Jana et al., J. Chem. Phys. 155, 024103 (2021)] resolves the order-of-limits problem in the original TM formulation while preserving its valuable essential behaviors. Those include performance on standard thermochemistry and solid data sets that is competitive with that of the most widely explored meta-generalized-gradient-approximation DFAs (SCAN and r2SCAN) while also providing superior performance on elemental solid magnetization. Puzzlingly however, rregTM proved to be intractable for de-orbitalization via the approach of Mejía-Rodríguez and Trickey [Phys. Rev. A 96, 052512 (2017)]. We report investigation that leads to diagnosis of how the regularization in rregTM of the z indicator functions (z = the ratio of the von-Weizsäcker and Kohn–Sham kinetic energy densities) leads to non-physical behavior. We propose a simpler regularization that eliminates those oddities and that can be calibrated to reproduce the good error patterns of rregTM. We denote this version as simplified, regularized Tao–Mo, sregTM. We also show that it is unnecessary to use rregTM correlation with sregTM exchange: Perdew–Burke–Ernzerhof correlation is sufficient. The subsequent paper shows how sregTM enables some progress on de-orbitalization. [ABSTRACT FROM AUTHOR]

The Pauli kinetic energy functional and its functional derivative, termed Pauli potential, play a crucial role in the successful implementation of orbital-free density functional theory for electronic structure calculations. However, the exact forms of these two quantities are not known. Therefore, perforce, one employs the approximate forms for the Pauli functional or Pauli potential for performing orbital-free density functional calculations. In the present study, we developed a feed-forward neural network-based representation for the Pauli potential using a 1-dimensional (1-D) model system. We expanded density in terms of basis functions, and the coefficients of the expansion were used as input to a feed-forward neural network. Using the neural network-based representation of the Pauli potential, we calculated the ground-state densities of the 1-D model system by solving the Euler equation. We calculated the Pauli kinetic energy using the neural network-based Pauli potential employing the exact relation between the Pauli kinetic energy functional and the potential. The sum of the neural network-based Pauli kinetic energy and the von Weizsäcker kinetic energy resulted in an accurate estimation of the total kinetic energy. The approach presented in this paper can be employed for the calculation of Pauli potential and Pauli kinetic energy, obviating the need for a functional derivative. The present study is an important step in the advancement of application of machine learning-based techniques toward the orbital-free density functional theory-based methods. [ABSTRACT FROM AUTHOR]

When an activated complex, as defined in transition state theory (TST), has a polyhedral shape, its kinetic energy is found to be diagonal in a system of spherical polar coordinates. If, in addition, the polyhedron is characterized by a high symmetry, then its dynamics considerably simplifies. An application of this approach to the most symmetrical TS known to date, i.e., that which controls the Cl− + CH3Cl → ClCH3 + Cl− SN2 nucleophilic substitution, is presented and an analytical expression of its potential energy surface is provided. In a substantial range around the saddle point, approximate equations of motion for the two components of the reaction coordinate, i.e., the antisymmetrical stretching motion of the ClCCl core and the wagging motion of the hydrogen triad, can be derived in an analytical form. During an extensive period of time, the main component of the reaction coordinate is governed by an unexpectedly simple equation of motion that depends on a single initial condition, irrespective of the other ones and of the internal energy. Reactive trajectories are observed to form a perfectly collimated bundle characterized by undetectable dispersion, thereby giving a spectacular example of regular dynamics in an anharmonic potential. Regularity and collimation are brought about by local symmetry, which is a widespread feature of potential energy surfaces. Anharmonicity is observed to influence the dynamics only at a late stage. As energy increases, trajectories tend to fan out and to deviate from the analytical equation. For the wagging motion, chaos sets in at much lower energies. [ABSTRACT FROM AUTHOR]

We investigate the electron emission from 3D chiral silver alloy nanohelices initiated by femtosecond laser pulses with a central photon energy of hν = 1.65 eV, well below the work function of the material. We find hot but thermally distributed electron spectra and a strong anisotropy in the electron yield with left- and right-circularly polarized light excitations, which invert in sign between left- and right-handed helices. We analyze the kinetic energy distribution and discuss the role of effective temperatures. Measurements of the reflectance and simulations of the absorbance of the helices based on retarded field calculations are compared to the anisotropy in photoemission. We find a significant enhancement of the anisotropy in the electron emission in comparison to the optical absorption. Neither simple thermionic nor a multiphoton photoemission can explain the experimentally observed asymmetries. Single photon deep-UV photoemission from these helices together with a change of the work function suggests a contribution of the chirally induced spin selectivity effect to the observed asymmetries. [ABSTRACT FROM AUTHOR]

We propose a stochastic Schrödinger equation in which the momentum is coupled to a white Gaussian noise. In the stochastic representation, the kinetic energy representing the self-interaction of momentum is reduced to a linear term of momentum. As such, the quantum evolution operator factorizes into two contributions due to the momentum and the potential, respectively. The exact quantum propagator thereby becomes an expectation of the stochastic one in which the amplitude results from the potential with a fluctuating position: The particle moves with a constant velocity, subjected to a complex Brownian motion. We demonstrate that the stochastic Schrödinger equation can be feasibly used to derive the quantum propagators for the linear potential and the harmonic oscillator system. Novel semiclassical and other approximations may be developed from the new representation of quantum mechanics. [ABSTRACT FROM AUTHOR]

Nucleation, the birth of a stable cluster from a disorder, is inherently stochastic. Yet up to date, there are no quantitative studies on NaCl nucleation that accounts for its stochastic nature. Here, we report the first stochastic treatment of NaCl-water nucleation kinetics. Using a recently developed microfluidic system and evaporation model, our measured interfacial energies extracted from a modified Poisson distribution of nucleation time show an excellent agreement with theoretical predictions. Furthermore, analysis of nucleation parameters in 0.5, 1.5, and 5.5 pl microdroplets reveals an interesting interplay between confinement effects and shifting of nucleation mechanisms. Overall, our findings highlight the need to treat nucleation stochastically rather than deterministically to bridge the gap between theory and experiment. [ABSTRACT FROM AUTHOR]

At the limit of an infinite confinement strength ω, the ground state of a system that comprises two fermions or bosons in harmonic confinement interacting through the Fermi–Huang pseudopotential remains strongly correlated. A detailed analysis of the one-particle description of this "contactium" reveals several peculiarities that are not encountered in conventional model systems (such as the two-electron harmonium atom, ballium, and spherium) involving Coulombic interparticle interactions. First of all, none of the natural orbitals (NOs) { ψ n (ω ; r) } of the contactium is unoccupied, which implies nonzero collective occupancies for all the angular momenta. Second, the NOs and their non-ascendingly ordered occupation numbers { ν n } turn out to be related to the eigenfunctions and eigenvalues of a zero-energy Schrödinger equation with an attractive Gaussian potential. This observation enables the derivation of their properties, such as the n − 4 / 3 asymptotic decay of ν n at the n → ∞ limit (which differs from that of n − 8 / 3 in the Coulombic systems), the independence of the confinement energy v n = ⟨ ψ n (ω ; r) | 1 2 ω 2 r 2 | ψ n (ω ; r) ⟩ of n , and the n − 2 / 3 asymptotic decay of the respective contribution ν n t n to the kinetic energy. Upon suitable scaling, the weakly occupied NOs of the contactium turn out to be virtually identical to those of the two-electron harmonium atom at the ω → ∞ limit, despite the entirely different interparticle interactions in these systems. [ABSTRACT FROM AUTHOR]

We study the linear energy response of the uniform electron gas to an external harmonic perturbation with a focus on resolving different contributions to the total energy. This has been achieved by carrying out highly accurate ab initio path integral Monte Carlo (PIMC) calculations for a variety of densities and temperatures. We report a number of physical insights into effects such as screening and the relative importance of kinetic and potential energies for different wave numbers. A particularly interesting finding is obtained from the observed non-monotonic behavior of the induced change in the interaction energy, which becomes negative for intermediate wave numbers. This effect is strongly dependent on the coupling strength and constitutes further direct evidence for the spatial alignment of electrons introduced in earlier works [T. Dornheim et al., Commun. Phys. 5, 304 (2022)]. The observed quadratic dependence on the perturbation amplitude in the limit of weak perturbations and the quartic dependence of perturbation amplitude corrections are consistent with linear and nonlinear versions of the density stiffness theorem. All PIMC simulation results are freely available online and can be used to benchmark new methods or as input for other calculations. [ABSTRACT FROM AUTHOR]

We study the electron-induced ion-pair dissociation to gas-phase oxygen molecules using a state-of-the-art velocity-map ion-imaging technique. The analysis is entirely based on the conical time-gated wedge-shaped velocity slice images of O−/O2 nascent anionic fragments, and the resulting observations are in favor of Van Brunt et al.'s report [R. J. Van Brunt and L. J. Kieffer, J. Chem. Phys. 60, 3057 (1974)]. A new image reconstruction method, Jacobian over parallel slicing, is introduced to overcome the drawback of ion exaggeration in determining the kinetic energy distribution from the time-gated parallel slicing technique, which offers an alternative approach to the wedge slicing method. Most importantly, the role of the quintet-heavy Rydberg state has been drawn out to the complex ion-pair formalism. The extracted kinetic energy and angular distributions from the wedge slice images reveal a high momentum transfer during the ion-pair dissociation process, which could be the finest rationale to observe the breakdown of dipole Born approximation driven by multipole moment associated with the incident electron beam. Three distinct dissociative momentum bands have been precisely identified for O− dissociation. However, radiationless Rydberg's predissociation continuum (≥15%) has become an inherent character of electron-induced ion-pair dissociation, which could be dealt with using the beyond Born–Oppenheimer treatment. The incoherent sum of Σ and Π symmetric-associated ion-pair final states has been precisely identified by modeling the angular distribution of O−/O2 for each of the kinetic energy bands. A negligibly small amount of forward–backward asymmetry is observed in the angular distribution of O−/O2, which might be explained by the dissociative state-specific quantum coherence mechanism as reported [Krishnakumar et al., Nat. Phys. 14, 149 (2018); Kumar et al., arXiv:2206.15024 (2022)] by Prabhudesai et al. [ABSTRACT FROM AUTHOR]

We study the rotational dynamics induced by the recoil effect in diatomic molecules using time-resolved two-color x-ray pump-probe spectroscopy. A short pump x-ray pulse ionizes a valence electron inducing the molecular rotational wave packet, whereas the second time-delayed x-ray pulse probes the dynamics. An accurate theoretical description is used for analytical discussions and numerical simulations. Our main attention is paid to the following two interference effects that influence the recoil-induced dynamics: (i) Cohen–Fano (CF) two-center interference between partial ionization channels in diatomics and (ii) interference between the recoil-excited rotational levels manifesting as the rotational revival structures in the time-dependent absorption of the probe pulse. The time-dependent x-ray absorption is computed for the heteronuclear CO and homonuclear N2 molecules as showcases. It is found that the effect of CF interference is comparable with the contribution from independent partial ionization channels, especially for the low photoelectron kinetic energy case. The amplitude of the recoil-induced revival structures for the individual ionization decreases monotonously with a decrease in the photoelectron energy, whereas the amplitude of the CF contribution remains sufficient even at the photoelectron kinetic energy below 1 eV. The profile and intensity of the CF interference depend on the phase difference between the individual ionization channels related to the parity of the molecular orbital emitting the photoelectron. This phenomenon provides a sensitive tool for the symmetry analysis of molecular orbitals. [ABSTRACT FROM AUTHOR]

The three-body breakup of [ C 2 H 2 ] 3 + formed in collision with Xe9+ moving at 0.5 atomic units of velocity is studied by using recoil ion momentum spectroscopy. Three-body breakup channels leading to (H+, C+, CH+) and (H+, H+, C 2 + ) fragments are observed in the experiment and their kinetic energy release is measured. The breakup into (H+, C+, CH+) occurs via concerted and sequential modes, whereas the breakup into (H+, H+, C 2 + ) shows only the concerted mode. By collecting events coming exclusively from the sequential breakup leading to (H+, C+, CH+), we have determined the kinetic energy release for the unimolecular fragmentation of the molecular intermediate, [ C 2 H ] 2 + . By using ab initio calculations, the potential energy surface for the lowest electronic state of [ C 2 H ] 2 + is generated, which shows the existence of a metastable state with two possible dissociation pathways. A discussion on the agreement between our experimental results and these ab initio calculations is presented. [ABSTRACT FROM AUTHOR]

Harriman showed that within finite basis sets of one-electron functions that form linearly independent products (LIP), differential and integral operators can be represented exactly and unambiguously by multiplicative (local) potentials. Although almost no standard basis sets of quantum chemistry form LIPs in a numerical sense, occupied self-consistent field (SCF) orbitals routinely do so. Using minimal LIP basis sets of occupied SCF orbitals, we construct multiplicative potentials for electronic kinetic energy and exact exchange that reproduce the Hartree–Fock and Kohn–Sham Hamiltonian matrices and electron densities for atoms and molecules. The results highlight fundamental differences between local and nonlocal operators and suggest a practical possibility of developing exact kinetic energy functionals within finite basis sets by using effective local potentials. [ABSTRACT FROM AUTHOR]

Machine learning techniques have received growing attention as an alternative strategy for developing general-purpose density functional approximations, augmenting the historically successful approach of human-designed functionals derived to obey mathematical constraints known for the exact exchange-correlation functional. More recently, efforts have been made to reconcile the two techniques, integrating machine learning and exact-constraint satisfaction. We continue this integrated approach, designing a deep neural network that exploits the exact constraint and appropriate norm philosophy to de-orbitalize the strongly constrained and appropriately normed (SCAN) functional. The deep neural network is trained to replicate the SCAN functional from only electron density and local derivative information, avoiding the use of the orbital-dependent kinetic energy density. The performance and transferability of the machine-learned functional are demonstrated for molecular and periodic systems. [ABSTRACT FROM AUTHOR]

The Hartree–Fock (HF) approximation has been an important tool for quantum-chemical calculations since its earliest appearance in the late 1920s and remains the starting point of most single-reference methods in use today. Intuition suggests that the HF kinetic energy should not exceed the exact kinetic energy; but no proof of this conjecture exists, despite a near century of development. Beginning from a generalized virial theorem derived from scaling considerations, we derive a general expression for the kinetic energy difference that applies to all systems. For any atom or ion, this trivially reduces to the well-known result that the total energy is the negative of the kinetic energy and, since correlation energies are never positive, proves the conjecture in this case. Similar considerations apply to molecules at their equilibrium bond lengths. We use highly precise calculations on Hooke's atom (two electrons in a parabolic well) to test the conjecture in a nontrivial case and to parameterize the difference between density functional and HF quantities, but find no violations of the conjecture. [ABSTRACT FROM AUTHOR]

O2 has a spin triplet multiplicity in its ground state, while the effect of its electron spin on O2-surface interaction is not well understood. In this study, the spin and/or alignment effects in O2 chemisorption on surfaces of Fe(110), Ni(111), and Co(0001) films grown on W(110) have been investigated with the use of a single spin-rotational state-selected O2 beam. The results indicate that the spin effects for the Fe and Ni films are similar in that the initial sticking probability (S0) of O2 is higher when the spins of O2 and the film are antiparallel, and the spin dependence in S0, which amounts to 30%–40% at thermal energy, decays with increasing the O2 kinetic energy (E0). In the case of the Fe/O2 system, however, considerable spin dependence was found to remain even at E0 > 0.2 eV and on the oxidized surface. It has also been shown that the barrier for O2 chemisorption increases in the order of Fe(110) < Ni(111) < Co(0001), and the difference in the observed alignment effect among the samples can be understood based on the difference in the barrier. [ABSTRACT FROM AUTHOR]

The asymmetric Hubbard dimer is a model that allows for explicit expressions of the Hartree–Fock (HF) and Kohn–Sham (KS) states as analytical functions of the external potential, Δv, and of the interaction strength, U. We use this unique circumstance to establish a rigorous comparison between the individual contributions to the correlation energies stemming from the two theories in the {U, Δv} parameter space. Within this analysis of the Hubbard dimer, we observe a change in the sign of the HF kinetic correlation energy, compare the indirect repulsion energies, and derive an expression for the "traditional" correlation energy, i.e., the one that corrects the HF estimate, in a pure site-occupation function theory spirit [Eq. (45)]. Next, we test the performances of the Liu–Burke and the Seidl–Perdew–Levy functionals, which model the correlation energy based on its weak- and strong-interaction limit expansions and can be used for both the traditional and the KS correlation energies. Our results show that, in the Hubbard dimer setting, they typically work better for the HF reference, despite having been originally devised for KS. These conclusions are somewhat in line with prior assessments of these functionals on various chemical datasets. However, the Hubbard dimer model allows us to show the extent of the error that may occur in using the strong-interaction ingredient for the KS reference in place of the one for the HF reference, as has been carried out in most of the prior assessments. [ABSTRACT FROM AUTHOR]

The semilocal form of an exchange hole is highly useful in developing non-local range-separated hybrid density functionals for finite and extended systems. The way to construct the conventional exact exchange hole model is based on either the Taylor series expansion or the reverse engineering technique from the corresponding exchange energy functional. Although the latter technique is quite popular in the context of generalized gradient approximation (GGA) functionals, the same for the meta-GGA functionals is not so much explored. Thus, in this study, we propose a reverse-engineered semilocal exchange hole of a meta-GGA functional, which only depends on the meta-GGA ingredient α (also known as the Pauli kinetic energy enhancement factor). The model is subsequently used to design the short-range-separated meta-GGA hybrid density functional. We show that the present method can be successfully applied for several challenging problems in the context of solids, especially for which the GGA based hybrid fails drastically. This assessment proves that the present functional is quite useful for materials sciences. Finally, we also use this method for several molecular test cases, where the results are also as comparative as its base semilocal functional. [ABSTRACT FROM AUTHOR]

The prediction of the thermodynamic and kinetic properties of chemical reactions is increasingly being addressed by machine-learning (ML) methods, such as artificial neural networks (ANNs). While a number of recent studies have reported success in predicting chemical reaction activation energies, less attention has been focused on how the accuracy of ML predictions filters through to predictions of macroscopic observables. Here, we consider the impact of the uncertainty associated with ML prediction of activation energies on observable properties of chemical reaction networks, as given by microkinetics simulations based on ML-predicted reaction rates. After training an ANN to predict activation energies, given standard molecular descriptors for reactants and products alone, we performed microkinetics simulations of three different prototypical reaction networks: formamide decomposition, aldol reactions, and decomposition of 3-hydroperoxypropanal. We find that the kinetic modeling predictions can be in excellent agreement with corresponding simulations performed with ab initio calculations, but this is dependent on the inherent energetic landscape of the networks. We use these simulations to suggest some guidelines for when ML-based activation energies can be reliable and when one should take more care in applications to kinetics modeling. [ABSTRACT FROM AUTHOR]

The kinetic energy dependences of the reactions of Pt+ (2D5/2) with SO2 were studied using a guided ion beam tandem mass spectrometer and theory. The observed cationic products are PtO+ and PtSO+, with small amounts of PtS+, all formed in endothermic reactions. Modeling the kinetic energy dependent product cross sections allows determination of the product bond dissociation energies (BDEs): D0(Pt+–O) = 3.14 ± 0.11 eV, D0(Pt+–S) = 3.68 ± 0.31 eV, and D0(Pt+–SO) = 3.03 ± 0.12 eV. The oxide BDE agrees well with more precise literature values, whereas the latter two results are the first such measurements. Quantum mechanical calculations were performed for PtO+, PtS+, PtO2+, and PtSO+ at the B3LYP and coupled-cluster with single, double, and perturbative triple [CCSD(T)] levels of theory using the def2-XZVPPD (X = T, Q) and aug-cc-pVXZ (X = T, Q, 5) basis sets and complete basis set extrapolations. These theoretical BDEs agree well with the experimental values. After including empirical spin–orbit corrections, the product ground states are determined as PtO+ (4Σ3/2), PtS+ (4Σ3/2), PtO2+ (2Σg+), and PtSO+ (2A′). Potential energy profiles including intermediates and transition states for each reaction were also calculated at the B3LYP/def2-TZVPPD level. Periodic trends in the thermochemistry of the group 9 metal chalcogenide cations are compared, and the formation of PtO+ from the Pt+ + SO2 reaction is compared with those from the Pt+ + O2, CO2, CO, and NO reactions. [ABSTRACT FROM AUTHOR]

A neural network-assisted molecular dynamics method is developed to reduce the computational cost of open boundary simulations. Particle influxes and neural network-derived forces are applied at the boundaries of an open domain consisting of explicitly modeled Lennard-Jones atoms in order to represent the effects of the unmodeled surrounding fluid. Canonical ensemble simulations with periodic boundaries are used to train the neural network and to sample boundary fluxes. The method, as implemented in the LAMMPS, yields temperature, kinetic energy, potential energy, and pressure values within 2.5% of those calculated using periodic molecular dynamics and runs two orders of magnitude faster than a comparable grand canonical molecular dynamics system. [ABSTRACT FROM AUTHOR]

A model based on the finite-basis representation of a vibrational Hamiltonian expressed in internal coordinates is developed. The model relies on a many-mode, low-order expansion of both the kinetic energy operator and the potential energy surface (PES). Polyad truncations and energy ceilings are used to control the size of the vibrational basis to facilitate accurate computations of the OH stretch and HOH bend intramolecular transitions of the water dimer ( H 2 16 O)2. Advantages and potential pitfalls of the applied approximations are highlighted. The importance of choices related to the treatment of the kinetic energy operator in reduced-dimensional calculations and the accuracy of different water dimer PESs are discussed. A range of different reduced-dimensional computations are performed to investigate the wavenumber shifts in the intramolecular transitions caused by the coupling between the intra- and intermolecular modes. With the use of symmetry, full 12-dimensional vibrational energy levels of the water dimer are calculated, predicting accurately the experimentally observed intramolecular fundamentals. It is found that one can also predict accurate intramolecular transition wavenumbers for the water dimer by combining a set of computationally inexpensive reduced-dimensional calculations, thereby guiding future effective-Hamiltonian treatments. [ABSTRACT FROM AUTHOR]

Protein–ligand interactions are crucial for a wide range of physiological processes. Many cellular functions result in these non-covalent "bonds" being mechanically strained, and this can be integral to proper cellular function. Broadly, two classes of force dependence have been observed—slip bonds, where the unbinding rate increases, and catch bonds, where the unbinding rate decreases. Despite much theoretical work, we cannot predict for which protein–ligand pairs, pulling coordinates, and forces a particular rate dependence will appear. Here, we assess the ability of MD simulations combined with enhanced sampling techniques to probe the force dependence of unbinding rates. We show that the infrequent metadynamics technique correctly produces both catch and slip bonding kinetics for model potentials. We then apply it to the well-studied case of a buckyball in a hydrophobic cavity, which appears to exhibit an ideal slip bond. Finally, we compute the force-dependent unbinding rate of biotin–streptavidin. Here, the complex nature of the unbinding process causes the infrequent metadynamics method to begin to break down due to the presence of unbinding intermediates, despite the use of a previously optimized sampling coordinate. Allowing for this limitation, a combination of kinetic and free energy computations predicts an overall slip bond for larger forces consistent with prior experimental results although there are substantial deviations at small forces that require further investigation. This work demonstrates the promise of predicting force-dependent unbinding rates using enhanced sampling MD techniques while also revealing the methodological barriers that must be overcome to tackle more complex targets in the future. [ABSTRACT FROM AUTHOR]

The unimolecular photodissociation dynamics of acetone spanning the entire S1 ← S0 absorption spectrum have been reinvestigated, with a focus on mechanisms that produce CO. At excitation wavelengths of λ > 305.8 nm, all photoproducts are formed on the S0 state after internal conversion. A roaming mechanism forming C2H6 + CO is active in the window λ = 311.2–305.8 nm. From λ = 305.8 to 262 nm, little or no CO is produced with the photochemistry dominated by the Norrish-type I C–C bond cleavage on the lowest excited triplet state, T1. At higher energy (λ < 262 nm), an increasing fraction of CH3CO radicals from the primary reaction have sufficient internal energy to spontaneously decompose to CH3 + CO. A new model is presented to account for the kinetic energy distribution of the secondary CH3 radical, allowing us to determine the height of the energetic barrier to CH3CO decomposition as 68 ± 4 kJ mol−1, which lies midway between previous measurements. The fraction of CO from triple fragmentation rises smoothly from 260 to 248 nm. We see no evidence of the return of roaming, or any other S0 reaction, in this higher energy region of the first electronic absorption band. [ABSTRACT FROM AUTHOR]

Most functional processes of biomolecules are rare events. Key to a rare event is the rare fluctuation that enables the energy activation process that precedes and powers crossing of the activation barrier. However, the physical nature of this rare fluctuation and how it enables energy activation and subsequently barrier crossing are unknown. We developed a novel metric, the reaction capacity pC, that rigorously defines the beginning and parameterizes the progress of energy activation. This enabled us to identify the rare fluctuation as a special phase-space condition that is necessary and sufficient for initiating systematic energy flow from the non-reaction coordinates into the reaction coordinates. The energy activation of a prototype biomolecular isomerization reaction is dominated by kinetic energy transferring into and accumulating in the reaction coordinates, administered by inertial forces alone. This mechanism for energy activation is fundamentally different from the mechanism suggested by Kramers theory. [ABSTRACT FROM AUTHOR]

The ion-induced fragmentation of C H 4 2 + into H+ and C H 3 + is studied using a cold target recoil ion momentum spectroscopy in coincidence with the charge state of the post-collision projectile. Using constant velocity Ar9+ and N3+, results from four different datasets are presented, with a selection on the final charge state of the projectile (Ar8+ or Ar7+ and N2+ or N+). Three distinct dissociation pathways (I, II, and III) are observed for each dataset, with the mean kinetic energy release values of around 4.7, 5.8, and 7.9 eV, respectively. The electronic states that are populated correspond to electronic configurations (1 t 2 ) − 2 and (2 a 1 ) − 1 (1 t 2 ) − 1 of the methane dication, C H 4 2 + . The relative branching ratios between the three pathways are discussed as a function of the charge state of the post-collision projectile, and a strong correlation with the specific nature of the ion–molecule interaction is found. The existing ab initio calculations have provided an explanation only for pathway II. In this article, we propose an explanation for pathway III, but pathway I still remains unexplained and requires further theoretical efforts. A discussion of the dependence of dissociation on the mode of excitation is presented. [ABSTRACT FROM AUTHOR]

The absence of a reliable formulation of the kinetic energy density functional has hindered the development of orbital free density functional theory. Using the data-aided learning paradigm, we propose a simple prescription to accurately model the kinetic energy density of any system. Our method relies on a dictionary of functional forms for local and nonlocal contributions, which have been proposed in the literature, and the appropriate coefficients are calculated via a linear regression framework. To model the nonlocal contributions, we explore two new nonlocal functionals—a functional that captures fluctuations in electronic density and a functional that incorporates gradient information. Since the analytical functional forms of the kernels present in these nonlocal terms are not known from theory, we propose a basis function expansion to model these seemingly difficult nonlocal quantities. This allows us to easily reconstruct kernels for any system using only a few structures. The proposed method is able to learn kinetic energy densities and total kinetic energies of molecular and periodic systems, such as H2, LiH, LiF, and a one-dimensional chain of eight hydrogens using data from Kohn–Sham density functional theory calculations for only a few structures. [ABSTRACT FROM AUTHOR]

Photodissociation of the van der Waals complex Ar–I2 after excitation into the Rydberg states of I2 has been investigated with velocity map imaging of photofragments. Formation of the translationally hot ions of argon Ar+ with three modes in kinetic energy distribution has been revealed. The measured dependence of the kinetic energy of Ar+ on the pumping photon energy indicates the appearance of Ar+ from three channels of the photodissociation of the linear intermediate Ar+–I–I− containing chemically bound argon. These channels are (1) dissociation into Ar++ I2−; (2) three-body dissociation into (Ar+)* + I* + I−, with (Ar+)* and I* being the 2P1/2 states of the species; and (3) two-body electron photodetachment, giving rise to Ar+ + I2 + e. Three indicated channels are similar to those established for the photodissociation of trihalide anions. This similarity confirms the conclusion on the formation of the Ar+–I–I− intermediate, which is isoelectronic to the trihalide anion Cl–I–I−. The mechanism of the Ar+–I–I− formation involves two-photon excitation of the complex Ar–I2 into the Rydberg state of I2 converted into the ion-pair state and further electron transfer from Ar to I+ of the ion-pair state. The self-assembling of the structure making the formation of the Ar+–I–I− intermediate energetically accessible is confirmed by modeling the dynamics in the excited linear complex Ar–I2. Photoexcitation of the van der Waals complexes of noble gases with halogens into the ion-pair states of halogen is supposed to be a promising approach for generating the new chemical compounds of noble gas atoms. [ABSTRACT FROM AUTHOR]