Ab initio chemical kinetics for the reactions of HNCN with O(P) and O2

The kinetics and mechanisms of the reactions of cyanomidyl radical (HNCN) with oxygen atoms and molecules have been investigated by ab initio calculations with rate constant prediction. The doublet and quartet state potential energy surfaces (PESs) of the two reactions have been calculated by single-point calculations at the CCSD(T)/6-311+G(3df, 2p) level based on geometries optimized at the CCSD/6-311++G(d, p) level. The rate constants for various product channels of the two reactions in the temperature range of 300-3000 K are predicted by variational transition state and RRKM theories. The predicted total rate constants of the O(³P) + HNCN reaction at 760 Torr Ar pressure can be represented by the expressions k_total (O + HNCN) = 3.12 × 10⁻¹⁰ × T^−0.05 exp (−37/T) cm³ molecule⁻¹ s⁻¹ at T = 300-3000 K. The branching ratios of primary channels of the O(³P) + HNCN are predicted: k₁ for producing the NO + CNH accounts for 0.72-0.64, k₂ + k₉ for producing the ³NH + NCO accounts for 0.27-0.32, and k₆ for producing the CN + HNO accounts for 0.01-0.07 in the temperature range studied. Meanwhile, the predicted total rate constants of the O₂ + HNCN reaction at 760 Torr Ar pressure can be represented by the expression, k_total(O₂ + HNCN) = 2.10 × 10⁻¹⁶ × T^1.28exp (−12200/T) cm³ molecule⁻¹ s⁻¹ at T = 300-3000 K. The predicted branching ratio for k₁₁ + k₁₃ producing HO₂ + ³NCN as the primary products accounts for 0.98-1.00 in the temperature range studied.     

High-pressure laminar flame speeds and kinetic modeling of carbon monoxide/hydrogen combustion

Laminar flame speeds were accurately measured for CO/H₂/air and CO/H₂/O₂/helium mixtures at different equivalence ratios and mixing ratios by the constant-pressure spherical flame technique for pressures up to 40 atmospheres. A kinetic mechanism based on recently published reaction rate constants is presented to model these measured laminar flame speeds as well as a limited set of other experimental data. The reaction rate constant of CO + HO₂ → CO₂ + OH was determined to be k = 1.15 × 10⁵T^2.278 exp(−17.55 kcal/RT) cm³ mol⁻¹ s⁻¹ at 300-2500 K by ab initio calculations. The kinetic model accurately predicts our measured flame speeds and the non-premixed counterflow ignition temperatures determined in our previous study, as well as homogeneous system data from literature, such as concentration profiles from flow reactor and ignition delay time from shock tube experiments.     

Theoretical study on the kinetics for OH reactions with CH3OH and C2H5OH

Kinetics and mechanisms for reactions of OH with methanol and ethanol have been investigated at the CCSD(T)/6-311 + G(3df, 2p)//MP2/6-311 + G(3df, 2p) level of theory. The total and individual rate constants, and product branching ratios for the reactions have been computed in the temperature range 200-3000 K with variational transition state theory by including the effects of multiple reflections above the wells of their pre-reaction complexes, quantum-mechanical tunneling and hindered internal rotations. The predicted results can be represented by the expressions k₁ = 4.65 × 10⁻²⁰ × T^2.68 exp(414/T) and k₂ = 9.11 × 10⁻²⁰ × T^2.58 exp(748/T) cm³ molecule⁻¹ s⁻¹ for the CH₃OH and C₂H₅OH reactions, respectively. These results are in reasonable agreements with available experimental data except that of OH + C₂H₅OH in the high temperature range. The former reaction produces 96-89% of the H₂O + CH₂OH products, whereas the latter process produces 98-70% of H₂O + CH₃CHOH and 2-21% of the H₂O + CH₂CH₂OH products in the temperature range computed (200-3000 K).     

Thermal decomposition of toluene: Overall rate and branching ratio

The two-channel thermal decomposition of toluene, C₆H₅CH₃ → C₆H₅CH₂ + H (1) and C₆H₅CH₃ → C₆H₅ + CH₃ (2), was investigated in shock tube experiments over the temperature range of 1400-1780 K at a pressure of 1.5 (±0.1) bar. Rate coefficients for reactions (1) and (2) were determined by monitoring benzyl radical (C₆H₅CH₂) absorption at 266 nm during the decomposition of toluene diluted in argon and modeling the temporal behavior of the benzyl concentration with a kinetic model. The first-order rate coefficients determined at a pressure of 1.5 bar are expressed by k₁(T) = 2.09 × 10¹⁵ exp (−87510 [cal/mol]/RT) [s⁻¹] and k₂(T) = 2.66 × 10¹⁶ exp (−97880 [cal/mol]/RT) [s⁻¹]. The resulting branching ratio, k₁/(k₁ + k₂), ranges from 0.8 at 1350 K to 0.6 at 1800 K.     

Metal ion-molecule kinetics at combustion temperatures: The reaction of Ca with O2

The high-temperature photochemistry (HTP) technique, previously used for reactions of neutral species, has been adapted to the study of atomic metal ion-molecule reactions. Ca⁺ ions were generated by 193 nm multi-photon photolysis of calcium acetyl acetonate and its pyrolysis fragments. The relative ion concentrations were monitored by laser-induced fluorescence at 393.4 nm. Ar was used as the bath gas. The data for the Ca⁺ + O₂ + M → CaO₂⁺ + M association reaction (1) are fitted by k₁(907-1425 K) = 3.5 × 10⁻³² exp(+3161 K/T) cm⁶ molecule⁻² s⁻¹. Combining with an approximate k₁(296 K) value in the literature leads to k₁(296-1425 K) = 5.8 × 10⁻²² (T/K)^−2.9 exp(−601 K/T) cm⁶ molecule⁻² s⁻¹. Over much of the observed temperature range reaction (1) has much smaller rate coefficients than the corresponding neutral Ca association reaction. Reaction (1) is shown to behave very similarly to the O₂ association reaction with neutral K atoms, with which Ca⁺ is iso-electronic. This suggests that the initial step is ion-pair complex formation of the superoxide Ca²⁺(O₂)⁻, which is also consistent with results from density functional calculations. The k₁ values are rationalized via Troe’s unimolecular formalism, which leads to good accord with the experiments.     

Optical characterization of fourfold (Td)- and sixfold (Oh)-transition-metal species in MgAl2O4:Co by time-resolved spectroscopy

This work investigates the origin of novel visible photoluminescence (PL) bands observed in the spinel MgAl₂O₄:Co²⁺. Besides the well-known fourfold-coordinated Co²⁺(T_d) PL at 670 nm [N.V. Kuleshov, V.P. Mikhailov, V.G. Scherbitsky, P.V. Prokoshin and K.V. Yumashev, J. Lumin. 55 (1993) 265.], a rich structured PL band at 686 nm was also observed that we associate with uncontrolled impurities of sixfold coordinated Cr³⁺(O_h) by time-resolved spectroscopy and lifetime measurements and their variation with temperature. We also show that the lifetime of the Co²⁺(T_d) emission at 670 nm varies from τ=6.7 μs to 780 ns on passing from T=10 to 290 K. This unexpected behaviour for T_d systems is related to the excited-state crossover (⁴T₁↔²E), making the emission band to transform from a narrow-like emission from ²E at low temperature to a broad structureless band from ⁴T₁ at room temperature.     

Kinetic modeling of the benzyl + HO2 reaction

Benzyl is a resonantly stabilized radical that commonly occurs as an intermediate in the combustion of aromatic compounds. The bimolecular reaction of benzyl with HO₂ is important in the oxidation of toluene, especially at low to moderate temperatures, where unimolecular decomposition of the benzyl radical is slow. We show that the addition of HO₂ to the methylene site in benzyl produces a vibrationally excited benzylhydroperoxide adduct, with over 60 kcal mol⁻¹ (251 kJ mol⁻¹) of excess energy above the ground state. RRKM simulations are performed on the benzyl + HO₂ reaction, using thermochemical and kinetic parameters obtained from ab initio calculations, with variational transition state theory (VTST) for treatment of barrierless radical + radical reaction kinetics. Our results reveal that the benzyl + HO₂ reaction proceeds predominantly to the benzoxyl radical + OH at temperatures of around 800 K and above, with the production of stabilized benzylhydroperoxide molecules dominating at lower temperatures. The heat of formation of the benzyl radical is calculated as 52.5 kcal mol⁻¹ (219.7 kJ mol⁻¹) at the G3B3 level of theory, in relative agreement with other recent determinations of this value.     

High-resolution FTIR measurement of the ν4 band of methylene fluoride-d2

A high-resolution (0.002 cm⁻¹) infrared absorption spectrum of methylene fluoride-d₂ (CD₂F₂) of the lowest fundamental mode ν₄ in the region from 460 to 610 cm⁻¹ has been measured on a Bruker IFS 120-HR Fourier transform infrared spectrometer. More than 3500 transitions have been assigned in this B-type band centered at 521.9 cm⁻¹. The data have been combined with upper state pure rotational measurements in a weighted least-squares fit to obtain molecular constants for the upper state resulting in an overall standard deviation of 0.00018 cm⁻¹. Accurate value for the band origin (521.9578036 cm⁻¹) has been obtained and inclusion of transitions with very high J (?60) and K_a (?34) values has resulted in improved precision for sextic centrifugal distortion constants, in particular D_K, H_KJ, and H_K.     

Structural and electrical properties of (Na0.85K0.15)0.5Bi0.5TiO3 thin films deposited on LaNiO3 and Pt bottom electrodes

(Na_0.85K_0.15)_0.5Bi_0.5TiO₃ thin films were deposited on LaNiO₃(LNO)/SiO₂/Si(1 0 0) and Pt/Ti/SiO₂/Si(1 0 0) substrates by metal-organic decomposition, and the effects of bottom electrodes LNO and Pt on the ferroelectric, dielectric and piezoelectric properties were investigated by ferroelectric tester, impedance analyzer and scanning probe microscopy, respectively. For the thin films deposited on LNO and Pt electrodes, the remnant polarization 2P_r are about 22.6 and 8.8 μC/cm² under 375 kV/cm, the dielectric constants 238 and 579 at 10 kHz, the dielectric losses 0.06 and 0.30 at 10 kHz, the statistic d_33eff values 95 and 81 pm/V. The improved piezoelectric properties could make (Na_1−xK_x)_0.5Bi_0.5TiO₃ thin film as a promising candidate for piezoelectric thin film devices.     

Dielectric properties and relaxation behavior of [TMA]2Zn0.5Cu0.5Cl4 compound

The [TMA]₂Zn_0.5Cu_0.5Cl₄ hybrid material was prepared and its dielectric spectra were measured in the 10⁻¹ Hz-10⁶ Hz frequency range and 200-305 K temperature interval. The dielectric permittivity showed a ferroelectric-paraelectric phase transition at 293 K. Double relaxation peaks are observed in the imaginary part of the electrical modulus, suggesting the presence of grain and grain boundary in the sample. The frequency dependent conductivity was interpreted in term of Jonscher's law: σ(ω)=σ_dc+Aωⁿ. The temperature dependent of the dc conductivity (σ_dc) was well described by the Arrhenius equation: σ_dcT=σ_o×exp(−E_a/kT).     

Intermolecular photoinduced electron-transfer processes between C60 and aniline derivatives in benzonitrile

The quenching behavior of the triplets of C₆₀ by various aniline derivatives (1a-d and 2a-e) was investigated by means of laser flash photolysis in benzonitrile at 293 K. Electron transfer process was proposed to be the main mechanism because of the direct detection of radical ions of aniline derivatives and C₆₀ in time-resolved transient absorption spectra. The quenching rate constants (k_q) of by different substrates determined at 740 nm approach or reach the diffusion-controlled limit. DFT method was employed to calculate the unknown oxidation potentials of substrates in solution. With these E_ox values, free energy changes (ΔG) were obtained through Rehm-Weller equation. Dependence of observed quenching rate constants on the free energy changes further indicates the photoinduced reactions between ³C₆₀^* and substrates proceed through an electron transfer mechanism. Obtained k_q values for the aniline derivatives are impacted obviously by ground-state configurations and the kinds substituents quantified by Hammett σ constant. Good correlation between log k_q and σ values conforms to the empirical Hammett equation. A more negative ρ value (−3.356) was gained for anilines (2a-e) than that of N,N-dimethylanilines (1a-d) (−1.382), which suggests a more susceptible reactivity for the former substrates. Charge density distribution of reaction center “N” originated from quantum calculation supports this suggestion. In addition, a relationship between quenching rate constants and solvent viscosity was gained from C₆₀/dimethyl-p-toluidine system in altered mixtures of acetonitrile and toluene.     

Interaction between titanium dioxide nanoparticles and human serum albumin revealed by fluorescence spectroscopy in the absence of photoactivation

Titanium dioxide (TiO₂) nanoparticles (NPs) are widely used as an important kind of biomaterials. The interaction between TiO₂ (P25) at 20 nm in diameter and human serum albumin (HSA) was studied by fluorescence spectroscopy in this work. Under the simulative physiological conditions, fluorescence data revealed the presence of a single class of binding site on HSA and its binding constants (K_a) were 2.18±0.04×10⁴, 0.87±0.05×10⁴, 0.68±0.06×10⁴ M⁻¹ at 298, 304 and 310 K, respectively. In addition, according to the Van’t Hoff equation, the thermodynamic functions standard enthalpy (ΔH⁰) and standard entropy (ΔS⁰) for the reaction were calculated to be −75.18±0.15 kJ mol⁻¹ and −170.11±0.38 J mol⁻¹ K⁻¹. These results indicated that TiO₂ NPs bond to HSA mainly by van der Waals force and hydrogen bonding formation in low dielectric media, and the electrostatic interactions cannot be excluded. Furthermore, the effects of common ions on the binding constant of TiO₂ NPs-HSA complex were discussed.     

A DFT study the solvent effects of H2 adsorption on Cu(h k l) surface

Density functional theory (DFT) combined with conductor-like solvent model (COSMO) have been performed to study the solvent effects of H₂ adsorption on Cu(h k l) surface. The result shows H₂ can not be parallel adsorbed on Cu(h k l) surface in gas phase and only vertical adsorbed. At this moment, the binding energies are small and H₂ orientation with respect to Cu(h k l) surfaces is not a determining parameter. In liquid paraffin, when H₂ adsorbs vertically on Cu(h k l) surface, solvent effects not only influences the adsorptive stability, but also improves the ability of H₂ activation; When H₂ vertical adsorption on Cu(h k l) surface at 1/4 and 1/2 coverage, H-H bond is broken by solvent effects. However, no stable structures at 3/4 and 1 ML coverage are found, indicating that it is impossible to get H₂ parallel adsorption on Cu(h k l) surfaces at 3/4 and 1 ML coverages due to the repulsion between adsorbed H₂ molecules.     

Perturbed angular correlation studies of Ta hyperfine interactions in Hf-Ni and Zr-Ni compounds

The hyperfine interaction experienced by ¹⁸¹Ta nuclei in the intermetallic compounds ZrNi₅, HfNi₅, and Hf₂Ni₇ has been investigated by perturbed angular correlation (PAC) spectroscopy. At temperatures T≥15 K the ¹⁸¹Ta angular correlation of appropriately annealed ZrNi₅ and HfNi₅ is unperturbed, indicating the absence of a magnetic hyperfine interaction. This observation rules out the possibility of spontaneous magnetic order of ZrNi₅ and HfNi₅ recently proposed in the literature. The temperature dependence of the electric quadrupole interaction of ¹⁸¹Ta in Hf₂Ni₇ suggests the existence of a reversible phase transformation at T≥500 K.     

On the kinetics of the C5H5 + C5H5 reaction

The possibility that the reaction between two cyclopentadienyl radicals (cC₅H₅) may lead to the production of naphthalene has been the subject of considerable theoretical and experimental studies. Though it has been proposed that this reaction may be the main channel for the formation of naphthalene in many combustion environments, the elementary mechanism leading from the initial adduct (C₅H₅_C₅H₅) to naphthalene is still not clear. In this study the portion of the C₁₀H₁₀ PES accessible to C₅H₅_C₅H₅ has been theoretically re-examined using density functional theory to locate stationary points and the CBS-QB3 computational protocol to determine energies. A new reaction pathway leading to the formation of a set of azulyl radicals was identified. Since it is known that azulyl radicals can easily decompose to naphthalene and atomic H, the proposed pathway provides an effective route for the formation of naphthalene. Channel specific kinetic constants were determined between 1100 K and 2000 K integrating the master equation for a PES comprising both this reaction pathway and the literature reaction pathway, which main product is the fulvalenyl radical. It was found that the main reaction channel is decomposition to reactants in the whole temperature range investigated and that the azulyl reaction channel is dominant over the fulvalenyl pathway up to 1450 K. The rate constants calculated at 1 bar for the azulyl and fulvalenyl reaction channels are 10^14.72T(K)^?0.853 exp(?3650/T(K)) and 10^10.30T(K)^0.951 exp(?7948/T(K)) cm³/mol/s, respectively. The rate constant for the formation of naphthalene through the azulyl channel is consistent with recent estimates based on the kinetic simulation of the pyrolysis and oxidation of cyclopentadiene.     

Kinetics for the reaction of phenyl radical with phenylacetylene and styrene

The kinetics for the reactions of C₆H₅ with phenylacetylene and styrene have been measured by CRDS in the temperature range 297-409 K under an Ar pressure of 3.6 Torr. The total rate constants can be given by the following Arrhenius expressions (in units of cm³ mol⁻¹ s⁻¹): k₁(C₆H₅ + C₆H₅C₂H) = 10^13.0±0.1exp [−(2430 ± 150)/RT] and k₂(C₆H₅ + C₆H₅C₂H₃) = 10^13.3±0.1 exp [−(2570 ± 180)/T]. Additional DFT and MP2 calculations have been carried out to assist our interpretation of the measured kinetic data. The addition of C₆H₅ to the terminal CH_x (x = 1 or 2) sites is predicted to be the dominant channel for both reactions. The calculated bimolecular rate constants are in reasonable agreement with experimental values for the temperature range studied.     

Short-range order analysis and some physical properties of InxSe1−x glasses

Bulk In_xSe_1−x (with x=5-25 at%) glasses were prepared using the melt-quench technique. Short range order(SRO) was examined by the X-ray diffraction using Cu(k_α) radiation in the wave vector interval 0.28≤k≤6.5 A⁰⁻¹.The SRO parameters have been obtained from the radial distribution function. The inter-atomic distance obtained from the first and second peak are r₁=0.263 and r₂=0.460 nm, which is equivalent In-Se and Se-Se bond length. The fundamental structural unit for the studied glasses is In₂Se₃ pyramid. Using the differential scanning calorimetry (DSC), the crystallization mechanism of In_xSe_1−x chalcogenide glass has been studied. The glass transition activation energy (E_g) is 289±0.3 kj/mol.There is a correlation amongst the glass forming ability, bond strength and the number of lone pair electrons. The utility of the Gibbs-Di Marzio relation was achieved by estimating T_g theoretically.     

Iron isotope effect on Tc in optimally-doped (Ba,K)Fe2As2 (Tc = 38 K) and SmFeAsO1−y (Tc = 54 K) superconductors

We report the iron isotope effect on a transition temperature (T_c) in an optimally-doped (Ba,K)Fe₂As₂ (T_c = 38 K) and SmFeAsO_1−y (T_c = 54 K) superconductors. In order to obtain the reliable isotope shift in T_c, twin samples with different iron isotope mass are synthesized in the same conditions (simultaneously) under high-pressure. We have found that (Ba,K)Fe₂As₂ shows an inverse iron isotope effect α_Fe = −0.18 ± 0.03 while SmFeAsO_1−y shows a small iron isotope effect α_Fe = −0.02 ± 0.01, where the isotope exponent α is defined by T_c ∼ M^−α (M is the isotopic mass). The results show that α_Fe changes in the iron-based superconductors depending on the system. The distinct iron isotope effects imply the exotic coupling mechanism in the iron-based superconductors.     

High-resolution infrared analysis of the ν7 band of cis-ethylene-d2 (cis-C2H2D2)

The spectrum of the ν₇ band of cis-ethylene-d₂ (cis-C₂H₂D₂) has been recorded with an unapodized resolution of 0.0063 cm⁻¹ in the 740-950 cm⁻¹ region using a Bruker IFS 125 HR Fourier transform infrared spectrometer. By fitting 2186 infrared transitions of ν₇ with a standard deviation of 0.00060 cm⁻¹ using a Watson’s A-reduced Hamiltonian in the I^r representation, accurate rovibrational constants for ν₇ = 1 state have been derived. The band center of ν₇ has been found to be 842.20957 ± 0.00004 cm⁻¹. In a simultaneous fit of 1331 infrared ground state combination differences from the present ν₇ transitions, together with 22 microwave frequencies, ground state constants have been improved. The rms deviation of the ground state fit was 0.00027 cm⁻¹.