Demonstrating the Connection between the Nonvalence Correlation-Bound Anions of Polyaromatic Hydrocarbons and the Image Potential States of Graphene Using a One-Electron Model Hamiltonian.

The ground and excited state nonvalence correlation-bound (NVCB) anion states of the C6n2H6n hexagonal polycylic aromatic hydrocarbons and of hexagonal C6n2 graphene nanoflakes are characterized using a one-electron model Hamiltonian which incorporates atomic electrostatic moments up to the quadrupole, coupled inducible charges and dipoles, and atom-centered Gaussians to describe the short-range repulsive interactions. Extrapolation of the calculated electron binding energies of the lowest energy symmetric and antisymmetric (with respect to the molecular plane) NVCB anions of both the polycylic aromatic hydrocarbons and the carbon nanoflakes to the n → ∞ limit yields binding energies that are in good agreement with those of the most stable symmetric and antisymmetric image potential states of freestanding graphene as determined from two-photon photoemission spectroscopy (2PPE) experiments.

Effect of N Atom Substitution on Electronic Resonances: A 2D Photoelectron Spectroscopic and Computational Study of Anthracene, Acridine, and Phenazine Anions.

The accommodation of an excess electron by polycyclic aromatic hydrocarbons (PAHs) has important chemical and technological implications ranging from molecular electronics to charge balance in interstellar molecular clouds. Here, we use two-dimensional photoelectron spectroscopy and equation-of-motion coupled-cluster calculations of the radical anions of acridine (C13H9N-) and phenazine (C12H8N2-) and compare our results for these species to those for the anthracene anion (C14H10-). The calculations predict the observed resonances and additionally find low-energy two-particle-one-hole states, which are not immediately apparent in the spectra, and offer a slightly revised interpretation of the resonances in anthracene. The study of acridine and phenazine allows us to understand how N atom substitution affects electron accommodation. While the electron affinity associated with the ground state anion undergoes a sizable increase with the successive substitution of N atoms, the two lowest energy excited anion states are not affected significantly by the substitution. The net result is that there is an increase in the energy gap between the two lowest energy resonances and the bound ground electronic state of the radical anion from anthracene to acridine to phenazine. Based on an energy gap law for the rate of internal conversion, this increased gap makes ground state formation progressively less likely, as evidenced by the photoelectron spectra.

Contemporary Prophylactic Antibiotic Practices and Adjunct Therapies in Autologous Fat Grafting Procedures: A Survey of The Aesthetic Society Members.

Background: Autologous fat grafting (AFG) is a widely used surgical technique that involves extracting a patient's own adipose tissue and transferring it to different areas of the body. This practice is still evolving. Guidelines for antibiotic prophylaxis and use of adjuncts in plastic surgery are currently limited, with a notable absence of standardized guidelines for AFG. Objectives: In this survey, we assess contemporary antibiotic practices and adjuncts in AFG procedures. Methods: A 52-question survey was emailed to 3106 active members of The Aesthetic Society. Two hundred and ninety-three responses were recorded, representing a 9% response rate. Results: We analyzed 288 responses. The most common AFG procedures were facial (38%), gluteal (34%), and breast (27%) augmentation. Preoperative antibiotics were used by 84.0% overall, with rates of 74.3%, 88.0%, and 92.7% in face, breast, and gluteal AFG, respectively. Lipoaspirate-antibiotic mixing was reported by 19.8%, mainly during gluteal AFG (46.9%), and less so in face (2.8%) and breast (8%) AFG. Notably, 46.9% of surgeons administered prolonged prophylaxis for 72 h or more. Tranexamic acid was utilized by 39.9% of the surveyed surgeons. Platelet-rich plasma was used by 5.6%. Doppler ultrasound was incorporated by 16.7% in AFG, with 21.5% in gluteal AFG, 14% in the face, and 19% in breast procedures. Conclusions: In this survey, we offer insights into antibiotic practices and adjunct therapies in AFG, especially intraoperative antibiotic mixing. Practices among members of The Aesthetic Society vary from guidelines. It is crucial to standardize practices and conduct further research to pave the way for evidence-based guidelines in AFG.

Use of bound state methods to calculate partial and total widths of shape resonances.

In this work we study the 2Π resonances of a two-site model system designed to mimic a smooth transition from the 2Πg temporary anion of N2 to the 2Π temporary anion of CO. The model system possesses the advantage that scattering and bound state (L2) methods can be directly compared without obfuscating electron-correlation effects. Specifically, we compare resonance parameters obtained with the complex Kohn variational (CKV) method with those from stabilization, complex absorbing potential, and regularized analytical continuation calculations. The CKV calculations provide p-wave and d-wave widths, the sum of which provides a good approximation of the total width. Then we demonstrate that the width obtained with modified bound state methods depends on the basis set employed: It can be the total width, a partial width, or an ill-defined sum of partial widths. Provided the basis set is chosen appropriately, widths from bound state methods agree well with the CKV results.

Application of a Fluctuating Charge Polarization Model to Large Polyaromatic Hydrocarbons and Graphene Nanoflakes.

We present a polarization model incorporating coupled fluctuating charges and point inducible dipoles that is able to accurately describe the dipole polarizabilities of small hydrocarbons and, for sufficiently large graphene nanoflakes, reproduce the classical image potential of an infinite conducting sheet. When our fluctuating charge model is applied to the hexagonal carbon nanoflake C60000 we attain excellent agreement with the image potential and induced charge distribution of a conducting sheet. With the inclusion of inducible dipole terms, the model predicts an image plane of zim = 1.3334 a0, which falls in line with prior estimates for graphene. We consider the case of two charges placed on opposite sides of C60000 and find that the fluctuating charge model reproduces classical electrostatics once again. By testing opposing and similar signs of the external charges, we conclude that an atomically thin molecule or extended system does not fully screen their interaction.

Isotope Effects in the Zundel-Eigen Isomerization of H+(H2O)6.

The isomerization pathway between the energetically low-lying Zundel and Eigen isomers of the protonated water hexamer was investigated using high-level ab initio calculations including a treatment of zero-point corrections. On the basis of these calculations, the Zundel-Eigen isomerization was found to proceed through a stable intermediate isomer, which consists of a four-membered ring with two single acceptor water molecules. The inclusion of vibrational zero-point energy is shown to be important for accurately establishing the relative energies of the three relevant isomers involved in the Zundel-Eigen isomerization. Diffusion Monte Carlo calculations including anharmonic vibrational effects show that all three isomers of H+(H2O)6 and D+(D2O)6 have well-defined structures. The energetic ordering of the three isomers changes upon deuteration. The implications of these results for the vibrational spectra of H+(H2O)6 and D+(D2O)6 are also discussed.

Crystallographic Mapping and Tuning of Water Adsorption in Metal-Organic Frameworks Featuring Distinct Open Metal Sites.

Crucial steps toward designing water sorption materials and fine-tuning their properties for specific applications include precise identification of adsorption sites and establishment of rigorous molecular-level insight into the water adsorption process. We report stepwise crystallographic mapping and density functional theory computations of adsorbed water molecules in ALP-MOF-1, a metal-organic framework decorated with distinct open metal sites and carbonyl functional groups that serve as water anchoring sites for seeding the nucleation of a complex water network. Identification of an unusual water adsorption step in ALP-MOF-1 motivated the tuning of metal ion composition to adjust water uptake. These studies provide direct evidence that the identity of the open metal sites in MOFs can dramatically affect water adsorption behavior between 0 and ∼20% RH and that multiple proximal water anchoring sites along the MOF skeleton facilitate water uptake which could be potentially useful for applications requiring rapid and energetically facile water sorption.

Non-Valence Anions of Pyridine and the Diazines.

The dipole-bound anions of pyridine, pyridazine, and pyrimidine are characterized using equation of motion coupled cluster singles and doubles calculations. These calculations predict that the anions of pyridine, pyrimidine and pyridazine are bound in the Born-Oppenheimer approximation by 0.05, 0.8, and 19.0 meV, respectively. The binding energies of pyrimidine and pyridazine are large enough that the anions will remain bound even when allowing for corrections to the Born-Oppenheimer approximation, while that of pyridine is a borderline case. We were unable to find a stable non-valence correlation-bound anion for pyrazine, which has a zero dipole moment.

Water Network Shape-Dependence of Local Interactions with the Microhydrated -NO2- and -CO2- Anionic Head Groups by Cold Ion Vibrational Spectroscopy.

We report the structural evolutions of water networks and solvatochromic response of the CH3NO2- radical anion in the OH and CH stretching regions by analysis of the vibrational spectra displayed by cryogenically cooled CH3NO2-·(H2O)n=1-6 clusters. The OH stretching bands evolve with a surprisingly large discontinuity at n = 6, which features the emergence of an intense, strongly red-shifted band along with a weaker feature that appears in the region assigned to a free OH fundamental. Very similar behavior is displayed by the perdeuterated carboxylate clusters, RCO2-·(H2O)n=5-7 (R = CD3CD2), indicating that this behavior is a general feature in the microhydration of the triatomic anionic domain and not associated with CH oscillators. Electronic structure calculations trace this behavior to the formation of a "book" isomer of the water hexamer that adopts a configuration in which one of the water molecules resides in an acceptor-acceptor-donor (AAD) (A = acceptor, D = donor) H-bonding site. Excitation of the bound OH in the AAD site explores the local network topology best suited to stabilize an incipient -XO2H-OH-(H2O)2 intracluster proton-transfer reaction. These systems thus provide particularly clear examples where the network shape controls the potential energy landscape that governs water network-mediated, intracluster proton transfer. The CH stretching bands of the CH3NO2-·(H2O)n=1-6 clusters also exhibit strong solvatochromic shifts, but in this case, they smoothly blue-shift with increasing hydration with no discontinuity at n = 6. This behavior is analyzed in the context of the solute-ion polarizability response and partial charge transfer to the water networks.

The binding of atomic hydrogen on graphene from density functional theory and diffusion Monte Carlo calculations.

In this work, density functional theory (DFT) and diffusion Monte Carlo (DMC) methods are used to calculate the binding energy of a H atom chemisorbed on the graphene surface. The DMC value of the binding energy is about 16% smaller in magnitude than the Perdew-Burke-Ernzerhof (PBE) result. The inclusion of exact exchange through the use of the Heyd-Scuseria-Ernzerhof functional brings the DFT value of the binding energy closer in line with the DMC result. It is also found that there are significant differences in the charge distributions determined using PBE and DMC approaches.

Vibrational Signatures of HNO3 Acidity When Complexed with Microhydrated Alkali Metal Ions, M+·(HNO3)(H2O)n=5 (M = Li, K, Na, Rb, Cs), at 20 K.

The speciation of strong acids like HNO3 under conditions of restricted hydration is an important factor in the rates of chemical reactions at the air-water interface. Here, we explore the trade-offs at play when HNO3 is attached to alkali ions (Li+-Cs+) with four water molecules in their primary hydration shells. This is achieved by analyzing the vibrational spectra of the M+·(HNO3)(H2O)5 clusters cooled to about 20 K in a cryogenic photofragmentation mass spectrometer. The local acidity of the acidic OH group is estimated by the extent of the red shift in its stretching frequency when attached to a single water molecule. The persistence of this local structural motif (HNO3-H2O) in all of these alkali metal clusters enables us to determine the competition between the effect of the direct complexation of the acid with the cation, which acts to enhance acidity, and the role of the water network in the first hydration shell around the ions, which acts to counter (screen) the intrinsic effect of the ion. Analysis of the vibrational features associated with the acid molecule, as well as those of the water network, reveals how cooperative interactions in the microhydration regime conspire to effectively offset the intrinsic enhancement of HNO3 acidity afforded by attachment to the smaller cations.

Real-Time Modulation of Hydrogen Evolution Activity of Graphene Electrodes Using Mechanical Strain.

This paper reports the effect of mechanically applied elastic strain on the hydrogen evolution reaction (HER) activity of graphene under acidic conditions. An applied tensile strain of 0.2% on a graphene electrode is shown to lead to a 1-3% increase in the HER current. The tensile strain increases HER activity, whereas compressive strain decreases it. Density functional theory (DFT) calculations using a periodic graphene slab model predict an increase in the adsorption energy of the H atom with growing tensile strain, consistent with an enhancement of the current density in HER, similar to that observed experimentally.

Temporary Anion Resonances of Pyrene: A 2D Photoelectron Imaging and Computational Study.

The low-energy electron-scattering resonances of pyrene were characterized using experimental and computational methods. Experimentally, a two-dimensional photoelectron imaging of the pyrene anion was used to probe the dynamics of resonances over the first 4 eV of the continuum. Computationally, the energies and character of the anion states were determined using equation-of-motion coupled cluster calculations, while taking specific care to avoid the collapse onto discretized continuum levels, and an application of the pairing theorem. Our results are in good agreement with the predictions of electron-scattering calculations that include an offset and with the pyrene anion absorption spectrum in a glass matrix. Taken together, we offer an assignment of the first five electronic resonances of pyrene. Some of the population in the lowest-energy 2B1u resonance was observed to decay to the ground electronic state of the anion, while all other resonances decay by a direct autodetachment. The astronomical relevance of a ground-state electron capture proceeding via a low-energy resonance in pyrene is discussed.

Frontiers of stochastic electronic structure calculations.

In recent years there has been a rapid growth in the development and application of new stochastic methods in electronic structure. These methods are quite diverse, from many-body wave function techniques in real space or determinant space to being used to sum perturbative expansions. This growth has been spurred by the more favorable scaling with the number of electrons and often better parallelization over large numbers of central processing unit (CPU) cores or graphical processing units (GPUs) than for high-end non-stochastic wave function based methods. This special issue of the Journal of Chemical Physics includes 33 papers that describe recent developments and applications in this area. As seen from the articles in the issue, stochastic electronic structure methods are applicable to both molecules and solids and can accurately describe systems with strong electron correlation. This issue was motivated, in part, by the 2019 Telluride Science Research Center workshop on Stochastic Electronic Structure Methods that we organized. Below we briefly describe each of the papers in the special issue, dividing the papers into six subtopics.

Two-Dimensional Adiabatic Model for Calculating Progressions Resulting from Stretch-Rock Coupling in Vibrational Spectra of Anion-Water Complexes.

Several anion-water dimers feature a distinct progression in the OH stretch region of their vibrational spectra. This progression arises from strong anharmonic couplings between the OH stretch and low-frequency intermolecular modes. In this work, we introduce a two-dimensional adiabatic model accounting explicitly for the water and anion rock degrees of freedom and use it to calculate the vibrational spectra of HCO2-·(H2O) and NO3-·(H2O). The spectra calculated by using this model are in excellent agreement with experiment, in terms of both peak spacings and lengths of the progressions, and represent a substantial improvement over earlier models.

Role of Overlap between the Discrete State and Pseudocontinuum States in Stabilization Calculations of Metastable States.

In a diabatic picture metastable states subject to decay by electron detachment can be viewed as arising from the coupling between a discrete state and a continuum. In treating such states with bound-state quantum chemical methods, the continuum is discretized. In this study, we elucidate the role of overlap in this interaction in the application of the stabilization method to temporary anion states. This is accomplished by use of a minimalist stabilization calculation on the lowest energy l=2 (D) resonance of the finite spherical well potential using two basis functions, one describing the diabatic discrete state and the other a diabatic discretized continuum state. We show that even such a simple treatment predicts a complex resonance energy in good agreement with the exact result. If the energy of the discrete state is assumed to be constant, which is tantamount to orthogonalizing the discretized continuum state to the discrete state, it is demonstrated that the square of the off-diagonal coupling has a maximum close to the crossing point of the orthogonalized diabatic curves and that the curvature in the coupling is responsible for the complex stationary point associated with the resonance. Moreover, this curvature is a consequence of the overlap between the two diabatic states.

A Fresh Look at the Role of the Coupling of a Discrete State with a Pseudocontinuum State in the Stabilization Method for Characterizing Metastable States.

The stabilization method is widely used to theoretically characterize temporary anions and other systems displaying resonances. In this approach, information about a metastable state is encoded in the interaction of a diabatic discrete state and discretized continuum solutions, the energy of which are varied by scaling the extent of the basis set. In this work, we identify the aspects of the coupling between the discrete state and the discretized continuum states that encode information about the existence of complex stationary points and, hence, complex resonance energies in stabilization graphs. This allows us to design a simple two-level model for extracting complex resonance energies from stabilization graphs. The resulting model is applied to the 2Πg anion state of N2.

The role of high-order electron correlation effects in a model system for non-valence correlation-bound anions.

The diffusion Monte Carlo (DMC), auxiliary field quantum Monte Carlo (AFQMC), and equation-of-motion coupled cluster (EOM-CC) methods are used to calculate the electron binding energy (EBE) of the non-valence anion state of a model (H2O)4 cluster. Two geometries are considered, one at which the anion is unbound and the other at which it is bound in the Hartree-Fock (HF) approximation. It is demonstrated that DMC calculations can recover from the use of a HF trial wave function that has collapsed onto a discretized continuum solution, although larger EBEs are obtained when using a trial wave function for the anion that provides a more realistic description of the charge distribution and, hence, of the nodal surface. For the geometry at which the cluster has a non-valence correlation-bound anion, both the inclusion of triples in the EOM-CC method and the inclusion of supplemental diffuse d functions in the basis set are important. DMC calculations with suitable trial wave functions give EBE values in good agreement with our best estimate EOM-CC result. AFQMC using a trial wave function for the anion with a realistic electron density gives a value of the EBE nearly identical to the EOM-CC result when using the same basis set. For the geometry at which the anion is bound in the HF approximation, the inclusion of triple excitations in the EOM-CC calculations is much less important. The best estimate EOM-CC EBE value is in good agreement with the results of DMC calculations with appropriate trial wave functions.

Toward a systematic improvement of the fixed-node approximation in diffusion Monte Carlo for solids-A case study in diamond.

While Diffusion Monte Carlo (DMC) is in principle an exact stochastic method for ab initio electronic structure calculations, in practice, the fermionic sign problem necessitates the use of the fixed-node approximation and trial wavefunctions with approximate nodes (or zeros). This approximation introduces a variational error in the energy that potentially can be tested and systematically improved. Here, we present a computational method that produces trial wavefunctions with systematically improvable nodes for DMC calculations of periodic solids. These trial wavefunctions are efficiently generated with the configuration interaction using a perturbative selection made iteratively (CIPSI) method. A simple protocol in which both exact and approximate results for finite supercells are used to extrapolate to the thermodynamic limit is introduced. This approach is illustrated in the case of the carbon diamond using Slater-Jastrow trial wavefunctions including up to one million Slater determinants. Fixed-node DMC energies obtained with such large expansions are much improved, and the fixed-node error is found to decrease monotonically and smoothly as a function of the number of determinants in the trial wavefunction, a property opening the way to a better control of this error. The cohesive energy extrapolated to the thermodynamic limit is in close agreement with the estimated experimental value. Interestingly, this is also the case at the single-determinant level, thus, indicating a very good error cancellation in carbon diamond between the bulk and atomic total fixed-node energies when using single-determinant nodes.