Search

The dispersive optical activity of two saturated cyclic amines, (R)-2-methylpyrrolidine (R-2MPY) and (S)-2-methylpiperidine (S-2MPI), has been interrogated under isolated and solvated conditions to elucidate the roles of large-amplitude motion associated with nitrogen-center inversion and ring-puckering dynamics. Experimental optical rotatory dispersion profiles were almost mirror images of one another and displayed parallel solvent dependencies. Quantum-chemical analyses built on density-functional and coupled-cluster methods revealed four low-lying conformers for each molecule, which are distinguished by axial/equatorial orientations of their amino hydrogens and methyl substituents. Chiroptical signatures predicted for these species were combined through an independent-conformer ansatz to simulate the ensemble-averaged response, with a polarizable continuum model (PCM) being used to treat implicit solute-solvent interactions. The intrinsic behavior observed for isolated (gaseous) R-2MPY and S-2MPI was reproduced best by merging coupled-cluster (CCSD) estimates of rotatory powers with thermal population fractions deduced from complete basis set (CBS-APNO) free-energy calculations. Although prior claims of sizable chiroptical contributions arising from helically twisted (chiral) heterocyclic frameworks could be discounted, less satisfactory agreement between experiment and theory was realized for solution phases. Response properties sustained modest isomer-dependent changes in the presence of PCM solvation, but the corresponding energy metrics showed systematic trends, whereby structures having larger electric-dipole moments were stabilized preferentially in media of high polarity. Despite the fact that R-2MPY conformations were predicted to undergo a progressive reordering of their relative energies across the six solvents of interest, S-2MPI was found to exhibit more pronounced solvent-induced perturbations at long wavelengths (viz., in regions far removed from electronic resonances). Experimental results are discussed in terms of the distinct ring-puckering mechanisms for R-2MPY and S-2MPI, which are expected to be dominated by hindered pseudorotation among envelope/twist motifs and semi-inversion between chairlike antipodes, respectively.

The ultraviolet (UV) photodissociation dynamics of the jet-cooled cyclohexyl (c-C6H11) radical is studied using the high-n Rydberg atom time-of-flight (HRTOF) technique. The cyclohexyl radical is produced by the 193 nm photodissociation of chlorocyclohexane and bromocyclohexane and is examined in the photolysis wavelength region of 232-262 nm. The H-atom photofragment yield (PFY) spectrum contains a broad peak centered at 250 nm, which is in good agreement with the UV absorption spectrum of the cyclohexyl radical and assigned to the 3p Rydberg states. The translational energy distributions of the H-atom loss product channel, P(ET)'s, are bimodal, with a slow (low ET) component peaking at ∼6 to 7 kcal/mol and a fast (high ET) component peaking at ∼44-48 kcal/mol. The fraction of the average translational energy in the total excess energy, ⟨fT⟩, is in the range of 0.16-0.25 in the photolysis wavelength region of 232-262 nm. The H-atom product angular distribution of the slow component is isotropic, while that of the fast component is anisotropic with an anisotropy parameter of β ≈ 0.5-0.7. The bimodal product translational energy and angular distributions indicate two dissociation pathways to the H + C6H10 products in cyclohexyl. The high-ET anisotropic component is from a repulsive, prompt dissociation on a repulsive potential energy surface coupling with the Rydberg excited states to produce H + cyclohexene. The low-ET isotropic component is consistent with the unimolecular dissociation of hot radical on the ground electronic state after internal conversion from the Rydberg states. The similarity of the photodissociation dynamics of the cyclohexyl radical to the previously studied small linear and branched alkyls expands on the understanding of the dissociation dynamics of alkyl radicals to include larger cyclic alkyl radicals.

Deciphering the exact electronic and geometric changes of photoexcited molecules is an important task not only to understand the fundamental atomistic mechanisms but also to rationally design molecular properties and functions. Here, we present a combined experimental and theoretical study of the twisted intramolecular charge transfer (TICT) process in hemithioindigo photoswitches. Using ultrafast transient IR spectroscopy as the main analytical method, a detailed understanding of the extent and direction of charge transfer within the excited molecule is obtained. At the same time, the geometrical distortion is monitored directly via changes of indicative vibrational modes over the time course of the photoreaction. These high-resolution data deliver a detailed molecular movie of the TICT process in this important class of chromophores with picosecond time resolution.

Due to the importance of predicting static and dynamic polarizabilities, the performance of various correlated linear response methods including random phase approximation (RPA), RPA(D), higher-order random phase approximation (HRPA), HRPA(D), second-order polarization propagator approximation (SOPPA), SOPPA(CC2), SOPPA(CCSD), CC2, and CCSD has been evaluated against CCSD(T) (static case) and CCSD (dynamic cases) for the T145 set of 145 organic molecules. The benchmark reveals that the HRPA(D) method has the best performance for both static and dynamic polarizabilities apart from CCSD. RPA(D) ranks second for the dynamic cases and third for the static case. Using coupled-cluster amplitudes in SOPPA(CCSD) and SOPPA(CC2), the SOPPA results are significantly improved. The HRPA method has the largest deviations from the reference values for both cases. In general, according to the performance and computational cost of the methods, the HRPA(D) and RPA(D) methods are proposed for calculations of static and dynamic polarizabilities of this and similar sets of molecules.

The challenge of assigning the absolute stereochemical configuration to a chiral compound can be overcome via accurate ab initio predictions of optical rotation, a sensitive molecular property that is further complicated by solvent effects. The solvent's "chiral imprint"-the transfer of the chirality from the solute to the surrounding achiral solvent-is explored here using conformational averaging and time-dependent density-functional theory. These complex solvent effects are taken into account via simple averaging over a molecular dynamics trajectory together with the explicit quantum mechanical consideration of the solvent molecules within the solute's cybotactic region and implicit modeling of the bulk solvent. We consider several axes along which the system's optical rotation varies, including the sampling of the dynamical trajectory, the quality of the one-electron basis set, and the use of continuum solvent models to account for bulk effects.

High-resolution photoelectron (PE) spectra of liquid methanol and ethanol were measured using a liquid microjet and He IIα radiation (40.813 eV). The vertical ionization energy and the ionization threshold were determined as 9.70 ± 0.07 and 8.69 ± 0.07 eV for methanol and 9.52 ± 0.07 and 8.52 ± 0.07 eV for ethanol, respectively. Individual photoemission bands observed for the liquids are well correlated with those in PE spectra of the gaseous samples also measured in the present study, except that the liquid band positions were shifted on average by −1.23 eV for methanol and −1.10 eV for ethanol as compared to the gas. The 5a′ and 7a′ bands of liquid methanol exhibit specifically larger broadening than other bands, for which we attempted spectral fitting with two components, similarly with the case of the 3a1 band of liquid water. PE spectra of both liquid and gaseous ethanol are congested partly due to the presence of the trans and gauche isomers; however, the overall band positions are generally in good agreement with predictions based on quantum chemical calculations. Comparison of the measured PE spectra with experimental and simulated X-ray emission spectra indicate that spectral differences in the lowest ionization band of both methanol and ethanol originate from involvement of nuclear dynamics in the X-ray emission process.

Multichannel thermal decomposition reactions of n-propyl radicals, 1-pentyl radicals, and toluene are investigated by solving a two-dimensional master equation formulated as a function of total energy (E) and angular momentum (J). The primary aim of this study is to elucidate the role of angular momentum in the kinetics of multichannel unimolecular reactions. The collisional transition processes of the reactants colliding with argon are characterized based on the classical trajectory calculations and implemented in the master equation. The rate constants calculated by using the two-dimensional master equation are compared with those of one-dimensional master equations. The consequence of the explicit treatment of angular momentum depends on the J dependence of the microscopic rate constants and is particularly emphasized in the thermal decomposition of toluene, for which the C-H and C-C bond fission channels are considered. The centrifugal effect is insignificant in the energetically favored C-H bond fission but is substantial in the energetically higher C-C bond fission, which causes rotational channel switching of the microscopic rate constants. The proper treatment of the J-dependent channel coupling effect, weak collisional transfer of J, and initial-J-dependent collisional energy transfer are found to be essential for predicting the branching fractions at low pressures.

An extended Debye-Huckel theory with fourth order gradient term is developed for electrolyte solutions; namely, the electric potential φ(r) of the bulk electrolyte solution can be described by ∇2φ(r) = κ2φ(r) + LQ2∇4φ(r), where the parameters κ and LQ are chosen to reproduce the first two roots of the dielectric response function of the bulk solution. Three boundary conditions for solving the electric potential problem are proposed based upon the continuity conditions of involving functions at the dielectric boundary, with which a boundary element method for the electric potential of a solute with a general geometrical shape and charge distribution is derived. Solutions for the electric potential of a spherical ion and a diatomic molecule are found and used to calculate their electrostatic solvation energies. The validity of the theory is successfully demonstrated when applied to binary as well as multicomponent primitive models of electrolyte solutions.

The multiconfiguration time-dependent Hartree (MCTDH) method and its generalization, the multilayer MCTDH (ML-MCTDH), result in equations of motion (EOMs) that are singular when there are virtual orbitals-the unoccupied single-particle functions-in the wave function expansion. For decades this singularity had been numerically removed by regularizing the reduced density matrix. In this Perspective we discuss our recent proposal to regularize the coefficient tensor instead, which has significant impact on both the efficiency and correctness of the EOMs in MCTDH and ML-MCTDH for challenging problems. We further demonstrate that when the system becomes large such that it is necessary to use ML-MCTDH with many layers, it is much more important to employ this new regularization scheme. We illustrate this point by studying a spin-boson model with a large bath that contains up to 100 000 modes. We show that even in the weak coupling regime the new regularization scheme is required to quickly rotate the virtual orbitals into the correct directions in Hilbert space. We argue that this situation can be common for applying a time-dependent tensor network approach to any large enough system.

Bright-state models are often applied to "deperturb" Fermi-coupled bands in molecular vibrational spectra, in cases where a harmonically forbidden transition "borrows" intensity from an energetically nearby allowed transition. However, forbidden transitions can also acquire intensity through anharmonic couplings on the potential energy surface ("mechanical anharmonicity") or dipole moment surface ("electrical anharmonicity") that are not accounted for within the bright-state model. In this work, we compare deperturbation shifts obtained by analysis of experimental data with those predicted using the bright-state model, for a series of discrete encapsulated chloride hydrate isotopomers. Predicted band center shifts and Fermi coupling matrix elements obtained using the bright-state model are larger than those estimated from experimental data.

Electron capture by the σ* LUMO of isoxazole triggers the dissociation of the O-N bond and the opening of the ring. Photodetachment of the resulting anion accesses a neutral structure, in which the O· and ·N bond fragments interact through the intact remainder of the molecular ring and via a 3 A gap created by the bond dissociation. These through-bond and through-space interactions result in a dense manifold of diradical states, including (in the order of increasing energy) a triplet, an open-shell singlet, a closed-shell singlet, and another triplet state. We report photoelectron spectra that reflect partially resolved signatures of these states. Remarkably, the structure of the isoxazole diradical manifold is qualitatively different from that of the analogous system in oxazole. The distinct properties of the two manifolds are explained by using a coupled-fragments molecular-orbital model. Consistent with the past conclusions [Culberson et al. Phys. Chem. Chem. Phys.2014, 16, 3964-3972], the lingering through-space interactions between the O· and ·C bond fragments in ring-open oxazole are responsible for the relative stabilization of the closed-shell singlet state, which correlates with the ground-state cyclic structure. In contrast, the placement of the N atom in the terminal position within the ring-open structure of isoxazole is the key factor leading to the near degeneracy of the π and σ* orbitals, favoring a triplet-state configuration.

In this study, an electrostatic potential (ESP) fitting method using constrained spatial electron density (cSED) expanded with preorthogonal natural atomic orbitals (pNAOs) was proposed. In this method, the electron density of a molecule is divided into spherical atom-centered electron densities and the expansion coefficient is determined to reproduce the ESP around the molecule. Our method was then applied to two systems: (i) a hydration reaction of cis-platin and (ii) a variety of organic/inorganic molecules. By evaluating the atomic charges along the hydration reaction, our method showed good conformational transferability, which cannot be obtained using conventional ESP fitting methods. Moreover, we successfully obtained the hydration structure along the reaction by coupling our method with a reference interaction site model (RISM). Reasonable data were obtained not only for organic molecules but also for inorganic molecules. This success came from the introduction of pNAOs as auxiliary basis sets in the charge fitting.

The spin-orbit charge transfer intersystem crossing (SOCT-ISC) and the formation of a long-lived charge transfer (CT) state were studied with a series of 4,4'-dimethoxy triphenylamine-BODIPY compact electron donor/acceptor dyads. Different torsion freedoms were applied in the dyads to tune the electronic coupling between the donor and acceptor, and a red-shifted CT absorption band was observed for one dyad. The dyads show solvent polarity-dependent singlet oxygen photosensitizing ability (quantum yields 3%-79%). Nanosecond transient absorption spectra of the dyad in nonpolar solvent confirm the formation of triplet states. The intrinsic triplet state lifetime is up to 383 μs (in fluid solution), which is much longer than that accessed with the heavy atom effect (276 μs). Intermolecular triplet photosensitizing of the dyads in a polar solvent produces a long-lived 3CT state (lifetime, τCT = 8.0 μs supported by the electron spin density surface analysis). The triplet state lifetime of the dyads doped in a Clear Flex 50 polymer film is exceptionally long (7.6-11.4 ms), and formation of a long-lived CT state (37 μs) was observed. Triplet-triplet annihilation upconversion was performed with the electron donor/acceptor dyads used as the triplet photosensitizer and perylene used as the triplet acceptor; the upconversion quantum yield is up to 15.8%.

Some organic pollutants in snowpacks undergo faster photodegradation than in solution. One possible explanation for such effect is that their UV-visible absorption spectra are shifted toward lower energy when the molecules are adsorbed at the air-ice interface. However, such bathochromic shift is difficult to measure experimentally. Here, we employ a multiscale/multimodel approach that combines classical and first-principles molecular dynamics, quantum chemical methods, and statistical learning to compute the light absorption spectra of two phenolic molecules in different solvation environments at the relevant thermodynamic conditions. Our calculations provide an accurate estimate of the bathochromic shift of the lowest-energy UV-visible absorption band when these molecules are adsorbed at the air-ice interface, and they shed light into its molecular origin.

The periodically oscillating electromagnetic potential of a photon can, in an electric-dipole transition, "shine" an electron from an anion's bound-state orbital directly into a continuum-state orbital. This occurs in photoelectron and photodetachment spectroscopy, both of which provide much information about the electronic structure of the anion. Alternatively, a molecular anion containing sufficient vibrational energy to "shake/rattle" an electron out of a bound-state orbital can induce electron detachment via a vibration-to-electronic nonadiabatic transition. In this case, the electron binding energy in the anion must be smaller than the vibrational energy-level spacing, so these processes involve anion states of low binding energy, and they eject electrons having low kinetic energy. If the anion's electron binding energy is even smaller, it is possible for a rotation-to-electronic energy transfer to "roll" an electron from the bound-state orbital into the continuum. For each of these mechanisms by which electron detachment can occur, there are different selection rules governing the angular distribution in which the electrons are ejected, and this manuscript discusses these rules, their origins, and their utility when using spectroscopic tools to probe the anion's electronic structure. Several examples of the shine-, shake/rattle-, and role-ejection of electrons from a range of experimental conditions are discussed as are similarities and differences among the corresponding selection rules. Of special novelty are the effects arising when electron ejection occurs from orbitals having very low electron binding energies and thus large radial extent.

Previously proposed theoretical schemes for estimating one-electron excitation energies using Kohn-Sham (KS) solutions with long-range corrected (LC) functionals are applied to the charge-transfer (CT) excitations of the ethylene···tetrafluoroethylene (C

Ab initio molecular dynamics simulations are presented to investigate the intramolecular vibrational energy redistribution (IVR) of an isolated nitromethane molecule. A number of IVR processes are simulated by monitoring the kinetic energy of vibrational modes under selective low-lying vibrational excitations from their ground states (Δν = 1 or 2). Evolution of the normal-mode kinetic energy gives the ultrafast energy transfer processes from parent modes to daughter modes intuitively. From the ultrafast vibrational transfer made by Fourier transformation of the time-dependent normal-mode kinetic energy, we can capture that the symmetry of the normal modes plays an important role in the anharmonic coupling between the vibrational modes. The results show three symmetry-dependent coupling mechanisms: direct symmetric coupling, overtone-assisted coupling, and rotation-assisted coupling. Furthermore, the calculated efficiencies of IVR also coincide with these mechanisms.

Al atoms generally adopt the +3 oxidation state and form stoichiometric oxides such as Al2O3 in the bulk phase. Among small cationic gas-phase clusters, near-stoichiometric clusters such as Al3O4+, Al3O5+, Al4O6+, Al4O7+, Al5O7+, and Al5O8+ have been readily generated in experimental studies. However, when a single Au atom was included in the clusters, oxygen-deficient clusters such as AuAl4O5+ were formed in high abundance; in these clusters, the Au atom accepted electron density from the Al atoms. The geometrical structures and atomic charges in the clusters suggest that a single Au atom can substitute for O atoms in Al oxide clusters. This propensity originates from the high electron and low oxygen affinities, which, together, constitute an unusual property of Au.

We have detected the H2O·DMA and H2O·TMA (DMA, dimethylamine; TMA, trimethylamine) bimolecular complexes at room temperature in the gas phase using Fourier transform infrared spectroscopy. For both complexes, five vibrational bands associated with the H2O molecule are observed and assigned. Within a reduced dimensional local mode framework, we set up a six-dimensional model, including the three H2O vibrational modes and three of the six intermolecular modes, all described with internal curvilinear coordinates. The single points on the potential energy surface and Eckart corrected dipole moment surface are calculated with the CCSD­(T)-F12a/cc-pVDZ-F12 method. Combining the measured and calculated transition intensities, we determine the Gibbs energy of complex formation of both complexes from each of the observed bands. The multiple determinations give similar Gibbs energies, for each complex, and increase the confidence in the combined experimental and theoretical approach, and improve the accuracy of the determined Gibbs energies. The average Gibbs energies of complex formation are found to be 5.0 ± 0.2 and 3.8 ± 0.2 kJ/mol for H2O·DMA and H2O·TMA, respectively. In addition to the experimental uncertainty, there is a potential error on the calculated intensities corresponding to 0.4 kJ/mol. However, the small spread among the four determinations suggests that this error is even less. The Gibbs energies of these complexes serve as accurate benchmarks for theoretical approaches that are prevalent in hydrogen bonding and nucleation studies.

We tested the influence of various parameters on the new particle formation rate predicted for the sulfuric acid–ammonia system using quantum chemistry and cluster distribution dynamics simulations, in our case, Atmospheric Cluster Dynamics Code (ACDC). We found that consistent consideration of the rotational symmetry number of monomers (sulfuric acid and ammonia molecules, and bisulfate and ammonium ions) leads to a significant rise in the predicted particle formation rate, whereas inclusion of the rotational symmetry number of the clusters only changes the results slightly, and only in conditions where charged clusters dominate the particle formation rate. This is because most of the clusters stable enough to participate in new particle formation have a rotational symmetry number of 1, and few exceptions to this rule are positively charged clusters. In contrast, the application of the quasi-harmonic correction for low-frequency vibrational modes tends to generally decrease predicted new particle formation rates and also significantly alters the slope of the formation rate curve plotted against the sulfuric acid concentration, which is a typical convention in atmospheric aerosol science. The impact of the maximum size of the clusters explicitly included in the simulations depends on the simulated conditions. The errors arising from a limited set of clusters are higher for higher evaporation rates, and thus tend to increase with temperature. Similarly, the errors tend to be higher for lower vapor concentrations. The boundary conditions for outgrowing clusters (that are counted as formed particles) have only a small influence on the results, provided that the definition is chemically reasonable and that the set of simulated clusters is sufficiently large. A comparison with data from the Cosmics Leaving OUtdoor Droplets (CLOUD) chamber and a cluster distribution dynamics model using older quantum chemistry input data shows improved agreement when using our new input data and the proposed combination of symmetry and quasi-harmonic corrections.

The vibrational partition function is calculated using the classical method of integration over the whole phase space. The calculations were done for the ground electronic state of a carbon monoxide molecule. The main focus is on temperature in the range 5000-20 000 K, which is common in hypersonic flows of gases and plasmas. The method presented here, because of the exclusion of the noninteracting part of canonical partition function according to the ideas of T.L. Hill, is applicable at temperatures of tens of thousands of Kelvins, where the standard expression for the vibrational partition function fails. At lower temperatures (here 1000-6000 K), the correct quantum results can be obtained with the help of Wigner-Kirkwood expansion. The influence of vibrations on the rotational partition function by bond-length elongation is examined, and the results are compared with the exact ro-vibrational partition function.

The spatial magnetic properties, through-space NMR shieldings (TSNMRS), of bent allene 1, the corresponding C-extended 1,3-butadiene derivative 2, and a number of related compounds 3 -20 have been calculated using the gauge-independent atomic orbital perturbation method, employing the nucleus-independent chemical shift concept and visualized as isochemical shielding surfaces of various sizes and directions. Prior to that, both structures and C-13 chemical shifts were calculated and compared with available experimental bond lengths and delta(C-13)/ppm values (also, as a quality criterion for the computed structures). Bond lengths, the delta(C-13)/ppm, and the TSNMRS values are employed to qualify and quantify the electronic structure of the studied compounds in terms of dative or classical electron-sharing bonds.

We investigate a few density functional theory-based reactivity indices of chemistry, with a view to arrive at an intercomparison and also consider their applications toward the problems of chemical significance. In particular, we propose to use the concepts of fugality and atom-atom polarizability to study the acidic strength of para-substituted benzoic acid derivatives. The nature of the variations and trends in the correlation of reactivity parameters and pKa values is shown to provide an insight into the applicability of these concepts to such reactions.

International audience; Calculating the spectroscopic properties of complex conjugated organic molecules in their relaxed state is far from simple. An additional complexity arises for flexible molecules in solution, where the rotational energy barriers are low enough so that non-minimum conformations may become dynamically populated. These metastable conformations quickly relax during the minimization procedures preliminary to DFT calculations, and so accounting for their contribution to the experimentally-observed properties is problematic. We describe a strategy for stabilising these non-minimum conformations in silico, allowing their properties to be calculated. Diadinoxanthin and alloxanthin present atypical vibrational properties in solution, indicating the presence of several conformations. Performing energy calculations in vacuo and polarizable continuum model calculations in different solvents, we found three different conformations with values for the δ dihedral angle of the end-ring ca. 0°, 180° and 90° with respect to the plane of the conjugated chain. The latter conformation, a non-global minimum, is not stable during the minimization necessary for modelling its spectroscopic properties. To circumvent this classical problem we used a Car-Parinello MD supermolecular approach, in which diadinoxanthin was solvated by water molecules so that metastable conformations were stabilized by hydrogen bonding interactions. We progressively removed the number of solvating waters to find the minimum required for this stabilisation. This strategy represents the first modelling of a carotenoid in a distorted conformation, and provides an accurate interpretation of the experimental data.

We employ and combine a number of recent developments in vibrational structure methods to push their current size limitations toward molecules with tens of modes and showcase their availability for the maleimide molecule. In particular, we assess the use of different rectilinear vibrational coordinates, namely, normal coordinates, hybrid optimized and localized coordinates, and flexible adaptation of local coordinates of nuclei coordinates. These different coordinate parameterizations are employed in conjunction with the adaptive density-guided approach to generate potential energy surfaces (PESs). A screening procedure is furthermore introduced, which provides estimates of the importance of individual terms in the PES, resulting in significant reductions in the computational cost of the PES construction. We find that all three sets of coordinates provide approximately the same level of accuracy in vibrational structure calculations and report fundamental excitation energies with a mean absolute deviation of less than 12 cm-1 when compared to experimental data. We expect that similar accuracy in vibrational structure calculations can be achieved for molecules of larger size using the proposed procedures.

Using density functional theory calculations, we have analyzed nonlinear optical properties of a series of T-graphene quantum dots differing in their shape and size. Electronic polarizability and first-order and second-order hyperpolarizability of these systems are investigated and shed light on their stability and electronic properties. Negative cohesive energy shows that they are energetically stable. The effect of size and incident frequency on their nonlinear responses are comprehensively discussed. Most of the systems exhibit a strong NLO response, and it is enhanced in the presence of an external field. All these systems show absorption maximum ranging from UV to visible window. Overall, this theoretical framework highlighted the nonlinear optical properties of T-graphene quantum dots that may provide valuable information in designing potential NLO materials.

By bridging graphene and benzene through a well-defined sequence of polycyclic aromatic hydrocarbons and their inherent shell structure, J. Phys. Chem. C, 2018, 122, 17526, it is shown that graphene is actually a coherent arrangement of interwoven benzene molecules, which are coordinated by aromaticity, shell structure and topology, all interrelated. The exotic properties of graphene are in fact macroscopic manifestations of aromaticity. At the molecular level this is revealed as a dynamical flipping of the atomic pz-orbitals belonging to different sublattices, leading to a dynamical interchange between aromatic and non-aromatic rings, associated with frontier occupied-unoccupied orbitals or valence-conduction bands interchange with opposite parities, in full agreement with many-body and topological insulators theory. Such dynamical interchange or coupling is driven by inversion symmetry-breaking, activated through geometrical frustration between sublattice and space group symmetry. This, for rectangular nanographenes, in contrast to hexagonal, leads to topological gapless edge states, which, contrary to opposite claims, are not spin polarized., 5 main figures 26 pages 6666 words

The reaction dynamics of excited-state intramolecular proton transfer (ESIPT) of 2,2'-dihydroxyazobenzene (2,2'-DHAB) was investigated by means of white-light supercontinuum femtosecond transient absorption spectroscopy. A coherent in-phase oscillation was observed in the entire wavelength range where stimulated emission of the photoproduct is dominant. This result indicates that the transition strength of the product state is dynamically modulated by a nuclear wavepacket motion (non-Condon effect). The observed vibration was assigned to the mode which modulates the distance between oxygen and hydrogen atoms. By integrating the result of time-dependent density functional theory calculation, the origin of the non-Condon effect was attributed to a dynamical change of configuration interaction between enol and keto characters along the vibrational coordinate, indicating that this vibration is strongly related to the reaction coordinate of ESIPT.

In this study we examine the accuracy of domain-based local pair natural orbital coupled cluster theory with single, double, and perturbative triple excitations (DLPNO-CCSD(T)) on a large benchmark data set. To this end, we use the recently published GMTKN55 superset of molecules that contains 1505 relative energies and 2462 single-point calculations. To our knowledge this is the most comprehensive benchmark evaluation of any highly correlated wave function based ab initio method to date. In the first part of the study, canonical CCSD(T) reference calculations were carried out on the entire test set in order to guarantee that the reference data are of uniform quality. Second, DLPNO-CCSD(T) calculations were carried out under identical conditions. The main finding is that with the exception of two data sets, all data sets have a MAD of 0.4 kcal/mol or less and the majority of sets have a MAD of less than 0.2 kcal/mol. For open shells, the accuracy of the DLPNO calculations was significantly improved through an iterative version of the triples correction.

The far-infrared spectrum of the weakly OH···π hydrogen-bonded benzene-water complex has been studied in neon and argon matrices, below 30 K. The in-plane water libration has been observed in both neon and argon for H2O and D2O complexed with C6H6 and C6D6 but not for the corresponding complexes involving HDO. Both H2O and D2O can tunnel between the two possible hydrogen bonds. This is not possible for HDO. The reported far-infrared observations have implications for the interpretation of the previously obtained molecular beam microwave spectrum of the benzene-water complex.

The effect of sulfur hydrogen bonding on the conformational equilibrium of methyl 3-mercaptopropionate was investigated using microwave spectroscopy in a supersonic jet expansion. The two most stable conformers (I and II) were assigned in the rotational spectra, and complex splitting patterns owing to the methyl internal rotation and SH tunneling motion were resolved and analyzed in detail. For both conformers, the experimental torsional barriers for the methyl top are similar and about 5.1 kJ mol(-1), revealing that their geometrical differences do not affect the methyl internal rotation. The experimentally derived rotational and centrifugal distortion constants, along with the methyl internal rotation barriers, are discussed and compared with results from density functional theory and ab initio calculations. Quantum theory of atoms in molecules, noncovalent interactions, and natural bond orbital analyses show that the global minimum geometry (I), which has the thiol hydrogen oriented toward the carbonyl of the ester, is stabilized by an SH center dot center dot center dot O=C hydrogen bond. The presence of a hydrogen bond is confirmed by the derivation of an accurate experimental geometry that reveals a hydrogen bond distance and S-H-O angle of 2.515(4) angstrom and 117.4(1)degrees, respectively. These results are key benchmarks to expand the current knowledge of sulfur hydrogen bonds and the relationship between internal motions and conformational preferences in esters.

(PMe3)2FeCl3 is an Fe(III) complex that exists in the intermediate-spin ground state in a distorted trigonal bipyramidal geometry. An electronic state with high-spin configuration lies in close vic...

Infrared spectra of the neutral dimethylamine-methanol cluster, DMA-CH3OH, were measured in the spectral range of 2800-3900 cm-1 using an infrared-vacuum ultraviolet (IR-VUV) scheme. Quantum chemical calculations and ab initio molecular dynamic (AIMD) simulations were carried out to understand the experimental spectral features. Experimental and theoretical results reveal the coexistence of N···HO and O···HN hydrogen-bonded structures. AIMD simulations show that the methyl group in methanol internally rotates around the N···O axis, addressing the dynamic effect of the fluctuation of hydrogen bonds on the vibrational features. The bonding analysis was performed to elucidate the nature of the intermolecular interaction between DMA and CH3OH. The present work provides the fundamental understanding of hydrogen-bonding networks in the amine-alcohol complexes.

International audience; We report on the photoionization of the anilino radical (C 6 H 5 NH) formed by H atom abstraction from aniline by F atoms in a flow-tube. The spectra were recorded from 7.8 to 9.7 eV by using a double-imaging photoelectron-photoion coincidence spectrometer with VUV radiation provided by the DESIRS beamline at the SOLEIL synchrotron. The vibrationally-resolved recorded threshold photoelectron spectrum (TPES) of the anilino radical, exhibiting an adiabatic ionization energy of 8.02 ± 0.02 eV, is accompanied by theoretically simulated spectra and geometry optimizations. The TPES of the aniline radical displays remarkable similarity with the threshold photoelectron spectrum of its structural isomer, 4-picolyl. This raises questions about the limitations to separating isomeric products in e.g. flame and combustion studies, with threshold photoelectron-photoion spectroscopy.

The coupling of a molecular ensemble to the confined electromagnetic modes of a microcavity can strongly modify the photophysics and photochemistry of the molecules upon photoexcitation. We investigate here how collective coupling effects lead to modifications of the mechanisms and rates of photochemical processes, in particular, photodissociation and nonradiative decay in NaI and pyrazine, respectively. We show that, after direct excitation into the lower polaritonic states, the lower-energy light-matter hybrid states, the dynamics of the molecular ensemble coupled to light is very similar to the dynamics of the corresponding isolated molecules. Conversely, excitation into the upper polaritonic states results in more complex dynamics that involve as a first step the population transfer toward the manifold of intermediate dark states. These dynamics differ substantially from those of the isolated molecules and may result in measurable time delays for nonradiative decay or excited-state reaction mechanisms. Similarly, we describe how addition of a buffer of nonreactive molecules coupled to the cavity mode can be used to delay the onset of the photochemical processes of the reactive part of the ensemble, where the buffer medium is more effective in inhibiting the reactive process than only reactive molecules in the cavity.

9 pags., 4 figs., 3 tabs., Complexes between metal cations and molecular hydrogen are systems quite amenable for precise spectroscopic and theoretical studies, and at the same time, they are relevant for applications in hydrogen storage and astrochemistry. In this work, we report new intermolecular potential energy surfaces and rovibrational states calculations for complexes involving molecular hydrogen and alkaline metal cations, M+−H2 (M+ = Na+, K+, Rb+, Cs+). The intermolecular potentials, formulated in an internally consistent way to emphasize differences in the properties of the systems, are represented by simple analytical expressions whose parameters have been optimized from comparison with accurate ab initio calculations. Properties of the low-lying bound statesbinding energies, frequencies, and rotational constantsare compared with previous measurements or computations and an overall good agreement is achieved, supporting the reliability of the present formulation. Variations of these properties as a function of the cation size and isotopic substitution, with a proper sequence of ortho and para rotational levels, are also discussed., The work has been funded by Spanish MINECO Grants FIS2017-84391-C2-2-P and FIS2017-83157-P. F.P. thanks the MIUR and the Universita degli Studi di Perugia for financial support to the project AMIS, through the program Dipartimenti di Eccellenza 2018−2022. Allocation of computing time by CESGA (Spain) is also acknowledged

A number of hydrogen-bonded complexes, formed between an alcohol donor and dimethylselenide, have been detected experimentally, at room temperature in the gas phase using FTIR spectroscopy. The Gibbs energy of complex formation has been determined from the measured integrated absorbance of the hydrogen-bonded OH stretching band and the calculated oscillator strength of the associated transition. The OH stretching frequency and Gibbs energy of the selenium hydrogen-bonded complexes are compared to those found in complexes with the same donor molecule and either dimethylether (O) or dimethylsulfide (S) as the acceptor molecule. For a given donor, we found a similar OH stretching frequency in the complexes for each of the three acceptors O, S, and Se. However, the Gibbs energies were found to be less positive (i.e., stronger bound) for the dimethylether complexes (OH·O), as compared to the dimethylsulfide (OH·S) and dimethylselenide (OH·Se) complexes, with the latter two having comparable Gibbs energies.

Aeration of common liquid solvents (e.g., cyclohexane, benzene, or water) gives rise to a unique absorption band associated with a transition from a ground state oxygen-solvent collision complex to an oxygen-solvent charge-transfer (CT) state. Although the absorption band is solvent-dependent, the energies involved invariably correspond to a one-photon transition in UV. We now show that this CT state can also be produced in a nonlinear two-photon process using visible light. Moreover, the CT state decays to produce an appreciable amount of singlet oxygen, O

Many water-anion complexes of the form X-·(H2O), where X- is a polyatomic anion, display a peak progression in the OH stretch region of the vibrational spectra with spacings of 65-85 cm-1. These progressions result from strong anharmonic coupling between the OH stretch and a low-frequency intermolecular rock vibration. In this study, we calculate these progressions in HCO2-·(H2O), NO3-·(H2O), and CS2-·(H2O) by use of a one-dimensional adiabatic model with rock potentials generated from ab initio energies and frequencies. The importance of using a geometry-dependent reduced mass in calculating the peak spacings is demonstrated. We find that the one-dimensional adiabatic model is more successful in predicting peak spacings in the spectrum of HCO2-·(H2O) than for NO3-·(H2O), for which the rock vibration is highly anharmonic.

The oxidation state, also called oxidation number, of atoms in molecules is a fundamental chemical concept. It is defined as the charge of an atom in a molecule after the ionic approximation of its heteronuclear bonds is applied. Even though for simple molecules the assignment of oxidation states is straightforward, redundancy and ambiguity do exist for others. In this work, we present a density-based framework to determine the oxidation state using the quantities from the information-theoretic approach. As a proof of concept, we examined six elements for a total of 49 molecules. Strong linear correlations were obtained with Shannon entropy, Ghosh-Berkowitz-Parr entropy, information gain, relative Renyi entropy of orders 2 and 3, and Hirshfeld charge. We also discovered that the crystal radius of elements plays the key role in rationalizing the system dependent nature of these strong linear correlations. The validity and effectiveness of our results were demonstrated by the examples of molecules containing elements with two or more oxidation states. Our results should be applicable to more complicated systems in assigning different oxidation states.

Experiments on the lowest lying singlet and triplet states of didehydropyrazine isomers (that focus on energy gaps, geometries, and vibrational frequencies) have been carried out computationally by implementing the improved virtual orbital-based multireference (MR) ab initio methods. Pyrazines possess a reasonable MR nature making their description challenging with the conventional quantum chemical approaches. Although wave functions of the diradicals usually have two dominant configurations, a larger reference space is warranted to consistently and accurately describe pyrazine diradicals indicating the complex nature of the systems. Present calculations predict a singlet ground state for ortho- and para-pyrazines, while a triplet ground state is suggested for the meta isomer. The adiabatic singlet-triplet energy gap for the para form is found to be notably higher (by 28.4 kcal/mol) compared to ortho (2.0 kcal/mol) and meta isomers (-5.1 kcal/mol). Accurate and reliable computations are imperative for forecasting the state-ordering in such diradicals. The structural properties obtained from the present calculations lend strong support toward future experimental explorations on these systems.

New molecular beam scattering experiments are reported for the water–hydrogen system. Integral cross sections of the first rotational excitations of para- and ortho-H2O by inelastic collisions with...

Based on the configurations generated by molecular dynamics (MD) simulations using the on-the-fly density-functional tight-bonding (DFTB) force field, we investigated performance of the discrete solvent reaction field (DRF) model coupled to time-dependent density functional theory (TD-DFT) for solvatochromic effect of one- and two-photon absorption phenomena. Dimethylaminonitrostilbene (DANS) molecule solvated in chloroform, dichloromethane, and dimethyl sulfoxide solvents was selected as a model system for our research purpose. For every selected MD/DFTB configuration, within the context of the DRF, solute molecule is represented by TD-DFT and solvent molecules are described by atomic charges and polarizabilities. The calculated one-photon absorption energies reproduce well the positive solvatochromic behavior of solvated DANS and are in good agreement with available experimental data. For the two-photon absorption cross section, even though our approach overshot the experimental data by about 20% in absolute magnitude, experimentally observed solvatochromic change was captured qualitatively in this work. At last, we examined the contributions of atomic charges and polarizabilities of solvent molecules to the solvatochromic shifts of properties of interest.

Curvilinear coordinate Monte Carlo phase space integration and a series of full-dimensional fitted potential energy surfaces are used to study the effectiveness of reduced-dimensional models for predicting rovibrational anharmonicity at high temperatures. Fully coupled and fully anharmonic, but classical, rovibrational partition functions Q are computed for 14 species with two or three fluxional modes (inversions or torsions) and as many as 30 degrees of freedom. These results are converted to semiclassical anharmonicity correction factors f and are analyzed alongside results obtained previously for 22 species with up to two fluxional modes. As expected, fluxional species show considerable variation in f at high temperatures; f is as small as 0.2 for acetone and is as large as 9 for methylene glycol at 2500 K. This set of full-dimensional results is used to test the accuracy of reduced-dimensional models where fluxional modes are treated as coupled to one another but as separable from the remaining nonflu...

Experimentalists working with diradicals are often facing the question of what kind of species among singlet or triplet diradicals or closed-shell molecules are observed. To treat large diradicals with a high density of electronic states, we propose a simplified version of the spin-flip time-dependent density functional theory (SF-TD-DFT) method for a fast computation of their state energies and absorption spectra with an accuracy similar to the nonsimplified scheme. An ultrafast tight-binding variant called SF-sTD-DFT-xTB is also developed to treat even larger systems. For a benchmark set of nine diradicals, good agreement between simplified and conventional SF excitation energies for standard functionals is found. This shows that the proposed parameterization is robust for a wide range of Fock exchange mixing values. With the asymptotically correct response integrals used in SF-sTD-DFT and a correction factor of

Two-photon absorption (2PA) spectra of liquid methanol and ethanol are reported for the energy range 7-10 eV from the first electronic excitation to close to the liquid-phase ionization potential. The spectra give detailed information on the electronic structures of these alcohols in the bulk liquid. The focus of this Article is to examine the electronic structure change compared with water on substitution of a hydrogen by an alkyl group. Continuous 2PA spectra are recorded in the broadband pump-probe fashion, with a fixed pump pulse in the UV region and a white-light continuum as a probe. Pump pulses of two different energies, 4.6 and 6.2 eV, are used to cover the spectral range up to 10 eV. In addition, theoretical 2PA cross sections for both molecules isolated in the gas phase are computed by the equation-of-motion coupled-cluster method with single and double substitutions (EOM-CCSD). These computational results are used to assign both the experimental 2PA and literature one-photon linear absorption spectra. The most intense spectral features are due to transitions to the Rydberg states, and the 2PA spectra are dominated by the totally symmetric 3pz ← 2pz transition in both alcohols. The experimental 2PA spectra are compared with the simulated 2PA spectra based on ab initio calculations that reveal a general blue shift of the excited transitions upon solvation. The effective 2PA thresholds in methanol and ethanol decrease to 6.9 eV compared with 7.8 eV for water. The analysis of the 2PA polarization ratio leads us to conclude that the excited states of ethanol deviate more markedly from water in the lower energy region compared with methanol. The polarization dependence of the 2PA spectra reveal the symmetries of the excited states within the measured energy range. Natural transition orbital calculations are performed to visualize the nature of the transitions and the orbitals participating during electronic excitation.

Solid-state density functional theory (DFT), molecular dynamics (MD), and terahertz (THz) spectroscopy were used to study the formation of enantiotropically related conformational Form I and Form II polymorphs of the pharmaceutical compound, probucol. DFT calculations were performed on the crystal systems to compare relative lattice energies and the solvent stabilization of the metastable Form II structure. The thermodynamics of solvent inclusion in the Form II·MeOH crystal system were determined from MD simulations, as was the favored conformation of molecular probucol in methanol and ethanol solutions. The findings from both solid-state DFT and MD calculations suggest that the preferred molecular orientations of the probucol molecule in solution and the probable inclusion of methanol in the crystal lattice during the crystallization process lead to the solvent selectivity of the probucol polymorph formation. The additional stabilization energy provided by the crystallization solvent facilitates the nucleation and growth of the Form II polymorph under conditions that favor this metastable crystal form over the thermodynamically stable Form I, despite the higher energy molecular and crystalline configurations of probucol Form II. We demonstrate the influence of solvent on the formation of pharmaceutical polymorphs and provide a molecular-level view of complex interactions leading to polymorphism using a combination of computational methods and THz spectral data.

The three-dimensional reference interaction site model molecular solvation theory with the Kovalenko–Hirata closure relation has been shown to produce excellent solvation characteristics for a large class of (bio)chemical systems in solution. Correct calculation of hydration free energy is central to successful application of any solvation model. In order to find out the best possible force-field parameters to be used for hydration free energy calculation with the aforementioned theory, we have developed an extended database containing a large number of experimental solvation free energies available in the current literature and used a plethora of theoretical models for assessment. The general Amber force field was found to perform satisfactorily, whereas special care should be taken in solute charge assignment with the universal force field.

Performance of the fixed-node diffusion quantum Monte Carlo method (FN-DMC) with a single Slater-Jastrow trial wavefunction using natural orbitals (NOs) from MP2, CCSD, CCSD(T), and CASSCF as well as density functional theory orbitals on atomization energies of the molecules in the Gaussian-2 set is investigated in this work. The effects of spin contamination and pseudopotentials (PPs) are also studied. Our results show that the DMC energy with NOs from MP2 or CCSD(T) is the lowest on average, whereas that from CASSCF is only lower than the DMC energy using Hartree-Fock orbitals. Atomization energies are generally underestimated with DMC, and mean absolute deviations of DMC atomization energies are about 2.8 kcal/mol with NOs from MP2 or CCSD(T) and about 2.7 kcal/mol with NOs from CASSCF or B3LYP orbitals. The better performance of the latter two orbitals is due to a more effective error cancellation. The accuracy of the present FN-DMC is similar to CCSD(T)/aug-cc-pVQZ on atomization energies. In addition, the error of DMC atomization energies tends to be larger for molecules with multiple bonds. DMC energies of the open-shell systems with spin-restricted orbitals are generally higher than those with spin-unrestricted orbitals if spin contamination is not serious and the atomization energies are improved to some extent. Furthermore, our results indicate that the error of the PPs is rather small at the CCSD(T) level, but could be more pronounced in DMC calculations and DMC atomization energies are improved with some newly developed PPs.

The accuracy of separable strategies for constructing full-dimensional potential energy surfaces for collisional energy transfer and collision rate calculations is studied systematically for three alcohols (A = methanol, ethanol, and butanol) and three bath gases (M = Ar, N2, and H2O). The fitting efficiency (defined as the number of ab initio data required to achieve parametrizations of a desired accuracy) is quantified for both pairwise (Buckingham or “exp6”) and nonpairwise (permutationally invariant polynomials, PIPs, of Morse variables) functional forms, for four sampling strategies, and as a function of the complexity and anisotropy of the interaction potential. We find that convergence with respect to the number of sampled ab initio data is largely independent of the choice of functional form but instead varies nearly linearly with the number of adjustable parameters and depends strongly on the sampling strategy. Specifically, the use of biased Sobol quasirandom sampling is ∼7× more efficient than ...