Graphene as Nanocarrier for Gold(I)-Monocarbene Complexes: Strength and Nature of Physisorption.

Gold(I) metal complexes are finding increasing applications as therapeutic agents against a variety of diseases. As their potential use as effective metallodrugs is continuously confirmed, the issue of their administration, distribution and delivery to desired biological targets emerges. Graphene and its derivatives possess attractive properties in terms of high affinity and low toxicity, suggesting that they can efficaciously be used as drug nanocarriers. In the present study, we computationally address the adsorption of a gold(I) N-heterocyclic monocarbene, namely, IMeAuCl (where IMe = 1,3-dimethylimidazol-2-ylidene), on graphene. The Au(I) N-heterocyclic carbene family has indeed shown promising anticancer activity and the N-heterocyclic ring could easily interact with planar graphene nanostructures. By means of high-level electronic structure approaches, we investigated the strength and nature of the involved interaction using small graphene prototypes, which allow us to benchmark the best-performing DFT functionals as well as assess the role of the different contributions to total interaction energies. Moreover, realistic adsorption enthalpies and free energy values are obtained by exploiting the optimal DFT method to describe the drug adsorption on larger graphene models. Such values (ΔHads = -18.4 kcal/mol and ΔGads= -7.20 kcal/mol for the largest C150H30 model) indicate a very favorable adsorption, mainly arising from the dispersion component of the interaction, with the electrostatic attraction also playing a non-negligible role.

Possible Formation and Destruction of the OD+ Ions in the Interstellar Medium.

The OH+ ion is an important constituent of the interstellar medium (ISM). It can be used as a probe of cosmic ray and X-ray ionization rates in molecular clouds as well as a tracer of oxygen chemistry. The deuterated variant of OH+, the OD+ ion, may also be present in the ISM despite the fact that it has not been detected yet. In this paper, we aim at providing quantitative insight into the OD+ chemistry and at accurately studying the possible formation and destruction processes of OD+ in the ISM. We study the formation and destruction of OD+ through the O+ + HD → OD+ + H and OH+ + D ↔ OD+ + H reactions that can occur in diffuse ISM. Reactive rate constants have been obtained from exact state-to-state quantum wave packet calculations for temperatures ranging from 10 to 1000 K. The new theoretical data are validated through a detailed comparison with available experimental data. The formation of OD+ is found to be less efficient than that of OH+. As a first application, the OD+/OH+ abundance ratio in ISM has been evaluated from a simple astrochemical model, and we found that this ratio can be larger than D/H abundance ratio only at low temperatures. These calculations may help in an astrochemical search of OD+ in cold ISM.

Quantum Effects on the D + H3+ → H2D+ + H Deuteration Reaction and Isotopic Variants.

The title reaction and its isotopic variants are studied using quasi-classical trajectory (QCT) (without taking into account corrections to account for the possible zero point energy breakdown) and ring polymer molecular dynamics (RPMD) methods with a full dimensional and accurate potential energy surface which presents an exchange barrier of approximately 0.144 eV. The QCT rate constant increases when the temperature decreases from 1500 to 10 K. On the contrary, the RPMD rate constant decreases with decreasing temperature, in semiquantitative agreement with recent experimental results. The present RPMD results are in between the thermal and translational experimental rate constants, extracted from the measured data to eliminate the initial vibrational excitation of H3+, obtained in an arc discharge. The difference between the present RPMD results and experimental values is attributed to the possible existence of non thermal vibrational excitation of H3+, not completely removed by the semiempirical model used for the analysis of the experimental results. Also, it is found that, below 200 K, the RPMD trajectories are trapped, forming long-lived collision complexes, with lifetimes longer than 1 ns. These collision complexes can fragment by either redissociating back to reactants or react to products, in the two cases tunneling through the centrifugal and reaction barriers, respectively. The contribution of the formation of the complex to the total deuteration rate should be calculated with more accurate quantum methods, as has been found recently for reactions of larger systems, and the present four atoms system is a good candidate to benchmark the adequacy of RPMD method at temperatures below 100 K.

Quantum and quasi-classical calculations for the S⁺ + H2(v,j) → SH⁺(v',j') + H reactive collisions.

State-to-state cross-sections for the S(+) + H2(v,j) → SH(+)(v',j') + H endothermic reaction are obtained using quantum wave packet (WP) and quasi-classical (QCT) methods for different initial ro-vibrational H2(v,j) over a wide range of translation energies. The final state distribution as a function of the initial quantum number is obtained and discussed. Additionally, the effect of the internal excitation of H2 on the reactivity is carefully studied. It appears that energy transfer among modes is very inefficient that vibrational energy is the most favorable for the reaction, and rotational excitation significantly enhances the reactivity when vibrational energy is sufficient to reach the product. Special attention is also paid to an unusual discrepancy between classical and quantum dynamics for low rotational levels while agreement improves with rotational excitation of H2. An interesting resonant behaviour found in WP calculations is also discussed and associated with the existence of roaming classical trajectories that enhance the reactivity of the title reaction. Finally, a comparison with the experimental results of Stowe et al. for S(+) + HD and S(+) + D2 reactions exhibits a reasonably good agreement with those results.

Ab initio studies of the Rg-NO(+)(X(1)Σ(+)) van der Waals complexes (Rg = He, Ne, Ar, Kr, and Xe).

We used the explicitly correlated variant of the coupled clusters method with single, double, and noniterative triple excitations [CCSD(T)-F12] to compute two-dimensional potential energy surfaces of van der Waals complexes formed by rare gas atoms (Rg) and NO(+)(X(1)Σ(+)) cations. We used the correlation-consistent, triple-zeta (cc-pVTZ-F12) atomic basis sets, and for Kr and Xe rare gases, we employed corresponding pseudopotential cc-pVTZ-PP-F12 atomic basis sets. These basis sets were additionally augmented with mid-bond functions. The complexes are all of skewed T-shape type with Rg atom being closer to the N-side. Using analytical representation of the potentials, we have estimated zero-point energy corrected dissociation energies from anharmonic calculations with BOUND program and also from the harmonic approximation. The binding energies increase with the polarization of the Rg atom in series from He to Xe and are 196 cm(-1), 360 cm(-1), 1024 cm(-1), 1434 cm(-1), and 2141 cm(-1), respectively. Their corresponding dissociation energies are 132 cm(-1), 300 cm(-1), 927 cm(-1), 1320 cm(-1), and 1994 cm(-1) for the complexes with He to Xe, respectively. We find good agreement with previous theoretical and experimental results. The harmonic vibrational frequencies were calculated for the bending and stretching modes of the Rg-NO(+) complexes.

Exchange and Inelastic OH(+) + H Collisions on the Doublet and Quartet Electronic States.

The exchange and inelastic state-to-state cross sections for the OH(+) + H collisions are computed from wave packet calculations using the doublet and quartet ground electronic potential energy surface (PES) correlating to the open shell reactants, for collision energies in the range of 1 meV to 0.7 eV. The doublet PES presents a deep insertion well, of ≈6 eV, but the exchange reaction has a rather low probability, showing that the mechanism is not statistical. This well is also responsible of a rather high rotational energy transfer, which makes the rigid-rotor approach overestimate the cross section for low Δj transitions and for high collisonal energies. The quartet PES, with a much shallower well, also presents a low exchange reaction cross section, but the inelastic state-to-state cross sections are very well reproduced by rigid-rotor calculations. When the electronic partition is used to obtain the total state-to-state cross section, the contribution of the doublet state becomes small, and the resulting total cross sections become close to those obtained for the quartet state. Thus, the total (quartet and doublet) cross sections for this open shell system can be reproduced rather satisfactorily by those obtained with the rigid-rotor approximation on the quartet state. Finally, we compare the new OH(+)-H cross sections with OH(+)-He ones recently computed. We found significant differences, especially for transitions with large Δj showing that specific OH(+)-H calculations had to be performed to accurately analyze the OH(+) emission from interstellar molecular clouds.

Quantum mechanical calculations of state-to-state cross sections and rate constants for the F + DCl → Cl + DF reaction.

We present accurate state-to-state quantum wave packet calculations of integral cross sections and rate constants for the title reaction. Calculations are carried out on the best available ground 1(2)A' global adiabatic potential energy surface of Deskevich et al. [J. Chem. Phys. 124, 224303 (2006)]. Converged state-to-state reaction cross sections have been calculated for collision energies up to 0.5 eV and different initial rotational and vibrational excitations, DCl(v = 0, j = 0 - 1; v = 1, j = 0). Also, initial-state resolved rate constants of the title reaction have been calculated in a temperature range of 100-400 K. It is found that the initial rotational excitation of the DCl molecule does not enhance reactivity, in contract to the reaction with the isotopologue HCl in which initial rotational excitation produces an important enhancement. These differences between the isotopologue reactions are analyzed in detail and attributed to the presence of resonances for HCl(v = 0, j), absent in the case of DCl(v = 0, j). For vibrational excited DCl(v = 1, j), however, the reaction cross section increases noticeably, what is also explained by another resonance.

Structures and optical properties of zinc oxide nanoclusters: a combined experimental and theoretical approach.

CONTEXT: In this study, theoretical and experimental analysis of the electrical, optical, and structural properties of a wurtzite-like zinc oxide (ZnO) nanostructure has been done. To investigate how quantum confinement affects the optical characteristics, two distinct ZnO clusters in nanowire structures have been investigated. The [(ZnO)55(H2O)4] system's HOMO-LUMO band gap (BG) was calculated to be 2.99 eV, which is quite close to the experimental measurement. It was found that the BG decreases with the increase in the number of atoms in the cluster in connection with the quantum confinement in nanoclusters. In addition, the lowest excitation energy in TD-DFT calculations of the identical system is in fairly good agreement with the experimental value with a difference of 0.1 eV. We conclude that the CAM-B3LYP functional has highly successful in reproducing the experimental data reported in the present study and previously reported experimental data. METHODS: The geometrical optimization of two different sizes of ZnO clusters ([(ZnO)25(H2O)4] and [ZnO)55(H2O)4]) was performed using the CAM-B3LYP functional with no symmetry constraints applied in the gas phase. LANL2DZ basis sets were used for the Zinc (Zn) atom and 6-31G* basis sets for the O and H atoms. To determine their optical and electronic properties, excited state calculations of the pre-optimized structures were performed using the Time-Dependent DFT (TD-DFT) method. Multiwfn, Gaussum 3.0, and GaussView 5.0 programs were used to visualize the results.

Computational Study of Doping in Dopamine with Halogens to Control Optical and Spectroscopic Properties.

In this research, a comprehensive study of dopamine was conducted using the theoretical first principles method due to its crucial importance as a hormone for the neurotransmission process in the animal body. Many basis sets and functionals were used for optimization of the compound to attain stability and find the appropriate energy point for the overall calculations. Then, the compound was doped with the first three members of the halogen family (fluorine, chlorine, and bromine) to analyze the effect of their presence in terms of change in their electronic properties, such as band gap and density of states, and spectroscopic parameters, such as nuclear magnetic resonance and Fourier transform infrared. It was found that the band gap of the system changes depending on the doping of halogens.

Accurate quantum wave packet calculations for the F + HCl → Cl + HF reaction on the ground 1(2)A' potential energy surface.

We present converged exact quantum wave packet calculations of reaction probabilities, integral cross sections, and thermal rate coefficients for the title reaction. Calculations have been carried out on the ground 1(2)A' global adiabatic potential energy surface of Deskevich et al. [J. Chem. Phys. 124, 224303 (2006)]. Converged wave packet reaction probabilities at selected values of the total angular momentum up to a partial wave of J = 140 with the HCl reagent initially selected in the v = 0, j = 0-16 rovibrational states have been obtained for the collision energy range from threshold up to 0.8 eV. The present calculations confirm an important enhancement of reactivity with rotational excitation of the HCl molecule. First, accurate integral cross sections and rate constants have been calculated and compared with the available experimental data.

Ab Initio Study on Dopant Relaxation Mechanism in Ti and Ce Cationically Substituted in Wurtzite Gallium Nitride.

The changes in properties of materials upon introduction of impurities is well documented but less is known about the location of foreign atoms in different hosts. This study is carried out with the motivation to explore dopant location in hexagonal GaN using density functional theory based calculations. The dopant site location of the individual dopants Ti, Ce, and Ti-Ce codoped wurtzite GaN was investigated by placing the dopants at cationic lattice sites as well as off-cationic sites along the c-axis. The geometry optimization relaxed individual dopants on cationic Ga sites but in the case of codoping Ce settled at site 7.8% away along [0001 ¯] and Ti adjusted itself at site 14% away along [0001] from regular cationic sites. The analysis of the results indicates that optimized geometry is sensitive to the starting position of the dopants. The magnetic exchange interactions between Ti and Ce ions are responsible for their structural relaxation in the matrix.

Structural and Electronic Properties of SnO Downscaled to Monolayer.

Two-dimensional (2D) SnO is a p-type semiconductor that has received research and industrial attention for device-grade applications due to its bipolar conductivity and transparent semiconductor nature. The first-principles investigations based on the generalized gradient approximation (GGA) level of theory often failed to accurately model its structure due to interlayer Van der Waals interactions. This study is carried out to calculate structural and electronic properties of bulk and layered structures of SnO using dispersion correction scheme DFT+D3 with GGA-PBE to deal with the interactions which revealed good agreement of the results with reported data. The material in three-dimensional bulk happened to be an indirect gap semiconductor with a band gap of 0.6 eV which is increased to 2.85 eV for a two-dimensional monolayer structure. The detailed analysis of the properties demonstrated that the SnO monolayer is a promising candidate for future optoelectronics and spintronics devices, especially thin film transistors.

Halogen Doping to Control the Band Gap of Ascorbic Acid: A Theoretical Study.

Ascorbic acid is an important antioxidant agent that acts as an electron donor and is involved in many physiological processes. Structural modification in ascorbic acid is a subject of extensive biochemical research due to its involvement in a variety of relevant phenomena including electron transport, complex redox reactions, neurochemical reactions, enzymatic reactions, and chemotherapeutic potential. In this work, the structure of ascorbic acid is modified via doping with the first three members of the halogen group to investigate the changes in the electronic structure and spectroscopic parameters using first-principles methods. To obtain the lowest-energy structures, different basis sets in density functional theory (DFT) and Hartree-Fock approaches were employed in the geometry optimization process. The potential energy maps of the structures were computed to study the molecular orientations and their optical and electrical properties. The spectroscopic properties were computed via UV-vis and nuclear magnetic resonance (NMR) spectroscopies to study the effects of doping into the compound. To obtain further insights into the chemical structure, the Fourier transform infrared (FT-IR) spectra of the materials were theoretically investigated. It was found that the band gap is sensitive to doping as we moved from fluorine to chlorine and then to bromine.

Effects of Yttrium Doping on Erbium-Based Hydroxyapatites: Theoretical and Experimental Study.

This is the first investigation of yttrium (Y) and erbium (Er) co-doped hydroxyapatite (HAp) structures, conducted using theoretical and experimental procedures. By using a wet chemical method, the materials were synthesized by varying the concentration of Y amounts of 0.13, 0.26, 0.39, 0.52, 0.65, and 0.78 at.% every virtual 10 atoms of calcium, whereas Er was kept fixed at 0.39 at.%. Spectroscopic, thermal, and in vitro biocompatibility testing were performed on the generated samples. Theoretical calculations were carried out to compute the energy bandgap, density of states, and linear absorption coefficient. The effects of Y concentration on thermal, morphological, and structural parameters were investigated in detail. Raman and Infrared (FTIR) spectroscopies confirmed the formation of the HAp structure in the samples. Theoretical investigations indicated that the increasing amount of Y increased the density from 3.1724 g cm-3 to 3.1824 g cm-3 and decreased the bandgap energy from 4.196 eV to 4.156 eV, except for the sample containing 0.39 at. % of the dopant, which exhibited a decrease in the bandgap. The values of linear absorption appeared reduced with an increase in photon energy. The samples exhibited cell viability higher than 110%, which revealed excellent biocompatibility for biological applications of the prepared samples.

Accurate time dependent wave packet calculations for the N + OH reaction.

We present accurate quantum calculations of state-to-state cross sections for the N + OH → NO + H reaction performed on the ground (3)A'' global adiabatic potential energy surface of Guadagnini et al. [J. Chem. Phys. 102, 774 (1995)]. The OH reagent is initially considered in the rovibrational state ν = 0, j = 0 and wave packet calculations have been performed for selected total angular momentum, J = 0, 10, 20, 30, 40,...,120. Converged integral state-to-state cross sections are obtained up to a collision energy of 0.5 eV, considering a maximum number of eight helicity components, Ω = 0,...,7. Reaction probabilities for J = 0 obtained as a function of collision energy, using the wave packet method, are compared with the recently published time-independent quantum mechanical one. Total reaction cross sections, state-specific rate constants, opacity functions, and product state-resolved integral cross-sections have been obtained by means of the wave packet method for several collision energies and compared with recent quasi-classical trajectory results obtained with the same potential energy surface. The rate constant for OH(ν = 0, j = 0) is in good agreement with the previous theoretical values, but in disagreement with the experimental data, except at 300 K.

Time-dependent wave packet and quasiclassical trajectory study of the C(3P) + OH(X 2Pi) --> CO(X 1Sigma+) + H(2S) reaction at the state-to-state level.

The first calculations of state-to-state reaction probabilities and product state-resolved integral cross sections at selected collision energies (0.05, 0.1, 0.5, and 1.0 eV) for the title reaction on the ab initio potential energy surface of [Zanchet et al. J. Phys. Chem. A 110, 12017 (2006)] with the OH reagent in selected rovibrational states (v = 0-2, j = 0-5) have been carried out by means of the real wave packet (RWP) and quasiclassical trajectory (QCT) methods. State-selected total reaction probabilities have been calculated for total angular momentum J = 0 in a broad range of collision energies. Integral cross sections and state-specific rate coefficients have been obtained from the corresponding J = 0 RWP reaction probabilities for initially selected rovibrational states by means of a capture model. The calculated RWP and QCT state-selected rate coefficients are practically temperature independent. Both RWP and QCT reaction probabilities, integral cross sections, and rate coefficients are almost independent of the initial rotational excitation. The RWP results are found to be in an overall good agreement with the corresponding QCT results. The present results have been compared with earlier wave packet calculations carried out on the same potential energy surface.

The effect of simulating body fluid on the structural properties of hydroxyapatite synthesized in the presence of citric acid.

In present work, the effect of citric acid (CA) addition in different amounts (0, 1, 5 and 10 ml) on the structure of hydroxyapatite (HAp) was investigated using X-ray diffraction (XRD), Fourier transform infrared (FTIR) spectroscopy, scanning electron microscopy (SEM) and energy dispersive X-ray (EDX) spectroscopy techniques. The crystallite dimensions, lattice parameters, unit cell volume, crystallinity percentage and Ca/P molar ratio were found to be affected by the CA content. To investigate the influence of CA on the bioactive properties of the HAp samples and to determine the optimum amount of CA, in vitro soaking tests in simulated body fluid (SBF) were performed. Although the samples' morphology was found to be affected by neither the amount of CA nor the soaking time in SBF, the soaking results revealed that the maximum changes in the Ca/P ratio were found for the HAp samples prepared in the presence of the highest amounts of CA, which pointed out to the highest bioactivity of these samples.

On the dynamics of the H+ +D2(v=0,j=0)-->HD+D + reaction: a comparison between theory and experiment.

The H+ +D2(v=0,j=0)-->HD+D + reaction has been theoretically investigated by means of a time independent exact quantum mechanical approach, a quantum wave packet calculation within an adiabatic centrifugal sudden approximation, a statistical quantum model, and a quasiclassical trajectory calculation. Besides reaction probabilities as a function of collision energy at different values of the total angular momentum, J, special emphasis has been made at two specific collision energies, 0.1 and 0.524 eV. The occurrence of distinctive dynamical behavior at these two energies is analyzed in some detail. An extensive comparison with previous experimental measurements on the Rydberg H atom with D2 molecules has been carried out at the higher collision energy. In particular, the present theoretical results have been employed to perform simulations of the experimental kinetic energy spectra.

Time-Dependent Quantum Wave Packet Calculations of Three-Dimensional He - O2 Inelastic Scattering.

We have studied a three-dimensional time-dependent quantum dynamics of He - O2 inelastic scattering by using a recently published ab initio potential energy surface. The state-to-state transition probabilities at zero total angular momentum have been calculated in the energy range of 0.12-0.59 eV, and the product rotational distributions are extracted. J-shifting approximation is used to estimate the probabilities for J > 0. The integral cross sections and thermal rate constants are then calculated.

A detailed quantum mechanical and quasiclassical trajectory study on the dynamics of the H+ + H2 --> H2 + H+ exchange reaction.

The H+ + H2 exchange reaction has been studied theoretically by means of a different variety of methods as an exact time independent quantum mechanical, approximate quantum wave packet, statistical quantum, and quasiclassical trajectory approaches. Total and state-to-state reaction probabilities in terms of the collision energy for different values of the total angular momentum obtained with these methods are compared. The dynamics of the reaction is extensively studied at the collision energy of E(coll)=0.44 eV. Integral and differential cross sections and opacity functions at this collision energy have been calculated. In particular, the fairly good description of the exact quantum results provided by the statistical quantum method suggests that the dynamics of the process is governed by an insertion mechanism with the formation of a long-lived collision complex.