Quantum mechanical estimation of Abraham hydrogen bond parameters using 1:1 donor–acceptor complexes

Hydrogen bond donor strength () and acceptor strength () have been successfully used in models of many environmental and chemical systems, and a number of computational methods have been developed to predict them. In this work, a quantum chemical Møller–Plesset perturbation (MP2) method is applied to estimate the binding free energies (ΔG_hbond) of several 1:1 hydrogen‐bonded complexes. A correlation between the binding free energies and hydrogen bond strength is established. This relationship can be used to develop an accurate computational model for predicting and using binding free energies. The accuracy of the method in predicting Abraham (root mean squared deviation (RMSD) = 0.0693) and (RMSD = 0.0677) are comparable to the empirical, fragment‐based ABSOLV method (RMSD = 0.1144 and 0.1281 for and , respectively). The binding free energy has been decomposed into its thermodynamic components, and it is demonstrated that the linear relationships in the dataset and the existence of magic point can be attributed to a constant entropy of reaction. Copyright © 2014 John Wiley & Sons, Ltd.     

Synthesis and thermochemical study of quinoxaline‐N‐oxides: enthalpies of dissociation of the N–O bond

The synthesis of three new quinoxaline mono‐N‐oxides derivatives, namely, 2‐tert‐butoxycarbonyl‐3‐methylquinoxaline‐N‐oxide, 2‐phenylcarbamoyl‐3‐ethylquinoxaline‐N‐oxide, and 2‐carbamoyl‐3‐methylquinoxaline‐N‐oxide, from their corresponding 1,4‐di‐N‐oxides is reported. Samples of these compounds were used for a thermochemical study, which allowed derivation of their gaseous standard molar enthalpies of formation, , from their enthalpies of formation in the condensed phase, , determined by static bomb combustion calorimetry, and from their enthalpies of sublimation, , determined by Calvet microcalorimetry. Finally, combining the for the quinoxaline‐N‐oxides derived in this work with literature values for the corresponding 1,4‐di‐N‐oxides and atomic oxygen, the bond dissociation enthalpies for cleavage of the first N?O bond in the di‐N‐oxides, DH₁(N–O), were obtained and compared with existing data. Copyright © 2011 John Wiley & Sons, Ltd.     

Kinetics and mechanisms of gas‐phase decarbonylation of α‐methyl‐trans‐cinamaldehyde and E‐2‐methyl‐2‐pentenal under homogeneous catalysis of hydrogen chloride

The kinetics of the gas‐phase elimination of α‐methyl‐trans‐cinamaldehyde catalyzed by HCl in the temperature range of 399.0–438.7 °C, and the pressure range of 38–165 Torr is a homogeneous, molecular, pseudo first‐order process and undergoing a parallel reaction to produce via (A) α‐methylstyrene and CO gas and via (B) β‐methylstyrene and CO gas. The decomposition of substrate E‐2‐methyl‐2‐pentenal was performed in the temperature range of 370.0–410.0 °C and the pressure range of 44–150 Torr also undergoing a molecular, pseudo first‐order reaction gives E‐2‐pentene and CO gas. These reactions were carried out in a static system seasoned reactions vessels and in the presence of toluene free radical inhibitor. The rate coefficients are given by the following Arrhenius expressions:

• Products formation from α‐methyl‐trans‐cinamaldehyde
• α‐methylstyrene :
• β‐methylstyrene :
• Products formation from E‐2‐methyl‐2‐pentenal
• E‐2‐pentene :

The kinetic and thermodynamic parameters for the thermal decomposition of α‐methyl‐trans‐cinamaldehyde suggest that via (A) proceeds through a bicyclic transition state type of mechanism to yield α‐methylstyrene and carbon monoxide, whereas via (B) through a five‐membered cyclic transition state to give β‐methylstyrene and carbon monoxide. However, the elimination of E‐2‐methyl‐2‐pentenal occurs by way of a concerted cyclic five‐membered transition state mechanism producing E‐2‐pentene and carbon monoxide. The present results support that uncatalyzed α‐β‐unsaturated aldehydes decarbonylate through a three‐membered cyclic transition state type of mechanism. Copyright © 2014 John Wiley & Sons, Ltd.     

Determining the excited‐state substituent constants of furyl and thienyl groups

Six series of styrene derivatives XCH═CHArY (total of 65) containing the styrene parent molecular skeleton were synthesized (here, Y is OMe, Me, H, F, Cl, CF₃, CN, and NO₂, and X is 2‐furyl, 3‐furyl, 2′‐methyl‐2‐furyl, 2‐thienyl, 3‐thienyl, and 2′‐methyl‐2‐theniyl). Their ultraviolet absorption spectra were measured in anhydrous ethanol, and their wavelength of absorption maximum λ_max was recorded. For the wavenumber ν_max (cm^?1, ν_max = 1/λ_max) of the obtained λ_max, a quantitative correlation analysis was performed, and 6 excited‐state substituent constants of groups X were obtained by means of curve‐fitting method. Taking the ν_max values of total 90 compounds of styrene derivatives as a data set (including 25 compounds from reference and 65 compounds of this work), a quantitative correlation analysis was performed, and the reliability of the obtained was verified. In addition, 12 samples of disubstituted Schiff bases (XCH═NArY) involving the above groups X were synthesized, and their ν_max values were recorded. Using these 12 ν_max together with the 14 ν_max values of Schiff bases taken from reference (total of 26 compounds), it was further verified that the values are reliable by means of quantitative correlation method.     

The substituent effect on the UV energy of 4,4′‐disubstituted benzylideneanilines

In this paper, 72 samples of disubstituted benzylideneanilines were all synthesized, and their UV data were measured in anhydrous ethanol. In the study on the UV energy of the titled compounds with single substituent changed, for the effect of the aniline substituent Y on the UV wavenumbers, its UV data can be correlated with a dual‐parameter equation; for the effect of benzylidene substituent X on the UV data, its UV energy can be correlated with a single‐parameter equation (Y is an electron‐withdrawing group and H) or a dual‐parameter equation (Y is an electron‐donating group). In the study on the UV energy of model compounds with double substituents changed, a correlation equation between the UV absorption wavenumbers and substituent constants and σ_p, was obtained. For 72 samples of 4,4′‐disubstituted benzylideneanilines, the correlation coefficient was 0.9876, and the standard deviation s was only 358.46 cm^–1. The equation can be used to predict well the UV energy of BA derivatives. It was found that Δσ² is a better parameter than σ_XY to scale the substituent cross‐interaction effect on the UV wavenumbers of benzylideneanilines molecules. The results implied that the law of substituent effect on the UV energy of titled compounds is different from that of substituted stilbenes, and it is helpful to understand the effect of substituent effects on the chemical and physical properties of conjugated compounds with an imine bridging group (C = N) or a nonplanar parent. Copyright © 2011 John Wiley & Sons, Ltd.     

Density functional theory study of methoxide promoted and Zn(II)‐complexed methoxide promoted cleavages of aryl‐ and alkyl acetates in methanol. Transition from concerted to stepwise processes as a function of leaving group ability

Density functional theory calculations were performed for the methanolysis reactions of a set of aryloxy and alkoxy acetates ( 1a , 1b , 1c , 1d , 1e , 1f , 1g , 1h , 1i , 1j , 1k , 1l , 1m ) promoted by methoxide and a 1,5,9‐triazacyclododecane‐complexed Zn^(II)‐methoxide [2(OCH₃)]⁺ in order to give free energies and structural data for the various intermediates and transition states along the reaction pathway. The methoxide‐promoted reactions experience a transition of pathways from enforced‐concerted addition of CH₃O^? to the C = O unit for substrates having a good aryloxy leaving groups (LGs) with strong electron withdrawers ( 1a , 1b , 1c , 1d , 1e ) to a two step process with rate‐limiting CH₃O^‐ attack on aryloxy acetates having higher (the pK_a of the parent phenol of the LG in methanol) values. Only in the case of the substrates 1i‐m having alkoxy LGs is there an observed change in rate‐limiting step that occurs at the quasi‐symmetrical point where the . The methanolysis process for the 2,4‐dinitrophenoxy substrate ( 1a ) promoted by [2(OCH₃)]⁺ involves transient binding of the substrate to the metal complex followed by a rate‐limiting, enforced‐concerted attack of Zn^(II)‐coordinated ^–OCH₃, with fast breakdown of an addition intermediate that does not have a significant lifetime. For substrates 1b,c having slightly less electron withdrawing substituents, the reaction has two steps with rate‐limiting attack and an unassisted LG departure. As the increases, the reaction still has two steps with rate‐limiting attack, but departure of the LG is now assisted by its coordination to the metal ion. For alkoxy containing substrates, a change in rate‐limiting step occurs centered at methoxy acetate, 1j , (when ) for which the second step of metal ion assisted departure of methoxide becomes partially rate‐limiting. The Brønsted plots computed for the methoxide‐promoted and [2(OCH₃)]⁺‐promoted methanolyses are compared with the previously determined experimental data and are analyzed as arising not from a common line attributable to all substrates but rather in terms of separate, but intersecting, plots for aryl‐ and alkyl acetates. Copyright © 2014 John Wiley & Sons, Ltd.     

Theoretical kinetic study on the decomposition of 1,5‐diaminotetrazole

1,5‐Daminotetrazole (DAT) is of much interest because of the practical significance and the diversity of characteristics. The study on the decomposition pathway and the kinetics of DAT has been performed based on the quantum chemistry theory. The minimum energy path (MEP) calculation has shown that NH₂N₃ and NH₂CN are the initially detected products of DAT. And the structures of reactant, products and transition state were optimized with MP2 methods using 6‐311G** basis sets, and the energies were refined using CCSD(T)/6‐311G** levels of theory. The calculated rate constants were obtained using the conventional transition‐state theory (TST) and the canonical variational transition‐state theory (CVT) methods. The calculation results indicated that the energy barrier of decomposition reaction is 47.98 kcal mol^?1 and the variational effect is small. In addition, the rate constants and the Arrhenius experience formula of DAT decomposition have been obtained between 200 and 2500 K temperature regions. The fitted three‐parameter expressions calculated using the TST and CVT methods are (TST) and (CVT). This work may provide the theoretical support for further experimental synthesis and testing. Copyright © 2015 John Wiley & Sons, Ltd.     

Gas‐phase elimination kinetics of selected aliphatic α,β‐unsaturated aldehydes catalyzed by hydrogen chloride

The gas‐phase elimination of 2‐methyl‐2‐propenal catalyzed by HCl yields propene and CO gas, while E‐2‐pentenal with the same catalyst gives butene and CO gas. The kinetics determinations were carried out in a static system with the reaction vessels deactivated with allyl bromide and the presence of the free radical inhibitor toluene. Temperature and pressure ranges were 350.0–410.0 °C and 34–76 Torr. The elimination reactions are homogeneous and unimolecular, and follow a first‐order rate law. The rate coefficients for the reactions are expressible by the following Arrhenius equations: Data from the kinetic and thermodynamic parameters of these catalyzed elimination reactions implies a mechanism of a concerted five‐membered cyclic transition state structure for the formation of the corresponding olefin and carbon monoxide. Copyright © 2015 John Wiley & Sons, Ltd.     

Solution properties of phenothiazine drug promazine hydrochloride with cationic hydrotropes in aqueous/electrolyte solution at different temperature

The current work deals with the mixed micellization phenomena of surface active promazine hydrochloride (PMZ) drug with cationic hydrotropes (para‐toluidine hydrochloride and ortho‐toluidine hydrochloride) in absence and occurrence of 50 mmol kg^?1 NaCl at five different temperature (293.15–313.15 K). PMZ is an amphiphilic phenothiazine drug and employed for the cure of mania and schizophrenia. Conductometry measurement was employed to gain a detailed picture of the interactions between drug and hydrotrope molecules. The experimental data were analyzed according to different mixing models within the outline of the pseudophase separation model. The evaluated values of critical micelle concentration (cmc) were found to be inferior than cmc^id values signifying attractive interactions involving the both components in the solutions. NaCl further reduces the cmc of pure amphiphiles and their mixed systems as a result of screening of the electrostatic repulsion between the polar head groups. The micellar mole fractions (X₁) of hydrotropes evaluated by various proposed models were constantly more than ideal values ( ) signifying high involvements of hydrotrope in mixed micelles. Activity coefficients ( and ) were always below one in all cases signifying synergism in mixed micelles. Thermodynamic parameters favor the process of micellization which is found to be entropy driven. The negative values of free energies of mixing demonstrated the stability of the mixed systems of drug and hydrotrope. Copyright © 2016 John Wiley & Sons, Ltd.     

Orientations of ZnO grown on GaN

On semipolar epitaxial ZnO grown by chemical vapor deposition consists of two distinct orientations as evidenced by transmission electron microscopy and X‐ray diffraction. The initially grown ZnO on GaN follows the GaN lattice with the epitaxial relationship of // and The other oriented ZnO domains then grow on faceted with and with good coherency with the ‐oriented grains. (© 2015 WILEY‐VCH Verlag GmbH &Co. KGaA, Weinheim)     

A generalized G‐function for the Quantum Rabi Model

The analytical solution of the quantum Rabi model is based on a transcendental function , the zeros of which determine the eigenenergies. is generalized here to a function , which allows a much better numerical control of the high‐energy part of the spectrum by an appropriate choice of the complex parameter z. Additionally, it is shown that all zeros of correspond to eigenvalues of the Hamiltonian as well as the zeros of for imaginary z.     

Theory of the Robin quantum wall in a linear potential. II. Thermodynamic properties

A theoretical analysis of the thermodynamic properties of the Robin wall characterized by the extrapolation length Λ in the electric field that pushes the particle to the surface is presented both in the canonical and two grand canonical representations and in the whole range of the Robin distance with the emphasis on its negative values which for the voltage‐free configuration support negative‐energy bound state. For the canonical ensemble, the heat capacity at exhibits a nonmonotonic behavior as a function of the temperature T with its pronounced maximum unrestrictedly increasing for the decreasing fields as and its location being proportional to . For the Fermi‐Dirac distribution, the specific heat per particle is a nonmonotonic function of the temperature too with the conspicuous extremum being preceded on the T axis by the plateau whose magnitude at the vanishing is defined as , with N being a number of the particles. The maximum of is the largest for and, similar to the canonical ensemble, grows to infinity as the field goes to zero. For the Bose‐Einstein ensemble, a formation of the sharp asymmetric feature on the ‐T dependence with the increase of N is shown to be more prominent at the lower voltages. This cusp‐like dependence of the heat capacity on the temperature, which for the infinite number of bosons transforms into the discontinuity of , is an indication of the phase transition to the condensate state. Some other physical characteristics such as the critical temperature and ground‐level population of the Bose‐Einstein condensate are calculated and analyzed as a function of the field and extrapolation length. Qualitative and quantitative explanation of these physical phenomena is based on the variation of the energy spectrum by the electric field.     

Conformation/Tautomerization effect on the pKa values of lumazine and 6‐thienyllumazine

Density functional theory [B3LYP/6‐311+G(d,p)] was used in combination with the conductor‐like polarizable continuum model (CPCM) solvation model to investigate the relative stability and site‐specific values of neutral and ionized tautomers of lumazine (LM) and 6‐thienylLM (TLM). Two types of populations should be taken into consideration when calculating the , tautomers, and conformers. The major tautomer of neutral LM in aqueous solution is 13‐LM (the 13 notation refers to the acidic protons being in positions 1 and 3 of LM) TLM has decreased acidity at N₈ relative to LM. Further, the trans conformer of TLM is more acidic than cis. Similar to the case of LM, for TLM, N₁ is more acidic than N₃ in the uracil part. However, N₈ is predicted to be a stronger acid than N₁ for TLM. This acidity enhancement is essentially because of a specific stabilization of the anion when the thienyl group replaces H. Two factors are responsible for the acidity strength of N₈: The thienyl ring upon deprotonation acts inductively as an electron‐withdrawing group, and the excess electron density is dispersed better when the system is trans and contains second‐row atoms. Accurate pK_a calculation requires that all conformers/tautomers be included into the calculation. Copyright © 2014 John Wiley & Sons, Ltd.     

CTED: the new aromaticity index based on corrected total electron density at bond critical points

In this work, the corrected total electron density based on ellipticity (ε) at C–C bond critical points in a given ring and bond length alternations was introduced to estimate π‐electron density distributions in the ring. Then, to evaluate aromaticity of rings with any number of members, the was normalized relative to a system assumed as a full aromatic, which is named as the corrected total electron density (CTED) aromaticity index. For a wide range of aromatic, nonaromatic and antiaromatic compounds, we have compared CTED index with the other commonly used aromatic indices, such as HOMA, PDI, FLU, NICS and recently introduced EL. CTED index was seen to be in agreement with the defined indices, and with general expectations. Hence, as similar to the other indices except from PDI, we have proposed that CTED index could be applied to study the aromaticity of rings without any restriction in the number of members of rings and used to analyze both the local and global aromatic character of rings as a new aromaticity index. Copyright © 2014 John Wiley & Sons, Ltd.     

Theory of the Robin quantum wall in a linear potential. I. Energy spectrum,polarization and quantum‐information measures

Information‐theoretical concepts are employed for the analysis of the interplay between a transverse electric field applied to a one‐dimensional surface and Robin boundary condition (BC), which with the help of the extrapolation length Λ zeroes at the interface a linear combination of the quantum mechanical wave function and its spatial derivative, and its influence on the properties of the structure. For doing this, exact analytical solutions of the corresponding Schrödinger equation are derived and used for calculating energies, dipole moments, position and momentum quantum information entropies and their Fisher information and and Onicescu information energies and counterparts. It is shown that the weak (strong) electric field changes the Robin wall into the Dirichlet, (Neumann, ), surface. This transformation of the energy spectrum and associated waveforms in the growing field defines an evolution of the quantum‐information measures; for example, it is proved that for the Dirichlet and Neumann BCs the position (momentum) quantum information entropy varies as a positive (negative) natural logarithm of the electric intensity what results in their field‐independent sum . Analogously, at and the position and momentum Fisher informations (Onicescu energies) depend on the applied voltage as () and its inverse, respectively, leading to the field‐independent product (). Peculiarities of their transformations at the finite nonzero Λ are discussed and similarities and differences between the three quantum‐information measures in the electric field are highlighted with the special attention being paid to the configuration with the negative extrapolation length.

     

Refined comparison theorems for the Dirac equation in d dimensions

A single spin‐1/2 particle obeys the Dirac equation in spatial dimension and is bound by an attractive central monotone potential which vanishes at infinity (in one dimension the potential is even). This work refines the relativistic comparison theorems which were derived by Hall 1 . The new theorems allow the graphs of the two comparison potentials and to crossover in a controlled way and still imply the spectral ordering for the eigenvalues at the bottom of each angular momentum subspace. More specifically in a simplest case we have: in dimension , if , then ; and in dimensions, if , where and , then .

     

Nodal theorems for the Dirac equation in d ≥ 1 dimensions

A single particle obeys the Dirac equation in spatial dimensions and is bound by an attractive central monotone potential that vanishes at infinity. In one dimension, the potential is even, and monotone for The asymptotic behavior of the wave functions near the origin and at infinity are discussed. Nodal theorems are proven for the cases and , which specify the relationship between the numbers of nodes n₁ and n₂ in the upper and lower components of the Dirac spinor. For , whereas for if and if where and This work generalizes the classic results of Rose and Newton in 1951 for the case Specific examples are presented with graphs, including Dirac spinor orbits      

Gas‐phase basicities of acetophenones toward lanthanum cation [La(OMe)]

The relative free energy changes (lanthanum cation basicity, LaCB[L₂]) for the reaction [La(OMe)₂]L ? La(OMe) + 2L were determined in the gas phase for m‐ and p‐substituted acetophenones based on the measurement of ligand exchange equilibria using an FT‐ICR mass spectrometer. The substituent effect on ΔLaCB[L₂] of acetophenone is described in terms of the Yukawa–Tsuno equation, ΔG = ρ(σ° + r⁺ Δ σ ), with a ρ value of ?11.2 and an r⁺ value of 0.49. From this result, a ρ value of ?7.0 and an r⁺ value of 0.49 were estimated for the monomeric complex [LLa(OMe)] with the aid of theoretical calculations. This ρ value was found to be significantly smaller than that for protonation, and even smaller than Li⁺ basicity. Such a small ρ value has been attributed to the largely ionic (ion–dipole interaction) nature of the bonding interaction between La(OMe) and the carbonyl oxygen atom and, in part, to the long distance between La(OMe) and the substituent. Contrary to the ρ value, the r⁺ value is identical in both La(OMe) and Li⁺ basicities, suggesting that the r⁺ value of 0.49 can be regarded as a limiting one in a series of Lewis cation basicities of the acetophenone system, H⁺ (0.86) > Me₃Si⁺ (0.75) > Me₃Ge⁺ (0.71) > Cu⁺ (0.60) > Li⁺ = La(OMe) (0.49). Since the binding interaction between La(OMe) or Li⁺ and a neutral ligand is mostly electrostatic, the moderate r⁺ was interpreted to result from the redistribution of the induced positive charge within the acetophenone moiety upon binding with a metal ion rather than transfer of positive charge from a metal ion to the aromatic moiety. Copyright © 2010 John Wiley & Sons, Ltd.