Lignin with controlled structural properties by N-heterocycle-based deep eutectic solvent extraction.

The complex and heterogeneous nature of the lignin macromolecule has presented a lasting barrier to its utilization. To achieve high lignin yield, the technical lignin extraction process usually severely modifies and condenses the native structure of lignin, which is a critical drawback for its utilization in conversion processes. In addition, there is no method capable of separating lignin from plant biomass with controlled structural properties. Here, we developed an N-heterocycle-based deep eutectic solvent formed between lactic acid and pyrazole (La-Py DES) with a binary hydrogen bonding functionality resulting in a high affinity toward lignin. Up to 93.7% of lignin was extracted from wheat straw biomass at varying conditions from 90 °C to 145 °C. Through careful selection of treatment conditions as well as lactic acid to pyrazole ratios, lignin with controlled levels of ether linkage content, hydroxyl group content, and average molecular weight can be generated. Under mild extraction conditions (90 °C to 120 °C), light-colored native-like lignin can be produced with up to 80% yield, whereas ether linkage-free lignin with low polydispersity can be obtained at 145 °C. Overall, this study offers a new strategy for native lignin extraction and generating lignin with controlled structural properties.

A critical comparison of CH⋯π versus π⋯π interactions in the benzene dimer: obtaining benchmarks at the CCSD(T) level and assessing the accuracy of lower scaling methods.

We have established CCSD(T)/CBS (Complete Basis Set) limits for 3 stationary points on the benzene dimer potential energy surface, corresponding to the π⋯π (parallel displaced or PD(C2h), minimum) and CH⋯π (T-shaped or T(C2v), transition state) and tilted T-shaped (or TT(Cs), minimum) bonding scenarios considering both the structure and binding energy. The CCSD(T)/CBS binding energies are -2.65 ± 0.02 (PD), -2.74 ± 0.03 (T), and -2.83 ± 0.01 kcal mol-1 (TT). To this end, the CH⋯π is ∼0.2 kcal mol-1 stronger than the π⋯π interaction, whereas the tilting of the CH donating benzene molecule with respect to the other benzene is worth 0.1 kcal mol-1. As previously discussed in the literature, the MP2 level of theory does not provide a close match for either the energy or structure, yet the SCS-MP2 yields structures in excellent agreement with respect to the CCSD(T) result. It is found that the SCS-MI-MP2 also gives optimized structures very close to SCS-MP2 (within ∼0.01 Å of the benchmark). Despite the closer match in structure, the spin-biased MP2 methods (SCS-, SCS-MI-, and SOS-MP2) incorrectly predict the relative stabilities of the isomers. That said, none of the spin biased MP2 methods offers a good compromise between energy and structure for the systems examined. Finally, the CCSD(T)/CBS benchmarks were used to assess the performance of 13 DFT functionals selected from different rungs of Jacob's ladder. Several functionals such as TPSS-D3, B3LYP-D3, B97-D, B97-D3, and B2PLYP-D3 provided a good description of the binding energies for both CH⋯π and π⋯π interactions, yielding values within 6% of the CCSD(T)/CBS benchmark values. Unlike the MP2 methods, these functionals correctly predict the relative stability of the PD(C2h) and T(C2v) dimers. Further, we find that there is no systematic improvement as Jacob's ladder is ascended (increased complexity of functional). The best functionals that result in a good compromise between structure and energy accuracy are B97-D3 and B2PLYP-D3 for both the CH⋯π and π⋯π interaction. Despite the impressive performance of these functionals, a challenge that remains is ensuring the transferability of these density functionals in accurately describing the interaction between dimers of larger aromatic molecules, the latter requiring high-level benchmarks for these systems.

G0W0 Ionization Potentials of First-Row Transition Metal Aqua Ions.

We report computations of the vertical ionization potentials within the GW approximation of the near-complete series of first-row transition metal (V-Cu) aqua ions in their most common oxidation states, i.e., V3+, Cr3+, Cr2+, Mn2+, Fe3+, Fe2+, Co2+, Ni2+, and Cu2+. The d-orbital occupancy of these systems spans a broad range from d2 to d9. All of the structures were first optimized at the density functional theory level using a large cluster of explicit water molecules that are embedded in a continuum solvation model. Vertical ionization potentials were computed with the one-shot G0W0 approach on a range of transition metal ion clusters (6, 18, 40, and 60 explicit water molecules), wherein the convergence with respect to the basis set size was evaluated using the systems with 40 water molecules. We assess the results using three different density functional approximations as starting points for the vertical ionization potential calculations, namely, G0W0@PBE, G0W0@PBE0, and G0W0@r2SCAN. While the predicted ground-state structures are similar to all three exchange-correlation functionals, the vertical ionization potentials were in closer agreement with experiment when using the G0W0@PBE0 and G0W0@r2SCAN approaches, with the r2SCAN-based calculations being significantly less expensive. Computed bond distances and vertical ionization potentials for all structures are in good agreement with available experimental data.

Direct functionalization of C-H bonds by electrophilic anions.

Alkanes and [B12X12]2- (X = Cl, Br) are both stable compounds which are difficult to functionalize. Here we demonstrate the formation of a boron-carbon bond between these substances in a two-step process. Fragmentation of [B12X12]2- in the gas phase generates highly reactive [B12X11]- ions which spontaneously react with alkanes. The reaction mechanism was investigated using tandem mass spectrometry and gas-phase vibrational spectroscopy combined with electronic structure calculations. [B12X11]- reacts by an electrophilic substitution of a proton in an alkane resulting in a B-C bond formation. The product is a dianionic [B12X11CnH2n+1]2- species, to which H+ is electrostatically bound. High-flux ion soft landing was performed to codeposit [B12X11]- and complex organic molecules (phthalates) in thin layers on surfaces. Molecular structure analysis of the product films revealed that C-H functionalization by [B12X11]- occurred in the presence of other more reactive functional groups. This observation demonstrates the utility of highly reactive fragment ions for selective bond formation processes and may pave the way for the use of gas-phase ion chemistry for the generation of complex molecular structures in the condensed phase.

Binding of saturated and unsaturated C6-hydrocarbons to the electrophilic anion [B12Br11]-: a systematic mechanistic study.

The highly reactive gaseous ion [B12Br11]- is a metal-free closed-shell anion which spontaneously forms covalent bonds with hydrocarbon molecules, including alkanes. Herein, we systematically investigate the reaction mechanism for binding of [B12Br11]- to the five hexane isomers yielding [B12Br11(C6H14)]-, as well as to cyclohexane and several hexene isomers (yielding [B12Br11(C6H12)]-) using collision-induced dissociation (CID), infrared photodissociation spectroscopy (IRPD) and computational methods. CID of the different [B12Br11(C6H14)]- ions results in distinct fragmentation patterns dependent on the structure of the hexane isomer. The observed fragmentation reactions provide insights into the addition mechanism of [B12Br11]- to hexane. Based on the observed CID patterns, we identified that either B-C bond formation through heterolytic C-C or C-H bond cleavages or B-H bond formation through heterolytic C-H cleavage occur dependent on the structure of the hexane isomer. Meanwhile, we observe identical CID spectra of adducts originating from isomers of C6H12. Spectroscopic investigations of adducts of 1-hexene and cyclohexane indicate the same product structure with an open C6 chain. Computational investigations evidenced that low lying transition states are present, which enable a ring opening reaction of cyclohexane when binding to [B12Br11]-.

Rational design of an argon-binding superelectrophilic anion.

Chemically binding to argon (Ar) at room temperature has remained the privilege of the most reactive electrophiles, all of which are cationic (or even dicationic) in nature. Herein, we report a concept for the rational design of anionic superelectrophiles that are composed of a strong electrophilic center firmly embedded in a negatively charged framework of exceptional stability. To validate our concept, we synthesized the percyano-dodecoborate [B12(CN)12]2-, the electronically most stable dianion ever investigated experimentally. It serves as a precursor for the generation of the monoanion [B12(CN)11]-, which indeed spontaneously binds Ar at 298 K. Our mass spectrometric and spectroscopic studies are accompanied by high-level computational investigations including a bonding analysis of the exceptional B-Ar bond. The detection and characterization of this highly reactive, structurally stable anionic superelectrophile starts another chapter in the metal-free activation of particularly inert compounds and elements.

Characterization of the alkali metal oxalates (MC2O4-) and their formation by CO2 reduction via the alkali metal carbonites (MCO2-).

The reduction of carbon dioxide to oxalate has been studied by experimental Collisionally Induced Dissociation (CID) and vibrational characterization of the alkali metal oxalates, supplemented by theoretical electronic structure calculations. The critical step in the reductive process is the coordination of CO2 to an alkali metal anion, forming a metal carbonite MCO2- able to subsequently receive a second CO2 molecule. While the energetic demand for these reactions is generally low, we find that the degree of activation of CO2 in terms of charge transfer and transition state energies is the highest for lithium and systematically decreases down the group (M = Li-Cs). This is correlated to the strength of the binding interaction between the alkali metal and CO2, which can be related to the structure of the oxalate moiety within the product metal complexes evolving from a planar to a staggered conformer with increasing atomic number of the interacting metal. Similar structural changes are observed for crystalline alkali metal oxalates, although the C2O42- moiety is in general more planar in these, a fact that is attributed to the increased number of interacting alkali metal cations compared to the gas-phase ions.

Properties of gaseous closo-[B6X6]2- dianions (X = Cl, Br, I).

Electronic structure, collision-induced dissociation (CID) and bond properties of closo-[B6X6]2- (X = Cl-I) are investigated in direct comparison with their closo-[B12X12]2- analogues. Photoelectron spectroscopy (PES) and theoretical investigations reveal that [B6X6]2- dianions are electronically significantly less stable than the corresponding [B12X12]2- species. Although [B6Cl6]2- is slightly electronically unstable, [B6Br6]2- and [B6I6]2- are intrinsically stable dianions. Consistent with the trend in the electron detachment energy, loss of an electron (e- loss) is observed in CID of [B6X6]2- (X = Cl, Br) but not for [B6I6]2-. Halogenide loss (X- loss) is common for [B6X6]2- (X = Br, I) and [B12X12]2- (X = Cl, Br, I). Meanwhile, X˙ loss is only observed for [B12X12]2- (X = Br, I) species. The calculated reaction enthalpies of the three competing dissociation pathways (e-, X- and X˙ loss) indicated a strong influence of kinetic factors on the observed fragmentation patterns. The repulsive Coulomb barrier (RCB) determines the transition state for the e- and X- losses. A significantly lower RCB for X- loss than for e- loss was found in both experimental and theoretical investigations and can be rationalized by the recently introduced concept of electrophilic anions. The positive reaction enthalpies for X- losses are significantly lower for [B6X6]2- than for [B12X12]2-, while enthalpies for X˙ losses are higher. These observations are consistent with a difference in bond character of the B-X bonds in [B6X6]2- and [B12X12]2-. A complementary bonding analysis using QTAIM, NPA and ELI-D based methods suggests that B-X bonds in [B12X12]2- have a stronger covalent character than in [B6X6]2-, in which X has a stronger halide character.

Properties of perhalogenated {closo-B10} and {closo-B11} multiply charged anions and a critical comparison with {closo-B12} in the gas and the condensed phase.

closo-Borate anions [closo-BnXn]2- are part of the most famous textbook examples of polyhedral compounds. Substantial differences in their reactivity and interactions with other compounds depending on the substituent X and cluster size n have been recognized, which favor specific closo-borates for different applications in cancer treatment, chemical synthesis, and materials science. Surprisingly, a fundamental understanding of the molecular properties underlying these differences is lacking. Here, we report our study comparing the electronic structure and reactivity of closo-borate anions [closo-BnXn]2- (X = Cl, Br, I, n = 10, 11, 12 in all combinations) in the gas phase and in solution. We investigated the free dianions and the ion pairs [nBu4N]+[closo-BnXn]2- by gas phase anion photoelectron spectroscopy accompanied by theoretical investigations. Strong similarities in electronic structures for n = 10 and 11 were observed, while n = 12 clusters were different. A systematic picture of the development in electronic stability along the dimension X is derived. Collision induced dissociation shows that fragmentation of the free dianions is mainly dependent on the substituent X and gives access to a large variety of boron-rich molecular ions. Fragmentation of the ion pair depends strongly on n. The results reflect the high chemical stability of clusters with n = 10 and 12, while those with n = 11 are much more prone to dissociation. We bridge our study to the condensed phase by performing comparative electrochemistry and reactivity studies on closo-borates in solution. The trends found at the molecular level are also reflected in the condensed-phase properties. We discuss how the gas phase values allow evaluation of the influence of the condensed phase on the electronic stability of closo-borates. A synthetic method via an oxidation/chlorination reaction yielding [closo-B10Cl10]2- from highly chlorinated {closo-B11} clusters is introduced, which underlines the intrinsically high reactivity of the {closo-B11} cage.

Gauging Molecular Orientation through Time Domain Simulations of Surface-Enhanced Raman Scattering.

Molecular reorientation dynamics modulate the physical and chemical properties of molecules at interfaces. This is particularly the case for organic molecules interacting with metallic surfaces-structural motifs of current interest to ultrasensitive chemical/biological detection/imaging. In this context, surface-enhanced Raman scattering (SERS) can be used to gauge the orientation of molecules in the immediate vicinity of plasmonic metals. Herein, we analyze SERS spectra of two aromatic thiols (4-mercaptobenzonitrile and thiophenol) chemisorbed onto silver substrates using ab initio molecular dynamics-based Raman spectral simulations that account for the optical response of molecules in the bulk and oriented molecules at the surface. Through a comparison with prior works aimed at gauging the orientation of our two model systems on metallic substrates, we describe the advantages and discuss the limitations of our described models and approach to gauging the orientation of molecules on metals through time domain simulations of their SERS signatures.

A benchmark photoelectron spectroscopic and theoretical study of the electronic stability of [B12H12]2.

We report a joint benchmark study on the electronic stability of closo-dodecaborate [B12H12]2- employing negative ion photoelectron spectroscopy and high level electronic structure methods. The photoelectron spectra of [B12H12]2-, measured at 266, 193, and 157 nm, yield the Adiabatic and Vertical Detachment Energies (ADE and VDE) of this dianion at 0.93 ± 0.05 eV and 1.15 ± 0.05 eV, respectively, along with a ∼3 eV Repulsive Coulomb Barrier (RCB) against electron detachment. Theoretical calculations at various levels of electronic structure theory confirm the high stability of this dianion. The ADE and VDE values calculated at the coupled cluster with single, double and a perturbative estimate of triple excitations/aug-cc-PVQZ level are 0.92 and 1.16 eV, in excellent agreement with the experimental benchmark values. The comparison between the experimental and the theoretical values obtained at different levels of theory indicate that the PBE0 density functional represents a cost-effective method of sufficient accuracy to describe the molecular properties of this dianion and associated compounds. The theoretical RCB was modeled after the electrostatic potential (ESP) and point charge method (PCM) along three different detachment pathways, viz., along the B-H bond, perpendicular to a B-B bond, and normal to a B-B-B triangle. It was found that detachment of the electron along the B-H bond is preferred, as this pathway is associated with RCBs between 2.3 eV (PCM) and 3.3 eV (ESP), values that bracket the experimental estimate of ∼3 eV.

Photoelectron spectroscopy of [Mo6X14]2- dianions (X = Cl-I).

Photoelectron spectroscopy and theoretical investigations have been performed to systematically probe the intrinsic electronic properties of [Mo6X14]2- (X = halogen). All three PE spectra of gaseous [Mo6X14]2- (X = Cl, Br, I) dianions, which were generated by electrospray ionization, exhibit multiple resolved peaks in the recorded binding energy range. Theoretical investigations on the orbital structure and charge distribution were performed to support interpretation of the observed spectra and were further extended onto [Mo6F14]2-, a dianion that was not available for the experimental study. The measured adiabatic (ADE) and vertical detachment energies (VDE) for X = Cl-I were well reproduced by density functional theory calculations (accuracy ∼0.1 eV). Corresponding ADE/VDE values for the dianions were found to be 1.48/2.13 (calc.) and 2.30/2.65, 2.30/2.62, and 2.20/2.42 eV (all expt.) for X = F, Cl, Br, and I, respectively, showing an interesting buckled trend of electron binding energy (EBE) along the halogen series, i.e., EBE (F) ≪ EBE (Cl) ∼ EBE (Br) > EBE (I). Molecular orbital analyses indicate different mixing of metal and halogen atomic orbitals, which is strongly dependent on the nature of X, and suggest that the most loosely bound electrons are detached mainly from the metal core for X = F and Cl, but from halide ligands for X = Br and I. The repulsive Coulomb barrier (RCB), estimated from the photon energy dependent spectra, decreases with increasing halogen size, from 1.8 eV for X = Cl to 1.6 eV for X = I. Electrostatic potential modeling confirms the experimental RCB values and predicts that the most favorable electron detaching pathway should lie via the face-bridging halide ligands.

Time Domain Simulations of Single Molecule Raman Scattering.

Nonequilibrium chemical phenomena are known to play an important role in single molecule microscopy and spectroscopy. Herein, we explore these effects through ab initio molecular dynamics (AIMD)-based Raman spectral simulations. We target an isolated aromatic thiol (thiobenzonitrile, TBN) as a prototypical molecular system. We first show that the essential features contained in the ensemble-averaged Raman spectrum of TBN can be reproduced by averaging over 18 short AIMD trajectories spanning a total simulation time of ∼60 ps. This involved more than 90 000 polarizability calculations at the B3LYP/def2-TZVP level of theory. We then illustrate that the short trajectories (∼3.3 ps total simulation time), where the accessible phase space is not fully sampled, provide a starting point for understanding key features that are often observed in measurements targeting single molecules. Our results suggest that a complete understanding of single molecule Raman scattering needs to account for molecular conformational flexibility and nonequilibrium chemical phenomena in addition to local optical fields and modified selection rules. The former effects are well-captured using the described AIMD-based single molecule Raman spectral simulations.

Electronic Structure and Stability of [B12X12]2- (X = F-At): A Combined Photoelectron Spectroscopic and Theoretical Study.

The stability and electron loss process of numerous multiply charged anions (MCAs) have been traditionally explained in terms of the classical Coulomb interaction between spatially separated charged groups. An understanding of these processes in MCAs with not well-separated excess charges is still lacking. We report the surprising properties and physical behavior of [B12X12]2-, X = F, Cl, Br, I, At, which are MCAs with not well-separated excess charges and cannot be described by the prevailing classical picture. In this series of MCAs, comprising a "boron core" surrounded by a "halogen shell", the sign of the total charge in these two regions changes along the halogen series from X = F-At. With the aid of experimental photoelectron spectroscopy and highly correlated ab initio electronic structure calculations, we demonstrate that the trend in the electronic stability of these MCAs is determined by the interplay between the Coulomb (de)stabilization originating from the "boron core" and "halogen shell" and the extension of the overlap between the orbitals from both regions. The second excess electron is always taken from the most positively charged region, viz., the "boron core" for X = F, Cl, and the surrounding "halogen shell" for X = I, At. This change in the physical behavior is attributed to the position of the highest occupied molecular orbital, which dwells in a region that is spatially separated from the one containing the excess negative charge. The unusual intrinsic electronic structure of the [B12X12]2- MCAs provides the basis for a molecular level understanding of their observed unique physical and chemical properties and a new paradigm for understanding the properties of these MCAs with not well-separated charges that departs from the prevailing model used for spatially separated charges that is based on their classical Coulomb interaction.

Coupled Cluster Studies of Ionization Potentials and Electron Affinities of Single-Walled Carbon Nanotubes.

In this paper, we apply equation-of-motion coupled cluster (EOM-CC) methods in the studies of the vertical ionization potentials (IPs) and electron affinities (EAs) for a series of single-walled carbon nanotubes (SWCNT). The EOM-CC formulations for IPs and EAs employing excitation manifolds spanned by single and double excitations (IP/EA-EOM-CCSD) are used to study the IPs and EAs of the SWCNTs as a function of the nanotube length. Several armchair nanotubes corresponding to C20nH20 models with n = 2-6 have been used in benchmark calculations. In agreement with previous studies, we demonstrate that the electronegativity of C20nH20 systems remains, to a large extent, independent of the nanotube length. We also compare IP/EA-EOM-CCSD results with those obtained with coupled cluster models with single and double excitations corrected by perturbative triples, CCSD(T), and density functional theory (DFT) using global and range-separated hybrid exchange-correlation functionals.

Nonequilibrium Chemical Effects in Single-Molecule SERS Revealed by Ab Initio Molecular Dynamics Simulations.

Recent developments in nanophotonics have paved the way for achieving significant advances in the realm of single-molecule chemical detection, imaging, and dynamics. In particular, surface-enhanced Raman scattering (SERS) is a powerful analytical technique that is now routinely used to identify the chemical identity of single molecules. Understanding how nanoscale physical and chemical processes affect single-molecule SERS spectra and selection rules is a challenging task and is still actively debated. Herein, we explore underappreciated chemical phenomena in ultrasensitive SERS. We observe a fluctuating excited electronic state manifold, governed by the conformational dynamics of a molecule (4,4'-dimercaptostilbene, DMS) interacting with a metallic cluster (Ag20). This affects our simulated single-molecule SERS spectra; the time trajectories of a molecule interacting with its unique local environment dictates the relative intensities of the observable Raman-active vibrational states. Ab initio molecular dynamics of a model Ag20-DMS system are used to illustrate both concepts in light of recent experimental results.

Intense fluorescence of Au20.

Ligand-protected Au clusters are non-bleaching fluorescence markers in bio- and medical applications. Here we show that their fluorescence can be an intrinsic property of the Au cluster itself. We find a very intense and sharp fluorescence peak located at λ=739.2 nm (1.68 eV) for Au20 clusters in a Ne matrix held at 6 K. The fluorescence reflects the Highest Occupied Molecular Orbital-Lowest Unoccupied Molecular Orbital (HOMO-LUMO) diabatic bandgap of the cluster. Au20 shows a very rich absorption fine structure reminiscent of well defined molecule-like quantum levels. These levels are resolved since Au20 has only one stable isomer (tetrahedral); therefore our sample is mono-disperse in cluster size and conformation. Density-functional theory (DFT) and time-dependent DFT calculations clarify the nature of optical absorption and predict both main absorption peaks and intrinsic fluorescence in fair agreement with experiment.

Time-Domain Simulations of Transient Species in Experimentally Relevant Environments.

Simulating the spectroscopic properties of short-lived thermal and photochemical reaction intermediates and products is a challenging task, as these species often feature atypical molecular and electronic structures. The complex environments in which such species typically reside in practice add further complexity to the problem. Herein, we tackle this problem in silico using ab initio molecular dynamics (AIMD) simulations, employing iso-CHBr3, namely H(Br)C-Br-Br, as a prototypical system. This species was chosen because it features both a nonconventional C-Br-Br bonding pattern, as well as a strong dependence of its spectral features on the local environment in which it resides, as illustrated in recent experimental reports. We simulate the UV-vis and IR spectra of iso-CHBr3 in the gas phase, as well as in a Ne cluster (64 atoms) and in a methylcyclohexane cage (14 solvent molecules) representative of the previously characterized matrix isolated and solvated iso-CHBr3 species. We exclusively perform fully quantum mechanical static and dynamic simulations. By comparing our condensed phase simulations to their experimental analogues, we stress the importance of (i) conformational sampling, even at cryogenic temperatures, and (ii) using a fully quantum mechanical description of both solute and bath to properly account for the experimental observables.

Benchmark theoretical study of the π-π binding energy in the benzene dimer.

We establish a new estimate for the binding energy between two benzene molecules in the parallel-displaced (PD) conformation by systematically converging (i) the intra- and intermolecular geometry at the minimum, (ii) the expansion of the orbital basis set, and (iii) the level of electron correlation. The calculations were performed at the second-order Møller-Plesset perturbation (MP2) and the coupled cluster including singles, doubles, and a perturbative estimate of triples replacement [CCSD(T)] levels of electronic structure theory. At both levels of theory, by including results corrected for basis set superposition error (BSSE), we have estimated the complete basis set (CBS) limit by employing the family of Dunning's correlation-consistent polarized valence basis sets. The largest MP2 calculation was performed with the cc-pV6Z basis set (2772 basis functions), whereas the largest CCSD(T) calculation was with the cc-pV5Z basis set (1752 basis functions). The cluster geometries were optimized with basis sets up to quadruple-ζ quality, observing that both its intra- and intermolecular parts have practically converged with the triple-ζ quality sets. The use of converged geometries was found to play an important role for obtaining accurate estimates for the CBS limits. Our results demonstrate that the binding energies with the families of the plain (cc-pVnZ) and augmented (aug-cc-pVnZ) sets converge [within <0.01 kcal/mol for MP2 and <0.15 kcal/mol for CCSD(T)] to the same CBS limit. In addition, the average of the uncorrected and BSSE-corrected binding energies was found to converge to the same CBS limit much faster than either of the two constituents (uncorrected or BSSE-corrected binding energies). Due to the fact that the family of augmented basis sets (especially for the larger sets) causes serious linear dependency problems, the plain basis sets (for which no linear dependencies were found) are deemed as a more efficient and straightforward path for obtaining an accurate CBS limit. We considered extrapolations of the uncorrected (ΔE) and BSSE-corrected (ΔEcp) binding energies, their average value (ΔEave), as well as the average of the latter over the plain and augmented sets (ΔEave) with the cardinal number of the basis set n. Our best estimate of the CCSD(T)/CBS limit for the π-π binding energy in the PD benzene dimer is De = -2.65 ± 0.02 kcal/mol. The best CCSD(T)/cc-pV5Z calculated value is -2.62 kcal/mol, just 0.03 kcal/mol away from the CBS limit. For comparison, the MP2/CBS limit estimate is -5.00 ± 0.01 kcal/mol, demonstrating a 90% overbinding with respect to CCSD(T). The spin-component-scaled (SCS) MP2 variant was found to closely reproduce the CCSD(T) results for each basis set, while scaled opposite spin (SOS) MP2 yielded results that are too low when compared to CCSD(T).

Toward enabling large-scale open-shell equation-of-motion coupled cluster calculations: triplet states of ß-carotene.

In this paper we discuss the application of novel parallel implementation of the coupled cluster (CC) and equation-of-motion coupled cluster methods (EOMCC) in calculations of excitation energies of triplet states in ß-carotene. Calculated excitation energies are compared with experimental data, where available. We also provide a detailed description of the new parallel algorithms for iterative CC and EOMCC models involving singles and doubles excitations.