Prosocial motives underlie scientific censorship by scientists: A perspective and research agenda.

Science is among humanity's greatest achievements, yet scientific censorship is rarely studied empirically. We explore the social, psychological, and institutional causes and consequences of scientific censorship (defined as actions aimed at obstructing particular scientific ideas from reaching an audience for reasons other than low scientific quality). Popular narratives suggest that scientific censorship is driven by authoritarian officials with dark motives, such as dogmatism and intolerance. Our analysis suggests that scientific censorship is often driven by scientists, who are primarily motivated by self-protection, benevolence toward peer scholars, and prosocial concerns for the well-being of human social groups. This perspective helps explain both recent findings on scientific censorship and recent changes to scientific institutions, such as the use of harm-based criteria to evaluate research. We discuss unknowns surrounding the consequences of censorship and provide recommendations for improving transparency and accountability in scientific decision-making to enable the exploration of these unknowns. The benefits of censorship may sometimes outweigh costs. However, until costs and benefits are examined empirically, scholars on opposing sides of ongoing debates are left to quarrel based on competing values, assumptions, and intuitions.

Two- and one-photon absorption spectra of aqueous thiocyanate anion highlight the role of symmetry in the condensed phase.

We present the two-photon absorption (2PA) spectrum of aqueous thiocyanate calculated using high-level quantum-chemistry methods. The 2PA spectrum is compared to the one-photon absorption (1PA) spectrum computed using the same computational protocol. Although the two spectra probe the same set of electronic states, the intensity patterns are different, leading to an apparent red-shift of the 2PA spectrum relative to the 1PA spectrum. The presented analysis explains the intensity patterns and attributes the differences between the 1PA and 2PA spectra to the native symmetry of isolated SCN - , which influences the spectra in the low-symmetry solvated environment. The native symmetry also manifests itself in variations of the polarization ratio (e.g., parallel vs. perpendicular cross sections) across the spectrum. The presented results highlight the potential of 2PA spectroscopy and high-level quantum-chemistry methods in studies of condensed-phase phenomena.

Core-ionization spectrum of liquid water.

We present state-of-the-art calculations of the core-ionization spectrum of water. Despite significant progress in procedures developed to mitigate various experimental complications and uncertainties, the experimental determination of ionization energies of solvated species involves several non-trivial steps such as assessing the effect of the surface potential, electrolytes, and finite escape depths of photoelectrons. This provides a motivation to obtain robust theoretical values of the intrinsic bulk ionization energy and the corresponding solvent-induced shift. Here we develop theoretical protocols based on coupled-cluster theory and electrostatic embedding. Our value of the intrinsic solvent-induced shift of the 1sO ionization energy of water is -1.79 eV. The computed absolute position and the width of the 1sO peak in photoelectron spectrum of water are 538.47 eV and 1.44 eV, respectively, agreeing well with the best experimental values.

Revisiting the benzene excimer using [2,2] paracyclophane model system: Experiment and theory.

We report high-level calculations of the excited states of [2,2]-paracyclophane (PCP), which was recently investigated experimentally by ultrafast pump-probe experiments on oriented single crystals [Haggag et al., ChemPhotoChem 6 e202200181 (2022)]. PCP, in which the orientation of the two benzene rings and their range of motion are constrained, serves as a model for studying benzene excimer formation. The character of the excimer state and the state responsible for the brightest transition are similar to those of the benzene dimer. The constrained structure of PCP allows one to focus on the most important degree of freedom, the inter-ring distance. The calculations explain the main features of the transient absorption spectral evolution. This brightest transition of the excimer is polarized along the inter-fragment axis. The absorption of the light polarized in the plane of the rings reveals the presence of other absorbing states of Rydberg character, with much weaker intensities. We also report new transient absorption data obtained by a broadband 8 fs pump, which time-resolve strong modulations of the excimer absorption. The combination of theory and experiment provides a detailed picture of the evolution of the electronic structure of the PCP excimer in the course of a single molecular vibration.

Signatures of the Bromine Atom and Open-Shell Spin Coupling in the X-ray Spectrum of the Bromobenzene Cation.

Tabletop X-ray spectroscopy measurements at the carbon K-edge complemented by ab initio calculations are used to investigate the influence of the bromine atom on the carbon core-valence transitions in the bromobenzene cation (BrBz+). The electronic ground state of the cation is prepared by resonance-enhanced two-photon ionization of neutral bromobenzene (BrBz) and probed by X-rays produced by high-harmonic generation (HHG). Replacing one of the hydrogen atoms in benzene with a bromine atom shifts the transition from the 1sC* orbital of the carbon atom (C*) bonded to bromine by ∼1 eV to higher energy in the X-ray spectrum compared to the other carbon atoms (C). Moreover, in BrBz+, the X-ray spectrum is dominated by two relatively intense transitions, 1sC→π* and 1sC*→σ*(C*-Br), where the second transition is enhanced relative to the neutral BrBz. In addition, a doublet peak shape for these two transitions is observed in the experiment. The 1sC→π* doublet peak shape arises due to the spin coupling of the unpaired electron in the partially vacant π orbital (from ionization) with the two other unpaired electrons resulting from the transition from the 1sC core orbital to the fully vacant π* orbitals. The 1sC*→σ* doublet peak shape results from several transitions involving σ* and vibrational C*-Br mode activations following the UV ionization, which demonstrates the impact of the C*-Br bond length on the core-valence transition as well as on the relaxation geometry of BrBz+.

Multiscale Simulations of the Covalent Inhibition of the SARS-CoV-2 Main Protease: Four Compounds and Three Reaction Mechanisms.

We report the results of computational modeling of the reactions of the SARS-CoV-2 main protease (MPro) with four potential covalent inhibitors. Two of them, carmofur and nirmatrelvir, have shown experimentally the ability to inhibit MPro. Two other compounds, X77A and X77C, were designed computationally in this work. They were derived from the structure of X77, a non-covalent inhibitor forming a tight surface complex with MPro. We modified the X77 structure by introducing warheads capable of reacting with the catalytic cysteine residue in the MPro active site. The reaction mechanisms of the four molecules with MPro were investigated by quantum mechanics/molecular mechanics (QM/MM) simulations. The results show that all four compounds form covalent adducts with the catalytic cysteine Cys 145 of MPro. From the chemical perspective, the reactions of these four molecules with MPro follow three distinct mechanisms. The reactions are initiated by a nucleophilic attack of the thiolate group of the deprotonated cysteine residue from the catalytic dyad Cys145-His41 of MPro. In the case of carmofur and X77A, the covalent binding of the thiolate to the ligand is accompanied by the formation of the fluoro-uracil leaving group. The reaction with X77C follows the nucleophilic aromatic substitution SNAr mechanism. The reaction of MPro with nirmatrelvir (which has a reactive nitrile group) leads to the formation of a covalent thioimidate adduct with the thiolate of the Cys145 residue in the enzyme active site. Our results contribute to the ongoing search for efficient inhibitors of the SARS-CoV-2 enzymes.

Magnetic exchange interactions in binuclear and tetranuclear iron(III) complexes described by spin-flip DFT and Heisenberg effective Hamiltonians.

Low-energy spectra of single-molecule magnets (SMMs) are often described by Heisenberg Hamiltonians. Within this formalism, exchange interactions between magnetic centers determine the ground-state multiplicity and energy separation between the ground and excited states. In this contribution, we extract exchange coupling constants (J) for a set of iron (III) binuclear and tetranuclear complexes from all-electron calculations using non-collinear spin-flip time-dependent density functional theory (NC-SF-TDDFT). For 12 binuclear complexes with J-values ranging from -6 to -132 cm-1 , our benchmark calculations using the short-range hybrid ωPBEh functional and 6-31G(d,p) basis set agree well with the experimentally derived values (mean absolute error of 4.7 cm-1 ). For the tetranuclear SMMs, the computed J constants are within 6 cm-1 from the experimentally derived values. We explore the range of applicability of the Heisenberg model by analyzing bonding patterns in these Fe(III) complexes using natural orbitals (NO), their occupations, and the number of effectively unpaired electrons. The results illustrate the efficiency of the spin-flip protocol for computing the exchange couplings and the utility of the NO analysis in assessing the validity of effective spin Hamiltonians.

Two Algorithms for Excited-State Quantum Solvers: Theory and Application to EOM-UCCSD.

Near-term quantum devices promise to revolutionize quantum chemistry, but simulations using the current noisy intermediate-scale quantum (NISQ) devices are not practical due to their high susceptibility to errors. This motivated the design of NISQ algorithms leveraging classical and quantum resources. While several developments have shown promising results for ground-state simulations, extending the algorithms to excited states remains challenging. This paper presents two cost-efficient excited-state algorithms inspired by the classical Davidson algorithm. We implemented the Davidson method into the quantum self-consistent equation-of-motion unitary coupled-cluster (q-sc-EOM-UCC) excited-state method adapted for quantum hardware. The circuit strategies for generating desired excited states are discussed, implemented, and tested. We demonstrate the performance and accuracy of the proposed algorithms (q-sc-EOM-UCC/Davidson and its variational variant) by simulations of H2, H4, LiH, and H2O molecules. Similar to the classical Davidson scheme, q-sc-EOM-UCC/Davidson algorithms are capable of targeting a small number of excited states of the desired character.

Influence of Transition Metal Electron Configuration on the Structure of Metal-EDTA Complexes.

The vibrational spectra of cold complexes of ethylenediaminetetraacetic acid (EDTA) with transition metal dications in vacuo show how the electronic structure of the metal provides a geometric template for interaction with the functional groups of the binding pocket. The OCO stretching modes of the carboxylate groups of EDTA serve as structural probes, informing on the spin state of the ion as well as the coordination number in the complex. The results highlight the flexibility of EDTA in accepting a large range of metal cations in its binding site.

State-Interaction Approach for Evaluating g-Tensors within EOM-CC and RAS-CI Frameworks: Theory and Benchmarks.

Among various techniques designed for studying open-shell species, electron paramagnetic resonance (EPR) spectroscopy plays an important role. The key quantity measured by EPR is the g-tensor, describing the coupling between an external magnetic field and molecular electronic spin. One theoretical framework for quantum chemistry calculations of g-tensors is based on response theory, which involves substantial developments that are specific to the underlying electronic structure models. A simplified and easier-to-implement approach is based on the state-interaction scheme, in which perturbation is included by considering a small number of states. We describe and benchmark the state-interaction approach using equation-of-motion coupled-cluster and restricted-active-space configuration interaction wave functions. The analysis confirms that this approach can deliver accurate results and highlights caveats of applying it, such as a choice of the reference state, convergence with respect to the number of states used in calculations, etc. The analysis also contributes toward a better understanding of challenges in calculations of higher-order properties using approximate wave functions.

Theory, implementation, and disappointing results for two-photon absorption cross sections within the doubly electron-attached equation-of-motion coupled-cluster framework.

The equation-of-motion coupled-cluster singles and doubles method with double electron attachment (EOM-DEA-CCSD) is capable of computing reliable energies, wave functions, and first-order properties of excited states in diradicals and polyenes that have a significant doubly excited character with respect to the ground state, without the need for including the computationally expensive triple excitations. Here, we extend the capabilities of the EOM-DEA-CCSD method to the calculations of a multiphoton property, two-photon absorption (2PA) cross sections. Closed-form expressions for the 2PA cross sections are derived within the expectation-value approach using response wave functions. We analyze the performance of this new implementation by comparing the EOM-DEA-CCSD energies and 2PA cross sections with those computed using the CC3 quadratic response theory approach. As benchmark systems, we consider transitions to the states with doubly excited character in twisted ethene and in polyenes, for which EOM-EE-CCSD (EOM-CCSD for excitation energies) performs poorly. The EOM-DEA-CCSD 2PA cross sections are comparable with the CC3 results for twisted ethene; however, the discrepancies between the two methods are large for hexatriene. The observed trends are explained by configurational analysis of the 2PA channels.

Simulating excited-state complex ensembles: Fluorescence and solvatochromism in amine-arene exciplexes.

Exciplexes are excited-state complexes formed as a result of partial charge transfer from the donor to the acceptor species when one moiety of the donor-acceptor pair is electronically excited. The arene-amine exciplex formed between oligo-(p-phenylene) (OPP) and triethylamine (TEA) is of interest in the catalytic photoreduction of CO2 because it can compete with complete electron transfer to the OPP catalyst. Therefore, formation of the exciplex can hinder the generation of a radical anion OPP·- necessary for subsequent CO2 reduction. We report an implementation of a workflow automating quantum-chemistry calculations that generate and characterize an ensemble of structures to represent this exciplex state. We use FireWorks, Pymatgen, and Custodian Python packages for high-throughput ensemble generation. The workflow includes time-dependent density functional theory optimization, verification of excited-state minima, and exciplex characterization with natural transition orbitals, exciton analysis, excited-state Mulliken charges, and energy decomposition analysis. Fluorescence spectra computed for these ensembles using Boltzmann-weighted contributions of each structure agree better with experiment than our previous calculations based on a single representative exciplex structure [Kron et al., J. Phys. Chem. A 126, 2319-2329 (2022)]. The ensemble description of the exciplex state also reproduces an experimentally observed red shift of the emission spectrum of [OPP-4-TEA]* relative to [OPP-3-TEA]*. The workflow developed here streamlines otherwise labor-intensive calculations that would require significant user involvement and intervention.

The Auger spectrum of benzene.

We present an ab initio computational study of the Auger electron spectrum of benzene. Auger electron spectroscopy exploits the Auger-Meitner effect, and although it is established as an analytic technique, the theoretical modeling of molecular Auger spectra from first principles remains challenging. Here, we use coupled-cluster theory and equation-of-motion coupled-cluster theory combined with two approaches to describe the decaying nature of core-ionized states: (i) Feshbach-Fano resonance theory and (ii) the method of complex basis functions. The spectra computed with these two approaches are in excellent agreement with each other and also agree well with experimental Auger spectra of benzene. The Auger spectrum of benzene features two well-resolved peaks at Auger electron energies above 260 eV, which correspond to final states with two electrons removed from the 1e1g and 3e2g highest occupied molecular orbitals. At lower Auger electron energies, the spectrum is less well resolved, and the peaks comprise multiple final states of the benzene dication. In line with theoretical considerations, singlet decay channels contribute more to the total Auger intensity than the corresponding triplet decay channels.

Dipole effects in the photoelectron angular distributions of the sulfur monoxide anion.

Photoelectron angular distributions (PADs) in SO- photodetachment using linearly polarized 355 nm (3.49 eV), 532 nm (2.33 eV), and 611 nm (2.03 eV) light were investigated via photoelectron imaging spectroscopy. The measurements at 532 and 611 nm access the X3Σ- and a1Δ electronic states of SO, whereas the measurements at 355 nm also access the b1Σ+ state. In aggregate, the photoelectron anisotropy parameter values follow the general trend with respect to electron kinetic energy (eKE) expected for π*-orbital photodetachment. The trend is similar to O2-, but the minimum of the SO- curve is shifted to smaller eKE. This shift is mainly attributed to the exit-channel interactions of the departing electron with the dipole moment of the neutral SO core, rather than the differing shapes of the SO- and O2- molecular orbitals. Of the several ab initio models considered, two approaches yield good agreement with the experiment: one representing the departing electron as a superposition of eigenfunctions of a point dipole-field Hamiltonian, and another describing the outgoing electron in terms of Coulomb waves originating from two separated charge centers, with a partial positive charge on the sulfur and an equal negative charge on the oxygen. These fundamentally related approaches support the conclusion that electron-dipole interactions in the exit channel of SO- photodetachment play an important role in shaping the PADs. While a similar conclusion was previously reached for photodetachment from σ orbitals of CN- (Hart, Lyle, Spellberg, Krylov, Mabbs, J. Phys. Chem. Lett., 2021, 12, 10086-10092), the present work includes the first extension of the dipole-field model to detachment from π* orbitals.

DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science.

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.

Spin-orbit couplings within spin-conserving and spin-flipping time-dependent density functional theory: Implementation and benchmark calculations.

We present a new implementation for computing spin-orbit couplings (SOCs) within a time-dependent density-functional theory (TD-DFT) framework in the standard spin-conserving formulation as well in the spin-flip variant (SF-TD-DFT). This approach employs the Breit-Pauli Hamiltonian and Wigner-Eckart's theorem applied to the reduced one-particle transition density matrices, together with the spin-orbit mean-field treatment of the two-electron contributions. We use a state-interaction procedure and compute the SOC matrix elements using zero-order non-relativistic states. Benchmark calculations using several closed-shell organic molecules, diradicals, and a single-molecule magnet illustrate the efficiency of the SOC protocol. The results for organic molecules (described by standard TD-DFT) show that SOCs are insensitive to the choice of the functional or basis sets, as long as the states of the same characters are compared. In contrast, the SF-TD-DFT results for small diradicals (CH2, NH2 +, SiH2, and PH2 +) show strong functional dependence. The spin-reversal energy barrier in a Fe(III) single-molecule magnet computed using non-collinear SF-TD-DFT (PBE0, ωPBEh/cc-pVDZ) agrees well with the experimental estimate.

Spectroscopic signatures of states in the continuum characterized by a joint experimental and theoretical study of pyrrole.

We report a combined experimental and theoretical investigation of electron-molecule interactions using pyrrole as a model system. Experimental two-dimensional electron energy loss spectra (EELS) encode information about the vibrational states of the molecule as well as the position and structure of electronic resonances. The calculations using complex-valued extensions of equation-of-motion coupled-cluster theory (based on non-Hermitian quantum mechanics) facilitate the assignment of all major EELS features. We confirm the two previously described π resonances at about 2.5 and 3.5 eV (the calculations place these two states at 2.92 and 3.53 eV vertically and 2.63 and 3.27 eV adiabatically). The calculations also predict a low-lying resonance at 0.46 eV, which has a mixed character-of a dipole-bound state and σ* type. This resonance becomes stabilized at one quanta of the NH excitation, giving rise to the sharp feature at 0.9 eV in the corresponding EELS. Calculations of Franck-Condon factors explain the observed variations in the vibrational excitation patterns. The ability of theory to describe EELS provides a concrete illustration of the utility of non-Hermitian quantum chemistry, which extends such important concepts as potential energy surfaces and molecular orbitals to states embedded in the continuum.

Feshbach-Fano approach for calculation of Auger decay rates using equation-of-motion coupled-cluster wave functions. II. Numerical examples and benchmarks.

X-ray photon absorption leads to the creation of highly excited species, which often decay through the Auger process. The theoretical treatment of Auger decay is challenging because of the resonance nature of the initial core-excited or core-ionized states and the continuous nature of the ejected electron. In Paper I [W. Skomorowski and A. I. Krylov, J. Chem. Phys. 154, 084124 (2021)], we have introduced a theoretical framework for computing Auger rates based on the Feshbach-Fano approach and the equation-of-motion coupled-cluster ansätze augmented with core-valence separation. The outgoing Auger electron is described with a continuum orbital. We considered two approximate descriptions-a plane wave and a Coulomb wave with an effective charge. Here, we use the developed methodology to calculate Auger transition rates in core-ionized and core-excited benchmark systems (Ne, H2O, CH4, and CO2). Comparison with the available experimental spectra shows that the proposed computational scheme provides reliable ab initio predictions of the Auger spectra. The reliability, cost efficiency, and robust computational setup of this methodology offer advantages in applications to a large variety of systems.

Feshbach-Fano approach for calculation of Auger decay rates using equation-of-motion coupled-cluster wave functions. I. Theory and implementation.

X-ray absorption creates electron vacancies in the core shell. These highly excited states often relax by Auger decay-an autoionization process in which one valence electron fills the core hole and another valence electron is ejected into the ionization continuum. Despite the important role of Auger processes in many experimental settings, their first-principles modeling is challenging, even for small systems. The difficulty stems from the need to describe many-electron continuum (unbound) states, which cannot be tackled with standard quantum-chemistry methods. We present a novel approach to calculate Auger decay rates by combining Feshbach-Fano resonance theory with the equation-of-motion coupled-cluster single double (EOM-CCSD) framework. We use the core-valence separation scheme to define projectors into the bound (square-integrable) and unbound (continuum) subspaces of the full function space. The continuum many-body decay states are represented by products of an appropriate EOM-CCSD state and a free-electron state, described by a continuum orbital. The Auger rates are expressed in terms of reduced quantities, two-body Dyson amplitudes (objects analogous to the two-particle transition density matrix), contracted with two-electron bound-continuum integrals. Here, we consider two approximate treatments of the free electron: a plane wave and a Coulomb wave with an effective charge, which allow us to evaluate all requisite integrals analytically; however, the theory can be extended to incorporate a more sophisticated description of the continuum orbital.