Fluxional Research Articles

Time-dependent density matrix renormalization group coupled with n-mode representation potentials for the excited state radiationless decay rate: Formalism and application to azulene

We propose a method for calculating the non-radiative decay rates for polyatomic molecules including anharmonic effects of the potential energy surface (PES) in the Franck-Condon region. The method combines the n-mode representation method to construct the ab initio PES and the nearly exact time-dependent density matrix renormalization group method (TD-DMRG) to simulate quantum dynamics. In addition, in the framework of TD-DMRG, we further develop an algorithm to calculate the final-state-resolved rate coefficient which is very useful to analyze the contribution from each vibrational mode to the transition process. We use this method to study the internal conversion (IC) process of azulene after taking into account the anharmonicity of the ground state PES. The results show that even for this semi-rigid molecule, the intramode anharmonicity enhances the IC rate significantly, and after considering the two-mode coupling effect, the rate increases even further. The reason is that the anharmonicity enables the C–H vibrations to receive electronic energy while C–H vibrations do not contribute on the harmonic PES as the Huang-Rhys factor is close to 0.

Reactivity of chloroacetamides toward sulfide + black carbon: Insights from structural analogues and dynamic NMR spectroscopy

Chloroacetamides are commonly used in herbicide formulations, and their occurrence has been reported in soils and groundwater. However, how their chemical structures affect transformation kinetics and pathways in the presence of environmental reagents such as hydrogen sulfide species and black carbon has not been investigated. In this work, we assessed the impact of increasing Cl substituents on reaction kinetics and pathways of six chloroacetamides. The contribution of individual pathways (reductive dechlorination vs. nucleophilic substitution) to the overall decay of selected chloroacetamides was differentiated using various experimental setups; both the transformation rates and product distributions were characterized. Our results suggest that the number of Cl substituents affected reaction pathways and kinetics: trichloroacetamides predominantly underwent reductive dechlorination whereas mono- and dichloroacetamides transformed via nucleophilic substitution. Furthermore, we synthesized eight dichloroacetamide analogs (Cl2CHC(=O)NRR’) with differing R groups and characterized their transformation kinetics. Dynamic NMR spectroscopy was employed to quantify the rotational energy barriers of dichloroacetamides. Our results suggest that adsorption of dichloroacetamides on black carbon constrained R groups from approaching the dichloromethyl carbon and subsequently favored nucleophilic attack. This study provides new insights to better predict the fate of chloroacetamides in subsurface environments by linking their structural characteristics to transformation kinetics and pathways.

Double Proton Tautomerism via Intra- or Intermolecular Pathways? The Case of Tetramethyl Reductic Acid Studied by Dynamic NMR: Hydrogen Bond Association, Solvent and Kinetic H/D Isotope Effects.

Using dynamic liquid-state NMR spectroscopy a degenerate double proton tautomerism was detected in tetramethyl reductic acid (TMRA) dissolved in toluene-d8 and in CD2Cl2. Similar to vitamin C, TMRA belongs to the class of reductones of biologically important compounds. The tautomerism involves an intramolecular HH transfer that interconverts the peripheric and the central positions of the two OH groups. It is slow in the NMR time scale around 200 K and fast at room temperature. Pseudo-first-order rate constants of the HH transfer and of the HD transfer after suitable deuteration were obtained by line shape analyses. Interestingly, the chemical shifts were found to be temperature dependent carrying information about an equilibrium between a hydrogen bonded dimer and a monomer forming two weak intramolecular hydrogen bonds. The structures of the monomer and the dimer are discussed. The latter may consist of several rapidly interconverting hydrogen-bonded associates. A way was found to obtain the enthalpies and entropies of dissociation, which allowed us to convert the pseudo-first-order rate constants of the reaction mixture into first-order rate constants of the tautomerization of the monomer. Surprisingly, these intrinsic rate constants were the same for toluene-d8 and CD2Cl2, but in the latter solvent more monomer is formed. This finding is attributed to the dipole moment of the TMRA monomer, compensated in the dimer, and to the larger dielectric constant of CD2Cl2. Within the margin of error, the kinetic HH/HD isotope effects were found to be of the order of 3 but independent of temperature. That finding indicates a stepwise HH transfer involving a tunnel mechanism along a double barrier pathway. The Arrhenius curves were described in terms of the Bell–Limbach tunneling model.

Revisiting conformations of methyl lactate in water and methanol.

The recently developed efficient protocols to implicit [Grimme et al., J. Phys. Chem. A 125, 4039-4054 (2021)] and explicit quantum mechanical modeling of non-rigid molecules in solution [Katsyuba et al., J. Phys. Chem. B 124, 6664-6670 (2020)] are applied to methyl lactate (ML). Building upon this work, a new combination scheme is proposed to incorporate solvation effects for the computation of infrared (IR) absorption spectra. Herein, Boltzmann populations calculated for implicitly solvated single conformers are used to weight the IR spectra of explicitly solvated clusters with a size of typically ten solvent molecules, i.e., accounting for the first solvation shell. It is found that in water and methanol, the most abundant conformers of ML are structurally modified relative to the gas phase, where the major form is ML1, in which the syn conformation of the -OH moiety is stabilized by a OH⋯O=C intramolecular hydrogen bond (HB). In solution, this syn conformation transforms to the gauche form because the intramolecular HB is disrupted by explicit water molecules that form intermolecular HBs with the hydroxyl and carbonyl groups. Similar changes induced by the gas-solution transition are observed for the minor conformers, ML2 and/or ML3, characterized by OH⋯OCH3 intramolecular HB in the gas phase. The relative abundance of ML1 is shown to decrease from ∼96% in gas to ∼51% in water and ∼92% in methanol. The solvent strongly influences frequencies, IR intensities, and normal modes, resulting in qualitatively different spectra compared to the gas phase. Some liquid-state conformational markers in the fingerprint region of IR spectra are revealed.

Efficient Quantum-Chemical Calculations of Acid Dissociation Constants from Free-Energy Relationships.

The calculation of acid dissociation constants (pKa) is an important task in computational chemistry and chemoinformatics. Theoretically and with minimal empiricism, this is possible from computed acid dissociation free energies via so-called linear free-energy relationships. In this study some modifications are introduced to the latter, providing a straightforward, broadly applicable protocol with an adjustable degree of sophistication for quantum chemistry-based calculations of pKa in water. It targets a wide pKa range (∼70 units) and medium-sized, flexible molecules. Herein, a focus is set on the recently published r2SCAN-3c and related efficient composite density functionals and the semiempirical GFN2-xTB method, including a newly introduced energy correction for heterolytic dissociation, both in combination with implicit solvation models. The performance is evaluated in comparison with experimental data, showing mean errors often smaller than a targeted 1 pKa unit accuracy. Larger deviations are observed only upon inclusion of challenging highly negative (<-5) or positive (>15) pKa values. Among all those tested, it is found that B97-3c is the best performing functional, although rather independently of the density functional theory (DFT) method used; low root-mean-square errors of 0.8-1.0 pKa units for typical drugs are obtained. For optimal performance, it is recommended to employ DFT functional specific free-energy relationship parameters. Additionally, a significant conformational dependence of the pKa values is revealed and quantified for some nonrigid drug molecules.

Revealing the Surface Structure of CdSe Nanocrystals by Dynamic Nuclear Polarization-Enhanced 77Se and 113Cd Solid-State NMR Spectroscopy.

Dynamic nuclear polarization (DNP) solid-state NMR (SSNMR) spectroscopy was used to obtain detailed surface structures of zinc blende CdSe nanocrystals (NCs) with plate or spheroidal morphologies which are capped by carboxylic acid ligands. 1D 113Cd and 77Se cross-polarization magic angle spinning (CPMAS) NMR spectra revealed distinct signals from Cd and Se atoms on the surface of the NCs, and those residing in bulk-like environments, below the surface. 113Cd cross-polarization magic-angle-turning (CP-MAT) experiments identified CdSe3O, CdSe2O2, and CdSeO3 Cd coordination environments on the surface of the NCs, where the oxygen atoms are presumably from coordinated carboxylate ligands. The sensitivity gain from DNP enabled natural isotopic abundance 2D homonuclear 113Cd-113Cd and 77Se-77Se and heteronuclear 113Cd-77Se scalar correlation solid-state NMR experiments which revealed the connectivity of the Cd and Se atoms. Importantly, 77Se{113Cd} scalar heteronuclear multiple quantum coherence (J-HMQC) experiments were used to selectively measure one-bond 77Se-113Cd scalar coupling constants (1J(77Se, 113Cd)). With knowledge of 1J(77Se, 113Cd), heteronuclear 77Se{113Cd} spin echo (J-resolved) NMR experiments were used to determine the number of Cd atoms bonded to Se atoms and vice versa. The J-resolved experiments directly confirmed that major Cd and Se surface species have CdSe2O2 and SeCd4 stoichiometries, respectively. Considering the crystal structure of zinc blende CdSe and the similarity of the solid-state NMR data for the platelets and spheroids, we conclude that the surface of the spheroidal CdSe NCs is primarily composed of {100} facets. The methods outlined here will generally be applicable to obtain detailed surface structures of various main group semiconductor nanoparticles.

Geometry of Pentaphenylantimony in Solution: Support for a Trigonal Bipyramidal Assignment from X-ray Absorption Spectroscopy and Vibrational Spectroscopic Data.

Pentaphenylantimony (SbPh5) has been previously crystallized in either a square pyramidal or trigonal bipyramidal geometry. Investigation of the solution-state structure of SbPh5 has been hampered by the extreme fluxionality of this compound, but previous vibrational spectroscopic studies concluded that it maintains a square pyramidal geometry in solution. This non-VSEPR-compliant geometry, which is also assumed by BiPh5 in the solid state, stands in contrast to the trigonal bipyramidal geometries of PPh5 and AsPh5. A range of phenomena have been invoked to explain this discrepancy, most notably, the increased importance of relativistic effects as group 15 is descended. We present crystallographic, spectroscopic, and computational data revealing that SbPh5 in fact assumes the VSEPR-compliant trigonal bipyramidal geometry in solution. In particular, Sb X-ray absorption spectroscopy (XAS) was used to obtain geometry-sensitive spectra that do not suffer from the slow spectroscopic time scale that has prevented NMR studies from elucidating the structure of this fluxional molecule. Sb K-edge and LIII-edge XAS spectra of crystalline solids featuring SbPh5 in either a square pyramidal (nonsolvate) or trigonal bipyramidal (cyclohexane hemisolvate or THF hemisolvate) form were compared to spectra of SbPh5 in solution. The solution-state spectra agree with those from solids containing trigonal bipyramidal SbPh5. The most diagnostic spectroscopic feature was the distribution of intensity in the Sb LIII pre-edge features. These distributions were rationalized using time-dependent density functional theory calculations that take into account spin-orbit coupling. Our use of Sb XAS not only resolves a long-standing physical inorganic question but also demonstrates more widely the utility of XAS in establishing the structures of fluxional main-group compounds. This conclusion was further supported by solid- and solution-state Raman data. Finally, we note that the present high-resolution diffraction data allow τ for nonsolvated SbPh5 to be revised to 0.216.

Efficient Quantum Chemical Calculation of Structure Ensembles and Free Energies for Nonrigid Molecules.

The application of quantum chemical, automatic multilevel modeling workflows for the determination of thermodynamic (e.g., conformation equilibria, partition coefficients, pKa values) and spectroscopic properties of relatively large, nonrigid molecules in solution is described. Key points are the computation of rather complete structure (conformer) ensembles with extremely fast but still reasonable GFN2-xTB or GFN-FF semiempirical methods in the CREST searching approach and subsequent refinement at a recently developed, accurate r2SCAN-3c DFT composite level. Solvation effects are included in all steps by accurate continuum solvation models (ALPB, (D)COSMO-RS). Consistent inclusion of thermostatistical contributions in the framework of the modified rigid-rotor-harmonic-oscillator approximation (mRRHO) based on xTB/FF computed PES is also recommended.

Distance degree vector and scalar sequences of corona and lexicographic products of graphs with applications to dynamic NMR and dynamics of nonrigid molecules and proteins

We have developed novel topological techniques to generate distance degree vector sequences (DDS) and distance degree scalar sequences (SS) of corona and lexicographic products of graphs both of which are of chemical and biochemical interest. Query ID="Q1" Text="Please confirm if the author names are presented accurately and in the correct sequence (given name, middle name/initial, family name). Author 1 Given name: [Medha Itagi] Last name [Huilgol]. Also, kindly confirm the details in the metadata are correct." We have obtained exact analytical expressions for such products for a variety of graphs and trees of chemical and biochemical interest. We have also outlined the applications of these graph products and their DDS to dynamic NMR spectroscopy and isomerization/conformational graphs of nonrigid molecules and proteins whose NMR and conformational dynamics graphs are expressible as lexicographic products of smaller graphs.

Floppy molecules—their internal dynamics, spectroscopy and applications

Abstract Floppy molecules can be defined as molecules performing large amplitude vibrations (LAVs). There are different types of LAVs among which the most common are inversion and internal rotation. Molecules with LAVs have been of great interest for a very long time since their dynamic, geometry and molecular spectra were very often considered as a challenge. In the review, we present an outline of the history and development of various theoretical approaches concerning molecules with LAVs. Different types of LAVs are described with the emphasis on inversion tunneling (wagging) and internal rotation (torsion). Furthermore, strategies for building explicit and effective Hamiltonians are given and explained in detail using a hydrazine molecule, which is an exemplary molecule performing three LAVs—two inversions and one internal rotation. Since floppy molecules play a significant role in numerous areas as chemistry, pharmacy, astrophysics, biology, agriculture etc., we also provide an overview of their applications.

Propane Dehydrogenation on Ga2O3-Based Catalysts: Contrasting Performance with Coordination Environment and Acidity of Surface Sites

α-Ga2O3, β-Ga2O3, and γ-Ga2O3 as well as the silica-supported catalysts γ-Ga2O3/SiO2, β-Ga2O3/SiO2, and Ga(NO3)3-derived Ga/SiO2 were prepared, characterized, and evaluated for propane dehydrogenation (PDH) at 550 °C. The coordination environment and acidity of surface sites in stand-alone and SiO2-supported Ga2O3 catalysts were studied using FTIR, 15N dynamic nuclear polarization surface-enhanced NMR spectroscopy (15N DNP SENS), and DFT modeling of the adsorbed pyridine probe molecule. The spectroscopic data suggest that the Lewis acidic surface Ga sites in γ-Ga2O3 and β-Ga2O3 (the latter obtained from colloidal nanocrystals of γ-Ga2O3 via thermal treatment at 750 °C) are similar, except that β-Ga2O3 contains a larger relative fraction of weak Ga3+ Lewis acid sites. In contrast, α-Ga2O3 features mostly strong Lewis acid sites. This difference in surface sites parallels their difference in catalytic activities: i.e., weak Lewis acid surface sites are more abundant in β-Ga2O3 relative to α-Ga2O3 and γ-Ga2O3 and the increased relative abundance of weak Lewis acidity correlates with a higher initial catalytic activity in PDH, 0.41 > 0.28 > 0.14 mmol C3H6 m–2 (Ga2O3) h–1 at 550 °C, for respectively β-, α-, and γ-Ga2O3 with initial propene selectivities of 86, 83, and 88%. Dispersion of γ-Ga2O3 or β-Ga2O3 on a silica support introduces strong as well as abundant weak Brønsted acidity to the catalysts, lowering the PDH selectivity. The γ-Ga2O3/SiO2 catalyst was slightly more active than β-Ga2O3/SiO2 in PDH (Ga normalized activity) with initial propene formation rates of 11 and 9 mol C3H6 mol Ga–1 h–1 (sel = 76 and 73%, respectively). However, these catalysts deactivated by ca. 55% within 100 min time on stream (TOS) due to coking. In contrast, Ga/SiO2, with mostly tetracoordinated surface Ga sites and abundant, strong Brønsted acid sites, gave a lower activity and selectivity in PDH (3.5 mol C3H6 mol Ga–1 h–1 and 49%, respectively) but showed no deactivation with TOS. DFT calculations using a fully dehydroxylated oxygen-deficient model β-Ga2O3 surface show that tetra- and pentacoordinated Ga Lewis acid sites bind pyridine more strongly than tricoordinated Ga sites and a higher relative fraction of strong Lewis acid sites correlates with increased coking. Overall, our results indicate that weakly Lewis acidic, tricoordinated Ga3+ sites are likely driving the superior PDH activity of β-Ga2O3.

Manipulating the Dynamics of Dark Excited States in Organic Materials for Phototheranostics.

Manipulating the dynamics of dark excited states (DESs), such as higher excited singlet or excited triplet states with no or small radiative decay, are of both fundamental and practical interests, an important application being photoactivated diagnosis and therapy (phototheranostics), which include photoacoustic (PA) imaging, photodynamic therapy (PDT), and photothermal therapy (PTT). However, the current understanding of DESs in organic structures is rather limited, thus making any rational manipulation of DES in organic materials very challenging.A DES decays primarily by radiationless transition through two pathways: (i) singlet-to-triplet intersystem crossing (ISC) and (ii) internal conversion (IC) relaxation. The deactivation of a DES via ISC can generate cytotoxic reactive oxygen species (ROS) for PDT, while IC could convert photons into heat for PA imaging and PTT. In this Account, we highlight our research on developing a fundamental understanding of structure-property relationships for manipulation of DESs in organic materials in relation to phototheranostic applications. We describe the application of femtosecond transient absorption (fs-TA) spectroscopy for obtaining valuable insights into the DES dynamics. Afterward, we present our work on DESs in nonrigid molecules that revealed greatly enhanced ISC through geometry twisting, which leads to an innovative pathway to develop organic materials exhibiting external stimuli-responsive reversible switching of ISC. We introduce the concept of smart PDT where highly efficient ISC imparted by geometry twisting in the acidic environment specific to tumors leads to very efficient and highly localized PDT, thus leaving surrounding healthy tissues at a different pH unaffected. This insightful understanding of ISC can lead to the development of more advanced photosensitizers for PDT. Two other emergent concepts from our work presented here are (1) significantly enhanced IC producing strong local heating by combining two-photon absorption with excited state absorption for cumulative multiphoton absorption, thus greatly increasing the strength of the PA signal for nonlinear PA imaging, and (2) shown by an example of an organic molecule, BODIPY, nanoscale charge-transfer state mediated strong IC in aggregate nanoparticles resulting in exceptionally high photothermal conversion efficiency of 61% for both PA and PTT. Some in vivo results of the phototheranostic studies using BODIPY are presented, providing an elegant example of nanoscale manipulation of the excited state dynamics.This Account concludes with a summary and discussion of future perspectives. We hope this Account will deepen the understanding of molecular and nanoscale control of excited state dynamics in organic materials, hopefully enticing a broad range of scientists within different disciplinary areas.

Accuracy Meets Interpretability for Computational Spectroscopy by Means of Hybrid and Double-Hybrid Functionals.

Accuracy and interpretability are often seen as the devil and holy grail in computational spectroscopy and their reconciliation remains a primary research goal. In the last few decades, density functional theory has revolutionized the situation, paving the way to reliable yet effective models for medium size molecules, which could also be profitably used by non-specialists. In this contribution we will compare the results of some widely used hybrid and double hybrid functionals with the aim of defining the most suitable recipe for all the spectroscopic parameters of interest in rotational and vibrational spectroscopy, going beyond the rigid rotor/harmonic oscillator model. We will show that last-generation hybrid and double hybrid functionals in conjunction with partially augmented double- and triple-zeta basis sets can offer, in the framework of second order vibrational perturbation theory, a general, robust, and user-friendly tool with unprecedented accuracy for medium-size semi-rigid molecules.

Circumambulatory Rearrangements of 5-Halo-1,2,3,4,5-pentaphenylcyclopentadienes

DFT quantum chemical calculations at the CAM-B3LYP/Def2TZVP level of theory showed that intramolecular migrations of halogens in 5-halo-1,2,3,4,5-pentaphenylcyclopentadienes (C5Ph5Hlg, Hlg = F, Cl, Br, I) involve chiral conformation of their molecules with a propeller arrangement of phenyl groups via 1,5-sigmatropic shifts around the five-membered ring through asymmetric transition states with energy barriers ΔE≠ZPE of 42.5 (F), 26.2 (Cl), 20.2 (Br), and 15.2 kcal/mol (I). The results were consistent with the data of dynamic NMR spectroscopy. The P and M enantiomers are readily interconvertible (ΔE≠ZPE = 1.7–3.7 kcal/mol) by way of synchronous flips of the phenyl groups. The calculated barriers to alternative 1,3-halogen shifts in C5Ph5Hlg are considerably higher than those for 1,5-shifts: ΔE≠ZPE = 60.7 (F), 38.6 (Cl), 32.0 (Br), and 27.9 kcal/mol (I).

Unsupervised search of low-lying conformers with spectroscopic accuracy: A two-step algorithm rooted into the island model evolutionary algorithm.

The fruitful interplay of high-resolution spectroscopy and quantum chemistry has a long history, especially in the field of small, semi-rigid molecules. However, in recent years, the targets of spectroscopic studies are shifting toward flexible molecules, characterized by a large number of closely spaced energy minima, all contributing to the overall spectrum. Here, artificial intelligence comes into play since it is at the basis of powerful unsupervised techniques for the exploration of soft degrees of freedom. Integration of such algorithms with a two-stage QM/QM' (Quantum Mechanical) exploration/refinement strategy driven by a user-friendly graphical interface is the topic of the present paper. We will address in particular: (i) the performances of different semi-empirical methods for the exploration step and (ii) the comparison between stochastic and meta-heuristic algorithms in achieving a cheap yet complete exploration of the conformational space for medium sized chromophores. As test cases, we choose three amino acids of increasing complexity, whose full conformer enumeration has been reached only very recently. Next, we show that systems in condensed phases can be treated at the same level and with the same efficiency when employing a polarizable continuum description of the solvent. Finally, the challenging issue represented by the vibrational circular dichroism spectra of some rhodium complexes with flexible ligands has been addressed, showing that our fully unsupervised approach leads to remarkable agreement with the experiment.

Deciphering High-Order Structural Correlations within Fluxional Molecules from Classical and Quantum Configurational Entropy.

We employ the kth nearest-neighbor estimator of configurational entropy in order to decode within a parameter-free numerical approach the complex high-order structural correlations in fluxional molecules going much beyond the usual linear, bivariate correlations. This generic entropy-based scheme for determining many-body correlations is applied to the complex configurational ensemble of protonated acetylene, a prototype for fluxional molecules featuring large-amplitude motion. After revealing the importance of high-order correlations beyond the simple two-coordinate picture for this molecule, we analyze in detail the evolution of the relevant correlations with temperature as well as the impact of nuclear quantum effects down to the ultralow temperature regime of 1 K. We find that quantum delocalization and zero-point vibrations significantly reduce all correlations in protonated acetylene in the deep quantum regime. Even at low temperatures up to about 100 K, most correlations are essentially absent in the quantum case and only gain importance at higher temperatures. In the high temperature regime, beyond roughly 800 K, the increasing thermal fluctuations are found to exert a destructive effect on the presence of correlations. At intermediate temperatures of approximately 100-800 K, a quantum-to-classical cross-over regime is found where classical mechanics starts to correctly describe trends in the correlations whereas it even qualitatively fails below 100 K. Finally, a classical description of the nuclei provides correlations that are in quantitative agreement with the quantum ones only at temperatures exceeding 1000 K. This data-intensive analysis has been made possible due to recent developments of machine learning techniques based on high-dimensional neural network potential energy surfaces in full dimensionality that allow us to exhaustively sample both the classical and quantum ensemble of protonated acetylene at essentially converged coupled cluster accuracy from 1 to more than 1000 K. The presented non-parametric analysis of correlations beyond usual linear two-coordinate terms is transferable to other system classes. The technique is also expected to complement and guide the analysis of experimental measurements, in particular multidimensional vibrational spectroscopy, by revealing the complex coupling between various degrees of freedom.

Facile Synthesis of Aryl-Substituted Cycloarenes via Bismuth(III) Triflate-Catalyzed Cyclization of Vinyl Ethers

Cycloarenes are an essential class of polycyclic aromatic hydrocarbons with unique electronic structure, but their synthesis is very challenging. Herein, we report a facile synthetic strategy prima...

Modified effect on properties of mesophase pitch prepared from various two-stage thermotreatments of FCC decant oil

The comparative study of different two-stage heat treatments was conducted to illustrate the modified effects on improving yield and spinnability of mesophase pitch. In such two-stage thermotreatments, a pressured condensation is employed to involve low-molecular-weight species into the formation of mesogens, thereby increasing pitch yields. Subsequently, vacuum polymerization could contribute to rapid development of mesophase and then reduce the residence time. The results indicate that the short duration could decrease occurrence of the intramolecular aromatization and intermolecular polymerization reactions that are readily induced by elimination of α-methyl and dehyroaromatization of naphthenic rings during the later stage of heat treatment. Consequently, the formed semi-rigid molecules in QS sub-fraction (i.e., fusible component) are seldom transformed into rigid ones mainly constituting QI sub-fraction (i.e., infusible component). Ultimately, sufficient QS sub-fraction, with a higher average Mw and a large amount of short alkyl and naphthenic substituents, is more likely to completely dissolve QI sub-fraction. Thus, the resultant mesophase exhibits an excellent fusible deformability. Besides, the structural feature of semi-rigid oligomer could somewhat loosen molecular association of layered stacking and contribute to the slippage between stacked mesogens due to its steric effects, thereby decreasing the softening points of mesophase pitch. Above modified effects are easily implementable during a suitable process of two-stage thermotreatment, such as the processes of preparing MP-4-0.01-1.0, MP-4-0.005-1.0 and MP-3-0.01-2.0 in this study. It is worthwhile to mention that the vacuum treatment is required to be conducted prior to large-scale coalescence of mesophase spheres in order to obtain spinnable mesophase pitch.

Derivation of ρ-dependent coordinate transformations for nonrigid molecules in the Hougen-Bunker-Johns formalism.

In this paper, we report a series of transformations for the construction of a Hamiltonian model for nonrigid polyatomic molecules in the framework of the Hougen-Bunker-Johns formalism (HBJ). This model is expressed in normal mode coordinates for small vibrations and in a specific coordinate ρ to describe the large amplitude motion. For the first time, a general procedure linking the "true" curvilinear coordinates to ρ is proposed, allowing the expression of the potential energy part in the same coordinate representation as the kinetic energy operator, whatever the number of atoms. A Lie group-based method is also proposed for the derivation of the reference configuration in the internal axis system. This work opens new perspectives for future high-resolution spectroscopy studies of nonrigid, medium-sized molecules using HBJ-type Hamiltonians. Illustrative examples and computation of vibrational energy levels on semirigid and nonrigid molecules are given to validate this method.

Quasi-symmetry effects in the threshold photoelectron spectrum of methyl isocyanate.

The vacuum-ultraviolet threshold photoelectron spectrum of methyl isocyanate CH3NCO has been recorded from 10.4 eV to 12 eV using synchrotron radiation and a coincidence technique allowing for a mass-discrimination of the photoelectron signal. A significant improvement is achieved over previous investigations as this experimental setup leads to a much more resolved spectrum. Ten sharp peaks and a broad feature spanning 1.2 eV were recorded. This spectrum consists of X̃+ 2A″←X̃ 1A' and Ã+ 2A'←X̃ 1A' ionizing transitions. For the former, the adiabatic ionization energy was determined experimentally to be 10.596(6) eV; for the latter, its value was estimated to be 10.759(50) eV. Seven sharp peaks could be assigned to vibrational modes of the cation X̃+ 2A″ and neutral X̃ 1A' ground electronic states involving only the NCO group atoms. Theoretical modeling of the threshold photoelectron spectrum has proven difficult as methyl isocyanate is a non-rigid molecule displaying large amplitude internal rotation of the methyl group and ∠CNC bending mode, leading to the quasi-symmetry. With the help of ab initio calculations, a theoretical model in which these two large amplitude motions are included in addition to the five small amplitude vibrational modes involving NCO group atoms is proposed. Comparison with the experimental spectrum shows that the broad feature and the strongest peak line positions are well accounted for; their intensities are also fairly well reproduced after adjusting a few parameters.