Gas phase basicities of polyfunctional molecules. Part 5: Non-aromatic sp2 nitrogen containing compounds.

This paper constitutes the fifth part of a general review of the gas-phase protonation thermochemistry of polyfunctional molecules (Part 1: Theory and methods, Mass Spectrom Rev 2007, 26:775-835, Part 2: Saturated basic sites, Mass Spectrom Rev 2012, 31:353-390, Part 3: Amino acids, Mass Spectrom Rev 2012, 31:391-435, Part 4: Carbonyl as basic site, Mass Spectrom Rev 2015, 34:493-534). This part is devoted to non-aromatic molecules characterized by a lone pair located on a sp2 nitrogen atom, it embraces functional groups such as imines, amidines, guanidines, diazenes, hydrazines, oximes, and phosphazenes. Specific examples are examined under five major chapters. In the first one, aliphatic and unsaturated (conjugated and cyclic) imines, hydrazones, and oximes are considered. A second chapter describes the protonation energetic of aliphatic, conjugated, or cyclic amidines. Guanidines, polyguanides, and biomolecules containing guanidine were examined in the third chapter. A fourth chapter describes the particular case of the phosphazene molecules. Finally, diazenes and azides were considered in the last chapter. Experimental data were re-evaluated according to the presently adopted basicity scale, i.e., PA(NH3 ) = 853.6 kJ/mol, GB (NH3 ) = 819 kJ/mol. Structural and energetic information given by G4MP2 quantum chemistry computations on typical systems are presented. © 2016 Wiley Periodicals, Inc. Mass Spec Rev 37:139-170, 2018.

The Antiferromagnetic Spin Coupling in Non-Kekulé Acenes-Impressive Polyradical Character Revealed by High-Level Multireference Methods.

Complete active space (CASSCF) and multireference (MR-CISD(Q) and MR-AQCC) calculations were performed for non-Kekulé analogues of acenes, dimethylenepolycyclobutadienes, with lengths of up to eight cyclobutadiene (CBD) units. Multireference calculations predict that the most stable energy state of the system is either triplet (if there is an odd number of CBD units) or singlet (if there is an even number of CBD units) due to antiferromagnetic spin coupling, which thus violates Hund's rule in larger molecules. We also show an impressive polyradical character in the system that increases with the size of the molecule, as witnessed by more than eleven unpaired electrons in the singlet state of the molecule with eight CBD units. Together with the small energy gap between singlet and higher multiplicity energy states even above the triplet state, this demonstrates the exceptional polyradical properties of these π-conjugated oligomeric chains.

QM/MM non-adiabatic decay dynamics of formamide in polar and non-polar solvents.

Non-adiabatic on-the-fly dynamics simulations of the photodynamics of formamide in water and n-hexane were performed using a QM/MM approach. It was shown that steric restrictions imposed by the solvent cage do not have an influence on the initial motion which leads to the lowest energy conical intersection seam. The initial deactivation in water is faster than in n-hexane and in the gas phase. However, most of the formamide molecules in water do not reach the ground state. The reason for the deactivation inefficiency in water is traced back to a decrease of close CO···HOH and NH···OH(2) contacts which fall in the range of hydrogen bonds. The energy deposition into H-bond breaking events leaves molecules with less energy for surmounting the CN dissociation barrier. In both solvents, after hopping to the ground state, the solvent cage keeps the HCO and NH(2) fragments or CO and NH(3) products in close proximity. Consequently, the number of trajectories where fast recombination happens is augmented with delayed recombinations that start when the dissociation fragments hit the cage wall and return back. The hot ground state formamide is formed in an internal conversion process identical to the path leading to CN photodissociation. In the case of aqueous formamide, good agreement with experimental results is achieved by combining dynamics simulations starting from the S(1) and the S(2) excited states collecting high and low energy trajectories, respectively.

Desymmetrisation of aromatic diamines and synthesis of non-symmetrical thiourea derivatives by click-mechanochemistry.

ortho- and para-Phenylenediamines were desymmetrised and quantitatively transformed into mono- and bis-(thio)ureas or mixed thiourea-ureas through a one-pot mechanochemical click reaction sequence; mechanochemical desymmetrisation proceeds quantitatively without excess reagents and allows the controlled extension of a molecular structure by combining normally competing reactions.

Click mechanochemistry: quantitative synthesis of "ready to use" chiral organocatalysts by efficient two-fold thiourea coupling to vicinal diamines.

Mechanochemical methods of neat grinding and liquid-assisted grinding have been applied to the synthesis of mono- and bis(thiourea)s by using the click coupling of aromatic and aliphatic diamines with aromatic isothiocyanates. The ability to modify the reaction conditions allowed the optimization of each reaction, leading to the quantitative formation of chiral bis(thiourea)s with known uses as organocatalysts or anion sensors. Quantitative reaction yields, combined with the fact that mechanochemical reaction conditions avoid the use of bulk solvents, enabled solution-based purification methods (such as chromatography or recrystallization) to be completely avoided. Importantly, by using selected model reactions, we also show that the described mechanochemical reaction procedures can be readily scaled up to at least the one-gram scale. In that way, mechanochemical synthesis provides a facile method to fully transform valuable enantiomerically pure reagents into useful products that can immediately be applied in their designed purpose. This was demonstrated by using some of the mechanochemically prepared reagents as organocatalysts in a model Morita-Baylis-Hillman reaction and as cyanide ion sensors in organic solvents. The use of electronically and sterically hindered ortho-phenylenediamine revealed that mechanochemical reaction conditions can be readily optimized to form either the 1:1 or the 1:2 click-coupling product, demonstrating that reaction stoichiometry can be more efficiently controlled under these conditions than in solution-based syntheses. In this way, it was shown that excellent stoichiometric control by mechanochemistry, previously established for mechanochemical syntheses of cocrystals and coordination polymers, can also be achieved in the context of covalent-bond formation.

Thermal reaction of [3,4]-benzo-8-substituted-3Z,5Z,7E-octatetraenes and quantum-chemical study of the (8π,6π)-electrocyclisation.

The first example of thermal (8π,6π)-electrocyclisation of 1,3,5,7-octatetraene with one double bond embedded in an aromatic moiety is described. By this process, [3,4]-benzo-8-substituted octatetraene derivatives, the cis,trans-1-(o-vinylphenyl)-4-(R = Me, Ph, 2-furyl)buta-1,3-dienes were transformed to a new endo-7-(R = Me, Ph, 2-furyl) and exo-7-(R = Me)-2,3-benzobicyclo[4.2.0]octa-2,4-dienes. Mechanism of reaction was also studied by DFT quantum-chemical calculations. The M06/6-311+G(d,p)//M06/6-31+G(d,p) calculations indicate that formation of the single endo-isomer in the case of phenyl and 2-furyl substituents is determined by higher activation barriers for exo-6π-electrocyclisation than for 8π-cycloreversion.

Photoinduced rearrangement of aromatic N-chloroamides to chloroaromatic amides in the solid state: inverted Π(N)-Σ(N) occupational stability of amidyl radicals.

We report a solid-state photochemical rearrangement reaction by which aromatic N-chloroamides exposed to UV light or sunlight are rapidly and efficiently converted to chloroaromatic amides. The course, the intermediate (nascent chlorine vs dichlorine) and the outcome of the reaction depend on the excitation (exposure time, wavelength, and intensity) and on inherent structural factors (the directing role of the substituents and, as demonstrated by the different reactivity of two polymorphs of N-chlorobenzanilide, the supramolecular structure). The photolysis of the chloroamides provides facile photochemical access to arylamidyl radicals as intermediates, which in the absence of strong hydrogen bond donors are stabilized in the reactant crystals by C-H/N-Cl···π interactions, thus, providing insight into their structure and chemistry. Thorough theoretical modeling of the factors determinant to the stability and the nature of the spin-hosting orbital evidenced that although the trans-Π(||) state (Np spin) of the amidyls is normally preferred over the trans-Σ(⊥) configuration (Nsp(2) spin), stabilization by aromatic conjugation, steric and geometry factors, as well as by electronic effects from the substituents can decrease the Π-Σ gap in these intermediates significantly, resulting in similar and, in the case of the orthogonal amide-phenyl disposition, even reversed population of the unpaired electron in the two orbitals. Quantitative correlation established that the inverted occupational spin stability and the Π(N)-Σ(N) crossover are collectively facilitated by the conformation, valence angle, and disposition of the amide group relative to the aromatic system. The stabilization and detection of a trans-Σ(⊥) radical was experimentally accomplished by steric locking of the orthogonal trans-amide conformation with double ortho-tert-butyl substitution at the phenyl ring. The effects of the single para-phenyl substituents on the relative occupational stability of the arylamidyl radical states point out to non-Hammett behavior. By including cumulative electronic effects from multiple substitutions, four distinct families of the aromatic amidyl radicals were identified. The Π(∥) state is the most stable structure of the N-phenylacetamidyl radical and of most of the substituted arylamidyls, although the Σ(⊥) and Π(⊥) states can also be stabilized by introducing tert-butyl and nitro groups, respectively.

Towards an environmentally-friendly laboratory: dimensionality and reactivity in the mechanosynthesis of metal-organic compounds.

We present a proof-of-principle study of an environmentally-friendly approach to laboratory research, in which the synthesis and structural characterisation of metal-organic complexes and frameworks are achieved without using bulk solvents; our study addresses the use of heteroditopic ligands for manipulating the dimensionality of metal-organic materials and describes how kinetic obstacles in such mechanosynthesis can be overcome.

Matrix-controlled photofragmentation of formamide: dynamics simulation in argon by nonadiabatic QM/MM method.

The short-time photodynamics (2 ps) of formamide embedded into an Ar matrix starting from the low-lying singlet excited S(1) (n(0)π*) and S(2) (ππ*) states were explored using a nonadiabatic photodynamics QM/MM approach. The interaction between formamide and the Ar matrix is taken into account at the MM level by means of Lennard-Jones potentials. This is the first example of exploring photodissociation of formamide with full nonadiabatic dynamics in a matrix and it nicely illustrates importance of considering environmental effects on photodissociation behavior of the peptide bond. It is shown that embedding of the formamide molecule in the argon matrix has strong impact on the outcome of the process. This is illustrated by formation of the 1:1 complex between ammonia and CO and prevention of full separation of the NH(2)˙ and HCO˙ subunits in the NH(2)˙ + HCO˙ radical pair. In addition, the argon matrix strongly influences the lifetime of the S(1) state, which increases by 211 fs relative to the gas phase.

The isomerization barrier in cyanocyclobutadienes: an ab initio multireference average quadratic coupled cluster study.

The energy profiles of the isomerization of mono, di-, and tetracyano-substituted cyclobutadienes (CBDs) are computed at the multireference average quadratic coupled cluster/complete active space self-consistent field level of theory. It was found that the energy barrier heights for the automerization reaction are 2.6 (tetracyano-CBD), 5.1 (1,3-dicyano-CBD), and 6.4 (cyano-CBD) kcal mol(-1), implying that they are lowered relative to that in the parent CBD (6.4 kcal mol(-1)), the monosubstituted derivative being an exception. Since the free CBD shuttles between two equivalent structures even at low temperature of 10 K, it follows that bond-stretch isomerism does not take place in cyanocyclobutadienes. Instead, these compounds exhibit rapid fluxional interconversion at room temperature between two bond-stretch isomers by the double bond flipping mechanism. The reason behind the decrease in the barrier heights is identified as a slightly enhanced resonance effect at the saddle points separating two (equivalent) bond-stretch isomers, compared to that in the equilibrium structures, predominantly due to the diradical character of the former. It is also shown that the energy gap between the singlet ground state saddle point structure and the first triplet equilibrium geometry decreases upon multiple substitution by the cyano groups. The splitting of the S and T energy is small being within the range of 6.5-8.2 kcal mol(-1).

Study of the mechanism of the N-CO photodissociation in N,N-dimethylformamide by direct trajectory surface hopping simulations.

Photodynamics of N,N-dimethylformamide (DMF) in its low-lying singlet excited n(O)-pi* and pi-pi* states have been explored by the direct trajectory surface hopping method based on multiconfigurational ab initio calculations. The dynamics simulations starting in the first excited singlet state (n(O)-pi*) showed that in 57% of trajectories, S(1) excited DMF decays to the ground via the crossing seam related to the N-CO bond stretching MXS1. In 41% of all trajectories, the relaxation process on the S(1) PES moves the molecule into the minimum, where it stays trapped until the end of simulation time. In simulations starting in the second excited state, all trajectories are found to deactivate through MXS5 (S(2)/S(1)) by very fast N-CO stretching. Because the N-CO dissociation process continues in the S(1) potential energy surface, most of the population overshoots the MXS1 and leads to fully dissociated electronically excited HCO and N(CH(3))(2) radicals. A mechanism for their deactivation by H-C-O and C-N-C bending modes is proposed.

Generation and dissociation pathways of singly and doubly protonated bisguanidines in the gas phase.

Para-bisguanidinyl benzene 1 and its N-permethylated derivative 2 are both sufficiently strong bases to afford not only the monocations [1+H]+ and [2+H]+, but also the doubly protonated ions, [1+2H]2+ and [2+2H]2+, in the gas phase. The title ions generated via electrospray ionization are probed by collision-induced dissociation experiments which inter alia reveal that the dicationic species [1+2H]2+ and [2+2H]2+ can even undergo fragmentation reactions with maintenance of the 2-fold charge. Complementary results from density functional theory predict PAs above 1000 kJ mol(-1) for the neutral compounds, i.e., PA(1) = 1025 kJ mol(-1) and PA(2) = 1067 kJ mol(-1). Due to the stabilization of the positive charge in the guanidinium ions and the para-phenylene spacer separating the basic sites, even the monocations bear sizable proton affinities, i.e., PA([1+H]+) = 740 kJ mol(-1) and PA([2+H]+) = 816 kJ mol(-1).

Simulation of the photodeactivation of formamide in the nO-pi* and pi-pi* states: an ab initio on-the-fly surface-hopping dynamics study.

The short-time photodynamics (1 ps) of formamide in its low-lying singlet excited n(O)-pi(*) and pi-pi(*) states have been investigated by the direct trajectory surface-hopping method based on multiconfigurational ab initio calculations. The simulations showed that in both states, the primary deactivation process is C-N bond dissociation. In the ground state, the energy is transferred to (a) translational motion of the HCO and NH(2) fragments, (b) additional C-H dissociation from the vibrationally hot HCO fragment, or (c) formation of NH(3) and CO. In addition to the C-N dissociation pathway, C-O bond fission is found to be an additional primary deactivation path in the pi-pi(*) dynamics. From fractional occupations of trajectories, lifetimes of formamide were estimated: tau(S(1))=441 fs and tau(S(2))=66 fs.

The nonadiabatic deactivation paths of pyrrole.

Multireference configuration interaction (MRCI) calculations have been performed for pyrrole with the aim of providing an explanation for the experimentally observed photochemical deactivation processes. Potential energy curves and minima on the crossing seam were determined using the analytic MRCI gradient and nonadiabatic coupling features of the COLUMBUS program system. A new deactivation mechanism based on an out-of-plane ring deformation is presented. This mechanism directly couples the charge transfer 1pipi* and ground states. It may be responsible for more than 50% of the observed photofragments of pipi*-excited pyrrole. The ring deformation mechanism should act complementary to the previously proposed NH-stretching mechanism, thus offering a more complete interpretation of the pyrrole photodynamics.

Automerization reaction of cyclobutadiene and its barrier height: an ab initio benchmark multireference average-quadratic coupled cluster study.

The problem of the double bond flipping interconversion of the two equivalent ground state structures of cyclobutadiene (CBD) is addressed at the multireference average-quadratic coupled cluster level of theory, which is capable of optimizing the structural parameters of the ground, transition, and excited states on an equal footing. The barrier height involving both the electronic and zero-point vibrational energy contributions is 6.3 kcal mol(-1), which is higher than the best earlier theoretical estimate of 4.0 kcal mol(-1). This result is confirmed by including into the reference space the orbitals of the CC sigma bonds beyond the standard pi orbital space. It places the present value into the middle of the range of the measured data (1.6-10 kcal mol(-1)). An adiabatic singlet-triplet energy gap of 7.4 kcal mol(-1) between the transition state (1)B(tg) and the first triplet (3)A(2g) state is obtained. A low barrier height for the CBD automerization and a small DeltaE((3)A(2g),(1)B(1g)) gap bear some relevance on the highly pronounced reactivity of CBD, which is briefly discussed.

Periodic trends in bond dissociation energies. A theoretical study.

Bond dissociation energies (BDEs) of all possible A-X single bonds involving the first- and second-row atoms, from Li to Cl, where the free valences are saturated by hydrogens, have been estimated through the use of the G3-theory and at the B3LYP/6-311+G(3df,2pd)//B3LYP/6-31G(2df,p) DFT level of theory. BDEs exhibit a periodical behavior. The A-X (A = Li, Be, B, Na, Mg, Al, and Si) BDEs show a steady increase along the first and the second row of the periodic table as a function of the atomic number Z(X). For A-X bonds involving electronegative atoms (A = C, N, O, F, P, S, and Cl) the bond energies achieve a maximum around Z(X) = 5. The same behavior is observed when BDEs are plotted against the electronegativity chi(X) of the atom X. Thus, for A-X bonds (A = Li, Be, B, Na, Mg, Al, Si), the BDEs for a fixed A increases, grosso modo, as the electronegativity differences between X and A increase, with some exceptions, which reflect the differences in the relaxation energies of the radicals produced upon the bond cleavage. A similar trend, albeit less pronounced, is found for single A-X bonds, where A = C, N, O, F, P, S, and Cl. However, there is an additional feature embodied in the enhancement of the strength of the A-boron bonds due to the ability of boron to act as a strong electron acceptor. The trends in bond lengths and charge densities at the bond critical points are in line with the aforementioned behavior.

On the bond-stretch isomerism in the benzo[1,2:4,5]dicyclobutadiene system--an ab initio MR-AQCC study.

Bond-stretch isomerism in benzo[1,2:4,5]dicyclobutadienle (BDCB) has been investigated using the MR-AQCC/6-31G(d) method, a high-level multireference ab initio approach including size-extensivity corrections. The applied theoretical approach includes both nondynamical and dynamical electron correlation effects. Full MR-AQCC geometry optimizations of localized (1) and delocalized (3) isomers as well as the transition structure (TS) have been determined using D2h, symmetry restriction. The calculations show that both isomers are approximately of equal stability separated by a barrier with a height of about 5 kcal mol(-1). Thus, the present results strongly indicate that benzof[1,3:4,5]dicyclobutadiene is a very good candidate for an organic compound exhibiting bond-stretch isomerism, since isomers 1 and 3 correspond to true minima on the double-well potential energy surface, which are separated by a sufficiently high barrier. It is particularly important to emphasize that isomer 3 represents a realization of the highly elusive quasi-[10]annulene.

Benzocyclopropenyl Anion: A Stable 8pi-Electron Species.

The conjugate base of benzocyclopropene has been generated in the gas phase. Its reactivity and thermodynamic stability were explored. The measured acidity is DeltaH degrees (acid)(benzocyclopropene) = 386 +/- 3 kcal/mol, and the electron affinity of benzocyclopropenyl radical is 0.51 eV < EA < 1.11 eV. Ab initio calculations satisfactorily reproduce the experimental results and provide additional insights. Benzocyclopropene is found to be 34.5 kcal/mol more acidic than the allylic position of cyclopropene and only 4 +/- 3 kcal/mol less acidic than toluene. These findings are explained in terms of the structure and electronic properties of benzocyclopropenyl anion.