分子旋光度的密度泛函计算及Soman绝对构型确定

手性化合物绝对构型的确定一直是有机化学特别是手性药物合成研究过程中经常遇到的问题之一. 利用密度泛函DFT/B3LYP量子化学计算方法, 通过计算两种已知旋光度分子的旋光度, 确认计算条件和验证计算结果的可靠性, 并将此方法应用于Soman旋光度的计算. 结合文献报道的Soman分子实验旋光度测定结果, 首次采用理论计算方法确定了梭曼(Soman)分子四个异构体的绝对构型, 分别为P(R),C(S)-Soman→P(＋)C(＋), P(R),C(R)-Soman→P(＋)C(－), P(S),C(S)-Soman→P(－)C(＋), P(S),C(R)-Soman→P(－)C(－).     

A Density Functional Theory Study of Optical Rotation in Some Aziridine and Oxirane Derivatives

We present time-dependent density functional theory (TDDFT) calculations of the electronic optical rotation (ORP) for seven oxirane and two aziridine derivatives in the gas phase and in solution and compare the results with the available experimental values. For seven of the studied molecules it is the first time that their optical rotation was studied theoretically and we have therefore investigated the influence of several settings in the TDDFT calculations on the results. This includes the choice of the one-electron basis set, the exchange-correlation functional or the particular polarizable continuum model (PCM). We can confirm that polarized quadruple zeta basis sets augmented with diffuse functions are necessary for converged results and find that the aug-pc-3 basis set is a viable alternative to the frequently employed aug-cc-pVQZ basis set. Based on our study, we cannot recommend the generalized gradient functional KT3 for calculations of the ORP in these compounds, whereas the hybrid functional PBE0 gives results quite similar to the long-range correct CAM−B3LYP functional. Finally, we observe large differences in the solvent effects predicted by the integral equation formalism of PCM and the SMD variant of PCM. For the majority of solute/solvent combinations in this study, we find that the SMD model in combination with the PBE0 functional and the aug-pc-3 basis set gives the best agreement with the experimental values.     

LC with Electrochemical Detection for Analysis of Caffeic Acid and Caffeic Acid Phenyl Ester in Propolis

A simple, sensitive, and accurate liquid chromatographic method with electrochemical detection (ECD) has been developed for analysis of caffeic acid and caffeic acid phenyl ester in propolis. ECD is more sensitive than diode-array detection for analysis of the two phenolic compounds. The calibration plots show linearity is good (r ≥ 0.9995) in the ranges tested. Recovery was from 93.5 to 97.4%.     

Structure and Internal Rotation of Molecules in Chlorodifluoroacetic Acid

The potential energy profile of an isolated CF₂ClCOOH molecule with a CF₂Cl group rotating around the C–C bond was determined by the Hartree–Fock method using the 6-31G(d) basis set. Barriers to internal rotation were estimated for this molecule; its geometrical parameters were found for the equilibrium and transition states that are due to the torsion potential with unequal wells. Crystal effect on CF₂Cl reorientations in solid chlorodifluoroacetic acid has been evaluated.     

丙酮酸分子结构与振动光谱的密度泛函理论研究

用密度泛函方法BLYP、B3LYP和从头算Hartree-Fock(HF)方法在6-31G*基组水平上对丙酮酸分子的几何结构（甲基的重叠式和交错式两种构象）和振动光谱分别进行了优化和计算,并给出了各种频率所对应的红外强度及拉曼活性,对光谱进行了指认。结果表明：在丙酮酸分子的两种构象中,重叠式比较稳定*B3LYP计算得到的构型参数与实验结果比较一致;在振动频率的计算中,BLYP未标度力场所计算的非CH3伸缩振动基频预测值和实验值的平均绝对偏差为10.4cm-1;而HF标度力场的平均绝对偏差为17.9cm-1。说明两者的结果与实验观测频率比较吻合,但B3LYP的频率计算值偏差（38.3cm-1）较大。根据振动频率的势能分布和红外光谱强度对此分子的振动基频进行了理论归属。     

Density Functional Theory Study of Infrared and Ultraviolet Spectra of Urea L-Malic Acid

During the measurement of atmospheric nitrate radical by long-path differential optical absorption spectroscopy, water vapor strong absorption could affect the measurement of nitrate radical and detection limits of system. Under the tropospheric condition, the optical density of water vapor absorption is non-linearly dependent on column density. An effective method was developed to eliminate the effect of water vapor absorption. Reference spectra of water vapor based on the daytime atmospheric absorption spectra, when fitted together with change of cross section with water vapor column densities, gave a more accurate fitting of water vapor absorptions, thus its effect on the measurements of nitrate radical could be restricted to a minimum and detection limits of system reached 3.6 ppt. The modified method was applied during an intensive field campaign in the Pearl River Delta, China. The NO3 concentration in polluted air masses varied from 3.6 ppt to 82.5 ppt with an average level of 23.6±1.8 ppt.     

巴比妥酸的质子迁移异构化反应及其芳香性的密度泛函研究

用密度泛函理论方法,在B3LYP/6-31G水平上研究了巴比妥酸的醇式和酮式结构的分子内质子迁移异构化反应,对反应势能面的研究发现,该化合物可能存在6种中间体和12种过渡态,对反应物、产物和过渡态进行了几何构型优化,在同一水平下计算了单点能量,并用频率振动模式和内稟坐标(IRC)确证了过渡态的存在,反应的活化能最小为107.50kJ·mol-1,最大为330.84kJ·mol-1。在GIAO-HF/6-31G//B3LYP/6-31G水平对反应物、产物,中间体和过渡态的非相关性化学位移(NICS)进行了计算,利用NICS值对反应过程中芳香性进行分析,发现如果分子中的C5与C6存在双键或者形成过渡态的四员环状结构,分子就有不同程度的芳香性,如果有C6上的羟基氢迁移到C5上,C5-C6以单键相结合,则分子就没有芳香性。     

Structural Investigation of DHICA Eumelanin Using Density Functional Theory and Classical Molecular Dynamics Simulations

Eumelanin is an important pigment, for example, in skin, hair, eyes, and the inner ear. It is a highly heterogeneous polymer with 5,6-dihydroxyindole-2-carboxylic acid (DHICA) and 5,6-dihydroxyindole (DHI) building blocks, of which DHICA is reported as the more abundant in natural eumelanin. The DHICA-eumelanin protomolecule consists of three building blocks, indole-2-carboxylic acid-5,6-quinone (ICAQ), DHICA and pyrrole-2,3,5-tricarboxylic acid (PTCA). Here, we focus on the self-assembly of DHICA-eumelanin using multi-microsecond molecular dynamics (MD) simulations at various concentrations in aqueous solutions. The molecule was first parameterized using density functional theory (DFT) calculations. Three types of systems were studied: (1) uncharged DHICA-eumelanin, (2) charged DHICA-eumelanin corresponding to physiological pH, and (3) a binary mixture of both of the above protomolecules. In the case of uncharged DHICA-eumelanin, spontaneous aggregation occurred and water molecules were present inside the aggregates. In the systems corresponding to physiological pH, all the carboxyl groups are negatively charged and the DHICA-eumelanin model has a net charge of −4. The effect of K+ ions as counterions was investigated. The results show high probability of binding to the deprotonated oxygens of the carboxylate anions in the PTCA moiety. Furthermore, the K+ counterions increased the solubility of DHICA-eumelanin in its charged form. A possible explanation is that the charged protomolecules favor binding to the K+ ions rather than aggregating and binding to other protomolecules. The binary mixtures show aggregation of uncharged DHICA-eumelanins; unlike the charged systems with no aggregation, a few charged DHICA-eumelanins are present on the surface of the uncharged aggregation, binding to the K+ ions.     

Computational Study of Vertical Ionization Potentials Using Density Functional Theory and Green's Function Methods

Over one hundred vertical ionization potentials (VIPs) were computed using density functional theory (DFT) and Green's function (GF) based methods. The DFT approaches include the unrestricted transition state (uTS) and unrestricted diffuse ionization (uDI) approximations using the Becke88‐Perdew86 exchange‐correlation functional. Green's function methods include the outer‐valence GF (OVGF) approach, the parametrized GF2 (pGF2), and the parametrized GF2 times screened interaction (pGW2) approximations. DFT computations of IPs using the uTS approximation was found to be nearly as accurate as those predicted using the elaborate OVGF method. The much more computationally efficient uDI approximation provides predictions of moderate accuracy and is recommended for computing IPs for larger molecules. We have observed that the average absolute deviations from a uDI calculation using poorer basis set (DZVP) and poorer geometry (AM1 optimization) is only slightly larger.     

The Density Functional Theory Account of Interplaying Long-Range Exchange and Dispersion Effects in Supramolecular Assemblies of Aromatic Hydrocarbons with Spin

Aromatic hydrocarbons with fused benzene rings and regular triangular shapes, called n-triangulenes according to the number of rings on one edge, form groundstates with n-1 unpaired spins because of topological reasons. Here, we focus on methodological aspects emerging from the density functional theory (DFT) treatments of dimer models of the n = 2 triangulene (called also phenalenyl), observing that it poses interesting new problems to the issue of long-range corrections. Namely, the interaction comprises simultaneous spincoupling and van der Waals effects, i.e., a technical conjuncture not considered explicitly in the benchmarks calibrating long-range corrections for the DFT account of supramolecular systems. The academic side of considering dimer models for calculations and related analysis is well mirrored in experimental aspects, and synthetic literature revealed many compounds consisting of stacked phenalenyl cores, with intriguing properties, assignable to their long-range spin coupling. Thus, one may speculate that a thorough study assessing the performance of state-of-the-art DFT procedures has relevance for potential applications in spintronics based on organic compounds.     

Accurate Quantification of Polycyclic Aromatic Hydrocarbons in Environmental Samples Using Deuterium-labeled Compounds as Internal Standards

For the accurate quantification of polycyclic aromatic hydrocarbons (PAHs) in environmental samples by isotope-dilution mass spectrometry (IDMS) and gas chromatography-mass spectrometry (GC-MS), we used deuterium-labeled PAHs (D-PAHs) as internal standards for microg g(-1)-level certified reference materials and corresponding calibration solutions. Although pressurized liquid extraction (PLE) causes significant biases in the analytical results for ng g(-1)-level samples (2.4 - 15%), biases in the analytical results by PLE (190 degrees C, 20 MPa, 2 cycles) were small (<2.3%) and negligible for microg g(-1)-level samples when expanded uncertainty (coverage factor k = 2) was considered.     

铜催化内炔的硅羧基化反应机理及区域选择性的密度泛函理论研究

用CO₂作为原料,在过渡金属催化下生成新的碳碳键是很重要的.在这类反应中,杂原子官能团和CO₂同时与不饱和的底物被催化生成高功能化的羧酸衍生物已经越来越受到人们的关注.本文采用密度泛函理论(DFT)方法,研究了金属铜催化剂催化内炔的硅羧基化反应机理.根据炔烃上甲基和苯环两个取代基的相对位置,提出了两条反应路径(path I:甲基和path II:苯环).计算结果表明炔烃插入Cu－Si键既是速率决速步骤也是区域选择决速步骤.在path I中,炔烃插入Cu－Si键的自由能为112.8 kJ·mol^-1,而在path II中为127.6kJ·mol^-1.显然, path I比path II在动力学上更有利,这与实验上两条路径对应产物的产率97: 3是一致的.分析表明区域选择性是由炔烃取代基甲基和苯环的电子效应决定的.     

运用密度泛函活性理论的信息论方法研究苯并富烯衍生物的芳香性

文献中已有越来越多的芳香性体系被发现,同时也有越来越多的芳香性指标被提出来,但是如何解释芳香化合物稳定性的起源以及理解芳香性的本质仍然是当今理论化学中一个悬而未决的难题。运用我们新近提出的密度泛函活性理论信息论方法,不久前我们曾对一系列富烯衍生物进行了系统研究并得到了一个全新的认识。本文进一步探讨苯并富烯衍生物的芳香性行为,目的在于考察一个或多个苯环与富烯连接之后其芳香性发生变化的情况。运用香农熵,费舍尔信息,Ghosh-Berkowitz-Parr熵,Onicescu信息能,信息增益,以及相对Rényi熵六个信息量,和四种芳香指标,ASE,HOMA,FLU和NICS,我们系统地研究了信息量和芳香性指标在单、双、三苯并富烯衍生物中的相关性。我们发现,不管是否有苯环与富烯相连,芳香指标和信息量的交叉相关性都是一样的。这表明,虽然苯环本身具有芳香性,但苯环与富烯相连并不能改变富烯的芳香性与反芳香性本质。苯并富烯衍生物与富烯衍生物的芳香性和反芳香性一致。苯并富烯衍生物的芳香性和反芳香性完全取决于富烯本身的芳香性和反芳香性。这些结果为认识和理解复杂体系芳香性和反芳香性起源和本质将提供有益的启示。     

Uranyl-Bound Tetra-Dentate Non-Innocent Ligands: Prediction of Structure and Redox Behaviour Using Density Functional Theory

Computational methods have been applied to understand the reduction potentials of [UO₂-salmnt-L] complexes (L=pyridine, DMSO, DMF and TPPO), and their redox behavior is compared with previous experiments in dichloromethane solution. Since the experimental results were inconclusive regarding the influence of the uranyl-bound tetra-dentate ‘salmnt’ ligand, here we will show that salmnt acts as a redox-active ligand and exhibits non-innocent behavior to interfere with the otherwise expected one-electron metal (U) reduction. We have employed two approaches to determine the uranyl (VI/V) reduction potentials, using a direct study of one-electron reduction processes and an estimation of the overall reduction using isodesmic reactions. Hybrid density functional theory (DFT) methods were combined with the Conductor-like Polarizable Continuum Model (CPCM) to account for solvation effects. The computationally predicted one-electron reduction potentials for the range of [UO₂-salmnt-L] complexes are in excellent agreement with shoulder peaks (∼1.4 eV) observed in the cyclic voltammetry experiments and clearly correlate with ligand reduction. Highly conjugated pi-bonds stabilize the ligand based delocalized orbital relative to the localized U f-orbitals, and as a consequence, the ligand traps the incoming electron. A second reduction step results in metal U(VI) to U(V) reduction, in good agreement with the experimentally assigned uranyl (VI/V) reduction potentials.     

Studying the Chemistry of Cationized Triacylglycerols Using Electrospray Ionization Mass Spectrometry and Density Functional Theory Computations

Analysis of triacylglycerols (TAGs), found as complex mixtures in living organisms, is typically accomplished using liquid chromatography, often coupled to mass spectrometry. TAGs, weak bases not protonated using electrospray ionization, are usually ionized by adduct formation with a cation, including those present in the solvent (e.g., Na⁺). There are relatively few reports on the binding of TAGs with cations or on the mechanisms by which cationized TAGs fragment. This work examines binding efficiencies, determined by mass spectrometry and computations, for the complexation of TAGs to a range of cations (Na⁺, Li⁺, K⁺, Ag⁺, NH₄ ⁺). While most cations bind to oxygen, Ag⁺ binding to unsaturation in the acid side chains is significant. The importance of dimer formation, [2TAG + M]⁺ was demonstrated using several different types of mass spectrometers. From breakdown curves, it became apparent that two or three acid side chains must be attached to glycerol for strong cationization. Possible mechanisms for fragmentation of lithiated TAGs were modeled by computations on tripropionylglycerol. Viable pathways were found for losses of neutral acids and lithium salts of acids from different positions on the glycerol moiety. Novel lactone structures were proposed for the loss of a neutral acid from one position of the glycerol moiety. These were studied further using triple-stage mass spectrometry (MS³). These lactones can account for all the major product ions in the MS³ spectra in both this work and the literature, which should allow for new insights into the challenging analytical methods needed for naturally occurring TAGs.

两种不同位置取代的对硝基苯乙烯基蒽的密度泛函研究

用量子化学从头算HF方法和DFT方法,对2-对硝基苯乙烯基蒽(2-pNSA)和9-对硝基苯乙烯基蒽(9-pNSA)进行计算.优化得到了平衡几何构型,并计算了简谐振动频率及在-180°～180°范围内旋转乙烯基和硝基苯基间的夹角得到了势能曲线,找到了2-pNSA的两种构象和9-pNSA的三种构象,讨论了分子的平面性及基态分子总能量.在PM3/CIS水平上计算了电子光谱,计算结果与实验相符合.     

基于多齿功能单体的三聚氰胺印迹材料的制备及应用

以三聚氰胺(MEL)为模板分子,采用沉淀聚合方法研究了以衣康酸(IA)为多齿功能单体,乙二醇二甲基丙烯酸酯(EGDMA)为交联剂时,组分之间摩尔比对三聚氰胺印迹聚合物(MIP)选择性的影响.研究发现,当MEL与IA摩尔比为1∶1.5时,所合成的MIP对三聚氰胺具有最佳选择性,选择性系数k=10.41(以三聚氰酸为对照物).此时IA的识别位点(羧基)与三聚氰胺的可识别位点(胺基)摩尔比为1∶1,因此,不会出现由过量功能单体引起的非特异性吸附.Scatchard分析结果表明,这种MIP能够在更宽的三聚氰胺浓度范围(1～100 mg/L)内只表现出特异性吸附.以灭蝇胺为对照物时,kNIP<1但kMIP>1,该结果证明MIP对分子的形状有良好的识别能力.最后,以制备的MIP为固相萃取吸附剂,结合高效液相色谱/紫外检测器(HPLC/UV)建立了三聚氰胺测定方法.当牛奶(奶粉)中MEL的加标浓度分别为15,600和4500μg/L(μg/kg)时检测的加标回收率为92.34%～109.4%.     

Hydrogen Bonding in the Dimer and Monohydrate of 2-Adamantanol: A Test Case for Dispersion-Corrected Density Functional Methods

Weakly-bound intermolecular clusters constitute reductionist physical models for non-covalent interactions. Here we report the observation of the monomer, the dimer and the monohydrate of 2-adamantanol, a secondary alcohol with a bulky ten-carbon aliphatic skeleton. The molecular species were generated in a supersonic jet expansion and characterized using broadband chirped-pulse microwave spectroscopy in the 2–8 GHz frequency region. Two different gauche-gauche O-H···O hydrogen-bonded isomers were observed for the dimer of 2-adamantanol, while a single isomer was observed for the monomer and the monohydrate. The experimental rotational parameters were compared with molecular orbital calculations using density functional theory (B3LYP-D3(BJ), B2PLYP-D3(BJ), CAM-B3LYP-D3(BJ), ωB97XD), additionally providing energetic and electron density characterization. The shallow potential energy surface makes the dimer an interesting case study to benchmark dispersion-corrected computational methods and conformational search procedures.