On the performance of density functional theory for symmetry-breaking problems

Structure of [18]Annulene Revisited: Challenges for Computing Benzenoid Systems

For cyclic conjugated structures, erratic computational results have been obtained with Hartree-Fock (HF) molecular orbital (MO) methods as well as density functional theory (DFT) with large HF-exchange contributions. In this work, the reasons for this unreliability are explored. Extensive computations on [18]annulene and related compounds highlight the pitfalls to be avoided and the due diligence required for such computational investigations. In particular, a careful examination of the MO singlet-stability eigenvalues is recommended. The appearance of negative eigenvalues is not (necessarily) problematic, but near-zero (positive or negative) eigenvalues can lead to dramatic errors in vibrational frequencies and related properties. DFT approaches with a lower HF admixture generally appear more robust in this regard for the description of benzenoid structures, although they may exaggerate the tendency toward planarity and C-C bond-equalization. For the iconic [18]annulene, the results support a nonplanar equilibrium structure. The density-fitted frozen natural orbital coupled-cluster singles and doubles with perturbative triples [DF-FNO CCSD(T)] method of electron correlation with an aug-pVQZ/aug-pVTZ basis set places the C2 global minimum 1.1 kcal mol-1 below the D6h stationary point.

Avoiding spin contamination and spatial symmetry breaking by exact-exchange-only optimized-effective-potential methods within the symmetrized Kohn-Sham framework

For open-shell atoms and molecules, Kohn-Sham (KS) methods typically resort to spin-polarized approaches that exhibit spin-contamination and often break spatial symmetries. As a result, the KS Hamiltonian operator and the KS orbitals do not exhibit the space and spin symmetry of the physical electron system. The KS formalism can be symmetrized in a rigorous way only in real space, only in spin space, or both in real and spin space. Within such symmetrized KS frameworks, we present exact-exchange-only optimized-effective-potential (OEP) methods that are free of spin contamination and/or spatial symmetry breaking. The effect of symmetrizations on the total energy and its parts and on the exchange potential is analyzed. The presented exact-exchange-only OEP methods may serve as a starting point for high-level symmetrized KS methods based, e.g., on the adiabatic-connection fluctuation-dissipation theorem.

Comment on "Theoretical study of the NO3 radical reaction with CH2ClBr, CH2ICl, CH2BrI, CHCl2Br, and CHClBr2" by I. Alkorta, J. M. C. Plane, J. Elguero, J. Z. Dávalos, A. U. Acuña and A. Saiz-Lopez, Phys. Chem. Chem. Phys. 2022, , 14365

This comment addresses a systematic error in the potential energy surfaces of the title reactions presented in the original article by Alkorta et al. The NO₃ radical has D_3h symmetry in the electronic ground state while the M08HX functional employed in the original article predicts an incorrect C_2v geometry and energy. By combining thermodynamic data for the OH + HNO₃ → H₂O + NO₃ reaction with spectroscopic data and results from M08HX calculations on HNO₃, H₂O and the OH radical, the ground state NO₃ radical energy is estimated to be 37 kJ mol^-1 lower than reported for the C_2v geometry.

Theoretical Investigations for Kinetics of the Chemical Reactions: H + SiClx (x = 1, 2, 3)

Hydrogen atoms and SiCl_x (x = 1, 2, 3) radicals coexist during the hydrogenation of silicon tetrachloride (STC, SiCl₄), an important process in the fabrication of industrial polysilicon. In this work, the mechanisms and kinetics of the reactions between H and SiCl_x (x = 1, 2, 3) were studied by theory. The structures and vibrational frequencies of reactants, products, intermediates, and transition states (TSs) were determined at the B2PLYP/may-cc-pVTZ level. The single-point energies of minima and saddle points were refined using the coupled-cluster single-double with triple perturbative (CCSD(T)) with the complete basis set extrapolation method. Some special treatments were designed to obtain reliable wave functions for unimolecular reactions without tight TSs by the density functional theory. Subsequently, Lennard-Jones (L-J) parameters between each intermediate (SiHCl_x) and bath gas (He) were obtained at the MP2/jul-cc-pVTZ level to derive reliable temperature- and pressure-dependent rate coefficients for unimolecular reactions according to the variational Rice-Ramsperger-Kassel-Marcus theory. For bimolecular reactions, rate coefficients were determined by the variational transition-state theory. The rate coefficients of barrierless reactions were derived based on the loose TSs with the maximum free energy. Finally, the master equation analysis was used to investigate the variation of the rate coefficients with pressure and temperature in the activated paths.

Symmetry Dilemma of Doubly Hybrid Density Functionals for Equilibrium Molecular Property Calculations

In electronic structure theory, when an approximate wavefunction tends to artifactually break the symmetry of the exact Hamiltonian, the corresponding method is referred to as having a "symmetry dilemma" problem. Such types of artifacts were often reported when Hartree-Fock (HF) and the low-level post-HF methods were used, while the traditional Kohn-Sham density functional theory (KS-DFT) methods were usually found to be more resistant to this breakdown. In this work, we present a systematic study on the reliability of the doubly hybrid (DH) DFT methods for several violable cases. Almost all the commonly used B2PLYP-type (bDH) functionals are shown to have a severe "symmetry dilemma" problem and yield dramatically unreliable molecular properties, such as dipole moment, vibrational frequency, and static polarizability at the equilibrium geometry. A one-parameter bDH functional model study demonstrates that such a problem is a combined effect of the inappropriate portion of the HF exchange (over 50%) for the self-consistent field (SCF) calculation and the augmentation of the second-order perturbative contribution. It is remarkable that the XYG3-type (xDH) functionals show a good capability to resist the artifactual symmetry breaking and yield reliable molecular properties when the same critical cases are calculated. In the xDH, there are two functionals of different purposes, namely, the SCF functional and the energy functional, which have different amounts of the HF exchange and different portions of the correlation contributions. The success of the xDHs can be attributed to this flexibility in xDH construction to avoid using an improperly large portion of the HF exchange in the SCF functional. The insights gained in this work are of significance for the development of an improved DH functional.

Dissociation of HCl in water nanoclusters: an energy decomposition analysis perspective

As known, small HCl-water nanoclusters display a particular dissociation behaviour, whereby at least four water molecules are required for the ionic dissociation of HCl. In this work, we examine how intermolecular interactions promote the ionic dissociation of such nanoclusters. To this end, a set of 45 HCl-water nanoclusters with up to four water molecules is introduced. Energy decomposition analysis based on absolutely localized molecular orbitals (ALMO-EDA) is employed in order to study the importance of frozen interaction, dispersion, polarization, and charge-transfer for the dissociation. The vertical ALMO-EDA scheme is applied to HCl-water clusters along a proton-transfer coordinate varying the amount of spectator water molecules. The corresponding ALMO-EDA results show a clear preference for the dissociated cluster only in the presence of four water molecules. Our analysis of adiabatic ALMO-EDA results reveals a push-pull mechanism for the destabilization of the HCl bond based on the synergy between forward and backward charge-transfer.

Probing Interfacial Effects on Ionization Energies: The Surprising Banality of Anion-Water Hydrogen Bonding at the Air/Water Interface

Liquid microjet photoelectron spectroscopy is an increasingly common technique to measure vertical ionization energies (VIEs) of aqueous solutes, but the interpretation of these experiments is subject to questions regarding sensitivity to bulk versus interfacial solvation environments. We have computed aqueous-phase VIEs for a set of inorganic anions, using a combination of molecular dynamics simulations and electronic structure calculations, with results that are in excellent agreement with experiment regardless of whether the simulation data are restricted to ions at the air/water interface or to those in bulk aqueous solution. Although the computed VIEs are sensitive to ion-water hydrogen bonding, we find that the short-range solvation structure is sufficiently similar in both environments that it proves impossible to discriminate between the two on the basis of the VIE, a conclusion that has important implications for the interpretation of liquid-phase photoelectron spectroscopy. More generally, analysis of the simulation data suggests that the surface activity of soft anions is largely a second or third solvation shell effect, arising from disruption of water-water hydrogen bonds and not from significant changes in first-shell anion-water hydrogen bonding.

Exploring spin symmetry-breaking effects for static field ionization of atoms: Is there an analog to the Coulson-Fischer point in bond dissociation?

Löwdin's symmetry dilemma is an ubiquitous issue in approximate quantum chemistry. In the context of Hartree-Fock (HF) theory, the use of Slater determinants with some imposed constraints to preserve symmetries of the exact problem may lead to physically unreasonable potential energy surfaces. On the other hand, lifting these constraints leads to the so-called broken symmetry solutions that usually provide better energetics, at the cost of losing information about good quantum numbers that describe the state of the system. This behavior has previously been extensively studied in the context of bond dissociation. This paper studies the behavior of different classes of HF spin polarized solutions (restricted, unrestricted, and generalized) in the context of ionization by strong static electric fields. We find that, for simple two electron systems, unrestricted Hartree-Fock (UHF) is able to provide a qualitatively good description of states involved during the ionization process (neutral, singly ionized, and doubly ionized states), whereas RHF fails to describe the singly ionized state. For more complex systems, even though UHF is able to capture some of the expected characteristics of the ionized states, it is constrained to a single M_s (diabatic) manifold in the energy surface as a function of field intensity. In this case, a better qualitative picture can be painted by using generalized Hartree-Fock as it is able to explore different spin manifolds and follow the lowest solution due to lack of collinearity constraints on the spin quantization axis.

Revealing the nature of electron correlation in transition metal complexes with symmetry breaking and chemical intuition

In this work, we provide a nuanced view of electron correlation in the context of transition metal complexes, reconciling computational characterization via spin and spatial symmetry breaking in single-reference methods with qualitative concepts from ligand-field and molecular orbital theories. These insights provide the tools to reliably diagnose the multi-reference character, and our analysis reveals that while strong (i.e., static) correlation can be found in linear molecules (e.g., diatomics) and weakly bound and antiferromagnetically coupled (monometal-noninnocent ligand or multi-metal) complexes, it is rarely found in the ground-states of mono-transition-metal complexes. This leads to a picture of static correlation that is no more complex for transition metals than it is, e.g., for organic biradicaloids. In contrast, the ability of organometallic species to form more complex interactions, involving both ligand-to-metal σ-donation and metal-to-ligand π-backdonation, places a larger burden on a theory's treatment of dynamic correlation. We hypothesize that chemical bonds in which inter-electron pair correlation is non-negligible cannot be adequately described by theories using MP2 correlation energies and indeed find large errors vs experiment for carbonyl-dissociation energies from double-hybrid density functionals. A theory's description of dynamic correlation (and to a less important extent, delocalization error), which affects relative spin-state energetics and thus spin symmetry breaking, is found to govern the efficacy of its use to diagnose static correlation.

The spin-polarized edge states of blue phosphorene nanoribbons induced by electric field and electron doping

Edge states of various two-dimensional materials such as graphene are intrinsically spin-polarized. In other materials, electric field and charge doping are required for introducing magnetism to their edges. In this work, by using first-principles calculations, we studied the effects of transverse electric field on the edge states of the armchair blue phosphorene nanoribbon (ABPNR), and found that a transverse electric field drives the edge electronic state occupied and at the same time spin-polarized. We also doped electrons to the ABPNR and found that these additional electrons occupy and spin-polarize the electronic states of both edges of the nanoribbon.

Accurate quantum dynamics simulation of the photodetachment spectrum of the nitrate anion (NO3 -) based on an artificial neural network diabatic potential model

The photodetachment spectrum of the nitrate anion (NO₃ ^-) is simulated from first principles using wavepacket quantum dynamics propagation and a newly developed accurate full-dimensional fully coupled five state diabatic potential model. This model utilizes the recently proposed complete nuclear permutation inversion invariant artificial neural network diabatization technique [D. M. G. Williams and W. Eisfeld, J. Phys. Chem. A 124, 7608 (2020)]. The quantum dynamics simulations are designed such that temperature effects and the impact of near threshold detachment are taken into account. Thus, the two available experiments at high temperature and at cryogenic temperature using the slow electron velocity-map imaging technique can be reproduced in very good agreement. These results clearly show the relevance of hot bands and vibronic coupling between the X̃ ²A₂ ^' ground state and the B̃ ²E' excited state of the neutral radical. This together with the recent experiment at low temperature gives further support for the proper assignment of the ν₃ fundamental, which has been debated for many years. An assignment of a not yet discussed hot band line is also proposed.

Polishing the Gold Standard: The Role of Orbital Choice in CCSD(T) Vibrational Frequency Prediction

While CCSD(T) with spin-restricted Hartree-Fock (RHF) orbitals has long been lauded for its ability to accurately describe closed-shell interactions, the performance of CCSD(T) on open-shell species is much more erratic, especially when using a spin-unrestricted HF (UHF) reference. Previous studies have shown improved treatment of open-shell systems when a non-HF set of molecular orbitals, like Brueckner or Kohn-Sham density functional theory (DFT) orbitals, is used as a reference. Inspired by the success of regularized orbital-optimized second-order Møller-Plesset perturbation theory (κ-OOMP2) orbitals as reference orbitals for MP3, we investigate the use of κ-OOMP2 orbitals and various DFT orbitals as reference orbitals for CCSD(T) calculations of the corrected ground-state harmonic vibrational frequencies of a set of 36 closed-shell (29 neutrals, 6 cations, 1 anion) and 59 open-shell diatomic species (38 neutrals, 15 cations, 6 anions). The aug-cc-pwCVTZ basis set is used for all calculations. The use of κ-OOMP2 orbitals in this context alleviates difficult cases observed for both UHF orbitals and OOMP2 orbitals. Removing two multireference systems and 12 systems with ambiguous experimental data leaves a pruned data set. Overall performance on the pruned data set highlights CCSD(T) with a B97 orbital reference (CCSD(T):B97), CCSD(T) with a κ-OOMP2 orbital reference (CCSD(T):κ-OOMP2), and CCSD(T) with a B97M-rV orbital reference (CCSD(T):B97M-rV) with RMSDs of 8.48 cm^-1, and 8.50 cm^-1, and 8.75 cm^-1 respectively, outperforming CCSD(T):UHF by nearly a factor of 5. Moreover, the performance on the closed- and open-shell subsets shows these methods are able to treat open-shell and closed-shell systems with comparable accuracy and robustness. CCSD(T) with RHF orbitals is seen to improve upon UHF for the closed-shell species, while spatial symmetry breaking in a number of restricted open-shell HF (ROHF) references leads CCSD(T) with ROHF reference orbitals to exhibit the poorest statistical performance of all methods surveyed for open-shell species. The use of κ-OOMP2 orbitals has also proven useful in diagnosing multireference character that can hinder the reliability of CCSD(T).

Third-Order Møller-Plesset Theory Made More Useful? The Role of Density Functional Theory Orbitals

The practical utility of Møller-Plesset (MP) perturbation theory is severely constrained by the use of Hartree-Fock (HF) orbitals. It has recently been shown that the use of regularized orbital-optimized MP2 orbitals and scaling of MP3 energy could lead to a significant reduction in MP3 error [Bertels, L. W.; J. Phys. Chem. Lett. 2019, 10, 4170 4176]. In this work, we examine whether density functional theory (DFT)-optimized orbitals can be similarly employed to improve the performance of MP theory at both the MP2 and MP3 levels. We find that the use of DFT orbitals leads to significantly improved performance for prediction of thermochemistry, barrier heights, noncovalent interactions, and dipole moments relative to the standard HF-based MP theory. Indeed, MP3 (with or without scaling) with DFT orbitals is found to surpass the accuracy of coupled-cluster singles and doubles (CCSD) for several data sets. We also found that the results are not particularly functional sensitive in most cases (although range-separated hybrid functionals with low delocalization error perform the best). MP3 based on DFT orbitals thus appears to be an efficient, noniterative O(N⁶) scaling wave-function approach for single-reference electronic structure computations. Scaled MP2 with DFT orbitals is also found to be quite accurate in many cases, although modern double hybrid functionals are likely to be considerably more accurate.

SAMPL7: Host-guest binding prediction by molecular dynamics and quantum mechanics

Statistical Assessment of Modeling of Proteins and Ligands (SAMPL) challenges provide routes to compare chemical quantities determined using computational chemistry approaches to experimental measurements that are shared after the competition. For this effort, several computational methods have been used to calculate the binding energies of Octa Acid (OA) and exo-Octa Acid (exoOA) host-guest systems for SAMPL7. The initial poses for molecular dynamics (MD) were generated by molecular docking. Binding free energy calculations were performed using molecular mechanics combined with Poisson-Boltzmann or generalized Born surface area solvation (MMPBSA/MMGBSA) approaches. The factors that affect the utility of the MMPBSA/MMGBSA approaches including solvation, partial charge, and solute entropy models were also analyzed. In addition to MD calculations, quantum mechanics (QM) calculations were performed using several different density functional theory (DFT) approaches. From SAMPL6 results, B3PW91-D3 was found to overestimate binding energies though it was effective for geometry optimizations, so it was considered for the DFT geometry optimizations in the current study, with single-point energy calculations carried out with B2PLYP-D3 with double-, triple-, and quadruple-ζ level basis sets. Accounting for dispersion effects, and solvation models was deemed essential for the predictions. MMGBSA and MMPBSA correlated better to experiment when used in conjunction with an empirical/linear correction.

The effects of strain and electric field on the half-metallicity of pristine and O-H/C-N-decorated zigzag graphene nanoribbons

In zigzag graphene nanoribbons (ZGNRs), the spin polarized edge states play a significant role in the electronic structure. The two ferromagnetically ordered edges anti-ferromagnetically coupled with each other, which would result in the half-metallicity under electric field. Given that the strain, external electric field, and edge decorations are the main means of tuning the magnetism and electronic property of one-dimentional materials. It motivates us to study the combine effects on ZGNRs of these methods. So, in present work, the corporate influences of the tensile strain, transverse electric field, and asymmetric edge decoration by  -OH and  -CN groups on the magnetism and electronic property of 8-ZGNR have been studied using the density functional theory. The calculational results indicate that the arising strain can modulate the response of electronic and magnetic properties to external electric field, improving the magnetism and extending the electric field range in which the ZGNR presents half-metallicity. In addition, the O-H/C-N groups decorated ZGNR possesses a lower critic electric field and a larger electric field range for realizing half-metallicity comparing with the unstrained pristine ZGNR.

Distinguishing artificial and essential symmetry breaking in a single determinant: approach and application to the C60, C36, and C20 fullerenes

The existence of a generalized Hartree–Fock solution in C₆₀ has led to controversy on whether C₆₀ is polyradicaloid (or strongly correlated). We attempt to end the controversy with κ-OOMP2 which removes the illusion of this artificial symmetry breaking. We conclude that C₆₀ is not strongly correlated.

Practical Density Functionals beyond the Overdelocalization-Underbinding Zero-Sum Game

Density functional theory (DFT) uses a density functional approximation (DFA) to add electron correlation to mean-field electronic structure calculations. Standard strategies (generalized gradient approximations GGAs, meta-GGAs, hybrids, etc.) for building DFAs, no matter whether based on exact constraints or empirical parametrization, all face a zero-sum game between overdelocalization (fractional charge error, FC) and underestimation of covalent bonding (fractional spin error, FS). This work presents an alternative strategy. Practical "Rung 3.5" ingredients are used to implement insights from hyper-GGA DFAs that reduce both FS and FC errors. Prototypes of this strategy qualitatively improve FS and FC error over 40 years of standard DFAs while maintaining low cost and practical evaluation of properties. Numerical results ranging from transition metal thermochemistry to absorbance peaks and excited-state geometry optimizations highlight this strategy's promise and indicate areas requiring further development.

Stable Radical Cations and Their π-Dimers Prepared from Ethylene- and Propylene-3,4-dioxythiophene Co-oligomers: Combined Experimental and Theoretical Investigations

Co-oligomers composed of two 3,4-ethylenedioxythiophene (EDOT) units and two or three 3,4-propylenedioxythiophene (ProDOT) units, i.e., 2E2P_Et and 2E3P_Et, were newly synthesized together with the ProDOT trimer 3P_Me. On the basis of cyclic voltammetry, the gaps between the first and second oxidation potentials (ΔE^1-2) of 2E2P_Et and 2E3P_Et were found to be larger than that of the previously synthesized ProDOT tetramer 4P_Hex. These co-oligomers gave the fairly stable radical cations 2E2P_Et^•+ and 2E3P_Et^•+ by chemical oxidation with AgSbF₆. The disproportionation of 2E2P_Et^•+ and 2E3P_Et^•+ into neutral and dicationic species, which was observed for 4P_Hex^•+, was inhibited in accord with the larger ΔE^1-2. Additionally, the formation of the π-dimers (3P_Me)₂²⁺, (2E2P_Et)₂²⁺, and (2E3P_E)₂²⁺ was clearly observed in dichloromethane solution at low temperatures with UV-vis-NIR spectroscopy. Furthermore, the π-dimerization enthalpies of 2E2P_Et^•+ and 2E3P_Et^•+ were greater than that of 3P_Me^•+, suggesting the formation of fully π-contacted structures. The structures of the π-dimers were optimized at the B97D3 method, and the calculated absorption spectra of the π-dimers obtained using TD-DFT methods were in reasonable agreement with the observed ones, supporting the reliability of the calculated structures.

Computational exploration of regioselectivity and atmospheric lifetime in NO3-initiated reactions of CH3OCH3 and CH3OCH2CH3

The NO₃-initiated reactions of CH₃OCH₃ and CH₃OCH₂CH₃ have been investigated by the BHandHLYP method in conjunction with the 6-311G(d,p) basis set. Thermodynamic and kinetic data are further refined using the comparatively accurate CCSD(T) method. According to the values of reaction enthalpies (ΔH_r,298^θ) and reaction Gibbs free energies (ΔG_r,298^θ) from CH₃OCH₂CH₃ with NO₃ system, we find that H-abstraction pathway from the α-CH₂ group is more exothermic. It is further confirmed by the calculated CH bond dissociation energy of CH₃OCH₂CH₃ molecule. All the rate constants, computed through means of canonical variational transition state with small-curvature tunneling correction, are fitted to the three-parameter expressions k₁=1.54×10^-23T^3.34exp(-1035.53/T) and k₂=3.55×10^-26T^4.31exp(-281.24/T)cm³molecule^-1s^-1 and branching ratios are computed over the temperature range 200-600K. The branching ratios are also discussed. The atmospheric lifetimes of CH₃OCH₃ and CH₃OCH₂CH₃ determined by the NO₃ radical are about 270 and 29days, respectively.

Thermoelectric Properties of Solution-Processed n-Doped Ladder-Type Conducting Polymers

Ladder-type "torsion-free" conducting polymers (e.g., polybenzimidazobenzophenanthroline (BBL)) can outperform "structurally distorted" donor-acceptor polymers (e.g., P(NDI2OD-T2)), in terms of conductivity and thermoelectric power factor. The polaron delocalization length is larger in BBL than in P(NDI2OD-T2), resulting in a higher measured polaron mobility. Structure-function relationships are drawn, setting material-design guidelines for the next generation of conducting thermoelectric polymers.

Polarons in Narrow Band Gap Polymers Probed over the Entire Infrared Range: A Joint Experimental and Theoretical Investigation

We investigate the photoinduced absorption (PIA) spectra of the prototypical donor-acceptor polymer [2,6-(4,4-bis(2-ethylhexyl)-4H-cyclopenta[2,1-b;3,4-b']dithiophene)-alt-4,7-(2,1,3-benzothiadiazole)] (C-PCPDTBT) and its silicon bridged variant Si-PCPDTBT over a spectral range from 0.07 to 1.5 eV. Comparison between time-dependent density functional theory simulations of the electronic and vibrational transitions of singlet excitons, triplet excitons, polarons, and bipolarons with the experimental results proves that the observed features are due to positive polarons delocalized on the polymer chains. We find that the more crystalline Si-bridged variant gives rise to a red-shift in the transition energies, especially in the mid-infrared (MIR) spectral range and furthermore observe that the pristine polymers' responses depend on the excitation energy. Blending with PCBM, on the other hand, leads to excitation-independent PIA spectra. By computing the response properties of molecular aggregates, we show that polarons are delocalized in not only the intra- but also the interchain direction, leading to intermolecular transitions which correspond well to experimental absorption features at the lowest energies.

CuNO2 and Cu(+)NO2 Revisited: A Comparative ab Initio and DFT Study

We have reinvestigated CuNO2 and Cu(+)NO2 at ab initio as well as at pure and hybrid DFT levels of approximation employing large ANO basis sets. The systems were fully optimized using the CCSD(T), QCISD(T), BPW91, PBE, PBE0, and B3LYP methods. Several stationary points (minima and transition structures) were found on the related potential energy surfaces (PES). The C2v bidentate η(2)-O,O isomer is calculated to be the most stable species on the CuNO2 PES, followed by two monodentate isomers [Formula: see text] the Cs η(1)-O and C2v η(1)-N species which are higher in energy by 12 and 14 kcal/mol, respectively, at CCSD(T)/Basis-II (where Basis-II is 21s15p10d6f4g/8s7p5d3f2g for Cu; 14s9p4d3f/5s4p3d2f for O and N). On the Cu(+)NO2 PES, the Cs monodentate η(1)-O trans (0 kcal/mol) and cis (+3 kcal/mol at CCSD(T)/Basis-II) isomers are found, followed by the C2v monodentate η(1)-N isomer (+14 kcal/mol at CCSD/Basis-II). In contrast to the pure DFT, the hybrid DFT methods perform reasonably well for predicting the relative stabilities (except for η(1)-N of CuNO2) and structures; however, their predictions of the bond dissociation energies are less reliable (for CuNO2 the difference is as much as 10 kcal/mol compared to the CCSD(T) values). The performance of the QCISD(T) method was analyzed, and, furthermore, the issue of symmetry breaking was investigated.

Even-odd product variation of the C(n)(+) + D(2) (n = 4-9) reaction: complexity of the linear carbon cation electronic states

We have studied reactions between linear Cn(+) (n = 4-9) and D2, using ion mobility mass spectrometry techniques and quantum chemical calculations in order to understand the complex reactivity of the linear cluster cations. Only linear CnD(+) products were observed for the odd (n = 5, 7, 9) linear clusters, while CnD2(+) was the main product for the even clusters. For the reaction rate constants determined for these two channels, we obtained the following two features: (1) the rate constant decreases with the size n, and (2) even-sized clusters have lower rate constants than neighboring odd-sized clusters. In the theoretical calculations using the CCSD(T) and B3LYP methods with the cc-pVTZ basis, we found that a low lying (2)Σ state in odd clusters may play an important role in these reactions. This opposes the previous interpretation that the (2)Πg/u state is the dominant electronic state for linear Cn(+) (n = 4-9) clusters. We showed that a barrierless radical abstraction forming CnD(+) occurs through a direct head on approach for the (2)Σ state Cn(+). In contrast, a carbene-like insertion forming CnD2(+) occurs through a sideways approach for the (2)Πg/u state Cn(+). We have concluded that the higher rate constants for the odd clusters come from the existence of symmetry broken (2)Σ states which are absent in even linear clusters.

Infrared vibrational and electronic transitions in the dibenzopolyacene family

We report experimental spectra in the mid-infrared (IR) and near-IR for a series of dibenzoacenes isolated in Ar matrices. The experiments are supported by Density Functional Theory (DFT) and Time-Dependent DFT (TD-DFT) calculations with both vibrational and electronic transitions studied. For the neutrals, we find good agreement between the experimental and B3LYP and BP86 results for all species studied. The band at about 1440 cm(-1) carries more intensity than in typical PAHs and increases in intensity with the size of the dibenzoacene molecule. For the ions the B3LYP approach fails to yield reasonable IR spectra for most systems and the BP86 approach is used. Electronic transitions dominate the vibrational bands in the mid-IR region for the large dibenzoacene ions. In spite of the very strong electronic transitions, there is still reasonable agreement between theory and experiment for the vibrational band positions. The experimental and theoretical results for the dibenzoacenes are also compared with those for the polyacenes.

Atmospheric formation of the NO3 radical from gas-phase reaction of HNO3 acid with the NH2 radical: proton-coupled electron-transfer versus hydrogen atom transfer mechanisms

The amidogen radical abstracts the hydrogen from nitric acid through a proton coupled electron transfer mechanism rather than by an hydrogen atom transfer process.