Normal Coordinate Analysis in Cartesian Space II. Planar Symmetric XYn Molecules

Force constants of planar symmetric XY_n type molecules in the Cartesian space have been derived in terms of the vibrational frequencies. The force constants and normal coordinates of some planar XY₃ molecules with available isotoptic data are worked out. Further study of infrared absorption spectra of square planar XY₄ complex ions will be needed for the present normal coordinate analysis.     

Ab initio Description of Nuclear Motion in Nonrigid M k XY n (k≥1) Molecules with Quasirigid XY n Fragments

For nonrigid M _k XY _n (k1) molecules, a model is constructed that describes the motion of k M nuclei relative to the quasirigid XY _n fragment taking into account 3^* k degrees of freedom. The parameters of the potential and kinetic terms of the model Hamiltonian are determined from results of ab initio calculations of the properties of a molecule and its fragments. Solutions of the corresponding Schrödinger equations are obtained by a variational method using bases constructed from products of spherical harmonics and harmonic–oscillator eigenfunctions. The form of model Hamiltonians for nonrigid MXY₄ and M₂XY₄ molecules with quasitetrahedral XY₄ fragments are discussed in detail. Group–theoretical analysis of the symmetry of the Hamiltonians is performed. It is shown that the molecular symmetry groups of nonrigid MXY₄ and M₂XY₄ molecules are the G ₂₄ and G ₄₈ groups, which are isomorphic to the T _d and O _h point groups, respectively.     

Nonempirical Study of Nuclear Dynamics in Nonrigid LiReO4 and K2SO4 Molecules

The previously developed nonempirical model for M _k XY _n (k 1) nonrigid molecules describing the motion of k M nuclei relative to a quasirigid XY _n fragment is used to study the dynamics of nuclei in LiReO₄ and K₂SO₄ molecules. The parameters of the kinetic and potential parts of the model Hamiltonian and dipole moment functions are determined from results of ab initio calculations of the molecules and their fragments by the Hartree—Fock and CISD + Q configurational interaction methods. The dynamic problem is solved by the variational method. Energy levels, transition frequencies, transition dipole moments, and mean geometrical parameters of the molecules are calculated. It is shown that the lowest energy levels of LiReO₄ and K₂SO₄ molecules can be described with high accuracy in a harmonic approximation using a quasirigid, single–minimum model, whereas for the energy levels located near and above the barriers of intramolecular rearrangements, these approximations are completely unsuitable.     

The NICS‐XY‐Scan: Identification of Local and Global Ring Currents in Multi‐Ring Systems

Nucleus‐independent chemical shift (NICS)‐based methods are very popular for the determination of the induced magnetic field under an external magnetic field. These methods are used mostly (but not only) for the determination of the aromaticity and antiaromaticity of molecules and ions, both qualitatively and quantitatively. The ghost atom that serves as the NICS probe senses the induced magnetic field and reports it in the form of an NMR chemical shift. However, the source of the field cannot be determined by NICS. Thus, in a multi‐ring system that may contain more than one induced current circuit (and therefore more than one source of the induced magnetic field) the NICS value may represent the sum of many induced magnetic fields. This may lead to wrong assignments of the aromaticity (and antiaromaticity) of the systems under study. In this paper, we present a NICS‐based method for the determination of local and global ring currents in conjugated multi‐ring systems. The method involves placing the NICS probes along the X axis, and if needed, along the Y axis, at a constant height above the system under study. Following the change in the induced field along these axes allows the identification of global and local induced currents. The best NICS type to use for these scans is NICS_πZZ, but it is shown that at a height of 1.7 Å above the molecular plane, NICS_ZZ provides the same qualitative picture. This method, namely the NICS‐XY‐scan, gives information equivalent to that obtained through current density analysis methods, and in some cases, provides even more details.     

A unified orbital model of delocalised and localised currents in monocycles, from annulenes to azabora-heterocycles

Why are some (4n+2)π systems aromatic, and some not? The ipsocentric approach to the calculation of the current density induced in a molecule by an external magnetic field predicts a four‐electron diatropic (aromatic) ring current for (4n+2)π carbocycles and a two‐electron paratropic (antiaromatic) current for (4n)π carbocycles. With the inclusion of an electronegativity parameter, an ipsocentric frontier‐orbital model also predicts the transition from delocalised currents in carbocycles to nitrogen‐localised currents in alternating azabora‐heterocycles, which rationalises the differences in (magnetic) aromaticity between these isoelectronic π‐conjugated systems. Ab initio valence‐bond calculations confirm the localisation predicted by the naïve model, and coupled‐Hartree–Fock calculations give current‐density maps that exhibit the predicted delocalised‐to‐localised/carbocycle–heterocycle transition.     

Quantum-chemical study of the formation mechanism of cucurbit[<Emphasis Type="Italic">n</Emphasis>]uril nanocavitands

The quantum-chemical study of the formation mechanism of cucurbit[n]urils CB[n] has been performed by the PBE density functional theory method. According to the thermodynamic characteristics of separate stages of CB[n] formation, the oligomer-bound water molecules formed in the course of synthesis are responsible for the predominant formation of CB[6] as compared with other CB[n] nanocavitands.     

Static and dynamic density functional investigation of hydrated beryllium dications

Using density functional theory with Becke's gradient correction for the exchange part and Vanderbilt's ultrasoft pseudopotentials, we investigated [Be(H₂O)_n]²⁺ clusters for n = 1 and 3. A new scheme implemented with the Car-Parrinello method is used, upon which no periodic boundary conditions are imposed so that isolated and possibly charged molecules can also be treated dynamically in a plane-wave basis. Harmonic vibrational frequencies are obtained via a fit of the molecular dynamics trajectory in terms of harmonic oscillators, for which we use a fragmentation scheme to analyze complex spectra. Based on this combination of techniques, we find good agreement with data from SCF calculations and correlated methods. © 1996 John Wiley & Sons, Inc.     

Korrelationen von „out-of-plane“-Kraftkonstanten planarerXYCZ-Moleküle mit Elektronegativitätswerten

Several linear relations of mean electronegativities with force constants for carbonyl and thiocarbonyl out-of-plane deformations of planarXYCZ molecules have been found. The influence of the correspondingG matrix element on these relations has been shown.

     

Perimeter Effects on Ring Currents in Polycyclic Aromatic Hydrocarbons: Circumcoronene and Two Hexabenzocoronenes

Current‐density maps are calculated at an ab initio level for the three symmetrical polycyclic aromatic hydrocarbons, circumcoronene [ 1 (D_6h)], hexabenzo[bc,ef,hi,kl,no,qr]coronene [ 2 a (D_6h) and 2 b (D_3d)], and hexabenzo[a,d,g,j,m,p]coronene [ 3 a (D_6h), 3 b (D₆) and 3 c (D_3d)], all of which can be formally derived by annelation of benzene rings to a coronene core. Whilst 1 is planar, 2 has a non‐planar minimum that is effectively isoenergetic with its planar form, and 3 has a well defined non‐planar structure. The shape of the molecular boundary rather than the planarity of the molecule plays the critical rôle in the character of the predicted currents. Formal deletion of outer hexagons from circumcoronene ( 1 ) in two different ways produces either hexabenzocoronene 2 with a prediction of disjoint local benzenoid diatropic currents linked by a global perimeter, or 3 with a giant diatropic perimeter current enclosing a weak paramagnetic circulation on the central hexagon. The current density map of 1 is effectively a superposition of those of 2 and 3 . Its strong diatropic perimeter current subsumes the six weaker diatropic benzenoid circulations evident in 2 , and bifurcates in the six outer benzenoid rings that form the corners of the giant hexagon; its benzene “hub” sustains a diatropic current, as would be expected from the partial cancellation of the strong diatropic hub current of 2 by the weaker paratropic hub current of 3 . The relationship between the three molecules is rationalised by considering orbital contributions to their current density maps.     

Bonding Character,Electron Delocalization,and Aromaticity of Cyclo[18]Carbon (C18) Precursors,C18-(CO)n (n=6, 4, and 2): Focusing on the Effect of Carbonyl (-CO) Groups**

The bonding character, electron delocalization, and aromaticity of the cyclo[18]carbon (C₁₈) precursors, C₁₈-(CO)_n (n=6, 4, and 2), have been studied by combining quantum chemical calculations and various electronic wavefunction analyses with different physical bases. It was found that C₁₈-(CO)_n (n=6, 4, and 2) molecules exhibit alternating long and short C−C bonds, and have out-of-plane and in-plane dual π systems (π^out and πⁱⁿ) perpendicular to each other, which are consistent with the relevant characteristics of C₁₈. However, the presence of carbonyl (-CO) groups significantly reduced the global electron conjugation of C₁₈-(CO)_n (n=6, 4, and 2) compared to C₁₈. Specifically, the -CO group largely breaks the extensive delocalization of πⁱⁿ system, and the π^out system is also affected by it but to a much lesser extent; as a consequence, C₁₈-(CO)_n (n=6, 4, and 2) with larger n shows weaker overall aromaticity. Mostly because of the decreased but still apparent π^out electron delocalization in the C₁₈-(CO)_n (n=6, 4, and 2), a notable diatropic induced ring current under the action of external magnetic field is observed, demonstrating the clear aromatic characteristic in the molecules. The correlation between C₁₈-(CO)_n (n=6, 4, and 2) and C₁₈ in terms of the gradual elimination of -CO from the precursors showed that the direct elimination of two CO molecules in C₁₈-(CO)_n (n=6, 4, and 2) has a synergistic mechanism, but it is kinetically infeasible under normal conditions due to the high energy barrier.     

Kramers' restricted hartree—fock method for polyatomic molecules using ab initio relativistic effective core potentials with spin—orbit operators

A two-component Kramers' restricted Hartree–Fock method (KRHF) has been developed for the polyatomic molecules with closed shell configurations. The present KRHF program utilizes the relativistic effective core potentials with spin–orbit operators at the Hartree–Fock (HF) level and produces molecular spinors obeying the double group symmetry. The KRHF program enables the variational calculation of spin–orbit interactions at the HF level. KRHF calculations have been performed for the HX, X₂, XY(X, Y = I, Br), and CH₃I molecules. It is demonstrated that the orbital energies from KRHF calculations are useful for the interpretation of spin-orbit splittings in photoelectron spectra. In all molecules studied, bond lengths are only slightly expanded, harmonic vibrational frequencies are reduced, and bond energies are significantly decreased by the spin–orbit interactions.     

Structural properties and the effect of 2,6‐diaminoanthraquinone on G‐tetrad,non‐G‐tetrads,and mixed tetrads—A density functional theory study

Telomerase inhibitor causes the attrition of telomere length and consequently leading to senescence which require a lag period for cancer cells to stop proliferating. Telomeric sequences form quadruplex structures stabilized by tetrads. The structural and electronic properties related with interaction of 2,6‐diaminoanthraquinone and tetrads are the key step to elucidate the anticancer activity. The present study has been focused on the stability of the isolated tetrads and the effect of interaction of 2,6‐diaminoanthraquinone with G‐tetrad, non‐G‐tetrads, and mixed tetrads using density functional theory method in both gas and aqueous phases. The solvent interaction with the molecular systems has increased the stability of the isolated tetrads and complexes. The sharing of electron density between the interacting molecules is shown through electron density difference maps. The atoms in molecules theory and natural bond orbital analysis have been performed to study the nature of hydrogen bonds in the inhibitor interacting complexes. The linear correlation is shown between electron density [ρ(r)], and its Laplacian [(²ρ(r)] at the bond critical points. The strong binding nature of 2,6‐diaminoanthraquinone with studied tetrads reveals that this inhibitor is suitable to stabilize the above tetrads and inhibit the telomerase activity. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Insights into the dynamics of evaporation and proton migration in protonated water clusters from Large‐scale Born–Oppenheimer direct dynamics

Large‐scale on‐the‐fly Born–Oppenheimer molecular dynamics simulations using recent advances in linear scaling electronic structure theory and trajectory integration techniques have been performed for protonated water clusters around the magic number (H₂O)_nH⁺, for n = 20 and 21. Besides demonstrating the feasibility and efficiency of the computational approach, the calculations reveal interesting dynamical details. Elimination of water molecules is found to be fast for both cluster sizes but rather insensitive to the initial geometry. The water molecules released acquire velocities compatible with thermal energies. The proton solvation shell changes between the well‐known Eigen and Zundel motifs and is characterized by specific low‐frequency vibrational modes, which have been quantified. The proton transfer mechanism largely resembles that of bulk water but one interesting variation was observed. © 2012 Wiley Periodicals, Inc.     

Theoretical and experimental study of the structure and stability of multiply Na-substituted glucose and 2,4,6-trihydroxyacetophenone derivatives

The equilibrium geometric parameters and the energetic and spectroscopic characteristics of low lying conformers for series of polyhydroxyl molecules and ions in which sodium atoms are successively substituted for the hydroxyl hydrogen atoms have been calculated by the density functional theory B3LYP method with the 6−31G* and 6−311+G** basis sets. The glucose derivatives [Glu − nH + nNa] and [Glu − nH + (n + 1)Na]⁺ (n = 1−5) and the 2,4,6-trihydroxyacetophenone derivatives [THAP − nH + nNa] and [THAP − nH + (n + 1)Na]⁺ (n = 1−4) have been considered. The affinities of the neutral [Glu − nH + nNa] and [THAP − nH + nNa] molecules for adding Na⁺ cations, as well as the energies of successive substitution of Na atoms for H atoms in the Glu and THAP molecules and the Glu⁺ and THAP⁺ ions in their reaction with sodium acetate molecules, have been estimated. Computations show that the first substitution of Na for H in ions is slightly exothermic and, presumably, can spontaneously occur under common conditions. Further substitutions are endothermic, but the required energy inputs are small. Therefore, successive substitutions for two, three, or more hydroxyl H atoms in the molecules and ions under consideration are possible at relatively low energy inputs. The computation results and conclusions are compared with the MALDI TOF mass spectral data for Na-substituted glucose and 2,4,6-trihydroxyacetophenone derivatives in the [glucose + CH₃COONa + THAP] system where, in addition to common Glu · Na⁺ and THAP · Na⁺ ion-molecular complexes, multiply substituted positive ions of the [Glu − nH + (n + 1)Na]⁺ (n = 1−4) and [THAP − nH + (n + 1)Na]⁺ (n = 1−3) type have been identified.     

Density functional studies of AnF6 (An=U,Np, and Pu) and UF6−nCln (n=1–6) using hybrid functionals: geometries and vibrational frequencies

We have compared the performance of widely used hybrid functionals for calculating the bond lengths and harmonic vibrational frequencies of AnF₆ (An=U, Np, and Pu) and UF_6?nCl_n (n=1–6) molecules using “small‐core” relativistic effective core potentials and extended basis sets. The calculated spectroscopic constants compare favorably with experimental data for the bond lengths (average error ≤ 0.01 Å) and vibrational frequencies (average error ≤ 7 cm^?1) of the AnF₆ molecules. The experimental vibrational frequencies of the stretching modes were available for most of the UF_6?nCl_n (n=1–6) molecules. The calculated vibrational frequencies are in good agreement with the experimental data to within 4.6 cm^?1 and 7.6 cm^?1 for selected Becke1 and Lee, Yang, Parr (B1LYP), and Becke3 and Perdew, Wang (B3PW91) functionals, respectively. We conclude that one can predict reliable geometries and vibrational frequencies for the unknown related systems using hybrid density functional calculations with the RECPs. The geometries and vibrational frequencies of the UF_6?nCl_n (n=1–6) molecules that have not been determined experimentally are also presented and discussed. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 2010–2017, 2001     

“Expanded” local mode approach for XY2 (C2v) molecules

For XY₂ (C_2v symmetry) molecules, the local mode method is modified on the basis of a careful analysis of the transformation coefficients l_Nαλ and the properties of the intramolecular potential function, and an expanded local mode approach is derived. It is shown that relationships between known spectroscopic parameters, found in the spectroscopic literature, can be improved. New relationships between the centrifugal distortion coefficients are obtained. It is shown that knowledge of the centrifugal distortion coefficients of the ground vibrational state allows the possibility of predicting values of both the band centres and the spectroscopic parameters of deformational bands, even under conditions of a total absence of initial information on these deformational bands. The prediction possibility of the derived expanded local mode approach is analyzed.     

Self-consistent Thomas–Fermi–Dirac theory,extended by Gell-Mann and Brueckner correlation,in terms of density n and its two reduced gradients ∇2n/n and ∇n/n

We prove the following results, relevant for the density functional theory: the Thomas–Fermi–Dirac theory, generalized to include the contribution due to the high electron density result of Gell-Mann and Brueckner for the correlation energy, is shown to lead to a differential equation for the self-consistent ground-state density n( r ) in atoms and molecules in the form F(n, { ∇ n/n}², ∇²n/n)=1, where the function F is given explicitly. A straightforward extension yields a similar result for the equation determining the Pauli plus exchange–correlation potential and for the divergence of the many-electron force. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 69: 145–149, 1998     

Tailoring the donor moieties in TPA-based organic dyes for efficient photovoltaic,optical and nonlinear optical response properties

The current study reports tailoring the electronic donor structures of organic dyes to modify their optical and nonlinear optical (NLO) response properties. Five (5) tri-phenyl amine (TPA) based Donor-π-Acceptor (D-π-A) organic dyes with the codes ICAA1 , ICAA2 , ICAA3 , ICAA4 , and ICAA5 were designed and investigated for their optical and NLO properties using quantum chemical methods. Optical and NLO properties of these dyes were studied by CAM-B3LYP method and 6-311G* basis set. The focus has been on the impact of adding secondary donors and shifting their substitutions at ortho (o), meta (m) and para (p) positions. Among all designed compounds, ICAA4 showed the highest amplitude of average third-order NLO polarizability <γ>, which is calculated to be 1316 × 10⁻³⁶ esu. Time-dependent Density Functional Theory (TD-DFT) method was used to determine how a change in the position of the donor affected the excitation energy (E_g) and NLO response properties. The findings showed that changing the position of the secondary donor results in a red shift among absorption spectra as well as the increase in their NLO responses. Complete process of intramolecular charge transfer (ICT) has been investigated in terms of different optical parameters such as frontier molecular orbitals (FMOs), molecular electrostatic potentials (MEPs), transition density matrix (TDMs), density of states (DOS), electron density difference (EDD), and natural bond orbital (NBO) analysis. Our calculations for study of ICT process indicate that p-position of methoxy group performs better among all other positions and even it has better NLO response properties than the compound with three collective methoxy groups. The calculated V_oc values of all designed molecules range from 1.09 to 1.30, all of them are positive while their ΔG_inject is found to be in the range of −0.87 to −1.79 eV indicating their decent potential for photovoltaic applications. The studied optical, NLO and photovoltaic parameters illustrated that ICAA1 to ICAA5 are appropriate molecules not only for NLO applications but also for efficient photovoltaic purposes.     

Computational investigation of the structures,properties, and host-guest chemistry of prism[n]arenes

The discovery of new and functional macrocyclic host compounds is an important part of supramolecular chemistry. Since the experimental synthesis, prism[n]arenes (Pr[n]As), a class of naphthol-based macrocyclic arenes, have attracted much attention. In this work, from the perspective of theoretical calculation and research, Pr[n]As (n = 4 ~ 7) were studied by density functional theory (DFT) calculations and molecular dynamics (MD) simulations. The prismatic configuration isomers, electronic structures, absorption spectra, and host-guest chemistry were discussed thoroughly. DFT calculation results showed that 1,5-, 3,7-, and “hybrid” 15,37-Pr[n]As were the most representative configurations with the rigid prismatic molecular skeleton. Based on time-dependent density functional theory (TD-DFT), the absorption spectra of Pr[n]As were all in the range of ultraviolet light which were mainly attributed to π-π* transitions. The molecular cavities of Pr[n]As were electron-rich and capable of accommodating a variety of cations or electron-conjugated molecules. MD simulation results showed that a Pr[n]A molecule was able to capture the guest molecule into its molecular cavity and maintain in the state of equilibrium in solvents.     

一系列水合和非水合铀酰卤族(F,Cl,Br)配合物在水溶液中的结构和紫外吸收光谱性质的密度泛函理论研究

本文考虑相对论效应并应用密度泛函理论(DFT)研究水溶液中UO2Xn(H2O)5-n(X=F,Cl,Br;n=1～4)和UO2Xn(X=F,Cl,Br;n=1～6)一系列水合和非水合铀酰化合物的结构和紫外吸收光谱性质。将这一系列物质命名为Xnm(X为F,Cl,和Br;n为卤素配体个数,m为水分子配体的个数)。在水溶液中,溶剂化效应采用类导体屏蔽模型(COSMO)并采用SAS溶剂接触曲面构造空穴模拟水溶剂对配合物的作用。配合物的紫外光谱性质采用考虑旋-轨耦合相对论效应的含时密度泛函(SO-TD-DFT)进行计算。U=O键随着F配体数目的增加而明显伸长,然而随Cl和Br配体数目的增加变化较小。随X配体数目的增加和水分子参与配位,铀与X的结合能逐渐减弱。配合物的紫外光谱计算表明铀酰氟的各种配合物并不出现特征吸收峰,而铀酰氯和铀酰溴的各种配合物均有特征吸收光谱。通过分子轨道分析可以很好解释光谱所体现的特征。