Journal of Chemical Physics最新文献

Low-energy (0-14 eV) resonance electron interaction and fragment species produced by dissociative electron attachment (DEA) for enantiomeric forms of glutamic acid (Glu) are studied under gas-phase conditions by means of DEA spectroscopy and density functional theory calculations. Contrary to a series of amino acids studied earlier employing the DEA technique, the most abundant species are not associated with the elimination of a hydrogen atom from the parent molecular negative ion. Besides this less intense closed-shell [Glu - H]- fragment, only two mass-selected negative ions, [Glu - 19]- and [Glu - 76]-, are detected within the same electron energy region, with the yield maximum observed at around 0.9 eV. This value matches well the energy of vertical electron attachment into the lowest normally empty π* COOH molecular orbital of Glu located at 0.88 eV according to the present B3LYP/6-31G(d) calculations. Although the detection of asymmetric DEA properties a priori is not accessible under the present experimental conditions, "chirality non-conservation" can be associated with some decay channels. Evidently, the measured spectra for the L- and D-forms are found to be identical, the results, nevertheless, being of interest for the forthcoming experiments utilizing spin-polarized electron beam as a chiral factor in the framework of conventional DEA technique.

The density-fitted equation-of-motion (EOM) orbital-optimized second-order perturbation theory (DF-EOM-OMP2) method is presented for the first time. In addition, κ-DF-EOM-MP2 and κ-DF-EOM-OMP2 methods are implemented with the addition of κ-regularization. The accuracy of the DF-EOM-OMP2, κ-DF-EOM-MP2, and κ-DF-EOM-OMP2 methods are compared to the density-fitted EOM-MP2 (DF-EOM-MP2), EOM coupled-cluster (CC) singles and doubles (DF-EOM-CCSD), and EOM-CCSD with the triples excitation correction model [EOM-CCSD(fT)] for excitation energies of many closed- and open-shell chemical systems. The excitation energies computed using different test cases and methods were compared to the EOM-CCSD(fT) method and mean absolute errors (MAEs) are presented. The MAE values of closed- and open-shell cases (closed-shell organic chromophores set, open-shell set, peptide radicals set, and radical set) according to the EOM-CCSD(fT) method show that the κ-regularization technique yields highly accurate results for the first excited states. Our results indicate that the κ-DF-EOM-MP2 and κ-DF-EOM-OMP2 methods perform noticeably better than the DF-EOM-MP2 and DF-EOM-OMP2 methods. They approach the EOM-CCSD quality, at a significantly reduced cost, for the computation of excitation energies. Especially, the κ-DF-EOM-MP2 method provides outstanding results for most test cases considered. Overall, we conclude that the κ-versions of DF-EOM-MP2 and DF-EOM-OMP2 emerge as a useful computational tool for the study of excited-state molecular properties.

The recent discovery of high-temperature superfluorescence in hybrid perovskite thin films has opened new possibilities for harnessing macroscopic quantum phenomena in nanotechnology. This study aimed to elucidate the mechanism that enables high-temperature superfluorescence in these systems. The proposed model describes a quasi-2D Wannier exciton in a thin film that interacts with phonons via the longitudinal optical phonon-exciton Fröhlich interaction. We show that the super-radiant properties of the coherent state in hybrid perovskites are stable against perturbations caused by the longitudinal optical phonon-exciton Fröhlich interaction. Using the multiconfiguration Hartree approach, we derive semiclassical equations of motion for a single-exciton wavefunction, where the vibrational degrees of freedom interact with the Wannier exciton through a mean-field Hartree term. Super-radiance is effectively described by a non-Hermitian term in the Hamiltonian. The analysis was then extended to multiple excited states using the semiclassical Hamiltonian as the basic model. We demonstrate that the ground state of the model exciton Hamiltonian with long-range interactions is a symmetric Dicke super-radiant state, where the Fröhlich interaction is nullified. The additional density matrix-based consideration draws an analogy between this system and stable systems, where the conservation laws determine the nullification of the constant (momentum-independent) decay rate part. In the exciton-phonon system, nullification is associated with the absence of a momentum-independent component in the Wannier exciton-phonon interaction coupling function.

We combine ab initio molecular dynamics (AIMD) simulations with an unsupervised machine learning approach to automate the search for possible configurations of CO oxidation reaction intermediates catalyzed by the atomically dispersed Pt1/TiO2 catalyst. Following the example of Roncoroni and co-workers [Phys. Chem. Chem. Phys. 25, 13741 (2023)], we employ t-distributed stochastic neighbor embedding and hierarchical density-based spatial clustering of applications with noise to reduce the dimensionality and cluster AIMD snapshots based on the local coordination environment of Pt. We identify new local minima, particularly in cases where CO2 is bound to the active site, because it can coordinate in various ways with both the metal and support. The new minima constitute additional elementary steps in some proposed pathways for CO oxidation, resulting in turnover frequencies that differ from prior estimates by several orders of magnitude. This work, therefore, demonstrates that configuration sampling is a necessary component of computational studies of catalytic cycles for atomically dispersed catalysts.

In single vibronic level (SVL) fluorescence experiments, the electronically excited initial state is also excited in one or several vibrational modes. Because computing such spectra by evaluating all contributing Franck-Condon factors becomes impractical (and unnecessary) in large systems, here we propose a time-dependent approach based on Hagedorn wavepacket dynamics. We use Hagedorn functions-products of a Gaussian and carefully generated polynomials-to represent SVL initial states because in systems whose potential is at most quadratic, Hagedorn functions are exact solutions to the time-dependent Schrödinger equation and can be propagated with the same equations of motion as a simple Gaussian wavepacket. Having developed an efficient recursive algorithm to compute the overlaps between two Hagedorn wavepackets, we can now evaluate emission spectra from arbitrary vibronic levels using a single trajectory. We validate the method in two-dimensional global harmonic models by comparing it with quantum split-operator calculations. In addition, we study the effects of displacement, distortion (squeezing), and Duschinsky rotation on SVL fluorescence spectra. Finally, we demonstrate the applicability of the Hagedorn approach to high-dimensional systems on a displaced, distorted, and Duschinsky-rotated harmonic model with 100 degrees of freedom.

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.

The recent discovery of cyano-substituted aromatic and two-ring polycyclic aromatic hydrocarbon molecules in Taurus Molecular Cloud-1 has prompted questions on how the electronic structure and excited-state dynamics of these molecules are linked with their existence and abundance. Here, we report a photodetachment and frequency- and angle-resolved photoelectron spectroscopy study of jet-cooled para-deprotonated benzonitrile (p-[Bzn-H]-). The adiabatic detachment energy was determined as 1.70 ± 0.01 eV, in good agreement with CCSD(T)/aug-cc-pVTZ calculations. The spectra across the first few electron-volts above threshold are dominated by prompt autodetachment processes associated with excitation of at least five short-lived (tens of femtoseconds) temporary anion shaped resonances since excitation cross sections are several orders of magnitude larger than direct photodetachment cross sections. The photoexcitation vibronic profile is dominated by a ≈640 cm-1 ring deformation mode. [Bzn-H]- lacks a valence-localized excited state situated below the detachment threshold and does not exhibit thermionic emission following excitation of the temporary anion resonances. Thus, [Bzn-H]- is unlikely to be stable in many interstellar environments.

Polarization-dependent photoemission electron microscopy (PD-PEEM) exploits spatial variation in the optical selection rules of materials to image domain formation and material organization on the nanoscale. In this Perspective, we discuss the mechanism of PD-PEEM that results in the observed image contrast in experiments and provide examples of a wide range of material domain structures that PD-PEEM has been able to elucidate, including molecular and polymer domains, local electronic structure and defect symmetry, (anti)ferroelectricity, and ferromagnetism. In the end, we discuss challenges and new directions that are possible with this tool for probing domain structure in materials, including investigating the formation of transient ordered states, multiferroics, and the influence of molecular and polymer order and disorder on excited state dynamics and charge transport.

The structures of metal nanoparticles (NPs) significantly influence their catalytic reactivities. Recent in situ experimental observations of dramatic structural changes in NPs underscore the need to establish a dynamic structure-property relationship that accounts for the reconstruction of NPs in reactive environments. Here, we present the MOSP, a free and open-source graphical user interface (GUI) package designed to simulate the structure and reactivity of metal NPs under operando conditions. MOSP integrates two models: the multiscale structure reconstruction model predicting equilibrium metal NP structures under specific reaction conditions and the kinetic Monte Carlo model simulating the reaction dynamics. This combination allows for the exploration of the dynamic structure-property relationships of NPs. MOSP enhances user accessibility through its intuitive GUI, facilitating easy input, post-processing, and visualization of simulation data. This article is the release note of MOSP, focusing on its implementation and functionality.

A newly parameterized combined density functional theory and multi-reference configuration interaction (DFT/MRCI) Hamiltonian, termed core-valence separation (CVS)-QE12, is defined for the computation of K-shell core-excitation and core-ionization energies. This CVS counterpart to the recently reported QE8 Hamiltonian [Costain et al., J. Chem. Phys, 160, 224106 (2024)] is parameterized by fitting to benchmark quality ab initio data. The definition of the CVS-QE12 and QE8 Hamiltonians differ from previous CVS-DFT/MRCI parameterizations in three primary ways: (i) the replacement of the BHLYP exchange-correlation functional with QTP17 to yield a balanced description of both core and valence excitation energies, (ii) the adoption of a new, three-parameter damping function, and (iii) the introduction of separate scaling of the core-valence and valence-valence Coulombic interactions. Crucially, the parameters of the CVS-QE12 Hamiltonian are obtained via fitting exclusively to highly accurate ab initio vertical core-excitation and ionization energies computed at the CVS-EOM-CCSDT level of theory. The CVS-QE12 Hamiltonian is validated against further benchmark computations and is found to furnish K-edge core vertical excitation and ionization energies exhibiting absolute errors ≤0.5 eV at low computational cost.