We continue the development of the in-house molecular dynamics software package SpaCIAL and test it for the computation of the collision-induced absorption coefficients for a neon (Ne) and krypton (Kr) gas mixture. An apodization procedure for the dipole autocorrelation function is implemented and tested. We also carry out a statistical study of the convergence rate with respect to ensemble size. The resulting absorption coefficients show a good accordance with quantum mechanical results. Comparison with laboratory measurements shows agreement within 10%–20% at T = 295 K. At T = 480 K, a larger difference of 40%–80% is observed, which can presumably be explained by experimental uncertainties. For the study, an empirical (Barker, Fisher, and Watts) interaction-potential [Mol. Phys. 21, 657 (1971)] for Ne–Kr has been developed. Ab initio {coupled cluster with singles and doubles (triples) [CCSD(T)]} potentials for Ne–Ne, Kr–Kr, and Ne–Kr have been computed, as well as the CCSD(T) interaction-induced Ne–Kr dipole moment curve.
We have implemented a scheme for classical molecular dynamics simulations of collision-induced absorption. The program has been applied to a gas mixture of argon (Ar) and krypton (Kr). The simulations are compared with accurate quantum dynamical calculations. The comparisons of the absorption coefficients show that classical molecular dynamics is correct within 10% for photon wave numbers up to 220 cm−1 at a temperature of 200 K for this system. At higher temperatures, the agreement is even better. Molecular dynamics accounts for many-body interactions, which, for example, give rise to continuous dimer formation and destruction in the gas. In this way, the method has an advantage compared with bimolecular classical (trajectory) treatments. The calculations are carried out with a new empirical Ar–Kr pair potential. This has been obtained through extensive analysis of experimental thermophysical and transport properties. We also present a new high level ab initio Ar–Kr potential curve for comparison, as well as ab initio interaction-induced dipole curves computed with different methods. In addition, the Ar–Kr polarizability and hyperpolarizability are reported. A comparison of the computed absorption spectra with an experiment taken at 300 K shows satisfactory agreement although a difference in absolute magnitude of 10%–15% persists. This discrepancy we attribute mainly to experimental uncertainty.
With the reformulation of the classical equations of collision-induced absorption, we present a method to perform the direct computation of the spectral density function. This way the absorption coefficient can be computed from classical molecular dynamics (MD) without the computationally demanding evaluation of the dipole autocorrelation function. In addition, we have developed an algorithm to extract the bound-to-bound dimer contribution to the MD simulated absorption. The method has been tested on the Ar–Xe rare gas system. Comparisons with quantum mechanical (QM) and conventional MD methods validate the approach. The obtained MD bound-to-bound spectra generally agree in shape and magnitude with QM results, including features stemming from rotations and vibrations of the Ar–Xe dimer.
In accordance with modern regulations and directives, the use of renewable biomass materials as precursors for the production of fuels for transportation purposes is to be strictly followed. Even though, there are problems related to processing, storage and handling in wide range of subsequent uses, since there must be a limit to the ratio of biofuels mixed with mineral raw materials. As a key factor with regards to these biomass sources pose a great risk of causing multiple forms of corrosion both to metallic and non-metallic structural materials. To assess the degree of corrosion risk to a variety of engineering alloys like low-carbon and stainless steels widely used as structural metals, this work is dedicated to investigating corrosion rates of economically reasonable engineering steel alloys in mixtures of raw gas oil and renewable biomass fuel sources under typical co-processing conditions. To model a desulphurising refining process, corrosion tests were carried out with raw mineral gasoline and its mixture with used cooking oil and animal waste lard in relative quantities of 10% (g/g). Co-processing was simulated by batch-reactor laboratory experiments. Experiments were performed at temperatures between 200 and 300ºC and a pressure in the gas phase of 90 bar containing 2% (m3/m3) hydrogen sulphide. The time span of individual tests were varied between 1 and 21 days so that we can conclude about changes in the reaction rates against time exposure of and extrapolate for longer periods of exposure. Initial and integral corrosion rates were defined by a weight loss method on standard size of coupons of all sorts of steel alloys. Corrosion rates of carbon steels indicated a linear increase with temperature and little variation with composition of the biomass fuel sources. Apparent activation energies over the first 24-hour period remained moderate, varying between 35.5 and 50.3 kJ mol−1. Scales developed on carbon steels at higher temperatures were less susceptible to spall and detach. Nonetheless, moderate deceleration of corrosion rates as a function of time are due to the less coherent, frequently spalling and low compactness, higher porosity of the scales evolved at lower and higher temperatures, respectively. On the surface of high alloy steels, sulphide scales of an enhanced barrier nature formed during the induction periods and the layer formation mechanism was found to be assisted by the increasing temperature as initial reaction rates considerably decreased over time. Nevertheless, corrosion-related sulphide conversion of metals and mass loss of the high alloys are strongly affected by the composition of the biomass fuel sources especially animal waste lard. Thermal activation in the first 24 hours decreased from 68.9 to 35.2 kJ mol−1. A greater degree of failure to develop protective sulphide scales was experienced by changing to composition of the biomass fuel sources than the impact of thermal activation between a narrow temperature range at around 100ºC. In accordance with the literature, high free fatty acid contents lead to high corrosion rates accounted for direct corrosion of high alloy steels and assisted solubilisation of corrosion products. In addition, the pronounced acceleration of sulphide corrosion is connected to the diminishing inhibition effect of the sulphide scales
Reactions of HBr with radicals are involved in atmospheric chemistry and in the mechanism of operation of bromine-containing flame retardants. The rate coefficients for two such reactions, HBr + OH and HBr + CH3, are available from earlier experiments at near or below room temperature, relevant for atmospheric chemistry, and in this domain, the activation energy for both has been found to be negative. However, no experimental data are available at combustion temperatures. In this work, to provide reliable data needed for modeling the action of brominated flame suppressants, we used the quasiclassical trajectory (QCT) method in combination with high-level ab initio potential energy surfaces to evaluate the rate coefficients of the two title reactions at combustion temperatures. The QCT calculations have been validated by reproducing the experimental rate coefficients at room temperature. At temperatures between 600 and 3200 K, the QCT rate coefficients display positive activation energies. We recommend the following extended Arrhenius expressions to describe the temperature dependence of the thermal rate coefficients: k6 = (9.86 ± 2.38) × 10–16T(1.23±0.03) exp[(5.93 ± 0.33) kJ mol–1/RT] cm3 molecule–1 s–1 for the HBr + OH → H2O + Br reaction, and k–2 = (4.06 ± 2.72) × 10–18T(1.83±0.08) exp[(7.53 ± 0.18) kJ mol–1/RT] cm3 molecule–1 s–1 for the HBr + CH3 → CH4 + Br reaction. The latter is in very good agreement with the formula proposed by Burgess et al. [Burgess, D. R., Jr.; Babushok, V. I.; Linteris, G. T.; Manion, J. A. A Chemical Kinetic Mechanism for 2-Bromo-3,3,3-trifluoropropene (2-BTP) Flame Inhibition. Int. J. Chem. Kinet. 2015, 47, 533−619, DOI: 10.1002/kin.20923]. The conventional transition state theory has been tested against the rate data obtained by the QCT method and was found to overestimate not only the rate coefficients but also the activation energies
Heterogeneous palladium catalysts were prepared by the immobilisation of palladium on supported ionic liquid phases (SILP) obtained by physisorption or grafting. They were characterised by different techniques such as 13C CP MAS NMR, 31P CP MAS NMR, FT-IR, XPS and ICP. The catalysts were used efficiently under solvent-free conditions. They led to the formation of amide products with good to excellent selectivity. Recyclability at 5–10 bar CO pressure has also been demonstrated. The same reusable catalysts were proved to lead to the selective formation of double carbonylation products under higher pressure in DMF.
Silica modified with organic dicationic moieties proved to be an excellent support for palladium catalysts used in the aminocarbonylation of aryl iodides. By an appropriate choice of the reaction conditions, the same catalyst could be used for selective mono- or double carbonylations leading to amide and α-ketoamide products, respectively. The best catalyst could be recycled for at least 10 consecutive runs with a loss of palladium below the detection limit. By the application of the new support, efficient catalyst recycling could be achieved under mild reaction conditions (under low pressure and in a short reaction time). Palladium-leaching data support a mechanism with dissolution—re-precipitation of the active palladium species.
A semiclassical method based on surface-hopping techniques is developed to model the dynamics of radiative association with electronic transitions. It can be proven that this method is an extension of the established semiclassical formula used in the characterization of diatomic molecule-formation. Our method is tested for diatomic molecules. It gives the same cross sections as the former semiclassical formula but, contrary to the former method, it allows us to follow the fate of the trajectories after the emission of a photon. This means that we can characterize the rovibrational states of the stabilized molecules. Using semiclassical quantization, we can obtain quantum state-resolved cross sections or emission spectra for the radiative association process. The calculated semiclassical state-resolved spectra show general agreement with the result of quantum mechanical perturbation theory. Furthermore, our surface-hopping model is not only applicable for the description of radiative association but it can be used for semiclassical characterization of any molecular process where spontaneous emission occurs
The rate coefficients and the cross-sections for the formation of imidogen (NH) molecule (and its isotopologues: 15NH and ND) through radiative association are determined by employing quantum mechanical perturbation theory, classical Larmor formula, and Breit–Wigner theory. We suggest the radiative association process as possible route for NH production in diffuse interstellar clouds.
We have developed the polyatomic extension of the established [M. Gustafsson, J. Chem. Phys. 138, 074308 (2013)] classical theory of radiative association in the absence of electronic transitions. The cross section and the emission spectrum of the process is calculated by a quasiclassical trajectory method combined with the classical Larmor formula which can provide the radiated power in collisions. We have also proposed a Monte Carlo scheme for efficient computation of ro-vibrationally quantum state resolved cross sections for radiative association. Besides the method development, the global potential energy and dipole surfaces for H + CN collisions have been calculated and fitted to test our polyatomic semiclassical method.
Cross sections and rate coefficients for the formation of BeH+ and BeD+ molecules in Be+ + H/D collisions through radiative association are calculated using quantum mechanical perturbation theory and Breit-Wigner theory. The local thermodynamic equilibrium limit of the molecule formation is also studied, since the process is also relevant in environments with high-density and/or strong radiation fields. The obtained rate coefficients may facilitate the kinetic modelling of BeH+/BeD+ production in astrochemical environments as well as the corrosion chemistry of thermonuclear fusion reactors.